Evolution of Shape Complementarity and Catalytic Efficiency from a Primordial Antibody Template
10.1126/science.286.5448.2345
Crossref journal-article
American Association for the Advancement of Science (AAAS)
Science (221)
Abstract

The crystal structure of an efficient Diels-Alder antibody catalyst at 1.9 angstrom resolution reveals almost perfect shape complementarity with its transition state analog. Comparison with highly related progesterone and Diels-Alderase antibodies that arose from the same primordial germ line template shows the relatively subtle mutational steps that were able to evolve both structural complementarity and catalytic efficiency.

Bibliography

Xu, J., Deng, Q., Chen, J., Houk, K. N., Bartek, J., Hilvert, D., & Wilson, I. A. (1999). Evolution of Shape Complementarity and Catalytic Efficiency from a Primordial Antibody Template. Science, 286(5448), 2345–2348.

Authors 7
  1. Jian Xu (first)
  2. Qiaolin Deng (additional)
  3. Jiangang Chen (additional)
  4. Kendall N. Houk (additional)
  5. Johannes Bartek (additional)
  6. Donald Hilvert (additional)
  7. Ian A. Wilson (additional)
References 55 Referenced 101
  1. Pauling L., Am. Sci. 36, 51 (1948); / Am. Sci. by Pauling L. (1948)
  2. ; W. P. Jencks Catalysis in Chemistry and Enzymology (McGraw-Hill New York 1969).
  3. 10.1126/science.7569920
  4. ; Acc. Chem. Res. 26 391 (1993); (10.1021/ar00032a001)
  5. Hilvert D., Top. Stereochem. 22, 83 (1999). (10.1002/9780470147313.ch2) / Top. Stereochem. by Hilvert D. (1999)
  6. Laschat S., Angew. Chem. Int. Ed. Engl. 35, 289 (1996). (10.1002/anie.199602891) / Angew. Chem. Int. Ed. Engl. by Laschat S. (1996)
  7. Haynes M. R., Stura E. A., Hilvert D., Wilson I. A., Science 263, 646 (1994); (10.1126/science.8303271) / Science by Haynes M. R. (1994)
  8. ; P. A. Patten et al. Science 271 1086 (1996); (10.1126/science.271.5252.1086)
  9. 10.1126/science.276.5319.1665
  10. ; J. B. Charbonnier et al. Science 275 1140 (1997); H. D. Ulrich et al. Nature 389 271 (1997). (10.1038/38470)
  11. Hilvert D., Hill K. W., Nared K. D., Auditor M.-T. M., J. Am. Chem. Soc. 111, 9261 (1989). (10.1021/ja00208a037) / J. Am. Chem. Soc. by Hilvert D. (1989)
  12. Quantum mechanical calculations at the RHF/6-31G(d) level were performed to optimize the geometries of reactants intermediates and transition states; energies were evaluated at the B3LYP/6-31G(d) level. Stationary points were characterized by RHF/3-21G frequency calculations (M. J. Frisch et al. Gaussian Gaussian Inc. Pittsburgh PA). The calculations show transition state 5 to be rate determining. The computed activation energy is 19.9 kcal/mol and increases to 22.1 kcal/mol in water (SCI-PCM) [M. Cossi et al. Chem. Phys. Lett. 255 327 (1996)]. A single hydrogen bond from a water molecule to NEM 2 lowers the activation energy by 0.6 kcal/mol. Calculations predict no effect on hydrogen bonding of water to the thiophene dioxide.
  13. Braisted A. C., Schultz P. G., J. Am. Chem. Soc. 112, 7430 (1990). (10.1021/ja00176a073) / J. Am. Chem. Soc. by Braisted A. C. (1990)
  14. V. E. Gouverneur et al. Science 275 1140 (1997)
  15. A. A. P. Meekel M. Resmini U. K. Pandit J. Chem. Soc. Chem. Commun. 571 (1995). (10.1039/C39950000571)
  16. J. T. Yli-Kauhaluoma et al. J. Am. Chem. Soc. 117 7041 (1995). (10.1021/ja00132a001)
  17. Pitt A. R., Stimson W. H., Suckling C. J., Marrero-Tellado J. J., Vazzana C., Israel J. Chem. 136, 171 (1996). (10.1002/ijch.199600023) / Israel J. Chem. by Pitt A. R. (1996)
  18. Seelig B., Jäschke A., Chem. Biol. 6, 167 (1999); (10.1016/S1074-5521(99)89008-5) / Chem. Biol. by Seelig B. (1999)
  19. Tarasow T. M., Tarasow S. L., Eaton B. E., Nature 389, 54 (1997). (10.1038/37950) / Nature by Tarasow T. M. (1997)
  20. Stura E. A., Fieser G. G., Wilson I. A., Immunomethods 3, 164 (1993); (10.1006/immu.1993.1051) / Immunomethods by Stura E. A. (1993)
  21. . The full length IgG1 was expressed from mouse myeloma cells and purified as described in (5). The 1E9 Fab fragment was generated by enzymatic digestion with 2% papain in the presence of cysteine.
  22. Haynes M. R., Lenz M., Taussig M. J., Wilson I. A., Hilvert D., Israel J. Chem. 36, 151 (1996). (10.1002/ijch.199600021) / Israel J. Chem. by Haynes M. R. (1996)
  23. Arevalo J. H., Taussig M. J., Wilson I. A., Nature 365, 859 (1993); (10.1038/365859a0) / Nature by Arevalo J. H. (1993)
  24. ; J. H. Arevalo et al. J. Mol. Biol. 231 103 (1993); J. H. Arevalo et al. J. Mol. Biol. 241 663 (1994).
  25. 10.1126/science.279.5358.1929
  26. Jones S., Thornton J. M., Proc. Natl. Acad. Sci. U.S.A. 93, 1320 (1996). (10.1073/pnas.93.3.1320) / Proc. Natl. Acad. Sci. U.S.A. by Jones S. (1996)
  27. G. M. Morris et al. J. Comp. Chem. 19 1639 (1998) (10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-B)
  28. G. M. Morris et al. J. Comput. Aided Mol. Des. 10 293 (1996); (10.1007/BF00124499)
  29. Goodsell D. S., Olson A. J., Proteins Struct. Funct. Genet. 8, 195 (1990). (10.1002/prot.340080302) / Proteins Struct. Funct. Genet. by Goodsell D. S. (1990)
  30. E. A. Kabat T. T. Wu H. M. Perry K. S. Gottesman C. Foeller Sequences of Proteins of Immunological Interest (U.S. Public Health Service National Institutes of Health Bethesda MD ed. 5 1991). See also the antibody database Web page
  31. The DJ regions of the V H gene of 1E9 come from its germ line precursor DFL 16.1 and J H 4 respectively with very common guanine-rich joints flanking both ends of the D gene whereas those of DB3 most likely come from SP25.8 and J H 1 with short VD and DJ junctions. For 39-A11 the DJ regions are most likely from DSP2.2 and J H 2 but with an uncommon nine-nucleotide VD junction. For 1E9 the use of the J H 4 minigene places a Met at position H100b providing a hydrophobic floor for the binding pocket that is highly adaptable to the ligand shape because of the flexibility of its side chain; DB3 and 39-A11 utilize J H 1 and J H 3 respectively resulting in a Phe at this position. Residue H100 which provides one of the walls of the cavity arises from the DJ junction; it is an Arg in 1E9 and 39-A11 and a Trp in DB3. The H97 residue (Thr in 1E9 Tyr in DB3 and Arg in 39-A11) interacts primarily with the linker moiety in the 1E9 and 39-A11 structures whereas in DB3 it provides a lid over one end of the hydrophobic cavity.
  32. Morea V., Tramontano A., Rustici M., Chothia C., Lesk A. M., J. Mol. Biol. 275, 269 (1998). (10.1006/jmbi.1997.1442) / J. Mol. Biol. by Morea V. (1998)
  33. Synthesis of 7 : 3 4 5-trichloro-2-thienylmethanol was reacted with ethyl iodoacetate under basic solid-liquid phase-transfer conditions to afford ethyl (3 4 5-trichlorothiophen-2-ylmethoxy) acetate in 69% yield based on recovered alcohol. The ester was quantitatively saponified and the resulting acid was oxidized with HOF-CH 3 CN to give 7 (73%). All new compounds gave satisfactory spectroscopic data.
  34. Reactions were carried out in 50 mM acetate buffer (pH 5.5) containing 20 mM NaCl and were monitored either by high-pressure liquid chromatography or by starch-I 2 bleaching at 606 nm. Initial velocities were determined at several concentrations of one substrate while the concentration of the other was kept constant. The uncatalyzed reaction was carried out under identical conditions and also in acetonitrile.
  35. Westheimer F. H., Adv. Enzymol. 24, 441 (1962). / Adv. Enzymol. by Westheimer F. H. (1962)
  36. 10.1073/pnas.68.8.1678
  37. Lewis C., Krämer T., Robinson S., Hilvert D., Science 253, 1019 (1991); (10.1126/science.1887215) / Science by Lewis C. (1991)
  38. Thorn S. N., Daniels R. G., Auditor M.-T. M., Hilvert D., Nature 373, 228 (1995); (10.1038/373228a0) / Nature by Thorn S. N. (1995)
  39. Hollenfelder F., Kirby A. J., Tawfik D., Nature 383, 60 (1996); (10.1038/383060a0) / Nature by Hollenfelder F. (1996)
  40. Kikuchi K., Hilvert D., J. Am. Chem. Soc. 118, 8184 (1996). (10.1021/ja9617395) / J. Am. Chem. Soc. by Kikuchi K. (1996)
  41. Warshel A., Florian J., Proc. Natl. Acad. Sci. U.S.A. 95, 5950 (1998). (10.1073/pnas.95.11.5950) / Proc. Natl. Acad. Sci. U.S.A. by Warshel A. (1998)
  42. Molecular dynamics simulations with the LIE method [
  43. Åqvist J., Medina C., Samuelsson J. E., Protein Eng. 7, 385 (1994); (10.1093/protein/7.3.385) / Protein Eng. by Åqvist J. (1994)
  44. ] derived from linear response theory were used to compute and analyze binding of reactants and transition states to the antibody. The free energy of binding is calculated as Δ G bind = 1/2 Δ V elec + αΔ V vdw where Δ V elec and Δ V vdw denote the difference of ensemble-averaged electrostatic and van der Waals interaction energies between the antibody-bound and aqueous-solvated states. The parameter α was obtained from fitting to known binding affinities of steroids to DB3 (14). Molecular dynamics simulations were computed with a modified version of AMBER4.1 (D. A. Pearlman et al. AMBER4.1 University of California San Francisco 1995). The free energy of binding of 1 plus 2 into a complex in the antibody is computed as −3.2 kcal/mol compared with experimental values of −5.2 kcal/mol from K m values. The free energy of binding of transition state 5 is −8.9 kcal/mol. Thus the predicted activation energy of the cycloaddition is lowered 5.7 kcal/mol relative to that in water (6).
  45. Romesberg F. E., Schultz P. G., Bioorg. Med. Chem. Lett. 9, 1741 (1999). (10.1016/S0960-894X(99)00281-4) / Bioorg. Med. Chem. Lett. by Romesberg F. E. (1999)
  46. A. Heine et al. Science 279 1934 (1998). (10.1126/science.279.5358.1934)
  47. 10.1016/S0076-6879(97)76066-X
  48. 10.1107/S0108767393007597
  49. A. T. Brünger X-PLOR Version 3.1: A System for X-Ray and NMR (Yale Univ. Press New Haven CT 1992)
  50. J. Badger et al. Proteins Struct. Funct. Genet. 35 25 (1999). (10.1002/(SICI)1097-0134(19990401)35:1<25::AID-PROT3>3.3.CO;2-M)
  51. 10.1107/S0021889883010985
  52. The molecular surfaces of the antibody antigen and antibody-antigen complex were constructed separately with the program GRASP [
  53. 10.1002/prot.340110407
  54. ] with a probe radius of 1.4 Å and the volume enclosed in each surface was then computed. The gap volume is taken as the difference between the volume of the complex and the sum of the volumes of the antibody and the antigen. Gap index is the gap volume divided by the buried surface area.
  55. Supported by NIH grants CA27489 (I.A.W.) and GM38273 (D.H.) and NSF grant CHE-9616772 (K.N.H.). We thank T. Horton N. Jiang and M. Auditor for technical assistance; the staff of SSRL beamline 9-1; R. A. Lerner A. Heine and D. J. Edwards for helpful discussions; the referees for helpful comments; N. R. Klinman for advice on germ line gene analysis; and K. Renner for computational assistance.
Dates
Type When
Created 23 years, 1 month ago (July 27, 2002, 5:37 a.m.)
Deposited 1 year, 7 months ago (Jan. 13, 2024, 5:34 a.m.)
Indexed 1 month, 3 weeks ago (July 2, 2025, 1:29 p.m.)
Issued 25 years, 8 months ago (Dec. 17, 1999)
Published 25 years, 8 months ago (Dec. 17, 1999)
Published Print 25 years, 8 months ago (Dec. 17, 1999)
Funders 0

None

@article{Xu_1999, title={Evolution of Shape Complementarity and Catalytic Efficiency from a Primordial Antibody Template}, volume={286}, ISSN={1095-9203}, url={http://dx.doi.org/10.1126/science.286.5448.2345}, DOI={10.1126/science.286.5448.2345}, number={5448}, journal={Science}, publisher={American Association for the Advancement of Science (AAAS)}, author={Xu, Jian and Deng, Qiaolin and Chen, Jiangang and Houk, Kendall N. and Bartek, Johannes and Hilvert, Donald and Wilson, Ian A.}, year={1999}, month=dec, pages={2345–2348} }