Teeter, M. M., Froimowitz, M., Stec, B. and. DuRand, C. J. (1994) J. Med. Chem. 37,2874-2888. Kyle, D. J., Chakravarty, S., Sinsko, J. A. and. Stormann, T. M. ...
Computer Modelling of Biological Molecules
3448-3462 3 Teeter, M. M., Froimowitz, M., Stec, B. and DuRand, C. J. (1994)J. Med. Chem. 37,2874-2888 4 Kyle, D. J., Chakravarty, S., Sinsko, J. A. and Stormann, T. M. (1994) J. Med. Chem. 37, 1347- 1345 5 Chakravarty, S., Connolly, M. and Kyle, D. J. (1995)Peptide Res. 8, 16-19 6 Zhang, D. Q. and Weinstein, H. (1993)J. Med. Chem. 36,934-938 7 Luo, X. C., Zhang, D. Q. and Weinstein, H. (1994) Protein Eng. 7, 1441-1448 8 Henderson, R., Baldwin, J. M., Ceska, T. A.,
9 10
11 12
13 14
Zemlin, F., Beckmann, E. and Downing, K. H. (1990)J. Mol. Biol. 213,899-929 Schertler, G. F. X., Villa, C. and Henderson, R. (1993)Nature (London) 362,770-772 Unger, V. M. and Schertler, G. F. X. (1995) Biophys. J. 68,1776-1786 Baldwin, J. M. (1993)EMBO J. 12,1693-1703 Strader, C. D., Fong, T. M., Tota, M. R. and Underwood, D. (1994) Annu. Rev. Biochem. 63, 101- 132 Fujinaga, M. and James, M. N. G. (1987)J. Mol. Biol. 195,373-396 Chen, Z. and Bode, W. (1983)J. Mol. Biol. 164, 283-3 1 1
15 Oliveira, L.,Paiva, A. C. M. and Vriend, G. (1993) J. Computer-Aided Mol. Design 7,649-658 16 Goodford, P. J. (1985)J. Med. Chem. 28,849-857 17 Weiner, S. J., Kollman, P. A., Case, D. A., Singh, U. C., Ghio, C., Alagona, G., Profeta, S., Jr, and Weiner, P. (1984)J. Am. Chem. SOC.106,765-784 18 Weiner, S. J., Kollman, P. A., Nguyen, D. T. and Case, D. A. (1986)J. Comput. Chem. 7,230-252 19 Singh, U. C., Weiner, P. K., Caldwell, J. W. and Kollman, P. A. (1988) AMBER, Version 4.0, Department of Pharmaceutical Chemistry, University of California, San Francisco 20 Gouldson, P. R., Winn, P. J. and Reynolds, C. A. (1995)J. Med. Chem. 38,4080-4086 21 Ferenczy, G. G., Reynolds, C. A. and Richards, W. G. (1990)J. Comput. Chem. 11, 159-169 22 Reynolds, C. A., Ferenczy, G. G. and Richards, W. G. (1992)J. Mol. Struct. THEOCHEM. 256, 249-269 23 Weinstein, H.,Mazurek, A. P., Osman, R. and Topiol, S. (1986)Mol. Pharmacol. 29,28-33 24 Alla, S.A., Buschko, J., Quitterer, U., Maidhof, A., Haasemann, M., Breipohl, G., Knolle, J. and Mulleresterl, W. (1993) J. Biol. Chem. 268, 17277- 17285 Received 7 August 1995
Quantum mechanical/molecular mechanical approaches to transition state structure: mechanism of sialidase action J. A. Barnes and I. H. Williams* School of Chemistry, University of Bath, Bath BA2 7AY, U.K.
Introduction Glycoside hydrolyses are very important biochemical processes [ 11 whose mechanisms pose a number of questions. Are they stepwise or concerted? Do they occur with retention or inversion of configuration at the anomeric centre? T h e s e issues demand a knowledge of transition state (TS) structure [2] which, from an experimental point of view, is best probed by means of kinetic isotope effects (KIEs). Primary KIEs indicate which atoms are directly involved in bond making or breaking in the rate-determining TS. Abbreviations used: TS, transition state; KIE, kinetic isotope effect; MM, molecular mechanics; QM, quantum mechanical; MO, molecular orbital; DANA, 2-deoxy-2,3-dehydro-N-acetylneuraminicacid; RMS; root mean square; FR, influenza with retention simulation; SR, S. typhimurzum with retention simulation; SI, S. typhzmurium with inversion simulation; MD, molecular dynamics. *To whom correspondence should be addressed.
Secondary KIEs, involving isotopic substitution at positions not directly involved in bond making or breaking, are often used to provide a measure of the location of the TS along the reaction coordinate between reactants and products. T h e most convincing examples of KIEs as probes of TS structure are those in which the effects of multiple isotopic substitutions are examined. From the viewpoint of theoretical modelling, questions of reactivity and TS structure require the use of a method capable of describing the idiosyncratic behaviour of electrons in making and breaking bonds. Molecular mechanics (MM) is inherently unsuited for this purpose, being an interpolation scheme for known properties of stable species. Quantum mechanical (QM) approaches are best able to provide the unknown information regarding TSs by means of molecular orbital (MO) theory. For large systems hybrid QM/MM methods may be applied, combining the merits of a QM treatment of bond making/break-
I996
263
Biochemical Society Transactions
264
ing with a MM description of the environment within which the chemical events occur, be it solvent or protein. In this paper we report use of the CHARMM program [31 to perform hybrid QM/MM calculations [4] upon enzyme-catalysed glycoside hydrolyses, combining an AM1 semiempirical MO theoretical description of the reacting system with a CHARMM22 empirical force field description of the surrounding protein and solvent molecules. We have successfully used a similar approach to study the (non-enzymic) acid-catalysed hydrolysis of a model for AMP embedded within a sphere of approx. 500 water molecules [5]; the reacting system was treated by AM1 and the solvent by the TIP3P empirical potential. The calculated energy surface for reaction suggested the possibility of concurrent stepwise and concerted mechanisms. KIEs may be computed by consideration of the shifts in vibrational frequencies resulting from isotopic substitution at various atoms of the reacting system, and comparison with experimental values provides a stern test of any TS structure predicted by theoretical modelling. We have modified CHARMM to compute second derivatives of the energy with respect to small displacements of the co-ordinates of the QM atoms within their MM environment, from which vibrational frequencies and KIEs may be obtained [6]. Four separate KIEs calculated for the stepwise mechanism on our QM/MM energy surface agreed within reported error of the experimentally measured values [7]. Computational modelling thus offers a tool, complementary with experiment, capable of providing valuable insight into complex reaction mechanisms. Sialidases (or neuraminidases, EC 3.2.1.18) are a family of glycohydrolases, found in many organisms, which catalyse the cleavage of terminal sialic acids (N-acetylneuraminic acid, NANA) which are linked a-ketosidically to glycoconjugates [8]. Only the a-anomers of sialic acid derivatives are substrates for these enzymes, and all known sialidases (with the possible exception of that from Salmonella typhimurium [9]) act with retention of configuration at the anomeric carbon. Many biological roles have been attributed to sialidases, including cell-cell recognition processes and the pathogenicity of some infections by some sialidase-bearing microorganisms [lo]. A sialidase is considered important for transporting the influenza virus through mucin [ 113 and for the elution of new virus from
Volume 24
infected cells [ 121. Sialidase inhibitors are therefore attractive targets for potential anti-influenza drugs [13]. Detailed knowledge of the mechanisms of sialidase action in enzymes from different organisms would be highly desirable. The three-dimensional structures of sialidases from human influenza B virus [14] and S. typhimurium [15] show considerable similarity in the general folding and spatial arrangement of catalytic residues, which are highly conserved, despite a sequence homology of only 15%. The main features of the mechanism proposed [16,17] for influenza sialidase are binding of the substrate in the (normally unfavourable) B2,5boat conformation, protonation of the aglycone by an unknown general acid followed by departure of the leaving group, stabilization of a sialosyl zwitterionic intermediate in the active site and general base catalysis of attack by water on the intermediate, yielding the sialic acid product. However, the sialidase from S. typhimurium is suggested to act by a single displacement mechanism involving a TS derived from the 'CS chair conformation [ 181.
MM modelling of sialidases T o X-ray crystallographic co-ordinates for sialidases from influenza B virus at 1.8 A resolution [ 141 and from S. typhimurium at 2.2 A resolution [15], each as bound inhibitor complexes with 2-deoxy-2,3-dehydro-N-acetylneuraminic acid (DANA), hydrogen atoms were added using the standard ionization state of each residue in each amino acid sequence. Each structure was subjected to 100 steps of steepest-descent MM energy minimization to remove bad contacts, followed by geometry optimization with the Powell algorithm to root mean square (RMS) gradient