Variable Regions of Antibodies and T-Cell Receptors May Not Be Sufficient in Molecular Simulations Investigating Binding

Immunoglobulin constant regions provide stabilization to the paratope and enforce epitope specificity

Constant domains in antibody molecules at the level of the Fab (CH1 and CL) have long been considered to be simple scaffolding elements that physically separate the paratope-defining variable (V) region from the effector function-mediating constant (C) regions. However, due to recent findings that C domains of different isotypes can modulate the fine specificity encoded in the V region, elucidating the role of C domains in shaping the paratope and influencing specificity is a critical area of interest. To dissect the relative contributions of each C domain to this phenomenon, we generated antibody fragments with different C regions omitted, using a set of antibodies targeting capsular polysaccharides from the fungal pathogen, Cryptococcus neoformans. Antigen specificity mapping and functional activity measurements revealed that V region-only antibody fragments exhibited poly-specificity to antigenic variants and extended to recognition of self-antigens, while measurable hydrolytic activity of the capsule was greatly attenuated. To better understand the mechanistic origins of the remarkable loss of specificity that accompanies the removal of C domains from identical paratopes, we performed molecular dynamics simulations which revealed increased paratope plasticity in the scFv relative to the corresponding Fab. Together, our results provide insight into how the remarkable specificity of immunoglobulins is governed and maintained at the level of the Fab through the enforcement of structural restrictions on the paratope by CH1 domains.

Contextualising the developability risk of antibodies with lambda light chains using enhanced therapeutic antibody profiling

Antibodies with lambda light chains (λ-antibodies) are generally considered to be less developable than those with kappa light chains (κ-antibodies). Though this hypothesis has not been formally established, it has led to substantial systematic biases in drug discovery pipelines and thus contributed to kappa dominance amongst clinical-stage therapeutics. However, the identification of increasing numbers of epitopes preferentially engaged by λ-antibodies shows there is a functional cost to neglecting to consider them as potential lead candidates. Here, we update our Therapeutic Antibody Profiler (TAP) tool to use the latest data and machine learning-based structure prediction, and apply it to evaluate developability risk profiles for κ-antibodies and λ-antibodies based on their surface physicochemical properties. We find that while human λ-antibodies on average have a higher risk of developability issues than κ-antibodies, a sizeable proportion are assigned lower-risk profiles by TAP and should represent more tractable candidates for therapeutic development. Through a comparative analysis of the low- and high-risk populations, we highlight opportunities for strategic design that TAP suggests would enrich for more developable λ-antibodies. Overall, we provide context to the differing developability of κ- and λ-antibodies, enabling a rational approach to incorporate more diversity into the initial pool of immunotherapeutic candidates.

Fast, Accurate, and System-Specific Variable-Resolution Modeling of Proteins

In recent years, a few multiple-resolution modeling strategies have been proposed, in which functionally relevant parts of a biomolecule are described with atomistic resolution, with the remainder of the system being concurrently treated using a coarse-grained model. In most cases, the parametrization of the latter requires lengthy reference all-atom simulations and/or the usage of off-shelf coarse-grained force fields, whose interactions have to be refined to fit the specific system under examination. Here, we overcome these limitations through a novel multiresolution modeling scheme for proteins, dubbed coarse-grained anisotropic network model for variable resolution simulations, or CANVAS. This scheme enables a user-defined modulation of the resolution level throughout the system structure; a fast parametrization of the potential without the necessity of reference simulations; and the straightforward usage of the model on the most commonly used molecular dynamics platforms. The method is presented and validated with two case studies, the enzyme adenylate kinase and the therapeutic antibody pembrolizumab, by comparing the results obtained with the CANVAS model against fully atomistic simulations. The modeling software, implemented in Python, is made freely available for the community on a collaborative github repository.

Embedding Dynamics in Intrinsic Physicochemical Profiles of Market-Stage Antibody-Based Biotherapeutics

