1
|
Rollins ZA, Widatalla T, Cheng AC, Metwally E. AbMelt: Learning antibody thermostability from molecular dynamics. Biophys J 2024:S0006-3495(24)00385-0. [PMID: 38851888 DOI: 10.1016/j.bpj.2024.06.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 03/16/2024] [Accepted: 06/04/2024] [Indexed: 06/10/2024] Open
Abstract
Antibody thermostability is challenging to predict from sequence and/or structure. This difficulty is likely due to the absence of direct entropic information. Herein, we present AbMelt where we model the inherent flexibility of homologous antibody structures using molecular dynamics simulations at three temperatures and learn the relevant descriptors to predict the temperatures of aggregation (Tagg), melt onset (Tm,on), and melt (Tm). We observed that the radius of gyration deviation of the complementarity determining regions at 400 K is the highest Pearson correlated descriptor with aggregation temperature (rp = -0.68 ± 0.23) and the deviation of internal molecular contacts at 350 K is the highest correlated descriptor with both Tm,on (rp = -0.74 ± 0.04) as well as Tm (rp = -0.69 ± 0.03). Moreover, after descriptor selection and machine learning regression, we predict on a held-out test set containing both internal and public data and achieve robust performance for all endpoints compared with baseline models (Tagg R2 = 0.57 ± 0.11, Tm,on R2 = 0.56 ± 0.01, and Tm R2 = 0.60 ± 0.06). In addition, the robustness of the AbMelt molecular dynamics methodology is demonstrated by only training on <5% of the data and outperforming more traditional machine learning models trained on the entire data set of more than 500 internal antibodies. Users can predict thermostability measurements for antibody variable fragments by collecting descriptors and using AbMelt, which has been made available.
Collapse
Affiliation(s)
- Zachary A Rollins
- Modeling and Informatics, Merck & Co., Inc., South San Francisco, California
| | - Talal Widatalla
- Modeling and Informatics, Merck & Co., Inc., South San Francisco, California
| | - Alan C Cheng
- Modeling and Informatics, Merck & Co., Inc., South San Francisco, California
| | - Essam Metwally
- Modeling and Informatics, Merck & Co., Inc., South San Francisco, California.
| |
Collapse
|
2
|
Park E, Izadi S. Molecular surface descriptors to predict antibody developability: sensitivity to parameters, structure models, and conformational sampling. MAbs 2024; 16:2362788. [PMID: 38853585 PMCID: PMC11168226 DOI: 10.1080/19420862.2024.2362788] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 05/29/2024] [Indexed: 06/11/2024] Open
Abstract
In silico assessment of antibody developability during early lead candidate selection and optimization is of paramount importance, offering a rapid and material-free screening approach. However, the predictive power and reproducibility of such methods depend heavily on the selection of molecular descriptors, model parameters, accuracy of predicted structure models, and conformational sampling techniques. Here, we present a set of molecular surface descriptors specifically designed for predicting antibody developability. We assess the performance of these descriptors by benchmarking their correlations with an extensive array of experimentally determined biophysical properties, including viscosity, aggregation, hydrophobic interaction chromatography, human pharmacokinetic clearance, heparin retention time, and polyspecificity. Further, we investigate the sensitivity of these surface descriptors to methodological nuances, such as the choice of interior dielectric constant, hydrophobicity scales, structure prediction methods, and the impact of conformational sampling. Notably, we observe systematic shifts in the distribution of surface descriptors depending on the structure prediction method used, driving weak correlations of surface descriptors across structure models. Averaging the descriptor values over conformational distributions from molecular dynamics mitigates the systematic shifts and improves the consistency across different structure prediction methods, albeit with inconsistent improvements in correlations with biophysical data. Based on our benchmarking analysis, we propose six in silico developability risk flags and assess their effectiveness in predicting potential developability issues for a set of case study molecules.
Collapse
Affiliation(s)
- Eliott Park
- Pharmaceutical Development, Genentech Inc, South San Francisco, CA, USA
| | - Saeed Izadi
- Pharmaceutical Development, Genentech Inc, South San Francisco, CA, USA
| |
Collapse
|
3
|
Math BA, Waibl F, Lamp LM, Fernández‐Quintero ML, Liedl KR. Cross-linking disulfide bonds govern solution structures of diabodies. Proteins 2023; 91:1316-1328. [PMID: 37376973 PMCID: PMC10952579 DOI: 10.1002/prot.26509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Accepted: 04/19/2023] [Indexed: 06/29/2023]
Abstract
In the last years, antibodies have emerged as a promising new class of therapeutics, due to their combination of high specificity with long serum half-life and low risk of side-effects. Diabodies are a popular novel antibody format, consisting of two Fv domains connected with short linkers. Like IgG antibodies, they simultaneously bind two target proteins. However, they offer altered properties, given their smaller size and higher rigidity. In this study, we conducted the-to our knowledge-first molecular dynamics (MD) simulations of diabodies and find a surprisingly high conformational flexibility in the relative orientation of the two Fv domains. We observe rigidifying effects through the introduction of disulfide bonds in the Fv -Fv interface and characterize the effect of different disulfide bond locations on the conformation. Additionally, we compare VH -VL orientations and paratope dynamics between diabodies and an antigen binding fragment (Fab) of the same sequence. We find mostly consistent structures and dynamics, indicating similar antigen binding properties. The most significant differences can be found within the CDR-H2 loop dynamics. Of all CDR loops, the CDR-H2 is located closest to the artificial Fv -Fv interface. All examined diabodies show similar VH -VL orientations, Fv -Fv packing and CDR loop conformations. However, the variant with a P14C-K64C disulfide bond differs most from the Fab in our measures, including the CDR-H3 loop conformational ensemble. This suggests altered antigen binding properties and underlines the need for careful validation of the disulfide bond locations in diabodies.
Collapse
Affiliation(s)
- Barbara A. Math
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Franz Waibl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Leonida M. Lamp
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Monica L. Fernández‐Quintero
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Klaus R. Liedl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| |
Collapse
|
4
|
Structural mechanism of Fab domain dissociation as a measure of interface stability. J Comput Aided Mol Des 2023; 37:201-215. [PMID: 36918473 PMCID: PMC10049950 DOI: 10.1007/s10822-023-00501-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2023] [Accepted: 02/23/2023] [Indexed: 03/16/2023]
Abstract
Therapeutic antibodies should not only recognize antigens specifically, but also need to be free from developability issues, such as poor stability. Thus, the mechanistic understanding and characterization of stability are critical determinants for rational antibody design. In this study, we use molecular dynamics simulations to investigate the melting process of 16 antigen binding fragments (Fabs). We describe the Fab dissociation mechanisms, showing a separation in the VH-VL and in the CH1-CL domains. We found that the depths of the minima in the free energy curve, corresponding to the bound states, correlate with the experimentally determined melting temperatures. Additionally, we provide a detailed structural description of the dissociation mechanism and identify key interactions in the CDR loops and in the CH1-CL interface that contribute to stabilization. The dissociation of the VH-VL or CH1-CL domains can be represented by conformational changes in the bend angles between the domains. Our findings elucidate the melting process of antigen binding fragments and highlight critical residues in both the variable and constant domains, which are also strongly germline dependent. Thus, our proposed mechanisms have broad implications in the development and design of new and more stable antigen binding fragments.
Collapse
|
5
|
Fernández-Quintero ML, Kokot J, Waibl F, Fischer ALM, Quoika PK, Deane CM, Liedl KR. Challenges in antibody structure prediction. MAbs 2023; 15:2175319. [PMID: 36775843 PMCID: PMC9928471 DOI: 10.1080/19420862.2023.2175319] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023] Open
Abstract
Advances in structural biology and the exponential increase in the amount of high-quality experimental structural data available in the Protein Data Bank has motivated numerous studies to tackle the grand challenge of predicting protein structures. In 2020 AlphaFold2 revolutionized the field using a combination of artificial intelligence and the evolutionary information contained in multiple sequence alignments. Antibodies are one of the most important classes of biotherapeutic proteins. Accurate structure models are a prerequisite to advance biophysical property predictions and consequently antibody design. Specialized tools used to predict antibody structures based on different principles have profited from current advances in protein structure prediction based on artificial intelligence. Here, we emphasize the importance of reliable protein structure models and highlight the enormous advances in the field, but we also aim to increase awareness that protein structure models, and in particular antibody models, may suffer from structural inaccuracies, namely incorrect cis-amide bonds, wrong stereochemistry or clashes. We show that these inaccuracies affect biophysical property predictions such as surface hydrophobicity. Thus, we stress the importance of carefully reviewing protein structure models before investing further computing power and setting up experiments. To facilitate the assessment of model quality, we provide a tool "TopModel" to validate structure models.
Collapse
Affiliation(s)
| | - Janik Kokot
- Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Franz Waibl
- Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Anna-Lena M. Fischer
- Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Patrick K. Quoika
- Center for Protein Assemblies (CPA), Physics Department, Chair of Theoretical Biophysics, Technical University of Munich, Garching, Germany
| | | | - Klaus R. Liedl
- CONTACT Klaus R. Liedl Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
6
|
Fernández-Quintero ML, Kroell KB, Grunewald LJ, Fischer ALM, Riccabona JR, Liedl KR. CDR loop interactions can determine heavy and light chain pairing preferences in bispecific antibodies. MAbs 2022; 14:2024118. [PMID: 35090383 PMCID: PMC8803122 DOI: 10.1080/19420862.2021.2024118] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022] Open
Abstract
As the current biotherapeutic market is dominated by antibodies, the design of different antibody formats, like bispecific antibodies, is critical to the advancement of the field. In contrast to monovalent antibodies, which consist of two identical antigen-binding sites, bispecific antibodies can target two different epitopes by containing two different antigen-binding sites. Thus, the rise of new formats as successful therapeutics has reignited the interest in advancing and facilitating the efficient production of bispecific antibodies. Here, we investigate the influence of point mutations in the antigen-binding site, the paratope, on heavy and light chain pairing preferences by using molecular dynamics simulations. In agreement with experiments, we find that specific residues in the antibody variable domain (Fv), i.e., the complementarity-determining region (CDR) L3 and H3 loops, determine heavy and light chain pairing preferences. Excitingly, we observe substantial population shifts in CDR-H3 and CDR-L3 loop conformations in solution accompanied by a decrease in bispecific IgG yield. These conformational changes in the CDR3 loops induced by point mutations also influence all other CDR loop conformations and consequentially result in different CDR loop states in solution. However, besides their effect on the obtained CDR loop ensembles, point mutations also lead to distinct interaction patterns in the VH-VL interface. By comparing the interaction patterns among all investigated variants, we observe specific contacts in the interface that drive heavy and light chain pairing. Thus, these findings have broad implications in the field of antibody engineering and design because they provide a mechanistic understanding of antibody interfaces, by identifying critical factors driving the pairing preferences, and thus can help to advance the design of bispecific antibodies.
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Katharina B Kroell
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Lukas J Grunewald
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Anna-Lena M Fischer
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Jakob R Riccabona
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Klaus R Liedl
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
7
|
Fernández-Quintero ML, Kroell KB, Bacher LM, Loeffler JR, Quoika PK, Georges G, Bujotzek A, Kettenberger H, Liedl KR. Germline-Dependent Antibody Paratope States and Pairing Specific V H-V L Interface Dynamics. Front Immunol 2021; 12:675655. [PMID: 34447370 PMCID: PMC8382685 DOI: 10.3389/fimmu.2021.675655] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Accepted: 06/24/2021] [Indexed: 11/13/2022] Open
Abstract
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 VH-VL 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 VH-VL 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 VH-VL 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.
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Katharina B Kroell
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Lisa M Bacher
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Johannes R Loeffler
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Patrick K Quoika
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Alexander Bujotzek
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Hubert Kettenberger
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Klaus R Liedl
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
8
|
Fernández-Quintero ML, Georges G, Varga JM, Liedl KR. Ensembles in solution as a new paradigm for antibody structure prediction and design. MAbs 2021; 13:1923122. [PMID: 34030577 PMCID: PMC8158028 DOI: 10.1080/19420862.2021.1923122] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
The rise of antibodies as a promising and rapidly growing class of biotherapeutic proteins has motivated numerous studies to characterize and understand antibody structures. In the past decades, the number of antibody crystal structures increased substantially, which revolutionized the atomistic understanding of antibody functions. Even though numerous static structures are known, various biophysical properties of antibodies (i.e., specificity, hydrophobicity and stability) are governed by their dynamic character. Additionally, the importance of high-quality structures in structure–function relationship studies has substantially increased. These structure–function relationship studies have also created a demand for precise homology models of antibody structures, which allow rational antibody design and engineering when no crystal structure is available. Here, we discuss various aspects and challenges in antibody design and extend the paradigm of describing antibodies with only a single static structure to characterizing them as dynamic ensembles in solution.
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Janos M Varga
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Klaus R Liedl
- Department of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
9
|
Fernández-Quintero ML, Kroell KB, Hofer F, Riccabona JR, Liedl KR. Mutation of Framework Residue H71 Results in Different Antibody Paratope States in Solution. Front Immunol 2021; 12:630034. [PMID: 33737932 PMCID: PMC7960778 DOI: 10.3389/fimmu.2021.630034] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Accepted: 01/20/2021] [Indexed: 12/13/2022] Open
Abstract
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 (VH-VL). 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 71H (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 71H, 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 71H on the paratope states and interface orientations has broad implications in the field of antibody engineering.
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Katharina B Kroell
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Florian Hofer
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Jakob R Riccabona
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Klaus R Liedl
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
10
|
Fernández-Quintero ML, Heiss MC, Pomarici ND, Math BA, Liedl KR. Antibody CDR loops as ensembles in solution vs. canonical clusters from X-ray structures. MAbs 2021; 12:1744328. [PMID: 32264741 PMCID: PMC7153821 DOI: 10.1080/19420862.2020.1744328] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
In the past decade, the relevance of antibodies as therapeutics has increased substantially. Therefore, structural and functional characterization, in particular of the complementarity-determining regions (CDRs), is crucial to the design and engineering of antibodies with unique binding properties. Various studies have focused on classifying the CDR loops into a small set of main-chain conformations to facilitate antibody design by assuming that certain sequences can only adopt a limited number of conformations. Here, we present a kinetic classification of CDR loop structures as ensembles in solution. Using molecular dynamics simulations in combination with strong experimental structural information, we observe conformational transitions between canonical clusters and additional dominant solution structures in the micro-to-millisecond timescale for all CDR loops, independent of length and sequence composition. Besides identifying all relevant conformations in solution, our results revealed that various canonical cluster medians actually belong to the same kinetic minimum. Additionally, we reconstruct the kinetics and probabilities of the conformational transitions between canonical clusters, and thereby extend the model of static canonical structures to reveal a dynamic conformational ensemble in solution as a new paradigm in the field of antibody structure design. Abbreviations: CDR: Complementary-determining region; Fv: Antibody variable fragment; PCCA: Perron cluster analysis; tICA: Time-lagged independent component analysis; VH: Heavy chain variable region; VL: Light chain variable region
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Martin C Heiss
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Nancy D Pomarici
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Barbara A Math
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| | - Klaus R Liedl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
11
|
Hofer F, Kamenik AS, Fernández-Quintero ML, Kraml J, Liedl KR. pH-Induced Local Unfolding of the Phl p 6 Pollen Allergen From cpH-MD. Front Mol Biosci 2021; 7:603644. [PMID: 33511157 PMCID: PMC7835895 DOI: 10.3389/fmolb.2020.603644] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 12/10/2020] [Indexed: 11/13/2022] Open
Abstract
Susceptibility to endosomal degradation is a decisive contribution to a protein's immunogenicity. It is assumed that the processing kinetics of structured proteins are inherently linked to their probability of local unfolding. In this study, we quantify the impact of endosomal acidification on the conformational stability of the major timothy grass pollen allergen Phl p 6. We use state of the art sampling approaches in combination with constant pH MD techniques to profile pH-dependent local unfolding events in atomistic detail. Integrating our findings into the current view on type 1 allergic sensitization, we characterize local protein dynamics in the context of proteolytic degradation at neutral and acidic pH for the wild type protein and point mutants with varying proteolytic stability. We analyze extensive simulation data using Markov state models and retrieve highly reliable thermodynamic and kinetic information at varying pH levels. Thereby we capture the impact of endolysosomal acidification on the structure and dynamics of the Phl p 6 mutants. We find that upon protonation at lower pH values, the conformational flexibilities in key areas of the wild type protein, i.e., T-cell epitopes and early proteolytic cleavage sites, increase significantly. A decrease of the pH even leads to local unfolding in otherwise stable secondary structure elements, which is a prerequisite for proteolytic cleavage. This effect is even more pronounced in the destabilized mutant, while no unfolding was observed for the stabilized mutant. In summary, we report detailed structural models which rationalize the experimentally observed cleavage pattern during endosomal acidification.
Collapse
|
12
|
Fernández-Quintero ML, Kroell KB, Heiss MC, Loeffler JR, Quoika PK, Waibl F, Bujotzek A, Moessner E, Georges G, Liedl KR. Surprisingly Fast Interface and Elbow Angle Dynamics of Antigen-Binding Fragments. Front Mol Biosci 2020; 7:609088. [PMID: 33330636 PMCID: PMC7732698 DOI: 10.3389/fmolb.2020.609088] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2020] [Accepted: 10/21/2020] [Indexed: 12/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Monica L. Fernández-Quintero
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Katharina B. Kroell
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Martin C. Heiss
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Johannes R. Loeffler
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Patrick K. Quoika
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Franz Waibl
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Alexander Bujotzek
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Ekkehard Moessner
- Roche Pharma Research and Early Development, Large Molecular Research, Roche Innovation Center Zurich, Schlieren, Switzerland
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Klaus R. Liedl
- Center for Molecular Biosciences Innsbruck (CMBI), Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
13
|
Waibl F, Fernández-Quintero ML, Kamenik AS, Kraml J, Hofer F, Kettenberger H, Georges G, Liedl KR. Conformational Ensembles of Antibodies Determine Their Hydrophobicity. Biophys J 2020; 120:143-157. [PMID: 33220303 PMCID: PMC7820740 DOI: 10.1016/j.bpj.2020.11.010] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Revised: 10/23/2020] [Accepted: 11/10/2020] [Indexed: 12/11/2022] Open
Abstract
A major challenge in the development of antibody biotherapeutics is their tendency to aggregate. One root cause for aggregation is exposure of hydrophobic surface regions to the solvent. Many current techniques predict the relative aggregation propensity of antibodies via precalculated scales for the hydrophobicity or aggregation propensity of single amino acids. However, those scales cannot describe the nonadditive effects of a residue’s surrounding on its hydrophobicity. Therefore, they are inherently limited in their ability to describe the impact of subtle differences in molecular structure on the overall hydrophobicity. Here, we introduce a physics-based approach to describe hydrophobicity in terms of the hydration free energy using grid inhomogeneous solvation theory (GIST). We apply this method to assess the effects of starting structures, conformational sampling, and protonation states on the hydrophobicity of antibodies. Our results reveal that high-quality starting structures, i.e., crystal structures, are crucial for the prediction of hydrophobicity and that conformational sampling can compensate errors introduced by the starting structure. On the other hand, sampling of protonation states only leads to good results when combined with high-quality structures, whereas it can even be detrimental otherwise. We conclude by pointing out that a single static homology model may not be adequate for predicting hydrophobicity.
Collapse
Affiliation(s)
- Franz Waibl
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck, Austria
| | - Monica L Fernández-Quintero
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck, Austria
| | - Anna S Kamenik
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck, Austria
| | - Johannes Kraml
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck, Austria
| | - Florian Hofer
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck, Austria
| | - Hubert Kettenberger
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Klaus R Liedl
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck, Austria.
| |
Collapse
|
14
|
Fernández-Quintero ML, Pomarici ND, Math BA, Kroell KB, Waibl F, Bujotzek A, Georges G, Liedl KR. Antibodies exhibit multiple paratope states influencing V H-V L domain orientations. Commun Biol 2020; 3:589. [PMID: 33082531 PMCID: PMC7576833 DOI: 10.1038/s42003-020-01319-z] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2020] [Accepted: 09/23/2020] [Indexed: 11/17/2022] Open
Abstract
In the last decades, antibodies have emerged as one of the most important and successful classes of biopharmaceuticals. The highest variability and diversity of an antibody is concentrated on six hypervariable loops, also known as complementarity determining regions (CDRs) shaping the antigen-binding site, the paratope. Whereas it was assumed that certain sequences can only adopt a limited set of backbone conformations, in this study we present a kinetic classification of several paratope states in solution. Using molecular dynamics simulations in combination with experimental structural information we capture the involved conformational transitions between different canonical clusters and additional dominant solution structures occurring in the micro-to-millisecond timescale. Furthermore, we observe a strong correlation of CDR loop movements. Another important aspect when characterizing different paratope states is the relative VH/VL orientation and the influence of the distinct CDR loop states on the VH/VL interface. Conformational rearrangements of the CDR loops do not only have an effect on the relative VH/VL orientations, but also influence in some cases the elbow-angle dynamics and shift the respective distributions. Thus, our results show that antibodies exist as several interconverting paratope states, each contributing to the antibody's properties.
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, A-6020, Innsbruck, Austria
| | - Nancy D Pomarici
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, A-6020, Innsbruck, Austria
| | - Barbara A Math
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, A-6020, Innsbruck, Austria
| | - Katharina B Kroell
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, A-6020, Innsbruck, Austria
| | - Franz Waibl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, A-6020, Innsbruck, Austria
| | - Alexander Bujotzek
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Klaus R Liedl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, A-6020, Innsbruck, Austria.
| |
Collapse
|
15
|
Pomés A, Mueller GA, Chruszcz M. Structural Aspects of the Allergen-Antibody Interaction. Front Immunol 2020; 11:2067. [PMID: 32983155 PMCID: PMC7492603 DOI: 10.3389/fimmu.2020.02067] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2020] [Accepted: 07/29/2020] [Indexed: 12/23/2022] Open
Abstract
The development of allergic disease involves the production of IgE antibodies upon allergen exposure in a process called sensitization. IgE binds to receptors on the surface of mast cells and basophils, and subsequent allergen exposure leads to cross-linking of IgE antibodies and release of cell mediators that cause allergy symptoms. Although this process is quite well-understood, very little is known about the epitopes on the allergen recognized by IgE, despite the importance of the allergen-antibody interaction for the allergic response to occur. This review discusses efforts to analyze allergen-antibody interactions, from the original epitope mapping studies using linear peptides or recombinant allergen fragments, to more sophisticated technologies, such as X-ray crystallography and nuclear magnetic resonance. These state-of-the-art approaches, combined with site-directed mutagenesis, have led to the identification of conformational IgE epitopes. The first structures of an allergen (egg lysozyme) in complex with Fab fragments from IgG antibodies were determined in the 1980s. Since then, IgG has been used as surrogate for IgE, due to the difficulty of obtaining monoclonal IgE antibodies. Technical developments including phage display libraries have contributed to progress in epitope mapping thanks to the isolation of IgE antibody constructs from combinatorial libraries made from peripheral blood mononuclear cells of allergic donors. Most recently, single B cell antibody sequencing and human hybridomas are new breakthrough technologies for finally obtaining human IgE monoclonal antibodies, ideal for epitope mapping. The information on antigenic determinants will facilitate the design of hypoallergens for immunotherapy and the investigation of the fundamental mechanisms of the IgE response.
Collapse
Affiliation(s)
- Anna Pomés
- Indoor Biotechnologies, Inc., Charlottesville, VA, United States
| | - Geoffrey A Mueller
- National Institute of Environmental Health Sciences, Durham, NC, United States
| | - Maksymilian Chruszcz
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, SC, United States
| |
Collapse
|
16
|
Fernández-Quintero ML, Pomarici ND, Loeffler JR, Seidler CA, Liedl KR. T-Cell Receptor CDR3 Loop Conformations in Solution Shift the Relative Vα-Vβ Domain Distributions. Front Immunol 2020; 11:1440. [PMID: 32733478 PMCID: PMC7360859 DOI: 10.3389/fimmu.2020.01440] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2020] [Accepted: 06/03/2020] [Indexed: 12/12/2022] Open
Abstract
T-cell receptors are an important part in the adaptive immune system as they are responsible for detecting foreign proteins presented by the major histocompatibility complex (MHC). The affinity is predominantly determined by structure and sequence of the complementarity determining regions (CDRs), of which the CDR3 loops are responsible for peptide recognition. We present a kinetic classification of T-cell receptor CDR3 loops with different loop lengths into canonical and non-canonical solution structures. Using molecular dynamics simulations, we do not only sample available X-ray structures, but we also observe a substantially broader CDR3 loop ensemble with various distinct kinetic minima in solution. Our results strongly imply, that for given CDR3 loop sequences several canonical structures have to be considered to characterize the conformational diversity of these loops. Our suggested dominant solution structures could extend the repertoire of available canonical clusters by including kinetic minimum structures present in solution. Thus, the CDR3 loops need to be characterized as conformational ensembles in solution. Furthermore, the conformational changes of the CDR3 loops follow the paradigm of conformational selection, because the experimentally determined binding competent state is present within this ensemble of pre-existing conformations without the presence of the antigen. We also identify strong correlations between the CDR3 loops and include combined state descriptions. Additionally, we observe a strong dependency of the CDR3 loop conformations on the relative Vα-Vβ interdomain orientations, revealing that certain CDR3 loop states favor specific interface orientations.
Collapse
MESH Headings
- Adaptive Immunity
- Animals
- Antigens/metabolism
- Complementarity Determining Regions/chemistry
- Complementarity Determining Regions/genetics
- Complementarity Determining Regions/metabolism
- Crystallography, X-Ray
- Histocompatibility Antigens/metabolism
- Humans
- Molecular Dynamics Simulation
- Peptides/metabolism
- Protein Binding
- Protein Conformation
- Protein Domains/genetics
- Receptors, Antigen, T-Cell, alpha-beta/chemistry
- Receptors, Antigen, T-Cell, alpha-beta/genetics
- Receptors, Antigen, T-Cell, alpha-beta/metabolism
- Structure-Activity Relationship
- T-Cell Antigen Receptor Specificity
Collapse
Affiliation(s)
| | | | | | | | - Klaus R. Liedl
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
17
|
Fernández-Quintero ML, Loeffler JR, Waibl F, Kamenik AS, Hofer F, Liedl KR. Conformational selection of allergen-antibody complexes-surface plasticity of paratopes and epitopes. Protein Eng Des Sel 2019; 32:513-523. [PMID: 32719844 PMCID: PMC7451023 DOI: 10.1093/protein/gzaa014] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Revised: 06/26/2020] [Accepted: 06/26/2020] [Indexed: 12/18/2022] Open
Abstract
Antibodies have the ability to bind various types of antigens and to recognize different antibody-binding sites (epitopes) of the same antigen with different binding affinities. Due to the conserved structural framework of antibodies, their specificity to antigens is mainly determined by their antigen-binding site (paratope). Therefore, characterization of epitopes in combination with describing the involved conformational changes of the paratope upon binding is crucial in understanding and predicting antibody-antigen binding. Using molecular dynamics simulations complemented with strong experimental structural information, we investigated the underlying binding mechanism and the resulting local and global surface plasticity in the binding interfaces of distinct antibody-antigen complexes. In all studied allergen-antibody complexes, we clearly observe that experimentally suggested epitopes reveal less plasticity, while non-epitope regions show high surface plasticity. Surprisingly, the paratope shows higher conformational diversity reflected in substantially higher surface plasticity, compared to the epitope. This work allows a visualization and characterization of antibody-antigen interfaces and might have strong implications for antibody-antigen docking and in the area of epitope prediction.
Collapse
Affiliation(s)
- Monica L Fernández-Quintero
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, A-6020 Innsbruck, Austria
| | - Johannes R Loeffler
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, A-6020 Innsbruck, Austria
| | - Franz Waibl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, A-6020 Innsbruck, Austria
| | - Anna S Kamenik
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, A-6020 Innsbruck, Austria
| | - Florian Hofer
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, A-6020 Innsbruck, Austria
| | - Klaus R Liedl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, A-6020 Innsbruck, Austria,To whom correspondence should be addressed. E-Mail:
| |
Collapse
|