|
|
||||||||
1 School of Pharmaceutical Sciences, Kitasato University, 5-9-1 Shirokane, Minato-ku, Tokyo 108-8641, Japan
2 School of Science, Kitasato University, Kanagawa 228-8555, Japan
Reprint requests to: Hideaki Umeyama, School of Pharmaceutical Sciences, Kitasato University, 5-9-1 Shirokane, Minato-ku, Tokyo 108-8641, Japan; e-mail: umeyamah{at}pharm.kitasato-u.ac.jp; fax: 81-3-3446-9553.
(RECEIVED March 30, 2003; FINAL REVISION June 13, 2003; ACCEPTED June 13, 2003)
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.03100803.
| Abstract |
|---|
|
|
|---|
fluctuations of HELFab differed from those of HELFv with normal mode calculation, and the dynamics can be thought to be related to the proteinprotein interactions. CLCH1 may have influence not only around local interfaces between CLCH1 and Fv, but also around the interacting regions between HEL and Fv, which are longitudinally distant. Eighteen water molecules were found in HELFv around the interface between HEL and Fv compared with one water molecule in HELFab. These solvent molecules may occupy the holes and channels, which may occur due to imperfect complementarity of the complex. Therefore, the suppression of atomic vibration around the interface between Fv and HEL can be thought to be related to favorable and compact interface formation by complete desolvation. It is suggested that the ability to control the antigenantibody affinity is obtained from modifying the CLCH1. The second upper loop in the constant domain of the light chain (UL2CL), which is a conserved gene in several light chains, showed the most remarkable fluctuation changes. UL2CL could play an important role and could be attractive for modification in protein engineering. Keywords: Immunoglobulin; constant domain; protein interaction; normal mode analysis; affinity
Abbreviations: Fab, antigen binding fragment of immunoglobulin Fv, variable domains of immunoglobulin CL, constant domain of immunoglobulin light chain CH1, first constant domain of immunoglobulin heavy chain HEL, hen egg white lysozyme VH, variable domain of immunoglobulin heavy chain VL, variable domain of immunoglobulin light chain CDR, complementary determining region CDR1-VL, first CDR of the variable domain in the light chain FR, lower loop as framework region UL, upper loop of constant domain
| Introduction |
|---|
|
|
|---|
NMA has been developed as another tool for the study of the dynamic behavior of biologic macromolecules (Levitt et al. 1985). The method includes some limitations, for instance, the harmonic approximation and conventional calculation in vacuo. Levitt et al. (1985) reported that the diffusing water molecules cannot be treated, as they do not move in a harmonic potential well, and that the presence of surrounding water would not be expected to have a major effect. Therefore, exclusion of water molecules in the calculation may not change the conclusion, and can avoid a cost of an additional six degrees of freedom per water molecule. NMA is still useful for understanding protein dynamics, especially for directly describing subtle fluctuations of rigid bodies (Soejima et al. 1999). Moreover, a more important usefulness of NMA is that it might functionally describe collective motions from the normal mode (Brooks and Karplus 1983; Go et al. 1983; Nishikawa and Go 1987; Hayward et al. 1995; Ishida et al. 1998; Jaaskelainen et al. 1998; Miller and Agard 1999; Nojima et al. 2002).
The three-dimensional structures of three monoclonal antibodies, D1.3, HyHEL-5, and HyHEL-10, have been studied during binding to different sites on the surface of hen egg white lysozyme (HEL). The antibodies bind to HEL specifically with six upper loops of variable domains (Fv), termed complementary determining regions (CDRs; Ward et al. 1989; Batra et al. 1990; Fishmann et al. 1991). In addition, the principles of the antigenantibody binding were explored with mutagenesis around the antigenantibody interface (Lavoie et al. 1992; Kam-Morgan et al. 1993; Tsumoto et al. 1994, 1996; Sundburg et al. 2000).
In the HyHEL-10 system, the crystal structure of HELFv (PDB code: 1C08 [PDB] ), which is a complex of HEL and Fv, had no significant differences from that of HELFab (PDB code: 3HFM [PDB] ), which is a complex of HEL and Fab. However, it was reported HELFv had a one order of magnitude lower dissociation constant than HELFab in the system: HELFab and HELFv had 450 x 108 M-1, 4.2 x 108 M-1, respectively (Smith-Gill et al. 1984; Lavoie et al. 1992; Kam-Morgan et al. 1993; Tsumoto et al. 1994). We focused on the dynamics of HELFab and HELFv, and attempted to explore the difference in the interaction between the antigen and the antibody.
Figure 1
shows the schematic structure of HELHyHEL-10 (Padlan et al. 1989). Fab is composed of the light and the heavy chains (L and H). The heavy chain cut between CH1 and CH2 has only one constant domain of CH1. Each chain folds into two domains: the variable and the constant domains. The variable domains (VL: residues 1107, VH: residues 1114) have four upper loops (CDRs and another) and three lower loops (FRs), respectively. The constant domains of each chain (CL: residues 108214, CH1: residues 115215) also have three upper (ULs) and three lower (no names) loops, respectively. It has been considered that six CDRs are placed between HEL and Fv for antigen binding. On the other hand, while six FRs and ULs in HELFab are placed between Fv and the constant domain of CLCH1, these functions are unclear.
|
| Results |
|---|
|
|
|---|
atoms were calculated. VL, VH, and HEL in HELFv show RMSDs of 0.62, 0.63, and 0.47 Å, respectively, for those in HELFab. These results suggest that the overall structure of HELFv is similar to that of HELFab. Interestingly, superimposition of VH around the gravity center of the self-domain and similar superimposition of VL gave RMSDs of 0.55 and 1.03 Å for VH, respectively. The difference of the latter RMSD may be due to the considerable conformational changes in the VH domain. In addition, 18 water molecules are found around the HELFv interface of the HELFv complex (1C08 [PDB] ), while only 1 is found around the interface of HELFab (3HFM [PDB] ).
Fluctuation changes of HELFab in the HyHEL-10 system due to Fab or HEL removal
Regarding all frequency modes, Figure 2
shows the fluctuation profiles of C
atoms on Fv and HEL in HyHEL-10 complexes. On HEL removal from the complex, the fluctuation changes appeared mainly around all CDRs and another upper loop, whose residue number is around 70, in each chain as shown in Figure 2A
. On the other hand, on Fab removal from the complex, the fluctuation changes appeared mainly around epitopes in HEL, as shown in Figure 2B
. The above results indicated that the fluctuation changes observed in this study were in the interacting regions between HEL and Fab (Padlan 1977; Padlan et al. 1989; Kondo et al. 1999). Namely, our results show that NMA calculation without water molecules can provide the proteinprotein interaction approximately. Hence, NMA may be useful for illustration of proteinprotein interactions between the antigen and the antibody, and further suppression of atomic vibration may express stronger proteinprotein interaction.
|
atom fluctuation profiles of Fv (Fig. 3A
|
|
|
| Discussion |
|---|
|
|
|---|
Antigenantibody binding is mainly thought to be a result of noncovalent interactions. It is compared to the relationship between "a key and a keyhole," and results in the almost complete desolvation of the interfaces in which hydrophobic clusters, hydrogen bonds, electrostatic interaction, and van der Waals interactions may contribute to the antigenantibody induced fit (Chothia and Janin 1975; Bhat et al. 1994). An antibody recognizes the antigenic surface with the variable domains, and a favorable interface formation may be necessary for the antigenantibody binding.
The parts of the schematic structure of HELFv, which are influenced by the CLCH1 domains, are illustrated in Figure 6
.
|
The RMSD results indicated that the three-dimensional structure of HEL-Fab is similar to that of HELFv, although the domain superimpositions are different. Namely, there are tiny but significant conformational changes in the HELFv complex. Considering the fluctuation differences around interfaces between HEL and Fv domains, the removal of constant domains from the HELFab complex may lead to the insertion of water molecules because of the imperfect complementarity.
It was reported that eight water molecules may mediate the stabilization of interaction in the hydrogen bonds formation, and the buried surface area increase was 100 Å2 at the interface in the HELFv complex in comparison with HELFab (Kondo et al. 1999). The increment of the buried surface area means that another induced fit may occur at the antigenantibody interface compared with the HELFab complex, and such exceptional conformational changes may bring about imperfect complementarity.
It is suggested that CLCH1 may play important roles in the perfect antigenantibody complementarity without some water molecules and the large suppression of induced local conformational changes at the antigenantibody interface.
Kabat et al. compared the sequences of the constant domains of human and mouse immunoglobulin light chains, and predicted that the conserved residues of 160175 may have unique functions (Kabat et al. 1975). The residues that Kabat et al. focused on correspond to UL2CL in this study. UL2CL shows considerable fluctuation changes on HELFv removal, and it is naturally considered that UL2CL is the most effective region in the CLCH1 domains at the interface between HEL and Fv. Therefore, UL2CL may play an important role for longitudinal HELFv interactions.
In the HELHyHEL-10 system, we can reach the following conclusions.
First, fluctuation changes may explain antigenantibody interactions, and the noncovalent interactions can be interpreted with the vibrational features. This is useful for illustrating the local interactions and distant interactions due to tiny dynamic changes.
Second, the constant domains of CLCH1 may have important functions not only for supporting the variable domains as a framework, but also for conserving the favorable interface between HEL and Fv. Incomplete desolvation around the binding interfaces and excess induced local conformational changes may bring about affinity loss of the proteinprotein interaction, and it may be compensated for by the energetically insufficient water-mediated hydrogen bonds formation.
Last, the second upper loop of the CL domain could play an important role in keeping the favorable binding interfaces as distant interactions.
These results have further implications. Regarding protein engineering, our results suggest that the antigenantibody affinity could be controlled by CLCH1 mutations without changing the antigenantibody binding specificity, and, UL2CL, especially, could be a most attractive region.
HELFab and HELFv complexes have an identical sequence of the Fv region and the RMSD around 0.6 with comparison between HELFv moieties: The binding specificity may be maintained. On the other hand, it was reported that HELFv had one order magnitude lower dissociation constant than HELFab. In the other words, the affinity of HELFv may differ without changing binding specificity, and could be controlled by constant domain mutation.
| Materials and methods |
|---|
|
|
|---|
First, energy minimization of the X-ray coordinates was performed with a slightly modified force field of AMBER, and NMA was done with torsion angles (Kamiya et al. 2003). We assumed that the complexes were in vacuo, and that the distance-dependent dielectric constant (r/Å) for electrostatic energy was maintained. The electrostatic potential and van del Waals potential were cut off at 9 Å and were smoothly switched to zero at 10 Å. A threshold of 0.04 kcal/mole Å for the maximal component of the gradients of atoms was used. Six energetically optimized structures on each were obtained under restriction conditions, and these were gradually relaxed near the experimental coordinates. The thermal fluctuations of the atoms were calculated under 300 K (Sumikawa et al. 1998; Soejima et al. 1999; Takeda-Shitaka et al. 1999). Eckarts condition was applied to obtain the fluctuation of each structure (Eckart 1935). The significance level was calculated on couples of six optimized conformations by Wilcoxons rank sum test (Nojima et al. 2002). Refinements of the protein coordinates and applying the statistics may reduce the unexpected effects such as the difference of the X-ray resolution as much as possible.
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Bhat, T.N., Bentley, G.A., Boulot, G., Greene, M.I., Tello, D., DallAcqua, W., Souchon, H., Schwarz, F.P., Mariuzza, R.A., and Poljak, R.J. 1994. Bound water molecules and conformational stabilization help mediate an antigen-antibody association. Proc. Natl. Acad. Sci. 91: 10891093.
Brooks, B. and Karplus, M. 1983. Harmonic dynamics of proteins: Normal modes and fluctuations in bovine pancreatic trypsin inhibitor. Proc. Natl. Acad. Sci. 80: 65716575.
Chothia, C. and Janin, J. 1975. Principles of proteinprotein recognition. Nature 256: 705708.[CrossRef][Medline]
Eckart, C. 1935. Some studies concerning rotating axes and polyatomic molecules. Phys. Rev. 47: 552558.[CrossRef]
Fishmann, T.O., Bentley, G.A., Bhat, T.N., Boulot, G., Mariuzza, R.A., Phillips, S.E., Tello, D., and Poljak, R.J. 1991. Crystallographic refinement of the three-dimensional structure of the FabD1.3lysozyme complex at 2.5 Å resolution. J. Biol. Chem. 266: 1291512920.
Go, N., Noguti, T., and Nishikawa, T. 1983. Dynamics of a small globular protein in terms of low-frequency vibrational modes. Proc. Natl. Acad. Sci. 80: 36963700.
Hayward, S., Kitao, A., and Go, N. 1995. Harmonicity and anharmonicity in protein dynamics: A normal mode analysis and principal component analysis. Proteins Struct. Funct. Genet. 23: 177186.[CrossRef][Medline]
Ishida, H., Jochi, Y., and Kidera, A. 1998. Dynamic structure of Subtilisin-Eglin c complex studied by normal mode analysis. Proteins Struct. Funct. Genet. 32: 324333.[CrossRef][Medline]
Jaaskelainen, S., Verma, C.S., Hubbard, R.E., Linko, P., and Caves, L.S.D. 1998. Conformational change in the activation of lipase: An analysis in terms of low-frequency normal modes. Protein Sci. 7: 13591367.[Abstract]
Kabat, E.A., Padlan, E.A., and Davies, D.R. 1975. Evolutionary and structural influences on light chain constant (CL) region of human and mouse immunoglobulins. Proc. Natl. Acad. Sci. 72: 27852788.
Kamiya, K., Sugawara, Y., and Umeyama, H. 2003. Algorithm for normal mode analysis with general internal coordinates. J. Comput. Chem. 24: 826841.[CrossRef][Medline]
Kam-Morgan, L.N.W., Smith-Gill, S.J., Taylor, M.G., Zhang, L., Wilson, A.C., and Kirsch, J.F. 1993. High-resolution mapping of the HyHEL-10 epitope of chicken lysozyme by site-directed mutagenesis. Proc. Natl. Acad. Sci. 90: 39583962.
Kondo, H., Shiroishi, M., Matsushima, M., Tsumoto, K., and Kumagai, I. 1999. Crystal structure of anti-hen egg white lysozyme antibody (HyHEL-10) Fv-antigen complex. J. Biol. Chem. 274: 2762327631.
Lavoie, T.B., Drohan, W.N., and Smith-Gill, S.J. 1992. Experimental analysis by site-directed mutagenesis of somatic mutation effects on affinity and fine specificity in antibodies specific for lysozyme. J. Immunol. 148: 503513.[Abstract]
Levitt, M., Sander, C., and Stern, P.S. 1985. Protein normal mode dynamics: Trypsin inhibitor, crambin, ribonuclease and lysozyme. J. Mol. Biol. 181: 423447.[CrossRef][Medline]
Miller, D.W. and Agard, D.A. 1999. Enzyme specificity under dynamic control: A normal mode analysis of
-lytic protease. J. Mol. Biol. 286: 267278.[CrossRef][Medline]
Nishikawa, T. and Go, N. 1987. Normal modes of vibration in bovine pancreatic trypsin inhibitor and its mechanical property. Proteins 2: 308329.[CrossRef][Medline]
Nojima, H., Takeda-Shitaka, M., Kurihara, Y., Adachi, M., Yoneda, S., Kamiya, K., and Umeyama, H. 2002. Dynamic characteristics of a peptide-binding groove of human HLA-A2 class 1 MHC molecules: Normal mode analysis of the antigen peptide-class 1 MHC complex. Chem. Pharm. Bull. 50: 12091214.
Padlan, E.A. 1977. Structural implications of sequence variability in immunoglobulins. Proc. Natl. Acad. Sci. 74: 25512555.
Padlan, E.A., Silverton, E.W., Sheriff, S., Cohen, G.H., Smith-Gill, S.J., and Davies, D.R. 1989. Structure of an antibodyantigen complex: Crystal structure of HyHEL-10 Fablysozyme complex. Proc. Natl. Acad. Sci. 86: 59385942.
Smith-Gill, S.J., Lavoie, T.B., and Mainhart, C.R. 1984. Antigenic regions defined by monoclonal antibodies correspond to structural domains of avian lysozyme. J. Immunol. 133: 384393.[Abstract]
Soejima, K., Kurihara, Y., Kamiya, K., and Umeyama, H. 1999. Dynamic character of the complex of human blood coagulation factor VIIa with the extracellular domain of human tissue factor: A normal mode analysis. FEBS Lett. 463: 1923.[CrossRef][Medline]
Sumikawa, H., Suzuki, E., Fukuhara, K., Nakajima, Y., Kamiya, K., and Umeyama, H. 1998. Dynamic structures of granulocyte colony-stimulating factor proteins studied by normal mode analysis: Two domain-type motions in low frequency modes. Chem. Pharm. Bull. 46: 10691077.
Sundberg, E.J., Urrutia, M., Braden, B.C., Isern, J., Tsuchiya, D., Fields, B.A., Malchiodi, E.L., Tormo, J., Schwarz, F.P., and Mariuzza, R.A. 2000. Estimation of hydrophobic effect in an antigen-antibody proteinprotein interface. Biochemistry 39: 1537515387.[CrossRef][Medline]
Takeda-Shitaka, M., Kamiya, K., Miyata, T., Ohkura, N., Madoiwa, S., Sakata, Y., and Umeyama, H. 1999. Structural studies of the interactions of normal and abnormal human plasmins with bovine basic pancreatic trypsin inhibitor. Chem. Pharm. Bull. 47: 322328.
Tsumoto, K., Ueda, Y., Maenaka, K., Watanabe, K., Ogasahara, K., Yutani, K., and Kumagai, I. 1994. Contribution to antibodyantigen interaction of structurally perturbed antigenic residues upon antibody binding. J. Biol. Chem. 269: 2877728782.
Tsumoto, K., Ogasahara, K., Ueda, Y., Watanabe, K., Yutani, K., and Kumagai, I. 1996. Role of salt bridge formation in antigenantibody interaction. J. Biol. Chem. 271: 3261232616.
Ward, E.S., Gussow, D., Griffiths, A.D., Jones, P.T., and Winter, G. 1989. Binding activities of a repertoire of single immunoglobulin variable domains secreted from Escherichia coli. Nature 341: 544546.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |