Next Article in Journal / Special Issue
Thinking outside the Laboratory: Analyses of Antibody Structure and Dynamics within Different Solvent Environments in Molecular Dynamics (MD) Simulations
Previous Article in Journal
Anti-Drug Antibodies: Emerging Approaches to Predict, Reduce or Reverse Biotherapeutic Immunogenicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), Singapore 138671, Singapore
2
p53 Laboratory, Agency for Science, Technology and Research (A*STAR), Singapore 138648, Singapore
*
Author to whom correspondence should be addressed.
Authors equally contributed to the work.
Submission received: 14 May 2018 / Revised: 2 June 2018 / Accepted: 5 June 2018 / Published: 7 June 2018
(This article belongs to the Collection Computational Antibody and Antigen Design)

Abstract

:
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.

Graphical Abstract

1. Introduction

Antibodies, called the “magic bullet” by Paul Erhlich [1,2,3], have shown great promise as therapeutic agents against numerous diseases [4], with many breakthroughs documented [5,6,7,8,9,10]. One promising isotype is IgA, whose predominant mucosal activation and localization can reduce systemic circulations and the associated side effects over the reigning IgG isotypes [9,10,11].
IgA (IgA1 and IgA2) is the major immunoglobulin isotype in adaptive mucosal immunity [12,13,14,15] and is responsible for several disease pathologies such as IgA nephropathy when they polymerize or self-aggregate [16,17]. Recently, a chimeric IgG-A antibody (with engineered CHγ2–CHα3 Fc region [18]) showed greater killing of Her2+ cancer cells by higher levels of complement-dependent cytotoxicity and activations of both neutrophils and macrophages. Although the chimeric IgG-A utilized only the CHα3 domain, this example clearly showed that the Fc of antibodies could be engineered towards various effector effects.
Fc manipulations have also been used to improve antibody half-lives [19] as well as to make bispecific antibodies [8] or create “sweeping” antibodies [20]. However, the overall effects of such constant region modifications on other antibody functions such as antigen bindings are not well established. Recently, there are increasing reports [21,22,23,24,25,26] of distant effects (likely allosteric communications) between the Fc region and the antigen-binding regions, with these studies typically based on IgG antibodies.
Similar to the case of IgG [22,26], our previous work [25] demonstrated that the heavy-chain constant regions can modulate antigen binding, most obviously for IgM and IgD, and to a lesser extent for IgA and its subtypes. To further investigate these effects, we generated mutations in the IgA constant regions and measured the antigen binding experimentally alongside computational analyses of allosteric communications between the constant and variable regions of these IgA antibodies.

2. Materials and Methods

2.1. Production of Recombinant Pertuzumab IgA Antibodies

Wild-type recombinant Pertuzumab IgA1 and IgA2 were synthesized and expressed as were previously described [25]. The mutations (C266Y/H317R for IgA1 and C253Y/H304R for IgA2): a conserved cysteine and the other residue randomly picked as a control were incorporated into the IgA constant regions by site-directed mutagenesis (Agilent Technologies, Santa Clara, CA, USA, Cat no. 200521). Produced IgA antibodies were quantified by spectrophotometric means using the extinction coefficient values determined from ProtParam [27]. Gel filtration figures were generated from Unicorn 6.0 software (GE Healthcare, Marlborough, MA, USA) with lines thickened using the GIMP 2.9.4 software. Purified antibody variants were analysed on reducing 10% SDS-PAGE gels and stained using Bio-Safe Coomassie stain (Bio-Rad, Hercules, CA, USA, Cat no. 1610786). Gel band sizes were determined using GelApp [28].

2.2. Binding Affinity Studies

Binding kinetics (using Blitz®, Fortebio, Fremont, CA, USA) of the antibodies to Her2 were carried out by pre-binding of HIS-tagged Her2 (Sino Biologicals Inc., Wayne, PA, USA, Cat no. H10004-H08H) onto the Ni-NTA (NTA) biosensors (Fortebio, Fremont, CA, USA, Cat no. 18-5101) as previously described and performed [25,29] using 1× kinetic buffer.

2.3. Modelling Full Antibody Structures of IgA1 and IgA2

Atomistic models of the two antibody variants IgA1 and IgA2 were constructed using two scattering-solved structures PDB: 2QTJ and PDB: 1R70 as templates for the Fc region, respectively. PDB: 1S78 was used as the template for the Pertuzumab Fab region. The resulting Cα-based backbones of the Fc regions were then used to construct the full-atom backbones and side chains using PULCHRA [30] and SCWRL4 [31], respectively. A standard procedure of energy minimization (5000 steps using steepest descent followed by conjugate gradient) was performed to remove possible clashes, using AMBER 14 [32]. Mutant IgA1 and IgA2 structures were modelled with corresponding mutations C266Y/H317R and C253Y/H304R, respectively.
The energy-minimized structures of the two variants (each including the wild type and mutant) were then subjected to coarse grain simulation (using Martini force field for proteins combined with ElNeDyn elastic network) to sample conformational changes of the whole antibody structures. The simulations were performed with time steps (dt) gradually increased from 15 fs to 22 fs during the equilibrium to accommodate ion wild motions, then fixed at dt = 22 fs during the production stages (3 × 1 μs) with the Verlet algorithm. Periodic boundary condition was also applied to avoid the finite size effects while simulating in explicit solvent (polarized water model; hence with PME). Temperature and pressure coupling schemes were used with the velocity rescale (V-rescale) and the Parrinello–Rahman barostat. Our analyses used the data from the last 600 ns (×3) of the simulated trajectories that reflected stable simulations, resulting in 3 × 1000 conformations.

2.4. Quantification of Allosteric Effects

We first used the minimized structures of the wild-type variants IgA1 and IgA2 to quantify the allosteric effects in both the Her2-binding and mutation events (as shown in Figure 1) using the server AlloSigMA [33], which have demonstrated successful quantification of allosteric effects in various benchmarked allosteric proteins [34,35,36,37]. The allosteric communications were estimated based on the responses of each residue (via residual free energy change Δgresidue) with respect to perturbations due to each of the events [33]. In this analysis, we simulated the mutations by initiating perturbations at the substituted positions (i.e., assigning “Up-mutation” in the AlloSigMA server to simulate larger residue substitutions). The resulting residue-wise allosteric free energies (with negative values indicating stabilizing and positive values indicating destabilizing effects) showed the quantified allosteric effect caused by the mutations. We then estimated the free energy change at the Her2 binding site (ΔgHer2site) and other corresponding regions (Δgregion) by averaging all Δgresidue values of the involving residues.
In addition, we estimated and clustered the distance ratio (i.e., distances between centre-of-mass of Fab regions versus those of both Fab and Fc as shown in Figure 2) of the wild-type structures resulting from the coarse grain simulation. Conformations nearest to the centroids were extracted and reverted to atomistic structures using Charmm36 force field and TIP3P water model, followed by short minimization and equilibrium. The structures were then used to study the spectrum of the allosteric effect driven by the domain motions.

2.5. Data Availability

The datasets generated and/or analysed during the current study are available upon reasonable request.

3. Results

