Hiding in plain sight: structure and sequence analysis reveals the importance of the antibody DE loop for antibody-antigen binding

Prospects for the computational humanization of antibodies and nanobodies

To be viable therapeutics, antibodies must be tolerated by the human immune system. Rational approaches to reduce the risk of unwanted immunogenicity involve maximizing the 'humanness' of the candidate drug. However, despite the emergence of new discovery technologies, many of which start from entirely human gene fragments, most antibody therapeutics continue to be derived from non-human sources with concomitant humanization to increase their human compatibility. Early experimental humanization strategies that focus on CDR loop grafting onto human frameworks have been critical to the dominance of this discovery route but do not consider the context of each antibody sequence, impacting their success rate. Other challenges include the simultaneous optimization of other drug-like properties alongside humanness and the humanization of fundamentally non-human modalities such as nanobodies. Significant efforts have been made to develop in silico methodologies able to address these issues, most recently incorporating machine learning techniques. Here, we outline these recent advancements in antibody and nanobody humanization, focusing on computational strategies that make use of the increasing volume of sequence and structural data available and the validation of these tools. We highlight that structural distinctions between antibodies and nanobodies make the application of antibody-focused in silico tools to nanobody humanization non-trivial. Furthermore, we discuss the effects of humanizing mutations on other essential drug-like properties such as binding affinity and developability, and methods that aim to tackle this multi-parameter optimization problem.

Unveiling the multifaceted landscape of N-glycosylation in antibody variable domains: Insights and implications

N-glycosylation at the antibody variable domain has emerged as an important modification influencing antibody function. Despite its significance, information regarding its role and regulation remains limited. To address this gap, we comprehensively explored antibody structures housing N-glycosylation within the Protein Data Bank, yielding fresh insights into this intricate landscape. Our findings revealed that among 208 structures, N-glycosylation was more prevalent in human and mouse antibodies containing IGHV1-8 and IGHV2-2 germline genes, respectively. Moreover, our research highlights the potential for somatic hypermutation to introduce N-glycosylation sites by substituting polar residues (Ser or Thr) in germline variable genes with asparagine. Notably, our study underscores the prevalence of N-glycosylation in antiviral antibodies, especially anti-HIV. Besides antigen-antibody interaction, our findings suggest that N-glycosylation may impact antibody specificity, affinity, and avidity by influencing Fab dimer formation and complementary-determining region orientation. We also identified different glycan structures in HIV and SARS-CoV-2 antibody proteomic datasets, highlighting disparities from the N-glycan structures between PDB antibodies and biological repertoires further highlighting the complexity of N-glycosylation patterns. Our findings significantly enrich our understanding of the N-glycosylation's multifaceted characteristics within the antibody variable domain. Additionally, they underscore the pressing imperative for a more comprehensive characterization of its impact on antibody function.

HLA3DB: comprehensive annotation of peptide/HLA complexes enables blind structure prediction of T cell epitopes

The class I proteins of the major histocompatibility complex (MHC-I) display epitopic peptides derived from endogenous proteins on the cell surface for immune surveillance. Accurate modeling of peptides bound to the human MHC, HLA, has been mired by conformational diversity of the central peptide residues, which are critical for recognition by T cell receptors. Here, analysis of X-ray crystal structures within our curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer pHLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work may be applied towards predicting antigen immunogenicity, and receptor cross-reactivity.

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.

A comparison of the binding sites of antibodies and single-domain antibodies

Antibodies are the largest class of biotherapeutics. However, in recent years, single-domain antibodies have gained traction due to their smaller size and comparable binding affinity. Antibodies (Abs) and single-domain antibodies (sdAbs) differ in the structures of their binding sites: most significantly, single-domain antibodies lack a light chain and so have just three CDR loops. Given this inherent structural difference, it is important to understand whether Abs and sdAbs are distinguishable in how they engage a binding partner and thus, whether they are suited to different types of epitopes. In this study, we use non-redundant sequence and structural datasets to compare the paratopes, epitopes and antigen interactions of Abs and sdAbs. We demonstrate that even though sdAbs have smaller paratopes, they target epitopes of equal size to those targeted by Abs. To achieve this, the paratopes of sdAbs contribute more interactions per residue than the paratopes of Abs. Additionally, we find that conserved framework residues are of increased importance in the paratopes of sdAbs, suggesting that they include non-specific interactions to achieve comparable affinity. Furthermore, the epitopes of sdAbs are only marginally less accessible than those of Abs: we posit that this may be explained by differences in the orientation and compaction of sdAb and Ab CDR-H3 loops. Overall, our results have important implications for the engineering and humanization of sdAbs, as well as the selection of the best modality for targeting a particular epitope.

