Antibody-Antigen Binding Interface Analysis in the Big Data Era

Structural trends in antibody-antigen binding interfaces: a computational analysis of 1833 experimentally determined 3D structures

Antibodies are attractive therapeutic candidates due to their ability to bind cognate antigens with high affinity and specificity. Still, the underlying molecular rules governing the antibody-antigen interface remain poorly understood, making in silico antibody design inherently difficult and keeping the discovery and design of novel antibodies a costly and laborious process. This study investigates the characteristics of antibody-antigen binding interfaces through a computational analysis of more than 850,000 atom-atom contacts from the largest reported set of antibody-antigen complexes with 1833 nonredundant, experimentally determined structures. The analysis compares binding characteristics of conventional antibodies and single-domain antibodies (sdAbs) targeting both protein- and peptide antigens. We find clear patterns in the number antibody-antigen contacts and amino acid frequencies in the paratope. The direct comparison of sdAbs and conventional antibodies helps elucidate the mechanisms employed by sdAbs to compensate for their smaller size and the fact that they harbor only half the number of complementarity-determining regions compared to conventional antibodies. Furthermore, we pinpoint antibody interface hotspot residues that are often found at the binding interface and the amino acid frequencies at these positions. These findings have direct potential applications in antibody engineering and the design of improved antibody libraries.

Understanding the Specific Implications of Amino Acids in the Antibody Development

As the demand for immunotherapy to treat and manage cancers, infectious diseases and other disorders grows, a comprehensive understanding of amino acids and their intricate role in antibody engineering has become a prime requirement. Naturally produced antibodies may not have the most suitable amino acids at the complementarity determining regions (CDR) and framework regions, for therapeutic purposes. Therefore, to enhance the binding affinity and therapeutic properties of an antibody, the specific impact of certain amino acids on the antibody's architecture must be thoroughly studied. In antibody engineering, it is crucial to identify the key amino acid residues that significantly contribute to improving antibody properties. Therapeutic antibodies with higher binding affinity and improved functionality can be achieved through modifications or substitutions with highly suitable amino acid residues. Here, we have indicated the frequency of amino acids and their association with the binding free energy in CDRs. The review also analyzes the experimental outcome of two studies that reveal the frequency of amino acids in CDRs and provides their significant correlation between the outcomes. Additionally, it discusses the various bond interactions within the antibody structure and antigen binding. A detailed understanding of these amino acid properties should assist in the analysis of antibody sequences and structures needed for designing and enhancing the overall performance of therapeutic antibodies.

A comprehensive immunoinformatics study to explore and characterize potential vaccine constructs against Ole e 9 allergen of Olea europaea

Immunoglobulin E (IgE)-mediated allergy, which affects more than 30% of the population, is the most prevalent hypersensitivity illness. In an atopic individual, even a small amount of allergen exposure can cause IgE antibodies to be produced. Due to the engagement of receptors that are highly selective for IgE, even tiny amounts of allergens can induce massive inflammation. This study focuses on the exploration and characterization of the allergen potential of Olea europaea allergen (Ole e 9) affecting the population in Saudi Arabia. A systematic computational approach was performed to identify potential epitopes of allergens and complementary determining regions of IgE. In support, physiochemical characterization and secondary structure analysis unravel the structural conformations of allergens and active sites. Epitope prediction uses a pool of computational algorithms to identify plausible epitopes. Furthermore, the vaccine construct was assessed for its binding efficiency using molecular docking and molecular dynamics simulation studies, which led to strong and stable interactions. This is because IgE is known to play a role in allergic responses, which facilitate host cell activation for an immune response. Overall, the immunoinformatics analysis advocates that the proposed vaccine candidate is safe and immunogenic and therefore can be pushed as a lead for in vitro and in vivo investigations.Communicated by Ramaswamy H. Sarma.

Prediction of Paratope-Epitope Pairs Using Convolutional Neural Networks