Adequate stability, manufacturability, and safety are crucial to bringing an antibody-based biotherapeutic to the market. Following the concept of holistic in silico developability, we introduce a physicochemical description of 91 market-stage antibody-based biotherapeutics based on orthogonal molecular properties of variable regions (Fvs) embedded in different simulation environments, mimicking conditions experienced by antibodies during manufacturing, formulation, and in vivo. In this work, the evaluation of molecular properties includes conformational flexibility of the Fvs using molecular dynamics (MD) simulations. The comparison between static homology models and simulations shows that MD significantly affects certain molecular descriptors like surface molecular patches. Moreover, the structural stability of a subset of Fv regions is linked to changes in their specific molecular interactions with ions in different experimental conditions. This is supported by the observation of differences in protein melting temperatures upon addition of NaCl. A DEvelopability Navigator In Silico (DENIS) is proposed to compare mAb candidates for their similarity with market-stage biotherapeutics in terms of physicochemical properties and conformational stability. Expanding on our previous developability guidelines (Ahmed et al. Proc. Natl. Acad. Sci. 2021, 118 (37), e2020577118), the hydrodynamic radius and the protein strand ratio are introduced as two additional descriptors that enable a more comprehensive in silico characterization of biotherapeutic drug candidates. Test cases show how this approach can facilitate identification and optimization of intrinsically developable lead candidates. DENIS represents an advanced computational tool to progress biotherapeutic drug candidates from discovery into early development by predicting drug properties in different aqueous environments.

Communication pathways bridge local and global conformations in an IgG4 antibody

The affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recent evidence, however, points toward a nontrivial relation between the CDR and distal sites: variations in the binding strengths have been observed upon mutating residues separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies, whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4 μs of MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR. These results support the hypothesis that pembrolizumab featues a multi-scale hierarchy of intertwined global and local conformational changes. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.

Application of phage display for T-cell receptor discovery

The immune system is tasked to keep our body unharmed and healthy. In the immune system, B- and T-lymphocytes are the two main components working together to stop and eliminate invading threats like virus particles, bacteria, fungi and parasite from attacking our healthy cells. The function of antibodies is relatively more direct in target recognition as compared to T-cell receptors (TCR) which recognizes antigenic peptides being presented on the major histocompatibility complex (MHC). Although phage display has been widely applied for antibody presentation, this is the opposite in the case of TCR. The cell surface TCR is a relatively large and complex molecule, making presentation on phage surfaces challenging. Even so, recombinant versions and modifications have been introduced to allow the growing development of TCR in phage display. In addition, the increasing application of TCR for immunotherapy has made it an important binding motif to be developed by phage display. This review will emphasize on the application of phage display for TCR discovery as well as the engineering aspect of TCR for improved characteristics.

Unraveling the Allosteric Communication Mechanisms in T-Cell Receptor-Peptide-Loaded Major Histocompatibility Complex Dynamics Using Molecular Dynamics Simulations: An Approach Based on Dynamic Cross Correlation Maps and Residue Interaction Energy Calculations

Antigen presentation by major histocompatibility complex (MHC) proteins to T-cell receptors (TCRs) plays a crucial role in triggering the adaptive immune response. Most of our knowledge on TCR-peptide-loaded major histocompatibility complex (pMHC) interaction stemmed from experiments yielding static structures, yet the dynamic aspects of this molecular interaction are equally important to understand the underlying molecular mechanisms and to develop treatment strategies against diseases such as cancer and autoimmune diseases. To this end, computational biophysics studies including all-atom molecular dynamics simulations have provided useful insights; however, we still lack a basic understanding of an overall allosteric mechanism that results in conformational changes in the TCR and subsequent T-cell activation. Previous hydrogen-deuterium exchange and nuclear magnetic resonance studies provided clues regarding these molecular mechanisms, including global rigidification and allosteric effects on the constant domain of TCRs away from the pMHC interaction site. Here, we show that molecular dynamics simulations can be used to identify how this overall rigidification may be related to the allosteric communication within TCRs upon pMHC interaction via essential dynamics and nonbonded residue-residue interaction energy analyses. The residues taking part in the rigidification effect are highlighted with an intricate analysis on residue interaction changes, which lead to a detailed outline of the complex formation event. Our results indicate that residues of the Cβ domain of TCRs show significant differences in their nonbonded interactions upon complex formation. Moreover, the dynamic cross correlations between these residues are also increased, in line with their nonbonded interaction energy changes. Altogether, our approach may be valuable for elucidating intramolecular allosteric changes in the TCR structure upon pMHC interaction in molecular dynamics simulations.

Surprisingly Fast Interface and Elbow Angle Dynamics of Antigen-Binding Fragments

Fab consist of a heavy and light chain and can be subdivided into a variable (V H and V L ) and a constant region (C H 1 and C L ). The variable region contains the complementarity-determining region (CDR), which is formed by six hypervariable loops, shaping the antigen binding site, the paratope. Apart from the CDR loops, both the elbow angle and the relative interdomain orientations of the V H -V L and the C H 1-C L domains influence the shape of the paratope. Thus, characterization of the interface and elbow angle dynamics is essential to antigen specificity. We studied nine antigen-binding fragments (Fab) to investigate the influence of affinity maturation, antibody humanization, and different light-chain types on the interface and elbow angle dynamics. While the CDR loops reveal conformational transitions in the micro-to-millisecond timescale, both the interface and elbow angle dynamics occur on the low nanosecond timescale. Upon affinity maturation, we observe a substantial rigidification of the V H and V L interdomain and elbow-angle flexibility, reflected in a narrower and more distinct distribution. Antibody humanization describes the process of grafting non-human CDR loops onto a representative human framework. As the antibody framework changes upon humanization, we investigated if both the interface and the elbow angle distributions are changed or shifted. The results clearly showed a substantial shift in the relative V H -V L distributions upon antibody humanization, indicating that different frameworks favor distinct interface orientations. Additionally, the interface and elbow angle dynamics of five antibody fragments with different light-chain types are included, because of their strong differences in elbow angles. For these five examples, we clearly see a high variability and flexibility in both interface and elbow angle dynamics, highlighting the fact that Fab interface orientations and elbow angles interconvert between each other in the low nanosecond timescale. Understanding how the relative interdomain orientations and the elbow angle influence antigen specificity, affinity, and stability has broad implications in the field of antibody modeling and engineering.

VH -VL interdomain dynamics observed by computer simulations and NMR

The relative orientation of the two variable domains, VH and VL , influences the shape of the antigen binding site, that is, the paratope, and is essential to understand antigen specificity. ABangle characterizes the VH -VL orientation by using five angles and a distance and compares it to other known structures. Molecular dynamics simulations of antibody variable domains (Fvs) reveal fluctuations in the relative domain orientations. The observed dynamics between these domains are confirmed by NMR experiments on a single-chain variable fragment antibody (scFv) in complex with IL-1β and an antigen-binding fragment (Fab). The variability of these relative domain orientations can be interpreted as a structural feature of antibodies, which increases the antibody repertoire significantly and can enlarge the number of possible binding partners substantially. The movements of the VH and VL domains are well sampled with molecular dynamics simulations and are in agreement with the NMR ensemble. Fast Fourier transformation of the ABangle metrics allows to assign timescales of 0.1-10 GHz to the fastest collective interdomain movements. The results clearly show the necessity of dynamics to understand and characterize the favorable orientations of the VH and VL domains implying a considerable binding interface flexibility and reveal in all antibody fragments (Fab, scFv, and Fv) very similar VH -VL interdomain variations comparable to the distributions observed for known X-ray structures of antibodies. SIGNIFICANCE STATEMENT: Antibodies have become key players as therapeutic agents. The binding ability of antibodies is determined by the antigen-binding fragment (Fab), in particular the variable fragment region (Fv). Antigen-binding is mediated by the complementarity-determining regions consisting of six loops, each three of the heavy and light chain variable domain VH and VL . The relative orientation of the VH and VL domains influences the shape of the antigen-binding site and is a major objective in antibody design. In agreement with NMR experiments and molecular dynamics simulations, we show a considerable binding site flexibility in the low nanosecond timescale. Thus we suggest that this flexibility and its implications for binding and specificity should be considered when designing and optimizing therapeutic antibodies.

Exploring the different states of wild-type T-cell receptor and mutant conformational changes towards understanding the antigen recognition

Recognition of proteolytic peptide fragments presented by major histocompatibility complex (MHC) on target cells by T-cell receptor (TCR) is among the most important interactions in the adaptive immune system. Several computational studies have been performed to investigate conformational and dynamical properties of TCRs for enhanced immunogenicity. Here, we present the large-scale molecular dynamics (MD) simulation studies of the two comprehensive systems consisting of the wild-type and mutant IG4 TCR in complex with the tumor epitope NY-ESO peptide (SLLMWITQC) and analyzed for mapping conformational changes of TCR in the states prior to antigen binding, upon antigen binding and after the antigen was released. All of the simulations were performed with different states of TCRs for each 1000 ns of simulation time, providing six simulations for time duration of 6000 ns (6µs). We show that rather than undergoing most critical conformational changes upon antigen binding, the high proportion of complementarity-determining region (CDR) loops change by comparatively small amount. The hypervariable CDRα3 and CDRβ3 loops showed significant structural changes. Interestingly, the TCR β chain loops showed the least changes, which is reliable with recent implications that β domain of TCR may propel antigen interaction. The mutant shows higher rigidity than wild-type even in released state; expose an induced fit mechanism occurring from the re-structuring of CDRα3 loop and can allow enhanced binding affinity of the peptide antigen. Additionally, we show that CDRα3 loop and peptide contacts are an adaptive feature of affinity enhanced mutant TCR.Communicated by Ramaswamy H. Sarma.

MHC binding affects the dynamics of different T-cell receptors in different ways

T cells use their T-cell receptors (TCRs) to scan other cells for antigenic peptides presented by MHC molecules (pMHC). If a TCR encounters a pMHC, it can trigger a signalling pathway that could lead to the activation of the T cell and the initiation of an immune response. It is currently not clear how the binding of pMHC to the TCR initiates signalling within the T cell. One hypothesis is that conformational changes in the TCR lead to further downstream signalling. Here we investigate four different TCRs in their free state as well as in their pMHC bound state using large scale molecular simulations totalling 26 000 ns. We find that the dynamical features within TCRs differ significantly between unbound TCR and TCR/pMHC simulations. However, apart from expected results such as reduced solvent accessibility and flexibility of the interface residues, these features are not conserved among different TCR types. The presence of a pMHC alone is not sufficient to cause cross-TCR-conserved dynamical features within a TCR. Our results argue against models of TCR triggering involving conserved allosteric conformational changes.

The interaction between T-cells and other cells is one of the most important interactions in the human immune system. If T-cells are not triggered major parts of the immune system cannot be activated or are not working effectively. Despite many years of research the exact mechanism of how a T-cell is initially triggered is not clear. One hypothesis is that conformational changes within the T-cell receptor (TCR) can cause further downstream signalling within the T-cell. In this study we computationally investigate the dynamics of four different TCRs in their free and bound configuration. Our large scale simulations show that all four TCRs react to binding in different ways. In some TCRs mainly the areas close to the binding region are affected while in other TCRs areas further apart from the binding region are also affected. Our results argue against a conserved structural activation mechanism across different types of TCRs.

NMR: an essential structural tool for integrative studies of T cell development, pMHC ligand recognition and TCR mechanobiology

Early studies of T cell structural biology using X-ray crystallography, surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) focused on a picture of the αβT cell receptor (αβTCR) component domains and their cognate ligands (peptides bound to MHC molecules, i.e. pMHCs) as static interaction partners. Moving forward requires integrating this corpus of data with dynamic technologies such as NMR, molecular dynamics (MD) simulations and real-time single molecule (SM) studies exemplified by optical tweezers (OT). NMR bridges relevant timescales and provides the potential for an all-atom dynamic description of αβTCR components prior to and during interactions with binding partners. SM techniques have opened up vistas in understanding the non-equilibrium nature of T cell signaling through the introduction of force-mediated binding measurements into the paradigm for T cell function. In this regard, bioforces consequent to T-lineage cell motility are now perceived as placing piconewton (pN)-level loads on single receptor-pMHC bonds to impact structural change and αβT-lineage biology, including peptide discrimination, cellular activation, and developmental progression. We discuss herein essential NMR technologies in illuminating the role of ligand binding in the preT cell receptor (preTCR), the αβTCR developmental precursor, and convergence of NMR, SM and MD data in advancing our comprehension of T cell development. More broadly we review the central hypothesis that the αβTCR is a mechanosensor, fostered by breakthrough NMR-based structural insights. Collectively, elucidating dynamic aspects through the integrative use of NMR, SM, and MD shall advance fundamental appreciation of the mechanism of T cell signaling as well as inform translational efforts in αβTCR and chimeric T cell (CAR-T) immunotherapies and T cell vaccinology.

How peptide/MHC presence affects the dynamics of the LC13 T-cell receptor

The interaction between T-cell receptors (TCRs) of T-cells and potentially immunogenic peptides presented by MHCs of antigen presenting cells is one of the most important mechanisms of the adaptive human immune system. A large number of structural simulations of the TCR/peptide/MHC system have been carried out. However, to date no study has investigated the differences of the dynamics between free TCRs and pMHC bound TCRs on a large scale. Here we present a study totalling 37 100 ns investigating the LC13 TCR in its free form as well as in complex with HLA-B*08:01 and different peptides. Our results show that the dynamics of the bound and unbound LC13 TCR differ significantly. This is reflected in (a) expected results such as an increased flexibility and increased solvent accessible surface of the CDRs of unbound TCR simulations but also in (b) less expected results such as lower CDR distances and compactness as well as alteration in the hydrogen bond network around CDR3α of unbound TCR simulations. Our study further emphasises the structural flexibility of TCRs and confirms the importance of the CDR3 loops for the adoption to MHC.

Avoiding False Positive Conclusions in Molecular Simulation: The Importance of Replicas

Molecular simulations are a computational technique used to investigate the dynamics of proteins and other molecules. The free energy landscape of these simulations is often rugged, and minor differences in the initial velocities, floating-point precision, or underlying hardware can cause identical simulations (replicas) to take different paths in the landscape. In this study we investigated the magnitude of these effects based on 310 000 ns of simulation time. We performed 100 identically parametrized replicas of 3000 ns each for a small 10 amino acid system as well as 100 identically parametrized replicas of 100 ns each for an 827 residue T-cell receptor/MHC system. Comparing randomly chosen subgroups within these replica sets, we estimated the reproducibility and reliability that can be achieved by a given number of replicas at a given simulation time. These results demonstrate that conclusions drawn from single simulations are often not reproducible and that conclusions drawn from multiple shorter replicas are more reliable than those from a single longer simulation. The actual number of replicas needed will always depend on the question asked and the level of reliability sought. On the basis of our data, it appears that a good rule of thumb is to perform a minimum of five to 10 replicas.

Conformational changes in antibody Fab fragments upon binding and their consequences on the performance of docking algorithms

BACKGROUND

The existence of conformational changes in antibodies upon binding has been previously established. However, existing analyses focus on individual cases and no quantitative study provides a more global view of potential moves and repacking, especially on recent data. The present study focuses on analyzing the conformational changes in various antibodies upon binding, providing quantitative observations to be exploited for antibody-related modeling.

METHODS

Cartesian and dihedral Root-Mean-Squared Deviations were calculated for different subparts of 27 different antibodies, for which X-ray structures in the bound and unbound states are available. Elbow angle variations were also calculated. Previously reported results of four docking algorithms were condensed into one score giving overall docking success for each of 16 antibody-antigen cases.

RESULTS

Very diverse movements are observed upon binding. While many loops stay very rigid, several others display side-chain repacking or backbone rearrangements, or both, at many different levels. Large conformational changes restricted to one or more antibody hypervariable loops were found to be a better indicator of docking difficulty than overall conformational variation at the antibody-antigen interface. However, the failure of docking algorithms on some almost-rigid cases shows that scoring is still a major bottleneck in docking pose prediction.

CONCLUSIONS

This study is aimed to help scientists working on antibody analysis and design by giving insights into the nature and the extent of conformational changes at different levels upon antigen binding.

Allosteric Effects between the Antibody Constant and Variable Regions: A Study of IgA Fc Mutations on Antigen Binding

Therapeutic antibodies have shifted the paradigm of disease treatments from small molecules to biologics, especially in cancer therapy. Despite the increasing number of antibody candidates, much remains unknown about the antibody and how its various regions interact. Recent findings showed that the antibody constant region can govern localization effects that are useful in reducing side effects due to systemic circulation by the commonly used IgG isotypes. Given their localized mucosal effects, IgA antibodies are increasingly promising therapeutic biologics. While the antibody Fc effector cell activity has been a focus point, recent research showed that the Fc could also influence antigen binding, challenging the conventional idea of region-specific antibody functions. To investigate this, we analysed the IgA antibody constant region and its distal effects on the antigen binding regions using recombinant Pertuzumab IgA1 and IgA2 variants. We found that mutations in the C-region reduced Her2 binding experimentally, and computational structural analysis showed that allosteric communications were highly dependent on the antibody hinge, providing strong evidence that we should consider antibodies as whole proteins rather than a sum of functional regions.