General Trends of the Camelidae Antibody VHHs Domain Dynamics

Conformational flexibility plays an essential role in antibodies' functional and structural stability. They facilitate and determine the strength of antigen-antibody interactions. Camelidae express an interesting subtype of single-chain antibody, named Heavy Chain only Antibody. They have only one N-terminal Variable domain (V_HH) per chain, composed of Frameworks (FRs) and Complementarity Determining regions (CDRs) like their VH and VL counterparts in IgG. Even when expressed independently, V_HH domains display excellent solubility and (thermo)stability, which helps them to retain their impressive interaction capabilities. Sequence and structural features of V_HH domains contributing to these abilities have already been studied compared to classical antibodies. To have the broadest view and understand the changes in dynamics of these macromolecules, large-scale molecular dynamics simulations for a large number of non-redundant V_HH structures have been performed for the first time. This analysis reveals the most prevalent movements in these domains. It reveals the four main classes of V_HHs dynamics. Diverse local changes were observed in CDRs with various intensities. Similarly, different types of constraints were observed in CDRs, while FRs close to CDRs were sometimes primarily impacted. This study sheds light on the changes in flexibility in different regions of V_HH that may impact their in silico design.

Structural Characterization of Nanobodies during Germline Maturation

Camelid heavy-chain antibody variable domains (VHH), nanobodies, are the smallest-known functional antibody fragments with high therapeutic potential. In this study, we investigate a VHH binding to hen egg-white lysozyme (HEL). We structurally and dynamically characterized the conformational diversity of four VHH variants to elucidate the antigen-binding process. For two of these antibodies, not only are the dissociation constants known, but also the experimentally determined crystal structures of the VHH in complex with HEL are available. We performed well-tempered metadynamics simulations in combination with molecular dynamics simulations to capture a broad conformational space and to reconstruct the thermodynamics and kinetics of conformational transitions in the antigen-binding site, the paratope. By kinetically characterizing the loop movements of the paratope, we found that, with an increase in affinity, the state populations shift towards the binding competent conformation. The contacts contributing to antigen binding, and those who contribute to the overall stability, show a clear trend towards less variable but more intense contacts. Additionally, these investigated nanobodies clearly follow the conformational selection paradigm, as the binding competent conformation pre-exists within the structural ensembles without the presence of the antigen.

Structural Classification of CDR-H3 in Single-Domain VHH Antibodies

The immune systems protect vertebrates from foreign molecules or antigens, and antibodies are important mediators of this system. The sequences and structural features of antibodies vary depending on species. Many of antibodies from vertebrates, including camelids, have both heavy and light chain variable domains, but camelids also have antibodies that lack the light chains. In antibodies that lack light chains, the C-terminal variable region is called the V_HH domain. Antibodies recognize antigens through six complementarity-determining regions (CDRs). The third CDR of the heavy chain (CDR-H3) is at the center of the antigen-binding site and is diverse in terms of sequence and structure. Due to the importance of antibodies in basic science as well as in medical applications, there have been many studies of CDR-H3s of antibodies that possess both light and heavy chains. However, nature of CDR-H3s of single-domain V_HH antibodies is less well studied. In this chapter, we describe current knowledge of sequence-structure-function correlations of single-domain V_HH antibodies with emphasis on CDR-H3. Based on the 370 crystal structures in the Protein Data Bank, we also attempt structural classification of CDR-H3 in single-domain V_HH antibodies and discuss lessons learned from the ever-increasing number of the structures.

Coevolved Canonical Loops Conformations of Single-Domain Antibodies: A Tale of Three Pockets Playing Musical Chairs

