Affinity improvement of a therapeutic antibody by structure-based computational design: generation of electrostatic interactions in the transition state stabilizes the antibody-antigen complex

Designing and bioengineering of CDRs with higher affinity against receptor-binding domain (RBD) of SARS-CoV-2 Omicron variant

The emergence of the SARS-CoV-2 Omicron variant highlights the need for innovative strategies to address evolving viral threats. This study bioengineered three nanobodies H11-H4, C5, and H3 originally targeting the Wuhan RBD, to bind more effectively to the Omicron RBD. A structure-based in silico affinity maturation pipeline was developed to enhance their binding affinities. The pipeline consists of three key steps: high-throughput in silico mutagenesis of complementarity-determining regions (CDRs), protein-protein docking for screening, and molecular dynamics (MD) simulations for assessment of the complex stability. A total of 741, 551, and 684 mutations were introduced in H11-H4, C5, and H3 nanobodies, respectively. Protein-protein docking and MD simulations shortlisted high-affinity mutants for H11-H4(6), C5(5), and H3(6). Further, recombinant production of H11-H4 mutants and Omicron RBD enabled experimental validation through Isothermal Titration Calorimetry (ITC). The H11-H4 mutants R27E, S57D, S107K, D108W, and A110I exhibited improved binding affinities with dissociation constant (KD) values ranging from ~8.8 to ~27 μM, compared to the H11-H4 nanobody KD of ~32 μM, representing a three-fold enhancement. This study demonstrates the potential of the developed in silico affinity maturation pipeline as a rapid, cost-effective method for repurposing nanobodies, aiding the development of robust prophylactic strategies against evolving SARS-CoV-2 variants and other pathogens.

Structure-guided mutations in CDRs for enhancing the affinity of neutralizing SARS-CoV-2 nanobody

The optimization of antibodies to attain the desired levels of affinity and specificity holds great promise for the development of next generation therapeutics. This study delves into the refinement and engineering of complementarity-determining regions (CDRs) through in silico affinity maturation followed by binding validation using isothermal titration calorimetry (ITC) and pseudovirus-based neutralization assays. Specifically, it focuses on engineering CDRs targeting the epitopes of receptor-binding domain (RBD) of the spike protein of SARS-CoV-2. A structure-guided virtual library of 112 single mutations in CDRs was generated and screened against RBD to select the potential affinity-enhancing mutations. Protein-protein docking analysis identified 32 single mutants of which nine mutants were selected for molecular dynamics (MD) simulations. Subsequently, biophysical ITC studies provided insights into binding affinity, and consistent with in silico findings, six mutations that demonstrated better binding affinity than native nanobody were further tested in vitro for neutralization activity against SARS-CoV-2 pseudovirus. Leu106Thr mutant was found to be most effective in virus-neutralization with IC50 values of ∼0.03 μM, as compared to the native nanobody (IC50 ∼0.77 μM). Thus, in this study, the developed computational pipeline guided by structure-aided interface profiles and thermodynamic analysis holds promise for the streamlined development of antibody-based therapeutic interventions against emerging variants of SARS-CoV-2 and other infectious pathogens.

Molecular Modeling Methods in the Development of Affine and Specific Protein-Binding Agents

High-affinity and specific agents are widely applied in various areas, including diagnostics, scientific research, and disease therapy (as drugs and drug delivery systems). It takes significant time to develop them. For this reason, development of high-affinity agents extensively utilizes computer methods at various stages for the analysis and modeling of these molecules. The review describes the main affinity and specific agents, such as monoclonal antibodies and their fragments, antibody mimetics, aptamers, and molecularly imprinted polymers. The methods of their obtaining as well as their main advantages and disadvantages are briefly described, with special attention focused on the molecular modeling methods used for their analysis and development.

High-throughput analysis system of interaction kinetics for data-driven antibody design

Surface plasmon resonance (SPR) is widely used for antigen-antibody interaction kinetics analysis. However, it has not been used in the screening phase because of the low throughput of measurement and analysis. Herein, we proposed a high-throughput SPR analysis system named "BreviA" using the Brevibacillus expression system. Brevibacillus was transformed using a plasmid library containing various antibody sequences, and single colonies were cultured in 96-well plates. Sequence analysis was performed using bacterial cells, and recombinant antibodies secreted in the supernatant were immobilized on a sensor chip to analyze their interactions with antigens using high-throughput SPR. Using this system, the process from the transformation to 384 interaction analyses can be performed within a week. This system utility was tested using an interspecies specificity design of an anti-human programmed cell death protein 1 (PD-1) antibody. A plasmid library containing alanine and tyrosine mutants of all complementarity-determining region residues was generated. A high-throughput SPR analysis was performed against human and mouse PD-1, showing that the mutation in the specific region enhanced the affinity for mouse PD-1. Furthermore, deep mutational scanning of the region revealed two mutants with > 100-fold increased affinity for mouse PD-1, demonstrating the potential efficacy of antibody design using data-driven approach.

In Vitro Affinity Maturation of Nanobodies against Mpox Virus A29 Protein Based on Computer-Aided Design

Mpox virus (MPXV), the most pathogenic zoonotic orthopoxvirus, caused worldwide concern during the SARS-CoV-2 epidemic. Growing evidence suggests that the MPXV surface protein A29 could be a specific diagnostic marker for immunological detection. In this study, a fully synthetic phage display library was screened, revealing two nanobodies (A1 and H8) that specifically recognize A29. Subsequently, an in vitro affinity maturation strategy based on computer-aided design was proposed by building and docking the A29 and A1 three-dimensional structures. Ligand-receptor binding and molecular dynamics simulations were performed to predict binding modes and key residues. Three mutant antibodies were predicted using the platform, increasing the affinity by approximately 10-fold compared with the parental form. These results will facilitate the application of computers in antibody optimization and reduce the cost of antibody development; moreover, the predicted antibodies provide a reference for establishing an immunological response against MPXV.

Arginine cluster introduction on framework region in anti-lysozyme antibody improved association rate constant by changing conformational diversity of CDR loops