Our previous work [25] suggested that the antibody heavy chain constant regions, but not light chain constant regions, influenced antigen binding beyond simple avidity effects, e.g., expected for IgM. To further investigate this phenomenon, we generated several disruptive mutations: C266Y/H317R in IgA1 and C253Y/H304R in IgA2 at the heavy chain constant region (CHα3) of our Pertuzumab IgA1 and IgA2 (i.e., intentionally substituting one of the disulphide-forming cysteine residues in the CH3 domain with the bulky residue tyrosine and another randomly selected histidine residue with similar positively charged arginine) to affect the heavy chain stability, and study the corresponding effects on the Her2-binding region.
Experimentally, we found the recombinant transient expression of both mutant variants to drop drastically (multiple folds) compared to the wild types (data not shown), implying the important role of the conserved disulphide-forming cysteine in the antibody stability [38]. A higher rate of aggregation in the isotype mutants was also observed (Figure S1). In addition, our binding kinetics measurements showed a significant decrease of the mutants by a log at 10−8 compared to the wild type at 10−9 for both IgA1 and IgA2. The major effects on IgA1 were at the association constant measurements for IgA1 and to a reduced extent, at the dissociation constant; for IgA2, the differences were less pronounced on both the association and dissociation constants (Figure 1A).
Our experimental measurements show that mutation-driven perturbations in the constant regions can affect Her2 binding, even with a few substitution mutations. This suggests clear allosteric communications between the two regions. We applied a structure-based statistical mechanical model [37] (using the AlloSigMA server [33]) to quantify these underlying allosteric effects. Results showed that the mutations caused stabilizing effects in the Fc region (Figure 1B and Table 1) that was compensated by increasing energy gain in distant regions (V-regions), thus destabilizing the Her2 binding site, which became more flexible (with ΔgHer2site > 0) due to contact losses between the Her2-interacting residues and their neighbours (Table 2). This destabilizing effect is more pronounced in the IgA1 mutant (Table 3 and Figure S2) than in IgA2.
We independently initiated computational perturbations at the Her2-binding site and studied the communication to the Fc in both the variants to simulate the bi-directional signal propagations between the two regions (Figure 1B,C). When comparing between the two events of mutations and Her2-binding, we noticed energy compensation to occur at the Fab regions of both IgA1 and IgA2. The IgA2-Fab domain compensated significantly more than that of IgA1 on the whole, distributing the energy compensation across the entire V-region, thus balancing the destabilization at the Her2 binding site (as Δ g H e r 2 s i t e I g A 2 < Δ g H e r 2 s i t e I g A 1 ) and retaining more binding to Her2. On the other hand, the energy compensation in the IgA1 Fab partially accommodated the changes at the Proline-rich hinge in terms of rigidity (Figure 1B). This might not be the case for the IgA2 short hinge, suggesting that the allosteric communication barrier was lifted with increasing hinge flexibility (Table 1 and Table 3).
Meanwhile, we hypothesized that the two mutation positions have functional selection constraints. To test this, we used the EVcoupling server [39] to investigate the residue couplings of the constant and variable regions (to demonstrate the local proximity as well as contacts in the antibody variants under the function dependent constraints). We found strong residue coupling networks forming independently within domains Fab or Fc, but weak links between the two domains (Figure 2). In the Fc domain, only the cysteine residue (C266 in IgA1 or C253 in IgA2) exhibited moderate functional dependence, which was expected for a conserved cysteine.
When additionally performing single perturbations at the individual positions, we found the allosteric signals to be elicited by the substitution of the histidine (H317 in IgA1 or H304 in IgA2) with a larger residue. Interestingly, the highly conserved disulphide-forming cysteine contributed to less destabilizing effects on the Her2 binding site (Table 3). These results indicate that the two domains Fab and Fc clearly communicate with each other, and that the histidine position might cause a bigger allosteric effect. We also found the destabilizing effect on the Her2-binding site to be a result of accumulative signalling, facilitated via the hinge. Therefore, the effect caused by the Fc mutations on the Her2 binding ability of the Fab domain (Figure 1A) is modulated by hinge flexibility.
The IgA1 isotype contains a longer hinge (connecting CH1 and the rest of the Fc region), the flexibility of which amplified with the mutation events (with Δ g h i n g e I g A 1 = 0.21 ± 0.20 whereas Δ g h i n g e I g A 2 = 0.08 ± 0.05). Scattering experiments on the wild-type IgA1 isotype (PDB: 2QTJ) showed that a rigid flanking hinge that separates Fab and Fc in an extended IgA1 conformation in solution is favoured. Our dynamics simulations of the wild-type IgA1 structure also showed that the two regions remained in constant proximity to each other. Nonetheless, diverse domain fluctuations between the two regions were observed in the conformational sampling of the IgA1 mutant (Figure 2A, top right), implying that mutations caused structural interferences and hence the hinge flexibility. On the other hand, the shorter IgA2 hinge did not allow for wild domain motions (Figure 2B, top right). Since the Her2-binding ability is mostly abolished in the IgA1 mutant but is retained in the IgA2 mutant, there is evidence that hinge flexibility modulates the propagation of the allosteric signals between the two distant domains.

4. Discussion