Single-domain antibodies (sdAbs) are a promising class of biotherapeutics with unique structural traits within their paratope region. The distribution of canonical conformations explored by their complementarity determining region (CDR) loops differs to some extent from conventional two-chain Fv fragments of monoclonal antibodies (mAbs). In this study, we explored in detail the canonical structures of sdAb CDR-H1 and CDR-H2 loops and compared those with mAbs from the IGHV3 and IGHV1 gene families. We surveyed the antibody structures catalogued in SAbDab and clustered the CDR canonical loops in Cartesian space. While most of the sdAb clusters were sub-populations of previously defined canonical Fv conformations of CDR-H1 and CDR-H2, our stricter clustering approach defined narrower clusters in sequence-space. Meticulous visual inspection of sub-populations allowed a clearer understanding of sequence-structure relationships. The packing densities within structural pockets contacted by CDR-H1 and CDR-H2 canonical conformations were analyzed on the premise that these pockets cannot be left vacant as they would leave exposed supportive hydrophobic residues. The fine resolution of the canonical clusters defined here revealed unique signatures within these pockets, including distinct structural complementarities between CDR-H1 and CDR-H2 canonical clusters, which could not be perceived with the previous coarser clusters. We highlight examples where a single residue change in CDR-H1 sequence is sufficient to induce a dramatic population shift in CDR-H2 conformation. This suggests that preferences in combining CDR-H1 and CDR-H2 emerged naturally during antibody evolution, leading to preferred sets of conserved amino acids at key positions in the framework as well as within the CDR loops. We outline a game of musical chairs that is necessary to maintain the integrity of the antibody structures that arose during evolution. Our study also provides refined CDR-H1 and CDR-H2 structural templates for sdAb homology modeling that could be leveraged for improved antibody design.

Conformational Entropy as a Potential Liability of Computationally Designed Antibodies

In silico antibody discovery is emerging as a viable alternative to traditional in vivo and in vitro approaches. Many challenges, however, remain open to enabling the properties of designed antibodies to match those produced by the immune system. A major question concerns the structural features of computer-designed complementarity determining regions (CDRs), including the role of conformational entropy in determining the stability and binding affinity of the designed antibodies. To address this problem, we used enhanced-sampling molecular dynamics simulations to compare the free energy landscapes of single-domain antibodies (sdAbs) designed using structure-based (DesAb-HSA-D3) and sequence-based approaches (DesAbO), with that of a nanobody derived from llama immunization (Nb10). Our results indicate that the CDR3 of DesAbO is more conformationally heterogeneous than those of both DesAb-HSA-D3 and Nb10, and the CDR3 of DesAb-HSA-D3 is slightly more dynamic than that of Nb10, which is the original scaffold used for the design of DesAb-HSA-D3. These differences underline the challenges in the rational design of antibodies by revealing the presence of conformational substates likely to have different binding properties and to generate a high entropic cost upon binding.

Delineating The RAS Conformational Landscape

Mutations in RAS isoforms (KRAS, NRAS, and HRAS) are among the most frequent oncogenic alterations in many cancers, making these proteins high priority therapeutic targets. Effectively targeting RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structural catalog of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. Here we present an expanded classification of RAS conformations based on analyses of the catalytic switch 1 (SW1) and switch 2 (SW2) loops. From 721 human KRAS, NRAS, and HRAS structures available in the Protein Data Bank (206 RAS-protein co-complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of Y32 in SW1 and Y71 in SW2. Clustering all well-modeled SW1 and SW2 loops using a density-based machine learning algorithm defined additional conformational subsets, some previously undescribed. Three SW1 conformations and nine SW2 conformations were identified, each associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation could be further subdivided based on the hydrogen bond type made between Y32 and the GTP γ-phosphate. Further analysis clarified the catalytic impact of G12D and G12V mutations and the inhibitor chemistries that bind to each druggable RAS conformation. Overall, this study has expanded our understanding of RAS structural biology, which could facilitate future RAS drug discovery.

Germline-Dependent Antibody Paratope States and Pairing Specific VH-VL Interface Dynamics

Antibodies have emerged as one of the fastest growing classes of biotherapeutic proteins. To improve the rational design of antibodies, we investigate the conformational diversity of 16 different germline combinations, which are composed of 4 different kappa light chains paired with 4 different heavy chains. In this study, we systematically show that different heavy and light chain pairings strongly influence the paratope, interdomain interaction patterns and the relative V_H-V_L interface orientations. We observe changes in conformational diversity and substantial population shifts of the complementarity determining region (CDR) loops, resulting in distinct dominant solution structures and differently favored canonical structures. Additionally, we identify conformational changes in the structural diversity of the CDR-H3 loop upon different heavy and light chain pairings, as well as upon changes in sequence and structure of the neighboring CDR loops, despite having an identical CDR-H3 loop amino acid sequence. These results can also be transferred to all CDR loops and to the relative V_H-V_L orientation, as certain paratope states favor distinct interface angle distributions. Furthermore, we directly compare the timescales of sidechain rearrangements with the well-described transition kinetics of conformational changes in the backbone of the CDR loops. We show that sidechain flexibilities are strongly affected by distinct heavy and light chain pairings and decipher germline-specific structural features co-determining stability. These findings reveal that all CDR loops are strongly correlated and that distinct heavy and light chain pairings can result in different paratope states in solution, defined by a characteristic combination of CDR loop conformations and V_H-V_L interface orientations. Thus, these results have broad implications in the field of antibody engineering, as they clearly show the importance of considering paired heavy and light chains to understand the antibody binding site, which is one of the key aspects in the design of therapeutics.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

Affinity maturation: highlights in the application of in vitro strategies for the directed evolution of antibodies

Affinity maturation is a key technique in protein engineering which is used to improve affinity and binding interactions in vitro, a process often required to fulfil the therapeutic potential of antibodies. There are many available display technologies and maturation methods developed over the years, which have been instrumental in the production of therapeutic antibodies. However, due to the inherent limitations in display capacity of these technologies, accommodation of expansive and complex library builds is still a challenge. In this article, we discuss our recent efforts in the affinity maturation of a difficult antibody lineage using an unbiased approach, which sought to explore a larger sequence space through the application of DNA recombination and shuffling techniques across the entire antibody region and selections using ribosome display. We also highlight the key features of several display technologies and diversification methods, and discuss the strategies devised by different groups in response to different challenges. Particular attention is drawn to examples which are aimed at the expansion of sequence, structural or experimental diversity through different means and approaches. Here, we provide our perspectives on these methodologies and the considerations involved in the design of effective strategies for the directed evolution of antibodies.

Mutation of Framework Residue H71 Results in Different Antibody Paratope States in Solution

Characterizing and understanding the antibody binding interface have become a pre-requisite for rational antibody design and engineering. The antigen-binding site is formed by six hypervariable loops, known as the complementarity determining regions (CDRs) and by the relative interdomain orientation (V_H-V_L). Antibody CDR loops with a certain sequence have been thought to be limited to a single static canonical conformation determining their binding properties. However, it has been shown that antibodies exist as ensembles of multiple paratope states, which are defined by a characteristic combination of CDR loop conformations and interdomain orientations. In this study, we thermodynamically and kinetically characterize the prominent role of residue 71^H (Chothia nomenclature), which does not only codetermine the canonical conformation of the CDR-H2 loop but also results in changes in conformational diversity and population shifts of the CDR-H1 and CDR-H3 loop. As all CDR loop movements are correlated, conformational rearrangements of the heavy chain CDR loops also induce conformational changes in the CDR-L1, CDR-L2, and CDR-L3 loop. These overall conformational changes of the CDR loops also influence the interface angle distributions, consequentially leading to different paratope states in solution. Thus, the type of residue of 71^H, either an alanine or an arginine, not only influences the CDR-H2 loop ensembles, but co-determines the paratope states in solution. Characterization of the functional consequences of mutations of residue 71^H on the paratope states and interface orientations has broad implications in the field of antibody engineering.

Anti-TNF Alpha Antibody Humira with pH-dependent Binding Characteristics: A constant-pH Molecular Dynamics, Gaussian Accelerated Molecular Dynamics, and In Vitro Study

Humira is a monoclonal antibody that binds to TNF alpha, inactivates TNF alpha receptors, and inhibits inflammation. Neonatal Fc receptors can mediate the transcytosis of Humira-TNF alpha complex structures and process them toward degradation pathways, which reduces the therapeutic effect of Humira. Allowing the Humira-TNF alpha complex structures to dissociate to Humira and soluble TNF alpha in the early endosome to enable Humira recycling is crucial. We used the cytoplasmic pH (7.4), the early endosomal pH (6.0), and pKa of histidine side chains (6.0-6.4) to mutate the residues of complementarity-determining regions with histidine. Our engineered Humira (W1-Humira) can bind to TNF alpha in plasma at neutral pH and dissociate from the TNF alpha in the endosome at acidic pH. We used the constant-pH molecular dynamics, Gaussian accelerated molecular dynamics, two-dimensional potential mean force profiles, and in vitro methods to investigate the characteristics of W1-Humira. Our results revealed that the proposed Humira can bind TNF alpha with pH-dependent affinity in vitro. The W1-Humira was weaker than wild-type Humira at neutral pH in vitro, and our prediction results were close to the in vitro results. Furthermore, our approach displayed a high accuracy in antibody pH-dependent binding characteristics prediction, which may facilitate antibody drug design. Advancements in computational methods and computing power may further aid in addressing the challenges in antibody drug design.