Antibodies are used for many therapeutic and biotechnological purposes. Because the affinity of an antibody to the antigen is critical for clinical efficacy of pharmaceuticals, many affinity maturation strategies have been developed. Although we previously reported an affinity maturation strategy in which the association rate of the antibody toward its antigen is improved by introducing a cluster of arginine residues into the framework region of the antibody, the detailed molecular mechanism responsible for this improvement has been unknown. In this study, we introduced five arginine residues into an anti-hen egg white lysozyme antibody (HyHEL10) Fab fragment to create the R5-mutant and comprehensively characterized the interaction between antibody and antigen using thermodynamic analysis, X-ray crystallography, and molecular dynamics (MD) simulations. Our results indicate that introduction of charged residues strongly enhanced the association rate, as previously reported, and the antibody-antigen complex structure was almost the same for the R5-mutant and wild-type Fabs. The MD simulations indicate that the mutation increased conformational diversity in complementarity-determining region loops and thereby enhanced the association rate. These observations provide the molecular basis of affinity maturation by R5 mutation.

Affinity maturation of antibody fragments: A review encompassing the development from random approaches to computational rational optimization

Routinely screened antibody fragments usually require further in vitro maturation to achieve the desired biophysical properties. Blind in vitro strategies can produce improved ligands by introducing random mutations into the original sequences and selecting the resulting clones under more and more stringent conditions. Rational approaches exploit an alternative perspective that aims first at identifying the specific residues potentially involved in the control of biophysical mechanisms, such as affinity or stability, and then to evaluate what mutations could improve those characteristics. The understanding of the antigen-antibody interactions is instrumental to develop this process the reliability of which, consequently, strongly depends on the quality and completeness of the structural information. Recently, methods based on deep learning approaches critically improved the speed and accuracy of model building and are promising tools for accelerating the docking step. Here, we review the features of the available bioinformatic instruments and analyze the reports illustrating the result obtained with their application to optimize antibody fragments, and nanobodies in particular. Finally, the emerging trends and open questions are summarized.

How can we discover developable antibody-based biotherapeutics?

Antibody-based biotherapeutics have emerged as a successful class of pharmaceuticals despite significant challenges and risks to their discovery and development. This review discusses the most frequently encountered hurdles in the research and development (R&D) of antibody-based biotherapeutics and proposes a conceptual framework called biopharmaceutical informatics. Our vision advocates for the syncretic use of computation and experimentation at every stage of biologic drug discovery, considering developability (manufacturability, safety, efficacy, and pharmacology) of potential drug candidates from the earliest stages of the drug discovery phase. The computational advances in recent years allow for more precise formulation of disease concepts, rapid identification, and validation of targets suitable for therapeutic intervention and discovery of potential biotherapeutics that can agonize or antagonize them. Furthermore, computational methods for de novo and epitope-specific antibody design are increasingly being developed, opening novel computationally driven opportunities for biologic drug discovery. Here, we review the opportunities and limitations of emerging computational approaches for optimizing antigens to generate robust immune responses, in silico generation of antibody sequences, discovery of potential antibody binders through virtual screening, assessment of hits, identification of lead drug candidates and their affinity maturation, and optimization for developability. The adoption of biopharmaceutical informatics across all aspects of drug discovery and development cycles should help bring affordable and effective biotherapeutics to patients more quickly.

The validation of artificial anti-monkeypox antibodies by in silico and experimental approaches

As a result of smallpox immunization programs that ended more than 40 years ago, a significant portion of the world's population is not immune. Moreover, due to the lack of anti-monkeypox drugs and vaccines against monkeypox, the spread of this virus may be the beginning of another challenge. In this study, novel antibodies against monkeypox virus were modeled based on a heavy chain of human antibody and a small peptide fragment. Docking of modeled antibodies with C19L protein showed the range of docking energy, and root-mean-square deviation (RMSD) was from -124 to -154 kcal/mL and 4-6 angstrom, respectively. Also, docking of modeled antibodies-C19L complex with gamma Fc receptor type I illustrated the range of docking energy, and RMSD was from -132 to -155 kcal/ml and 5-7 angstrom, respectively. Moreover, molecular dynamics simulation showed that antibody 62 had the highest stability with the lowest energy level and RMSD. Interestingly, no modeled antibodies had immunogenicity, allergenicity, and toxicity. Although all of them had good stability, only antibodies 25, 28, 54, and 62 had a half-life of >10 h. Moreover, the interaction between C19L protein and anti-C19L antibodies (wild-type and synthetic) was evaluated by the SPR method. We found that KD in synthetic antibodies was lower than wild antibody. In terms of δH°, TδS°, and δG°, the results were consistent with binding parameters. Here, the lowest value of thermodynamic parameters was obtained for antibody 62. These data show that the synthetic antibodies, especially antibody 62, had a higher affinity than the wild-type antibody.

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.

Structure-Based Affinity Maturation of Antibody Based on Double-Point Mutations

Structure-based site-directed affinity maturation of antibodies can be expanded by multiple-point mutations to obtain various mutants. However, selecting the appropriate number of promising mutants for experimental evaluation from the vast number of combinations of multiple-point mutations is challenging. In this report, we describe how to narrow candidate mutants using the so-called weak interaction analysis such as CH-π and CH-O in addition to widely recognized interactions such as hydrogen bonds.

Analytical Method for Experimental Validation of Computer-Designed Antibody

In the computational design of antibodies, the interaction analysis between target antigen and antibody is an essential process to obtain feedback for validation and optimization of the design. Kinetic and thermodynamic parameters as well as binding affinity (K_D) allow for a more detailed evaluation and understanding of the molecular recognition. In this chapter, we summarize the conventional experimental methods which can calculate K_D value (ELISA, FP), analyze a binding activity to actual cells (FCM), and evaluate the kinetic and thermodynamic parameters (ITC, SPR, BLI), including high-throughput analysis and a recently developed experimental technique.

In Silico Maturation of a Nanomolar Antibody against the Human CXCR2