Antibodies play a central role in the adaptive immune response of vertebrates through the specific recognition of exogenous or endogenous antigens. The rational design of antibodies has a wide range of biotechnological and medical applications, such as in disease diagnosis and treatment. However, there are currently no reliable methods for predicting the antibodies that recognize a specific antigen region (or epitope) and, conversely, epitopes that recognize the binding region of a given antibody (or paratope). To fill this gap, we developed ImaPEp, a machine learning-based tool for predicting the binding probability of paratope-epitope pairs, where the epitope and paratope patches were simplified into interacting two-dimensional patches, which were colored according to the values of selected features, and pixelated. The specific recognition of an epitope image by a paratope image was achieved by using a convolutional neural network-based model, which was trained on a set of two-dimensional paratope-epitope images derived from experimental structures of antibody-antigen complexes. Our method achieves good performances in terms of cross-validation with a balanced accuracy of 0.8. Finally, we showcase examples of application of ImaPep, including extensive screening of large libraries to identify paratope candidates that bind to a selected epitope, and rescoring and refining antibody-antigen docking poses.

A general approach for selection of epitope-directed binders to proteins

Directing antibodies to a particular epitope among many possible on a target protein is a significant challenge. Here, we present a simple and general method for epitope-directed selection (EDS) using a differential phage selection strategy. This involves engineering the protein of interest (POI) with the epitope of interest (EOI) mutated using a systematic bioinformatics algorithm to guide the local design of an EOI decoy variant. Using several alternating rounds of negative selection with the EOI decoy variant followed by positive selection on the wild-type POI, we were able to identify highly specific and potent antibodies to five different EOI antigens that bind and functionally block known sites of proteolysis. Among these, we developed highly specific antibodies that target the proteolytic site on the CUB domain containing protein 1 (CDCP1) to prevent its proteolysis allowing us to study the cellular maturation of this event that triggers malignancy. We generated antibodies that recognize the junction between the pro- and catalytic domains for three different matrix metalloproteases (MMPs), MMP1, MMP3, and MMP9, that selectively block activation of each of these enzymes and impair cell migration. We targeted a proteolytic epitope on the cell surface receptor, EPH Receptor A2 (EphA2), that is known to transform it from a tumor suppressor to an oncoprotein. We believe that the EDS method greatly facilitates the generation of antibodies to specific EOIs on a wide range of proteins and enzymes for broad therapeutic and diagnostic applications.

Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔG Predictions

The development of highly potent antibodies and antibody fragments as binding agents holds significant implications in fields such as biosensing and biotherapeutics. Their binding strength is intricately linked to the arrangement and composition of residues at the binding interface. Computational techniques offer a robust means to predict the three-dimensional structure of these complexes and to assess the affinity changes resulting from mutations. Given the interdependence of structure and affinity prediction, our objective here is to disentangle their roles. We aim to evaluate independently six side-chain reconstruction methods and ten binding affinity estimation techniques. This evaluation was pivotal in predicting affinity alterations due to single mutations, a key step in computational affinity maturation protocols. Our analysis focuses on a data set comprising 27 distinct antibody/hen egg white lysozyme complexes, each with crystal structures and experimentally determined binding affinities. Using six different side-chain reconstruction methods, we transformed each structure into its corresponding mutant via in silico single-point mutations. Subsequently, these structures undergo minimization and molecular dynamics simulation. We therefore estimate ΔΔG values based on the original crystal structure, its energy-minimized form, and the ensuing molecular dynamics trajectories. Our research underscores the critical importance of selecting reliable side-chain reconstruction methods and conducting thorough molecular dynamics simulations to accurately predict the impact of mutations. In summary, our study demonstrates that the integration of conformational sampling and scoring is a potent approach to precisely characterizing mutation processes in single-point mutagenesis protocols and crucial for computational antibody design.

Modelling the assembly and flexibility of antibody structures

Antibodies are large protein assemblies capable of both specifically recognising antigens and engaging with other proteins and receptors to coordinate immune action. Traditionally, structural studies have been dedicated to antibody variable regions, but efforts to determine and model full-length antibody structures are emerging. Here we review the current knowledge on modelling the structures of antibody assemblies, focusing on their conformational flexibility and the challenge this poses to obtaining and evaluating structural models. Integrative modelling approaches, combining experiments (cryo-electron microscopy, mass spectrometry, etc.) and computational methods (molecular dynamics simulations, deep-learning based approaches, etc.), hold the promise to map the complex conformational landscape of full-length antibody structures.

Predicting the immune escape of SARS-CoV-2 neutralizing antibodies upon mutation

COVID-19 has resulted in millions of deaths and severe impact on economies worldwide. Moreover, the emergence of SARS-CoV-2 variants presented significant challenges in controlling the pandemic, particularly their potential to avoid the immune system and evade vaccine immunity. This has led to a growing need for research to predict how mutations in SARS-CoV-2 reduces the ability of antibodies to neutralize the virus. In this study, we assembled a set of 1813 mutations from the interface of SARS-CoV-2 spike protein's receptor binding domain (RBD) and neutralizing antibody complexes and developed a machine learning model to classify high or low escape mutations using interaction energy, inter-residue contacts and predicted binding free energy change. Our approach achieved an Area under the Receiver Operating Characteristics (ROC) Curve (AUC) of 0.91 using the Random Forest classifier on the test dataset with 217 mutations. The model was further utilized to predict the escape mutations on a dataset of 29,165 mutations located at the interface of 83 RBD-neutralizing antibody complexes. A small subset of this dataset was also validated based on available experimental data. We found that top 10 % high escape mutations were dominated by charged to nonpolar mutations whereas low escape mutations were dominated by polar to nonpolar mutations. We believe that the present method will allow prioritization of high/low escape mutations in the context of neutralizing antibodies targeting SARS-CoV-2 RBD region and assist antibody design for current and emerging variants.

Structural analysis of light chain-driven bispecific antibodies targeting CD47 and PD-L1

In contrast to natural antibodies that rely mainly on the heavy chain to establish contacts with their cognate antigen, we have developed a bispecific antibody format in which the light chain (LC) drives antigen binding and specificity. To better understand epitope-paratope interactions in this context, we determined the X-ray crystallographic structures of an antigen binding fragment (Fab) in complex with human CD47 and another Fab in complex with human PD-L1. These Fabs contain a κ-LC and a λ-LC, respectively, which are paired with an identical heavy chain (HC). The structural analysis of these complexes revealed the dominant contribution of the LCs to antigen binding, but also that the common HC provides some contacts in both CD47 and PD-L1 Fab complexes. The anti-CD47 Fab was affinity optimized by diversifying complementary-determining regions of the LC followed by phage display selections. Using homology modeling, the contributions of the amino acid modification to the affinity increase were analyzed. Our results demonstrate that, despite a less prominent role in natural antibodies, the LC can mediate high affinity binding to different antigens and neutralize their biological function. Importantly, Fabs containing a common variable heavy (VH) domain enable the generation of bispecific antibodies retaining a truly native structure, maximizing their therapeutic potential.

SiteFerret: Beyond Simple Pocket Identification in Proteins

We present a novel method for the automatic detection of pockets on protein molecular surfaces. The algorithm is based on an ad hoc hierarchical clustering of virtual probe spheres obtained from the geometrical primitives used by the NanoShaper software to build the solvent-excluded molecular surface. The final ranking of putative pockets is based on the Isolation Forest method, an unsupervised learning approach originally developed for anomaly detection. A detailed importance analysis of pocket features provides insight into which geometrical (clustering) and chemical (amino acidic composition) properties characterize a good binding site. The method also provides a segmentation of pockets into smaller subpockets. We prove that subpockets are a convenient representation to pinpoint the binding site with great precision. SiteFerret is outstanding in its versatility, accurately predicting a wide range of binding sites, from those binding small molecules to those binding peptides, including difficult shallow sites.