We set out to investigate the mechanism of the allosteric effects that the antibody constant regions elicited on the antigen binding as suggested in our previous work [25]. Working on our IgA1 and IgA2 models, we sought to study such effects in greater details by analysing the effects of Her2 binding when two mutations were introduced, at a conserved cysteine (C266/C253) and a random selected control histidine (H317/H304). While we acknowledge that the starting structure dependency remains a challenge in our approach, we have agreement between the computational and experimental results, where the computational observations could provide us with insight into the allostery phenomenon between the variable and constant regions.
Experimentally, we found that the double mutants were produced at a significantly reduced rate in our transient expressions, and that there were compromises in terms of binding ability to Her2, albeit at different magnitudes, despite unmodified variable regions. This further demonstrated that the constant region had clear significant effects on antigen binding, agreeing with other such studies on IgGs [21,22,24,40], in the IgA context. To understand the mechanism underlying our experimental observations, we showed that the allosteric signalling propagated bi-directionally between the two distant regions, from Her2-binding region to the Fc region and vice versa, on the IgA isotypes. In fact, the domain-linking hinge mediated such communication signalling, demonstrating that the hinge flexibility modulated the level of energy compensations at the distal regions.
We postulated that the energy compensation at the Fabs of the IgAs variants to accommodate the flexibility changes upon balance at the Her2-binding site (Table 1 and Table 3) to maintain the Her2-binding ability. As a result, our findings are in partial agreement with several previous studies [41,42,43] that the antigen-binding region requires a certain level of rigidity. On the contrary, other studies [44,45] suggested no significant conformational differences found in these regions, perhaps due to the different antigen binding contacts. It should also be noted that many of these previous studies were performed on antibody fragments, e.g., Fab or Fv, or on IgG, and not on the whole IgA. Thus, together with our other findings [25,26,46,47] also involving other proteins, we have shown the need to examine proteins as a whole for comprehensive holistic investigations, especially on allostery. In the context of the whole antibody, we found that the hinge rigidity and the allosteric communications between the distant constant and variable regions both contribute to antigen binding.
Our results here are also consistent with other previous studies [21,23] that showed induced conformational changes at the Fc due to the binding of antigen, hence consequently regulating the Fc receptor binding. Since there might also be other residues in the Fc region that could allosterically affect the antigen-binding site, we believe there would be a network of allosteric residues that drive the allosteric signalling between the regions. Our future work will thereby explore hotspots or mutation boundaries that can improve antigen binding by taking the advantages of allosteric communications from constant regions towards a rational and targeted approach for Fc engineering.
In conclusion, using both experimental and computational approaches, we were able to show that the constant region is important in the ability of the variable regions in binding antigens. We showed this in the less studied IgA, and that such effects were computationally found to be mediated by the hinge region of the antibody. This further illustrates the need to consider antibodies as a whole, rather than merely a sum of its regions.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4468/7/2/20/s1, Figure S1: Pertuzumab isotype/subtype biophysical analysis. Figure S2: Number of contacts within 6 Å (0.6 nm) of the Her2-interacting residues of the two variants IgA1 and IgA2.

Author Contributions

C.T.-T.S. designed and performed the computational analyses. W.-H.L. and W.-L.L. performed wet-lab experiments. W.-H.L., C.T.-T.S. and S.K.-E.G. analysed the results. C.T.-T.S., W.-H.L.,and S.K.-E.G. wrote the manuscript. S.K.-E.G. supervised the whole study. All authors read and approved the manuscript.

Funding

This work was supported by the Agency for Science, Technology, and Research (A*STAR) Joint Council Office (JCO1334i00050) and the A*STAR BioMedical Research Council Young Investigator Grant (BMRC YIG 1510651021) in Singapore.

Acknowledgments