The steady increase in computational power in the last 50 years is opening unprecedented opportunities in biology, as computer simulations of biological systems have become more accessible and can reproduce experimental results more accurately. Here, we wanted to test the ability of computer simulations to replace experiments in the limited but practically useful scope of improving the biochemical characteristics of the abN48 antibody, a nanomolar antagonist of the CXC chemokine receptor 2 (CXCR2) that was initially selected from a combinatorial antibody library. Our results showed a good correlation between the computed binding energies of the antibody to the peptide target and the experimental binding affinities. Moreover, we showed that it is possible to design new antibody sequences in silico with a higher affinity to the desired target using a Monte Carlo Metropolis algorithm. The newly designed sequences had an affinity comparable to the best ones obtained using in vitro affinity maturation and could be obtained within a similar timeframe. The methodology proposed here could represent a valid alternative for improving antibodies in cases in which experiments are too expensive or technically tricky and could open an opportunity for designing antibodies for targets that have been elusive so far.

Electrostatic Features for the Receptor Binding Domain of SARS-COV-2 Wildtype and Its Variants. Compass to the Severity of the Future Variants with the Charge-Rule

Electrostatic intermolecular interactions are important in many aspects of biology. We have studied the main electrostatic features involved in the interaction of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with the human receptor Angiotensin-converting enzyme 2 (ACE2). As the principal computational tool, we have used the FORTE approach, capable to model proton fluctuations and computing free energies for a very large number of protein-protein systems under different physical-chemical conditions, here focusing on the RBD-ACE2 interactions. Both the wild-type and all critical variants are included in this study. From our large ensemble of extensive simulations, we obtain, as a function of pH, the binding affinities, charges of the proteins, their charge regulation capacities, and their dipole moments. In addition, we have calculated the pK_as for all ionizable residues and mapped the electrostatic coupling between them. We are able to present a simple predictor for the RBD-ACE2 binding based on the data obtained for Alpha, Beta, Gamma, Delta, and Omicron variants, as a linear correlation between the total charge of the RBD and the corresponding binding affinity. This "RBD charge rule" should work as a quick test of the degree of severity of the coming SARS-CoV-2 variants in the future.

RosettaSurf-A surface-centric computational design approach

Proteins are typically represented by discrete atomic coordinates providing an accessible framework to describe different conformations. However, in some fields proteins are more accurately represented as near-continuous surfaces, as these are imprinted with geometric (shape) and chemical (electrostatics) features of the underlying protein structure. Protein surfaces are dependent on their chemical composition and, ultimately determine protein function, acting as the interface that engages in interactions with other molecules. In the past, such representations were utilized to compare protein structures on global and local scales and have shed light on functional properties of proteins. Here we describe RosettaSurf, a surface-centric computational design protocol, that focuses on the molecular surface shape and electrostatic properties as means for protein engineering, offering a unique approach for the design of proteins and their functions. The RosettaSurf protocol combines the explicit optimization of molecular surface features with a global scoring function during the sequence design process, diverging from the typical design approaches that rely solely on an energy scoring function. With this computational approach, we attempt to address a fundamental problem in protein design related to the design of functional sites in proteins, even when structurally similar templates are absent in the characterized structural repertoire. Surface-centric design exploits the premise that molecular surfaces are, to a certain extent, independent of the underlying sequence and backbone configuration, meaning that different sequences in different proteins may present similar surfaces. We benchmarked RosettaSurf on various sequence recovery datasets and showcased its design capabilities by generating epitope mimics that were biochemically validated. Overall, our results indicate that the explicit optimization of surface features may lead to new routes for the design of functional proteins.

A glimpse on metazoan ZNFX1 helicases, ancient players of antiviral innate immunity

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.

Development of a potent high-affinity human therapeutic antibody via novel application of recombination signal sequence-based affinity maturation

Therapeutic antibody development requires discovery of an antibody molecule with desired specificities and drug-like properties. For toxicological studies, a therapeutic antibody must bind the ortholog antigen with a similar affinity to the human target to enable relevant dosing regimens, and antibodies falling short of this affinity design goal may not progress as therapeutic leads. Herein, we report the novel use of mammalian recombination signal sequence (RSS)–directed recombination for complementarity-determining region–targeted protein engineering combined with mammalian display to close the species affinity gap of human interleukin (IL)-13 antibody 731. This fully human antibody has not progressed as a therapeutic in part because of a 400-fold species affinity gap. Using this nonhypothesis-driven affinity maturation method, we generated multiple antibody variants with improved IL-13 affinity, including the highest affinity antibody reported to date (34 fM). Resolution of a cocrystal structure of the optimized antibody with the cynomolgus monkey (or nonhuman primate) IL-13 protein revealed that the RSS-derived mutations introduced multiple successive amino-acid substitutions resulting in a de novo formation of a π–π stacking–based protein–protein interaction between the affinity-matured antibody heavy chain and helix C on IL-13, as well as an introduction of an interface-distant residue, which enhanced the light chain–binding affinity to target. These mutations synergized binding of heavy and light chains to the target protein, resulting in a remarkably tight interaction, and providing a proof of concept for a new method of protein engineering, based on synergizing a mammalian display platform with novel RSS-mediated library generation.

An integrated computational pipeline for designing high-affinity nanobodies with expanded genetic codes

Protein engineering and design principles employing the 20 standard amino acids have been extensively used to achieve stable protein scaffolds and deliver their specific activities. Although this confers some advantages, it often restricts the sequence, chemical space, and ultimately the functional diversity of proteins. Moreover, although site-specific incorporation of non-natural amino acids (nnAAs) has been proven to be a valuable strategy in protein engineering and therapeutics development, its utility in the affinity-maturation of nanobodies is not fully explored. Besides, current experimental methods do not routinely employ nnAAs due to their enormous library size and infinite combinations. To address this, we have developed an integrated computational pipeline employing structure-based protein design methodologies, molecular dynamics simulations and free energy calculations, for the binding affinity prediction of an nnAA-incorporated nanobody toward its target and selection of potent binders. We show that by incorporating halogenated tyrosines, the affinity of 9G8 nanobody can be improved toward epidermal growth factor receptor (EGFR), a crucial cancer target. Surface plasmon resonance (SPR) assays showed that the binding of several 3-chloro-l-tyrosine (3MY)-incorporated nanobodies were improved up to 6-fold into a picomolar range, and the computationally estimated binding affinities shared a Pearson's r of 0.87 with SPR results. The improved affinity was found to be due to enhanced van der Waals interactions of key 3MY-proximate nanobody residues with EGFR, and an overall increase in the nanobody's structural stability. In conclusion, we show that our method can facilitate screening large libraries and predict potent site-specific nnAA-incorporated nanobody binders against crucial disease-targets.

Delicate balance among thermal stability, binding affinity, and conformational space explored by single-domain VHH antibodies

The high binding affinities and specificities of antibodies have led to their use as drugs and biosensors. Single-domain V_HH antibodies exhibit high specificity and affinity but have higher stability and solubility than conventional antibodies as they are single-domain proteins. In this work, based on physicochemical measurements and molecular dynamics (MD) simulations, we have gained insight that will facilitate rational design of single-chain V_HH antibodies. We first assessed two homologous V_HH antibodies by differential scanning calorimetry (DSC); one had a high (64.8 °C) and the other a low (58.6 °C) melting temperature. We then generated a series of the variants of the low stability antibody and analyzed their thermal stabilities by DSC and characterized their structures through MD simulations. We found that a single mutation that resulted in 8.2 °C improvement in melting temperature resulted in binding affinity an order of magnitude lower than the parent antibody, likely due to a shift of conformational space explored by the single-chain V_HH antibody. These results suggest that the delicate balance among conformational stability, binding capability, and conformational space explored by antibodies must be considered in design of fully functional single-chain V_HH antibodies.

Computational and Rational Design of Single-Chain Antibody against Tick-Borne Encephalitis Virus for Modifying Its Specificity

Tick-borne encephalitis virus (TBEV) causes 5−7 thousand cases of human meningitis and encephalitis annually. The neutralizing and protective antibody ch14D5 is a potential therapeutic agent. This antibody exhibits a high affinity for binding with the D3 domain of the glycoprotein E of the Far Eastern subtype of the virus, but a lower affinity for the D3 domains of the Siberian and European subtypes. In this study, a 2.2-fold increase in the affinity of single-chain antibody sc14D5 to D3 proteins of the Siberian and European subtypes of the virus was achieved using rational design and computational modeling. This improvement can be further enhanced in the case of the bivalent binding of the full-length chimeric antibody containing the identified mutation.

Analysis of Binding Modes of Antigen-Antibody Complexes by Molecular Mechanics Calculation

Antibodies are one of the most important protein molecules in biopharmaceutics. Due to the recent advance in technology for producing monoclonal antibodies, many structural data are available on the antigen-antibody complexes. To characterize the molecular interaction in antigen-antibody recognition, we computationally analyzed 500 complex structures by molecular mechanics calculations. The presence of Ser and Tyr is markedly large in the complementarity-determining regions (CDRs). Although Ser is abundant in CDRs, its contribution to the binding score is not large. Instead, Tyr, Asp, Glu, and Arg significantly contribute to the molecular interaction from the viewpoint of the binding score. The decomposition of the binding score suggests that the hydrophilic interaction is predominant in all CDRs compared with the hydrophobic one. The contribution of the heavy chain is larger than that of the light chain. In particular, H2 and H3 largely contribute to the binding interaction. Tyr is a main contributing residue both in H2 and H3. The positively charged residue Arg also significantly contributes to the binding score in H3, while the contribution of Lys is small. The appearance of Ser is remarkable in H2, and Asp is abundant in H3. The non-charged polar residues, Thr, Asn, and Gln, appear much in H2, compared to appearing in H3. The negatively charged residues Asp and Glu significantly contribute to the binding score in H3. The contributions of Phe and Trp are not large in spite that the aromatic residues are capable of making the π-π or CH-π interaction. Gly is commonly abundant both in H2 and H3. The average distance of the shortest direct hydrogen bond between the antigen and antibody is longer than that of the hydrogen bonds observed in the complexes between compounds and their target proteins. Therefore, the antigen-antibody interface is not so tight as the compound-target protein interface. The calculation of shape complementarity is consistent with the result of the hydrogen bonds in that the fitness of the antigen-antibody contact is not so high as that of the compound-target protein contact. There exist many water molecules at the antigen-antibody interface. These findings suggest that Tyr, Asp, Glu, and Arg are rich in H3 and work as major contributors for the interaction with the antigen. Ser, Thr, Asn, and Gln are rich in H2 and support the interaction with enhancing molecular fitness. Gly is helpful in increasing flexibility and geometrical diversity. Because the antigen-antibody binding is fundamentally hydrophilic-driven, the non-polar residues are unfavorable for mediating the contact even for the aromatic residues such as Phe and Trp.

Optimization of therapeutic antibodies

In this review, we have summarized the current landscape of therapeutic antibody optimization for successful development. By engineering antibodies with display technology, computer-aided design and site mutagenesis, various properties of the therapeutic antibody candidates can be improved with the purpose of enhancing their safety, efficacy and developability. These properties include antigen binding affinity and specificity, biological efficacy, pharmacokinetics and pharmacodynamics, immunogenicity and physicochemical developability features. A best-in-class strategy may require the optimization of all these properties to generate a good therapeutic antibody.

IsAb: a computational protocol for antibody design

The design of therapeutic antibodies has attracted a large amount of attention over the years. Antibodies are widely used to treat many diseases due to their high efficiency and low risk of adverse events. However, the experimental methods of antibody design are time-consuming and expensive. Although computational antibody design techniques have had significant advances in the past years, there are still some challenges that need to be solved, such as the flexibility of antigen structure, the lack of antibody structural data and the absence of standard antibody design protocol. In the present work, we elaborated on an in silico antibody design protocol for users to easily perform computer-aided antibody design. First, the Rosetta web server will be applied to generate the 3D structure of query antibodies if there is no structural information available. Then, two-step docking will be used to identify the binding pose of an antibody-antigen complex when the binding information is unknown. ClusPro is the first method to be used to conduct the global docking, and SnugDock is applied for the local docking. Sequentially, based on the predicted binding poses, in silico alanine scanning will be used to predict the potential hotspots (or key residues). Finally, computational affinity maturation protocol will be used to modify the structure of antibodies to theoretically increase their affinity and stability, which will be further validated by the bioassays in the future. As a proof of concept, we redesigned antibody D44.1 and compared it with previously reported data in order to validate IsAb protocol. To further illustrate our proposed protocol, we used cemiplimab antibody, a PD-1 checkpoint inhibitor, as an example to showcase a step-by-step tutorial.

64Cu-labeled minibody D2101 visualizes CDH17-positive gastric cancer xenografts with short waiting time

OBJECTIVE

We previously reported In-labeled anti-cadherin17 (CDH17) IgG visualized CDH17-positive gastric cancer xenografts. Unfortunately, a long waiting time was required to obtain high-contrast images due to long blood retention (blood half-life: 26 h). To accelerate blood clearance, we have developed anti-CDH17 minibody (D2101 minibody) and evaluated the pharmacokinetics in gastric cancer mouse models.

METHODS

Two different single chain Fvs (scFvs), D2101 mutant and D2111, were developed from each parental IgG. The binding ability to CDH17 and stability in plasma were evaluated. D2101 minibody, constructed based on D2101 mutant scFv, was labeled with Cu (Cu-D2101 minibody), and the in-vitro and in-vivo properties were evaluated by cell ELISA, biodistribution experiments, and PET imaging in mice bearing CDH17-positive AGS and CDH17-negative MKN74 tumors.

RESULTS

D2101 mutant and D2111 scFvs showed similar affinities to CDH17. D2101 mutant scFv was more stable than D2111 scFv in plasma. No loss of binding affinity of the D2101 minibody by chelate conjugation and radiolabeling procedures was observed. The biodistribution of Cu-D2101 minibody showed high uptake in AGS tumors and low uptake in MKN74. The blood half-life of Cu-D2101 minibody was 6.5 h. Improved blood clearance of Cu-D2101 minibody provided high tumor-to-blood ratios compared with the previous results of parental IgG in AGS xenograft mice. PET studies showed consistent results with biodistribution studies.

CONCLUSIONS

Cu-D2101 minibody exhibited higher tumor-to-blood ratios at earlier time points than those of the radiolabeled parental IgG. Cu-D2101 minibody has potential as an immunoimaging agent for CDH17-positive tumors.

Generating the Engineered Form of a Nanobody Against Placenta Growth Factor with High Antiangiogenesis Potential

Antibody engineering is a dynamic field in antibody industry. Over 30% of novel monoclonal antibodies (mAbs) in R&D and clinical trials are engineered forms. Affinity enhancement contributes to the development of new binders that are not only effective in low dose and cost but also improve some drawbacks of antibody production. After previous successful work on in silico affinity maturation of nanobody against placenta growth factor and finding the best engineered nanobodies (Mut2:S31D and Mut4:R45E), according to bioinformatic parameters and molecular dynamics (MD) simulation results, in this study we focused on experimental confirmation of affinity enhancement of a mutant form of nanobody. So, we cloned and expressed two of four mutant forms in pHEN6c vector. Affinity binding was assayed by enzyme-linked immunosorbent assay on purified mutants, with results showing that 10-time enhancement in affinity compared with the native form associated with MD simulation results. We checked the effectiveness of these mutant nanobodies in angiogenesis inhibition by human umbilical vein endothelial cell proliferation and 3D capillary tube formation. EC50 of mut2, mut4, and native in proliferation assay was 110, 140, and 190 ng/mL, respectively, and that in 3D capillary tube formation was 80, 83, and 100 ng/mL. The results of functional studies revealed strong effectiveness of Mut2 followed by Mut4 compared with the native form. Our study confirmed that in silico approach could facilitate development of novel versions of mAb with better characteristics, which could save cost and time.

Review: PET imaging with macro- and middle-sized molecular probes

Recent progress in radiolabeling of macro- and middle-sized molecular probes has been extending possibilities to use PET molecular imaging for dynamic application to drug development and therapeutic evaluation. Theranostics concept also accelerated the use of macro- and middle-sized molecular probes for sharpening the contrast of proper target recognition even the cellular types/subtypes and proper selection of the patients who should be treated by the same molecules recognition. Here, brief summary of the present status of immuno-PET, and then further development of advanced technologies related to immuno-PET, peptidic PET probes, and nucleic acids PET probes are described.

Progress in Nanorobotics for Advancing Biomedicine

Nanorobotics, which has long been a fantasy in the realm of science fiction, is now a reality due to the considerable developments in diverse fields including chemistry, materials, physics, information and nanotechnology in the past decades. Not only different prototypes of nanorobots whose sizes are nanoscale are invented for various biomedical applications, but also robotic nanomanipulators which are able to handle nano-objects obtain substantial achievements for applications in biomedicine. The outstanding achievements in nanorobotics have significantly expanded the field of medical robotics and yielded novel insights into the underlying mechanisms guiding life activities, remarkably showing an emerging and promising way for advancing the diagnosis & treatment level in the coming era of personalized precision medicine. In this review, the recent advances in nanorobotics (nanorobots, nanorobotic manipulations) for biomedical applications are summarized from several facets (including molecular machines, nanomotors, DNA nanorobotics, and robotic nanomanipulators), and the future perspectives are also presented.

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.

Structure-based design and discovery of novel anti-tissue factor antibodies with cooperative double-point mutations, using interaction analysis

The generation of a wide range of candidate antibodies is important for the successful development of drugs that simultaneously satisfy multiple requirements. To find cooperative mutations and increase the diversity of mutants, an in silico double-point mutation approach, in which 3D models of all possible double-point mutant/antigen complexes are constructed and evaluated using interaction analysis, was developed. Starting from an antibody with very high affinity, four double-point mutants were designed in silico. Two of the double-point mutants exhibited improved affinity or affinity comparable to that of the starting antibody. The successful identification of two active double-point mutants showed that a cooperative mutation could be found by utilizing information regarding the interactions. The individual single-point mutants of the two active double-point mutants showed decreased affinity or no expression. These results suggested that the two active double-point mutants cannot be obtained through the usual approach i.e. a combination of improved single-point mutants. In addition, a triple-point mutant, which combines the distantly located active double-point mutation and an active single-point mutation collaterally obtained in the process of the double-point mutation strategy, was designed. The triple-point mutant showed improved affinity. This finding suggested that the effects of distantly located mutations are independent and additive. The double-point mutation approach using the interaction analysis of 3D structures expands the design repertoire for mutants, and hopefully paves a way for the identification of cooperative multiple-point mutations.

mCSM-AB2: guiding rational antibody design using graph-based signatures

MOTIVATION

A lack of accurate computational tools to guide rational mutagenesis has made affinity maturation a recurrent challenge in antibody (Ab) development. We previously showed that graph-based signatures can be used to predict the effects of mutations on Ab binding affinity.

RESULTS

Here we present an updated and refined version of this approach, mCSM-AB2, capable of accurately modelling the effects of mutations on Ab-antigen binding affinity, through the inclusion of evolutionary and energetic terms. Using a new and expanded database of over 1800 mutations with experimental binding measurements and structural information, mCSM-AB2 achieved a Pearson's correlation of 0.73 and 0.77 across training and blind tests, respectively, outperforming available methods currently used for rational Ab engineering.

AVAILABILITY AND IMPLEMENTATION

mCSM-AB2 is available as a user-friendly and freely accessible web server providing rapid analysis of both individual mutations or the entire binding interface to guide rational antibody affinity maturation at http://biosig.unimelb.edu.au/mcsm_ab2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein electrostatics: From computational and structural analysis to discovery of functional fingerprints and biotechnological design

Computationally driven engineering of proteins aims to allow them to withstand an extended range of conditions and to mediate modified or novel functions. Therefore, it is crucial to the biotechnological industry, to biomedicine and to afford new challenges in environmental sciences, such as biocatalysis for green chemistry and bioremediation. In order to achieve these goals, it is important to clarify molecular mechanisms underlying proteins stability and modulating their interactions. So far, much attention has been given to hydrophobic and polar packing interactions and stability of the protein core. In contrast, the role of electrostatics and, in particular, of surface interactions has received less attention. However, electrostatics plays a pivotal role along the whole life cycle of a protein, since early folding steps to maturation, and it is involved in the regulation of protein localization and interactions with other cellular or artificial molecules. Short- and long-range electrostatic interactions, together with other forces, provide essential guidance cues in molecular and macromolecular assembly. We report here on methods for computing protein electrostatics and for individual or comparative analysis able to sort proteins by electrostatic similarity. Then, we provide examples of electrostatic analysis and fingerprints in natural protein evolution and in biotechnological design, in fields as diverse as biocatalysis, antibody and nanobody engineering, drug design and delivery, molecular virology, nanotechnology and regenerative medicine.

Antibody interface prediction with 3D Zernike descriptors and SVM

MOTIVATION

Antibodies are a class of proteins capable of specifically recognizing and binding to a virtually infinite number of antigens. This binding malleability makes them the most valuable category of biopharmaceuticals for both diagnostic and therapeutic applications. The correct identification of the antigen-binding residues in the antibody is crucial for all antibody design and engineering techniques and could also help to understand the complex antigen binding mechanisms. However, the antibody-binding interface prediction field appears to be still rather underdeveloped.

RESULTS

We present a novel method for antibody interface prediction from their experimentally solved structures based on 3D Zernike Descriptors. Roto-translationally invariant descriptors are computed from circular patches of the antibody surface enriched with a chosen subset of physico-chemical properties from the AAindex1 amino acid index set, and are used as samples for a binary classification problem. An SVM classifier is used to distinguish interface surface patches from non-interface ones. The proposed method was shown to outperform other antigen-binding interface prediction software.

AVAILABILITY AND IMPLEMENTATION

Linux binaries and Python scripts are available at https://github.com/sebastiandaberdaku/AntibodyInterfacePrediction. The datasets generated and/or analyzed during the current study are available at https://doi.org/10.6084/m9.figshare.5442229.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Directed evolution methods for overcoming trade-offs between protein activity and stability

Engineered proteins are being widely developed and employed in applications ranging from enzyme catalysts to therapeutic antibodies. Directed evolution, an iterative experimental process composed of mutagenesis and library screening, is a powerful technique for enhancing existing protein activities and generating entirely new ones not observed in nature. However, the process of accumulating mutations for enhanced protein activity requires chemical and structural changes that are often destabilizing, and low protein stability is a significant barrier to achieving large enhancements in activity during multiple rounds of directed evolution. Here we highlight advances in understanding the origins of protein activity/stability trade-offs for two important classes of proteins (enzymes and antibodies) as well as innovative experimental and computational methods for overcoming such trade-offs. These advances hold great potential for improving the generation of highly active and stable proteins that are needed to address key challenges related to human health, energy and the environment.

Computer-guided library generation applied to the optimization of single-domain antibodies

Computer-guided library generation is a plausible strategy to optimize antibodies. Herein, we report the improvement of the affinity of a single-domain camelid antibody for its antigen using such approach. We first conducted experimental and computational alanine scanning to describe the precise energetic profile of the antibody-antigen interaction surface. Based on this characterization, we hypothesized that in-silico mutagenesis could be employed to guide the development of a small library for phage display with the goal of improving the affinity of an antibody for its antigen. Optimized antibody mutants were identified after three rounds of selection, in which an alanine residue at the core of the antibody-antigen interface was substituted by residues with large side-chains, generating diverse kinetic responses, and resulting in greater affinity (>10-fold) for the antigen.

Antibody Structure and Function: The Basis for Engineering Therapeutics

Antibodies and antibody-derived macromolecules have established themselves as the mainstay in protein-based therapeutic molecules (biologics). Our knowledge of the structure–function relationships of antibodies provides a platform for protein engineering that has been exploited to generate a wide range of biologics for a host of therapeutic indications. In this review, our basic understanding of the antibody structure is described along with how that knowledge has leveraged the engineering of antibody and antibody-related therapeutics having the appropriate antigen affinity, effector function, and biophysical properties. The platforms examined include the development of antibodies, antibody fragments, bispecific antibody, and antibody fusion products, whose efficacy and manufacturability can be improved via humanization, affinity modulation, and stability enhancement. We also review the design and selection of binding arms, and avidity modulation. Different strategies of preparing bispecific and multispecific molecules for an array of therapeutic applications are included.

Structural and thermodynamic basis for the recognition of the substrate-binding cleft on hen egg lysozyme by a single-domain antibody

Single-domain antibodies (VHHs or nanobodies), developed from heavy chain-only antibodies of camelids, are gaining attention as next-generation therapeutic agents. Despite their small size, the high affinity and specificity displayed by VHHs for antigen molecules rival those of IgGs. How such small antibodies achieve that level of performance? Structural studies have revealed that VHHs tend to recognize concave surfaces of their antigens with high shape-complementarity. However, the energetic contribution of individual residues located at the binding interface has not been addressed in detail, obscuring the actual mechanism by which VHHs target the concave surfaces of proteins. Herein, we show that a VHH specific for hen egg lysozyme, D3-L11, not only displayed the characteristic binding of VHHs to a concave region of the surface of the antigen, but also exhibited a distribution of energetic hot-spots like those of IgGs and conventional protein-protein complexes. The highly preorganized and energetically compact interface of D3-L11 recognizes the concave epitope with high shape complementarity by the classical lock-and-key mechanism. Our results shed light on the fundamental basis by which a particular VHH accommodate to the concave surface of an antigens with high affinity in a specific manner, enriching the mechanistic landscape of VHHs.

Nanoscale integration of single cell biologics discovery processes using optofluidic manipulation and monitoring

The new and rapid advancement in the complexity of biologics drug discovery has been driven by a deeper understanding of biological systems combined with innovative new therapeutic modalities, paving the way to breakthrough therapies for previously intractable diseases. These exciting times in biomedical innovation require the development of novel technologies to facilitate the sophisticated, multifaceted, high-paced workflows necessary to support modern large molecule drug discovery. A high-level aspiration is a true integration of "lab-on-a-chip" methods that vastly miniaturize cellulmical experiments could transform the speed, cost, and success of multiple workstreams in biologics development. Several microscale bioprocess technologies have been established that incrementally address these needs, yet each is inflexibly designed for a very specific process thus limiting an integrated holistic application. A more fully integrated nanoscale approach that incorporates manipulation, culture, analytics, and traceable digital record keeping of thousands of single cells in a relevant nanoenvironment would be a transformative technology capable of keeping pace with today's rapid and complex drug discovery demands. The recent advent of optical manipulation of cells using light-induced electrokinetics with micro- and nanoscale cell culture is poised to revolutionize both fundamental and applied biological research. In this review, we summarize the current state of the art for optical manipulation techniques and discuss emerging biological applications of this technology. In particular, we focus on promising prospects for drug discovery workflows, including antibody discovery, bioassay development, antibody engineering, and cell line development, which are enabled by the automation and industrialization of an integrated optoelectronic single-cell manipulation and culture platform. Continued development of such platforms will be well positioned to overcome many of the challenges currently associated with fragmented, low-throughput bioprocess workflows in biopharma and life science research.

Homology Modeling-Based in Silico Affinity Maturation Improves the Affinity of a Nanobody

Affinity maturation and rational design have a raised importance in the application of nanobody (VHH), and its unique structure guaranteed these processes quickly done in vitro. An anti-CD47 nanobody, Nb02, was screened via a synthetic phage display library with 278 nM of K_D value. In this study, a new strategy based on homology modeling and Rational Mutation Hotspots Design Protocol (RMHDP) was presented for building a fast and efficient platform for nanobody affinity maturation. A three-dimensional analytical structural model of Nb02 was constructed and then docked with the antigen, the CD47 extracellular domain (CD47^ext). Mutants with high binding affinity are predicted by the scoring of nanobody-antigen complexes based on molecular dynamics trajectories and simulation. Ultimately, an improved mutant with an 87.4-fold affinity (3.2 nM) and 7.36 °C higher thermal stability was obtained. These findings might contribute to computational affinity maturation of nanobodies via homology modeling using the recent advancements in computational power. The add-in of aromatic residues which formed aromatic-aromatic interaction plays a pivotal role in affinity and thermostability improvement. In a word, the methods used in this study might provide a reference for rapid and efficient in vitro affinity maturation of nanobodies.

Biophysical characterization and single-chain Fv construction of a neutralizing antibody to measles virus

The measles virus (MV) is a major cause of childhood morbidity and mortality worldwide. We previously established a mouse monoclonal antibody, 2F4, which shows high neutralizing titers against eight different genotypes of MV. However, the molecular basis for the neutralizing activity of the 2F4 antibody remains incompletely understood. Here, we have evaluated the binding characteristics of a Fab fragment of the 2F4 antibody. Using the MV infectious assay, we demonstrated that 2F4 Fab inhibits viral entry via either of two cellular receptors, SLAM and Nectin4. Surface plasmon resonance (SPR) analysis of recombinant proteins indicated that 2F4 Fab interacts with MV hemagglutinin (MV-H) with a K_D value at the nm level. Furthermore, we designed a single-chain Fv fragment of 2F4 antibody as another potential biopharmaceutical to target measles. The stable 2F4 scFv was successfully prepared by the refolding method and shown to interact with MV-H at the μm level. Like 2F4 Fab, scFv inhibited receptor binding and viral entry. This indicates that 2F4 mAb uses the receptor-binding site and/or a neighboring region as an epitope with high affinity. These results provide insight into the neutralizing activity and potential therapeutic use of antibody fragments for MV infection.

Trends in therapeutic antibody affinity maturation: From in-vitro towards next-generation sequencing approaches

Current advances in antibody engineering driving the strongest growth area in biotherapeutic agents development. Affinity improvement that is mainly important for biological activity and clinical efficacy of therapeutic antibodies, has still remained a challenging task. In the human body, during a course of immune response affinity maturation increase antibody activity by several rounds of somatic hypermutation and clonal selection in the germinal center. The final outputs are antibodies representing higher affinity and specificity against a particular antigen. In the realm of biotechnology, exploring of mutations which improve antibody affinity while preserving its specificity and stability is an extremely time-consuming and laborious process. Recent advances in computational algorithms and DNA sequencing technologies help researchers to redesign antibody structure to achieve desired properties such as improved binding affinity. In this review, we briefly described the principle of affinity maturation and different corresponding in vitro techniques. Also, we recapitulated the most recent advancements in the field of antibody affinity maturation including computational approaches and next-generation sequencing (NGS).

Experimentally guided computational antibody affinity maturation with de novo docking, modelling and rational design

Antibodies are an important class of therapeutics that have significant clinical impact for the treatment of severe diseases. Computational tools to support antibody drug discovery have been developing at an increasing rate over the last decade and typically rely upon a predetermined co-crystal structure of the antibody bound to the antigen for structural predictions. Here, we show an example of successful in silico affinity maturation of a hybridoma derived antibody, AB1, using just a homology model of the antibody fragment variable region and a protein-protein docking model of the AB1 antibody bound to the antigen, murine CCL20 (muCCL20). In silico affinity maturation, together with alanine scanning, has allowed us to fine-tune the protein-protein docking model to subsequently enable the identification of two single-point mutations that increase the affinity of AB1 for muCCL20. To our knowledge, this is one of the first examples of the use of homology modelling and protein docking for affinity maturation and represents an approach that can be widely deployed.

Structure-guided design fine-tunes pharmacokinetics, tolerability, and antitumor profile of multispecific frizzled antibodies

Aberrant activation of Wnt/β-catenin signaling occurs frequently in cancer. However, therapeutic targeting of this pathway is complicated by the role of Wnt in stem cell maintenance and tissue homeostasis. Here, we evaluated antibodies blocking 6 of the 10 human Wnt/Frizzled (FZD) receptors as potential therapeutics. Crystal structures revealed a common binding site for these monoclonal antibodies (mAbs) on FZD, blocking the interaction with the Wnt palmitoleic acid moiety. However, these mAbs displayed gastrointestinal toxicity or poor plasma exposure in vivo. Structure-guided engineering was used to refine the binding of each mAb for FZD receptors, resulting in antibody variants with improved in vivo tolerability and developability. Importantly, the lead variant mAb significantly inhibited tumor growth in the HPAF-II pancreatic tumor xenograft model. Taken together, our data demonstrate that anti-FZD cancer therapeutic antibodies with broad specificity can be fine-tuned to navigate in vivo exposure and tolerability while driving therapeutic efficacy.

Exploring designability of electrostatic complementarity at an antigen-antibody interface directed by mutagenesis, biophysical analysis, and molecular dynamics simulations

Antibodies protect organisms from a huge variety of foreign antigens. Antibody diversity originates from both genetic and structural levels. Antigen recognition relies on complementarity between antigen-antibody interfaces. Recent methodological advances in structural biology and the accompanying rapid increase of the number of crystal structures of proteins have enabled atomic-level manipulation of protein structures to effect alterations in function. In this study, we explored the designability of electrostatic complementarity at an antigen-antibody interface on the basis of a crystal structure of the complex. We designed several variants with altered charged residues at the interface and characterized the designed variants by surface plasmon resonance, circular dichroism, differential scanning calorimetry, and molecular dynamics simulations. Both successes and failures of the structure-based design are discussed. The variants that compensate electrostatic interactions can restore the interface complementarity, enabling the cognate antigen-antibody binding. Retrospectively, we also show that these mutational effects could be predicted by the simulations. Our study demonstrates the importance of charged residues on the physical properties of this antigen-antibody interaction and suggests that computational approaches can facilitate design of antibodies that recognize a weakly immunogenic antigen.

Optimisation of human VH domain antibodies specific to Mycobacterium tuberculosis heat shock protein (HSP16.3)

Mycobacterium tuberculosis (Mtb) 16.3 kDa heat shock protein 16.3 (HSP16.3) is a latency-associated antigen that can be targeted for latent tuberculosis (TB) diagnostic and therapeutic development. We have previously developed human V_H domain antibodies (dAbs; clone E3 and F1) specific against HSP16.3. In this work, we applied computational methods to optimise and design the antibodies in order to improve the binding affinity with HSP16.3. The V_H domain antibodies were first docked to the dimer form of HSP16.3 and further sampled using molecular dynamics simulation. The calculated binding free energy of the HSP16.3-dAb complexes showed non-polar interactions were responsible for the antigen-antibody association. Per-residue free energy decomposition and computational alanine scanning have identified one hotspot residue for E3 (Y391) and 4 hotspot residues for F1 (M394, Y396, R397 and M398). These hotspot residues were then mutated and evaluated by binding free energy calculations. Phage ELISA assay was carried out on the potential mutants (E3_Y391W, F1_M394E, F1_R397N and F1_M398Y). The experimental assay showed improved binding affinities of E3_Y391W and F1_M394E against HSP16.3 compared with the wild type E3 and F1. This case study has thus showed in silico methods are able to assist in optimisation or improvement of antibody-antigen binding.

In Silico Methods in Antibody Design

Antibody therapies with high efficiency and low toxicity are becoming one of the major approaches in antibody therapeutics. Based on high-throughput sequencing and increasing experimental structures of antibodies/antibody-antigen complexes, computational approaches can predict antibody/antigen structures, engineering the function of antibodies and design antibody-antigen complexes with improved properties. This review summarizes recent progress in the field of in silico design of antibodies, including antibody structure modeling, antibody-antigen complex prediction, antibody stability evaluation, and allosteric effects in antibodies and functions. We listed the cases in which these methods have helped experimental studies to improve the affinities and physicochemical properties of antibodies. We emphasized how the molecular dynamics unveiled the allosteric effects during antibody-antigen recognition and antibody-effector recognition.