We thank Santhosh Kumar Karthikeyan, Phua Ser Xian, Jonathan Kuan-Wai Tham, Chua Hua Lun, Ryan Chang, Joshua Yi Yeo and Chan Kwok Fong for their help in generating AlloSigMA data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Strebhardt, K.; Ullrich, A. Paul Ehrlich’s magic bullet concept: 100 years of progress. Nat. Rev. Cancer 2008, 8, 473–480. [Google Scholar] [CrossRef] [PubMed]
  2. Kohler, G.; Milstein, C. Continuous cultures of fused cells secreting antibody of predefined specificity. Nature 1975, 256, 495–497. [Google Scholar] [CrossRef] [PubMed]
  3. Ehrlich, P. Über den jetzigen Stand der Karzinomforschung. In Beiträge zur Experimentellen Pathologie und Chemotherapie; Akademische Verlagsgesellschaft: Leipzig, Germany, 1909; pp. 117–164. [Google Scholar]
  4. Strohl, W.R. Current progress in innovative engineered antibodies. Protein Cell 2018, 9, 86–120. [Google Scholar] [CrossRef] [PubMed]
  5. Igawa, T.; Ishii, S.; Tachibana, T.; Maeda, A.; Higuchi, Y.; Shimaoka, S.; Moriyama, C.; Watanabe, T.; Takubo, R.; Doi, Y.; et al. Antibody recycling by engineered pH-dependent antigen binding improves the duration of antigen neutralization. Nat. Biotechnol. 2010, 28, 1203–1207. [Google Scholar] [CrossRef] [PubMed]
  6. Robinson, M.; Ke, N.; Lobstein, J.; Peterson, C.; Szkodny, A.; Mansell, T.; Tuckey, C.; Riggs, P.; Colussi, P.; Noren, C.; et al. Efficient expression of full-length antibodies in the cytoplasm of engineered bacteria. Nat. Commun. 2015, 6, 8072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Park, H.; Yoon, H.; Jung, S. The highly evolvable Antibody Fc domain. Trends Biotechnol. 2016, 34, 895–908. [Google Scholar] [CrossRef] [PubMed]
  8. Ha, J.; Kim, J.; Kim, Y. Immunoglobulin Fc heterodimer platform technology: From design to applications in therapeutic antibodies and proteins. Front. Immunol. 2016, 7, 394. [Google Scholar] [CrossRef] [PubMed]
  9. Leusen, J.H. IgA as therapeutic antibody. Mol. Immunol. 2015, 68, 35–39. [Google Scholar] [CrossRef] [PubMed]
  10. Valerius, T.; Stockmeyer, B.; van Spriel, A.; Graziano, R.; van den Herik-Oudijk, I.; Repp, R.; Deo, Y.; Lund, J.; Kalden, J.; Gramatzki, M.; et al. FcαRI (CD89) as a novel trigger molecule for bispecific antibody therapy. Blood 1997, 90, 4485–4492. [Google Scholar] [PubMed]
  11. Bakema, J.E.; Egmond, M.V. Immunoglobulin A. mAbs 2011, 3, 352–361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Schroeder, H.J.; Cavacini, L. Structure and function of immunoglobulins. J. Allergy Clin. Immunol. 2010, 125, S41–S52. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Woof, J.M.; Kerr, M.A. IgA function—Variations on a theme. Immunology 2004, 113, 175–177. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, K.; Cerutti, A. The function and regulation of Immunoglobulin D. Curr. Opin. Immunol. 2011, 23, 345–352. [Google Scholar] [CrossRef] [PubMed]
  15. Preud’homme, J.-L.; Petit, I.; Barra, A.; Morel, F.; Lecron, J.-C.; Lelievre, E. Structural and functional properties of membrane and secreted IgD. Mol. Immunol. 2000, 37, 871–887. [Google Scholar] [CrossRef]
  16. Bonner, A.; Perrier, C.; Corthesy, B.; Perkins, S.J. Solution structure of human secretory component and implications for biological function. J. Biol. Chem. 2007, 282, 16969–16980. [Google Scholar] [CrossRef] [PubMed]
  17. Bonner, A.; Furtado, P.; Almogren, A.; Kerr, M.; Perkins, S. Implications of the near-planar solution structure of human myeloma dimeric IgA1 for mucosal immunity and IgA nephropathy. J. Immunol. 2008, 180, 1008–1018. [Google Scholar] [CrossRef] [PubMed]
  18. Kelton, W.; Mehta, N.; Charab, W.; Lee, J.; Lee, C.-H.; Kojima, T.; Kang, T.H.; Georgiou, G. IgGA: A “Cross-Isotype” engineered human Fc antibody domain that displays both IgG-like and IgA-like effector functions. Chem. Biol. 2014, 21, 1603–1609. [Google Scholar] [CrossRef] [PubMed]
  19. Borrok, M.; Mody, N.; Lu, X.; Kuhn, M.; Wu, H.; Dall’Acqua, W.; Tsui, P. An “Fc-Silenced” IgG1 format with extended half-life designed for improved stability. J. Pharm. Sci. 2017, 106, 1008–1017. [Google Scholar] [CrossRef] [PubMed]
  20. Igawa, T.; Maeda, A.; Haraya, K.; Tachibana, T.; Iwayanagi, Y.; Mimoto, F.; Higuchi, Y.; Ishii, S.; Tamba, S.; Hironiwa, N.; et al. Engineered monoclonal antibody with novel antigen-sweeping activity in vivo. PLoS ONE 2013, 8, e63236. [Google Scholar] [CrossRef] [PubMed]
  21. Zhao, J.; Nussinov, R.; Ma, B. Antigen induced dynamic conformation changes of antibody to facilitate recognition of Fc receptors. Biophys. J. 2018, 114 (Suppl. 1), 233a. [Google Scholar] [CrossRef]
  22. Janda, A.; Bowen, A.; Greenspan, N.S.; Casadevall, A. Ig constant region effects on Variable region structure and function. Front. Microbiol. 2016, 7, 22. [Google Scholar] [CrossRef] [PubMed]
  23. Oda, M.; Kozono, H.; Morii, H.; Azuma, T. Evidence of allosteric conformational changes in the antibody constant region upon antigen binding. Int. Immunol. 2003, 15, 417–426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yang, D.; Kroe-Barrett, R.; Singh, S.; Roberts, C.J.; Laue, T.M. IgG cooperativity—Is there allostery? Implications for antibody functions and therapeutic antibody development. mAbs 2017, 9, 1231–1252. [Google Scholar] [CrossRef] [PubMed]
  25. Lua, W.; Ling, W.; Yeo, J.; Poh, J.; Lane, D.; Gan, S. The effects of antibody engineering CH and CL in Trastuzumab and Pertuzumab recombinant models: Impact on antibody production and antigen-binding. Sci. Rep. 2018, 8, 718. [Google Scholar] [CrossRef] [PubMed]
  26. Ling, W.-L.; Lua, W.-H.; Poh, J.-J.; Yeo, J.Y.; Lane, D.P.; Gan, S.K.-E. Effect of VH-VL families in Pertuzumab and Trastuzumab recombinant production, Her2 and FcγIIA binding. Front. Immunol. 2018, 9, 469. [Google Scholar] [CrossRef] [PubMed]
  27. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.; Appel, R.; Bairoch, A. Protein identification and analysis tools on the ExPASy server. In The Proteomics Protocols Handbook; Walker, J.M., Ed.; Humana Press: New York, NY, USA, 2005; pp. 571–607. [Google Scholar]
  28. Sim, J.Z.; Nguyen, P.V.; Lee, H.K.; Gan, S. GelApp: Mobile gel electrophoresis analyser. Nat. Methods Appl. Notes 2015. [Google Scholar] [CrossRef]
  29. Lua, W.H.; Gan, S.; Lane, D.; Verma, C. A search for synergy in the binding kinetics of Trastuzumab and Pertuzumab whole and F(ab) to Her2. NPJ Breast Cancer 2015, 11, 15012. [Google Scholar] [CrossRef] [PubMed]
  30. Rotkiewicz, P.; Skolnick, J. Fast procedure for reconstruction of full-atom protein models from reduced representations. J. Comput. Chem. 2008, 29, 1460–1465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Krivov, G.; Shapovalov, M.; Dunbrack, R.J. Improved prediction of protein side-chain conformations with SCWRL4. Proteins 2009, 77, 778–795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Case, D.A.; Berryman, J.T.; Betz, R.M.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; et al. AMBER 14; University of California: San Francisco, CA, USA, 2015. [Google Scholar]
  33. Guarnera, E.; Tan, Z.W.; Zheng, Z.; Berezovsky, I.N. AlloSigMA: Allosteric Signalling and Mutation Analysis server. Bioinformatics 2017, 33, 3996–3998. [Google Scholar] [CrossRef] [PubMed]
  34. Chiang, R.; Gan, S.; Su, C. A computational study for rational HIV-1 non-nucleoside Reverse Transcriptase inhibitor selection and discovery of novel allosteric pockets for inhibitor design. Biosci. Rep. 2018. [Google Scholar] [CrossRef] [PubMed]
  35. Su, C.T.T.; Kwoh, C.K.; Verma, C.S.; Gan, S.K.E. Modeling the full length HIV-1 Gag polyprotein reveals the role of its p6 subunit in viral maturation and the effect of non-cleavage site mutations in protease drug resistance. J. Biomol. Struct. Dyn. 2017. [Google Scholar] [CrossRef] [PubMed]
  36. Kurochkin, I.V.; Guarnera, E.; Wong, J.H.; Eisenhaber, F.; Berezovsky, I.N. Toward allosterically increased catalytic activity of insulin-degrading enzyme against amyloid peptides. Biochemistry 2017, 56, 228–239. [Google Scholar] [CrossRef] [PubMed]
  37. Guarnera, E.; Berezovsky, I.N. Structure-based statistical mechanical model accounts for the causality and energetics of allosteric communication. PLoS Comput. Biol. 2016, 12, e1004678. [Google Scholar] [CrossRef] [PubMed]
  38. Feige, M.J.; Hendershot, L.M.; Buchner, J. How antibodies fold. Trends Biochem. Sci. 2009, 35, 189–198. [Google Scholar] [CrossRef] [PubMed]
  39. Hopf, T.A.; Colwell, L.J.; Sheridan, R.; Rost, B.; Sander, C.; Marks, D.S. Three-dimensional structures of membrane proteins from genomic sequencing. Cell 2012, 149, 1607–1621. [Google Scholar] [CrossRef] [PubMed]
  40. Knapp, B.; Dunbar, J.; Alcala, M.; Deane, C.M. Variable regions of antibodies and T-cell Receptors may not be sufficient in molecular simulations investigating binding. J. Chem. Theory Comput. 2017, 13, 3097–3105. [Google Scholar] [CrossRef] [PubMed]
  41. Zimmermann, J.; Oakman, E.L.; Thorpe, I.F.; Shi, X.; Abbyad, C.L.B.; Boxer, S.G.; Romesberg, F.E. Antibody evolution constrains conformational heterogeneity by tailoring protein dynamics. Proc. Natl. Acad. Sci. USA 2006, 103, 13722–13727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Chong, L.T.; Duan, Y.; Wang, L.; Massova, I.; Kollman, P.A. Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7. Proc. Natl. Acad. Sci. USA 1999, 96, 14330–14335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Li, T.; Tracka, M.B.; Uddin, S.; Casas-Finet, J.; Jacobs, D.J.; Livesay, D.R. Rigidity emerges during Antibody evolution in three distinct antibody systems: Evidence from QSFR analysis of Fab fragments. PLoS Comput. Biol. 2015, 11, e1004327. [Google Scholar] [CrossRef] [PubMed]
  44. Midelfort, K.; Hernandez, H.; Lippow, S.; Tidor, B.; Drennan, C.; Wittrup, K. Substantial energetic improvement with minimal structural perturbation in a high affinity mutant antibody. J. Mol. Biol. 2004, 343, 685–701. [Google Scholar] [CrossRef] [PubMed]
  45. Jeliazkov, J.R.; Sljoka, A.; Kuroda, D.; Tsuchimura, N.; Katoh, N.; Tsumoto, K.; Gray, J.J. Repertoire analysis of Antibody CDR-H3 loops suggests affinity maturation does not typically result in rigidfication. Front. Immunol. 2018, 9, 413. [Google Scholar] [CrossRef] [PubMed]
  46. Su, C.T.-T.; Ling, W.-L.; Lua, W.-H.; Poh, J.-J.; Gan, S.K.-E. The role of antibody Vκ Framework 3 region towards antigen binding: Effects on recombinant production and Protein L binding. Sci. Rep. 2017, 7, 3766. [Google Scholar] [CrossRef] [PubMed]
  47. Phua, S.-X.; Lua, W.-H.; Gan, S.K.-E. Role of FcaR EC2 region in extracellular membrane localization. Cell Cycle 2018. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Synergistic allosteric effects by the two IgA1 and IgA2 constant region mutations on the Her2-binding variable regions. (A) Binding kinetics analysis of the isotype variants IgA1 and IgA2 to Her2, using the antibodies at 200 nM to 25 nM to pre-loaded Her2 on NTA biosensors. The binding kinetics was measured using Blitz®. All experiments were performed in triplicates independently. The binding kinetics values of the wild-type IgA1 and IgA2 shown were obtained from our previous work [25]. (B) Surface presentations of the quantified allosteric communications (presented by residual allosteric free energy change Δgresidue) demonstrate destabilizing effects on the Her2 binding region caused by the mutations (black dots) in both the mutant constant region variants. (C) The quantified allosteric effects shown were based on the event of Her2 bindings. In (B,C), the effects were estimated using the minimized structures of both the wild-type IgA subtypes for perturbations with respect to mutating or binding events.
Figure 1. Synergistic allosteric effects by the two IgA1 and IgA2 constant region mutations on the Her2-binding variable regions. (A) Binding kinetics analysis of the isotype variants IgA1 and IgA2 to Her2, using the antibodies at 200 nM to 25 nM to pre-loaded Her2 on NTA biosensors. The binding kinetics was measured using Blitz®. All experiments were performed in triplicates independently. The binding kinetics values of the wild-type IgA1 and IgA2 shown were obtained from our previous work [25]. (B) Surface presentations of the quantified allosteric communications (presented by residual allosteric free energy change Δgresidue) demonstrate destabilizing effects on the Her2 binding region caused by the mutations (black dots) in both the mutant constant region variants. (C) The quantified allosteric effects shown were based on the event of Her2 bindings. In (B,C), the effects were estimated using the minimized structures of both the wild-type IgA subtypes for perturbations with respect to mutating or binding events.
Antibodies 07 00020 g001
Figure 2. Residue couplings and hinge-dependent domain proximity of the two variants IgA1 (in A) and IgA2 (in B). The sequence alignments of the mutation regions were performed using the EVcoupling server [39] and the coupling values were mapped back to the minimized structures of the variants. Distances were estimated between the centre of mass of the Fc region and that of the Her2-binding sites. Note that in the distance distributions of the IgA2 mutant (in B) showed the plotted values obtained from two replicates that successfully reached equilibrium in the given time scales.
Figure 2. Residue couplings and hinge-dependent domain proximity of the two variants IgA1 (in A) and IgA2 (in B). The sequence alignments of the mutation regions were performed using the EVcoupling server [39] and the coupling values were mapped back to the minimized structures of the variants. Distances were estimated between the centre of mass of the Fc region and that of the Her2-binding sites. Note that in the distance distributions of the IgA2 mutant (in B) showed the plotted values obtained from two replicates that successfully reached equilibrium in the given time scales.
Antibodies 07 00020 g002
Table 1. Accumulative allosteric effect on various regions (represented by Δgregion) using the minimized structures of IgA1 and IgA2, when mutated or when bound to Her2.
Table 1. Accumulative allosteric effect on various regions (represented by Δgregion) using the minimized structures of IgA1 and IgA2, when mutated or when bound to Her2.
Δgregion (kcal/mol)
Mutating EventHer2-Binding Event
RegionIgA1IgA2IgA1IgA2
Her2 binding site0.180.120.06−0.05
Fab0.150.380.080.02
Fc (CH2–CH3)−0.46−0.080.03−0.007
Hinge0.050.04−0.020.03
Table 2. Percentages * of native contacts (%) in the minimized structures of IgA1 and IgA2 mutants, showing the contact loss as compared to the wild type.
Table 2. Percentages * of native contacts (%) in the minimized structures of IgA1 and IgA2 mutants, showing the contact loss as compared to the wild type.
All Heavy Atoms/Cα Atoms
IgA1IgA2
Wild type100/100100/100
Mutant (replicate 1)69.7/58.357.1/61.8
Mutant (replicate 2)66.5/56.755.7/61.8
Mutant (replicate 3)66.5/56.754.7/61.8
* Percentage = (number of native contacts in mutant)/(number of native contacts in wild type).
Table 3. Allosteric free energy (Δgregion) estimated in mutation events (accumulative or single mutations) using different conformations of the wild-type variants IgA1 and IgA2.
Table 3. Allosteric free energy (Δgregion) estimated in mutation events (accumulative or single mutations) using different conformations of the wild-type variants IgA1 and IgA2.
Δgregion (kcal/mol) in the Event of Mutations
IgA1 (C266/H317)IgA2 (C253/H304)IgA1 (C266)IgA2 (C253)IgA1 (H317)IgA2 (H304)
Her2 site0.38 ± 0.320.12 ± 0.060.18 ± 0.160.06 ± 0.070.28 ± 0.230.09 ± 0.05
Fab0.62 ± 0.540.31 ± 0.120.28 ± 0.260.23 ± 0.090.5 ± 0.360.13 ± 0.08
Fc (CH2–CH3)−0.52 ± 0.05−0.16 ± 0.27−0.4 ± 0.07−0.06 ± 0.29−0.25 ± 0.09−0.1 ± 0.03
Hinge0.21 ± 0.20.08 ± 0.050.08 ± 0.09−0.003 ± 0.0320.20 ± 0.170.07 ± 0.04

Share and Cite

MDPI and ACS Style

Su, C.T.-T.; Lua, W.-H.; Ling, W.-L.; Gan, S.K.-E. Allosteric Effects between the Antibody Constant and Variable Regions: A Study of IgA Fc Mutations on Antigen Binding. Antibodies 2018, 7, 20. https://0-doi-org.brum.beds.ac.uk/10.3390/antib7020020

AMA Style

Su CT-T, Lua W-H, Ling W-L, Gan SK-E. Allosteric Effects between the Antibody Constant and Variable Regions: A Study of IgA Fc Mutations on Antigen Binding. Antibodies. 2018; 7(2):20. https://0-doi-org.brum.beds.ac.uk/10.3390/antib7020020

Chicago/Turabian Style

Su, Chinh Tran-To, Wai-Heng Lua, Wei-Li Ling, and Samuel Ken-En Gan. 2018. "Allosteric Effects between the Antibody Constant and Variable Regions: A Study of IgA Fc Mutations on Antigen Binding" Antibodies 7, no. 2: 20. https://0-doi-org.brum.beds.ac.uk/10.3390/antib7020020

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop