computer programs

Journal of Applied MetLigDB: a web-based database for the Crystallography identification of chemical groups to design ISSN 0021-8898 metalloprotein inhibitors

Received 12 April 2011 Accepted 10 June 2011 Hwanho Choi,‡ Hongsuk Kang‡ and Hwangseo Park*

Department of Bioscience and Biotechnology, Sejong University, 98 Kunja-dong, Kwangjin-ku, Seoul 143-747, Korea. Correspondence e-mail: [email protected]

MetLigDB (http://silver.sejong.ac.kr/MetLigDB) is a publicly accessible web- based database through which the interactions between a variety of chelating groups and various central metal ions in the of metalloproteins can be explored in detail. Additional information can also be retrieved, including protein and inhibitor names, the amino acid residues coordinated to the central metal ion, and the binding affinity of the inhibitor for the target metalloprotein. Although many metalloproteins have been considered promising targets for drug discovery, it is difficult to discover new inhibitors because of the difficulty in designing a suitable chelating moiety to impair the catalytic activity of the central metal ion. Because both common and specific chelating groups can be identified for varying metal ions and the associated coordination environments, # 2011 International Union of Crystallography MetLigDB is expected to give users insight into designing new inhibitors of Printed in Singapore – all rights reserved metalloproteins for drug discovery.

1. Introduction 2008), the coordination geometry around the metal ion (Castagnetto The importance of metal ions in biological systems is increasingly et al., 2002; Andreini, Bertini, Cavallaro, Holliday & Thornton, 2009), drawing attention, as indicated by the recent emergence of terms such catalytic mechanisms (Andreini, Bertini, Cavallaro, Najmanovich & as metallome and metallomics (Shi & Chance, 2008). This is not Thornton, 2009) and interactions in protein structures (Hemavathi et surprising because about 40% of proteins require metal ions in their al., 2010). Therefore, we have constructed MetLigDB to provide active sites for biological activity (Andreini et al., 2008). Throughout further insight when designing new inhibitors of metalloproteins. To this paper we will use the term ‘active site’ for all kinds of proteins, this end, users are enabled to examine the interactions between a although it is normally referred to . Most metalloproteins known inhibitor and the active-site metal ion cluster of a metallo- play a key role in biological processes and therefore have been protein of interest in the Protein Data Bank (PDB; Berman et al., considered to be promising targets for drug discovery. Metal ions in 2000). This allows the identification of the chemical moieties appro- the active site can participate in the enzymatic reaction as a receptor priate for the chelative inactivation of a specific metal ion with a given for lone-pair electrons, which has the effect of weakening the coordination environment in the metalloprotein. The coordination chemical bonds in the (Armstrong, 2000). Metal ions may patterns of a variety of chemical moieties can be found and visualized also play the role of a Lewis acid catalyst, stabilizing the transition in MetLigDB. In order for the relative binding affinities of the inhi- state or an intermediate in which the negative charge is developed as bitors to be estimated, their Ki or IC50 (binding affinity and 50% a consequence of (Kleifeld et al., 2003). The central inhibitory concentration) values against known metalloproteins are metal ion is thus indispensable for the biological activity of metallo- also provided. This paper summarizes the data content and the basic proteins. technical aspects employed in the construction of MetLigDB. In order for a small-molecule inhibitor to regulate the activity of a metalloprotein, the inhibitor must contain a chemical group that can bind to the central metal ion. In the case of virtual screening for 2. Methodology inhibitors, the covalent-like interactions between metal centers and their ligands, coupled with their large electrostatic potentials, make MetLigDB is a publicly accessible web-based database designed by metalloproteins a challenge for many docking scoring functions taking advantage of well known web technologies such as a stable (Irwin et al., 2005). Despite such difficulties, a number of chemical web server (Apache), a fast database engine (MySQL) and a groups that chelate and inactivate metal ions at the active site of powerful web scripting language (PHP). By virtue of these open- metalloproteins have been reported, including carboxylate, phos- source tools, we have been able to construct interactive web query phorate, hydroxamate, -keto acid and diol moieties. interfaces and functions to allow remote viewing and searching. Although some databases and web servers for metalloproteins Structural information about the metalloprotein–inhibitor interac- have already been reported, all the biochemical information has been tions was collected from the PDB, which contains a total of 71 635 limited to the identification of metal-binding residues (Hsin et al., experimentally determined structures made public before 8 March 2011. In this study, the metalloprotein and the inhibitor refer to a protein that has at least one metal ion in its physiological site and a ‡ These authors contributed equally to this work. small molecule that has been shown to impair the biological activity

878 doi:10.1107/S0021889811022503 J. Appl. Cryst. (2011). 44, 878–881 computer programs of a protein, respectively. Starting from these structures, the data Table 1 content of MetLigDB for metalloprotein–inhibitor interaction was Number of metal ions in the simplified model for metalloprotein–inhibitor complexes in MetLigDB. obtained in a few steps of data filtration and simplification, as depicted in Fig. 1. The first step involved the identification of the Metal ion Number of protein–ligand complexes Number of metalloproteins structures in the PDB that contain at least one metal ion and Zn2+ 436 142 simultaneously a small-molecule inhibitor in the active site. Among Mn2+ 188 86 2+ 3+ the 3045 structures that satisfied this criterion, those in which the Fe /Fe 61 28 Ni2+ 21 8 central metal ion forms at least one coordination bond with the Mg2+ 14 12 inhibitor were selected for further processing. As the distance limit Cu2+ 88 ˚ Co2+ 22 for defining a coordination bond, we used R0 + 0.5 A as suggested by Mo2+ 21 Shi et al. (2007), where R0 represents the average interatomic distance Total 732 287 between the central metal ion and a ligand atom as given in the PDB. The structural data file in PDB format was then prepared to show the detailed coordination pattern by selecting the central metal ion, the With respect to the ligand-binding affinities, data from multiple small-molecule inhibitor, and the amino acid residues and the research groups were collected to obtain the average values of Ki or structural water molecules coordinated to the metal ion. The other IC50 of the inhibitor for the target protein. The journals from which parts of the original PDB file of the protein–inhibitor complex were the binding affinity data were drawn include Journal of Medicinal removed for simplicity and visual clarity. On the basis of this proce- Chemistry, Journal of Biological Chemistry, Bioorganic and Medic- dure, 732 simplified structural models for metalloprotein–inhibitor inal Chemistry Letters, and Biochemistry. Despite an extensive search complexes were obtained for 250 proteins and eight types of metal for binding affinity data, Ki or IC50 values for only 41% of the ion. metalloprotein–inhibitor complexes contained in MetLigDB were available in the literature. This is in contrast to the case of nonmetal enzymes, for which the binding affinity data can be found for most of the inhibitors present in the PDB (Liu et al., 2007). The majority of the binding affinity data in MetLigDB were obtained from inhibition studies, while a small portion of data was obtained from the more informative methods such as isothermal titration calorimetry (Pierce et al., 1999).

3. Contents of database The web page of MetLigDB is organized in a hierarchical fashion. As can be seen in Fig. 2, users can select the kind of metal ion or protein on the menu bar. In MetLigDB, the numbers of metalloproteins that belong to the categories , , , , and are 87, 108, 306, 110, 27 and 8, respectively. There are also 86 metalloproteins in the database whose enzyme Figure 1 The flowchart used for the identification of a metalloprotein–inhibitor complex commission (EC) numbers have not been assigned. Each entry involving a direct metal–ligand coordination in the PDB. contains the PDB code, the protein name, the EC number, the protein classification, the kind of metal ion in the active site, the IUPAC name of the small-molecule inhibitor, the binding affi- nity and the binding mode of the inhibitor to the active-site metal ion cluster. The coordination pattern for a chelating group of the inhibitor with the central metal ion can be displayed with Jmol (http:// www.jmol.org), an open-source Java viewer for chemical structures in three dimensions. Users may be enabled to detect a chemical moiety suitable for the inactivation of the catalytic action of the metalloprotein of interest by visual inspection of the simpli- fied structural model for metalloprotein– inhibitor complexes. Table 1 lists the number of entries in MetLigDB for a given metal ion. The zinc– ligand complexes appear to be predomi- nant, occupying 59.6% of the total data Figure 2 A screenshot of the main window of MetLigDB. In total, seven interface windows are presented to users, each of content. Such an abundance of zinc metallo- which can be shown by clicking on the tabs at the top of the main window. proteins is consistent with the fact that zinc

J. Appl. Cryst. (2011). 44, 878–881 Choi, Kang and Park  MetLigDB 879 computer programs

Table 2 Representative ligand moieties for the central metal ions in metalloproteins. Indicated in red are the atoms coordinated to the metal ion.

Ligand group Metal ion Target protein

Zn2+,Mn2+,Ni2+ Matrix metalloproteinase Histone acetylase Thermolysin Formylmethionine deformylase Peptide deformylase

Zn2+,Mn2+,Fe2+ 6-Pyruvoyltetrahydropterin synthase Glucose-6-phosphate isomerase Xylose isomerase d-Psicose 3-epimerase !-Amidase Nitrobenzene dioxygenase

Zn2+,Mn2+,Fe2+ 2-Isopropylmalate synthase Isocitrate dehydrogenase Pyruvate kinase Fumarylacetoacetase Sialic acid synthase Factor inhibiting hif-1 alpha Leucocyanidin oxygenase

Zn2+ Carbonic anhydrase 2

Mn2+ Methionine aminopeptidase-2 Figure 3 Visualization of information retrieved from selecting an entry in MetLigDB. This example illustrates that the hydroxamate group should be an effective chelator for a zinc ion in the active site of matrix metalloproteinases. Fe2+ Isopenicillin N synthase

value of 3 nM. This indicates that the hydroxamate moiety can be an effective chelator that should be included in the molecular structure Ni2+ 1-Aminocyclopropane-1-carboxylate synthase of matrix metalloproteinase inhibitors. Indeed, a vast majority of matrix metalloproteinase inhibitors are known to possess the Co2+ Trypsin hydroxamate group as the specific zinc-binding moiety (Sang et al., 2006). Some representative chemical moieties suitable for the chelative inhibition of metal ions in the active sites of metalloproteins are Mo2+ Nitrate reductase summarized in Table 2. The full list of chelators is given in the supplementary material1 and is readily available under the ‘Chelating Group’ tab on the MetLigDB web site. First of all, the carboxylate moiety is found to be the most popular chelating group for metal ions (see supplementary material for details). It is observed in the active is the trace metal observed most frequently in biological systems: sites of more than 20 metalloproteins in coordination with various more than 300 enzymes have been identified with the Zn2+ ion in the metal ions including Zn2+,Mn2+,Fe2+,Ni2+ and Cu2+ ions. The catalytic site as a (Auld, 2001; Cox & McLendon, 2000). The phosphorate group also appears to be a common chelating group, proportions of Mn2+,Fe2+/Fe3+,Ni2+,Mg2+ and Cu2+ metalloproteins which can be coordinated to almost all kinds of metal ions in the contained in MetLigDB are 25.7, 8.3, 2.9, 1.9 and 1.1%, respectively. active site of metalloproteins. Because they have such a wide range of There are only two structures of Co2+ and Mo2+ metalloproteins in inhibitory activity, the carboxylate and phosphorate groups seem to complex with an inhibitor with a chelating group. The rarity of such be inadequate as a chemical group in a small-molecule drug because structures is not surprising because the coordination spheres of Co2+ nonspecific inhibition of various metalloproteins could cause serious and Mo2+ ions can be formed in metalloproteins only under specific side effects. In the case of Mg2+ metalloproteins, only the phos- conditions (Williams, 2003). MetLigDB is thus likely to provide much phorate ion and its analogs are found in the active site, which indi- information on identifying a suitable chelating group for zinc or cates a difficulty in designing a specific inhibitor of Mg2+ manganese metalloproteins, whereas only limited information may be metalloproteins. available for cobalt and molybdenum metalloproteins. Hydroxamate, vicinal diol and -keto acid moieties are also found Fig. 3 shows an example of the results of a data search for a suitable to be effective chelating groups that can be coordinated to at least chelating group for the Zn2+ ion in the active site of matrix met- three kinds of metal ions in the active site of metalloproteins. Besides alloproteinases. In this case, the two O atoms of the hydroxamate the common ligands, some specific chelating groups are also identified group appear to be coordinated to the zinc ion to form a trigonal- bipyramidal coordination geometry with three histidine residues at 1 Supplementary material discussed in this paper is available from the IUCr the bottom of the active site. This strong coordination of the electronic archives (Reference: WF5093). Services for accessing this material hydroxamate group makes the inhibitor very potent, with an IC50 are described at the back of the journal.

880 Choi, Kang and Park  MetLigDB J. Appl. Cryst. (2011). 44, 878–881 computer programs for various metal ions. For example, sulfonamide, 3-amino-2-hydroxy- This work was supported by grant No. K11061 from the Korea propionamide, thiazolidin-5-ylmethanol, 2-aminoethnaol, bis(1H- Institute of Oriental Medicine. imidazol-2-yl)methanone and vicinal dithiol groups appear to play 2+ 2+ 2+ 2+ 2+ 2+ the role of specific chelator for Zn ,Mn ,Fe ,Ni ,Co and Mo References ions, respectively. Therefore, these chemical groups can be consid- Andreini, C., Bertini, I., Cavallaro, G., Holliday, G. L. & Thornton, J. M. ered as key structural components in designing specific inhibitors of (2008). J. Biol. Inorg. Chem. 13, 1205–1208. metalloproteins owing to the selectivity in the chelative inactivation Andreini, C., Bertini, I., Cavallaro, G., Holliday, G. L. & Thornton, J. M. of the metal ion in the active site. MetLigDB can thus make it (2009). Bioinformatics, 25, 2088–2089. possible for a user to select a suitable chelating group that should be Andreini, C., Bertini, I., Cavallaro, G., Najmanovich, R. J. & Thornton, J. M. (2009). J. Mol. Biol. 388, 356–380. included in the inhibitor of a metalloprotein of interest, which is a Armstrong, R. N. (2000). Biochemistry, 39, 13625–13632. starting point in designing new metalloprotein inhibitors. Auld, D. S. (2001). Biometals, 14, 271–313. Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., Shindyalov, I. N. & Bourne, P. E. (2000). Nucleic Acids Res. 28, 235–242. 4. Conclusions Castagnetto, J. M., Hennessy, S. W., Roberts, V. A., Getzoff, E. D., Tainer, J. A. & Pique, M. E. (2002). Nucleic Acids Res. 30, 379–382. We have summarized the chemical moieties that have been shown to Cox, E. H. & McLendon, G. L. (2000). Curr. Opin. Chem. Biol. 4, 162–165. play the role of chelator for a metal ion in the active site of Hemavathi, K., Kalaivani, M., Udayakumar, A., Sowmiya, G., Jeyakanthan, J. metalloproteins, and thereby to inhibit the biological activity of the & Sekar, K. (2010). J. Appl. Cryst. 43, 196–199. Hsin, K., Sheng, Y., Harding, M. M., Taylor, P. & Walkinshaw, M. D. (2008). J. target metalloprotein. Both common and specific chelating groups Appl. Cryst. 41, 963–968. are found for various metal ions. Previously it has been difficult to Irwin, J. J., Raushel, F. M. & Shoichet, B. K. (2005). Biochemistry, 44, 12316– design inhibitors of metalloproteins because of the difficulty in 12328. identifying chemical moieties suitable for the chelative inactivation of Kleifeld, O., Frenkel, A., Martin, J. M. & Sagi, I. (2003). Nat. Struct. Biol. 10, the central metal ion. In this regard, MetLigDB is expected to be 98–103. Liu, T., Lin, Y., Wen, X., Jorissen, R. N. & Gilson, M. K. (2007). Nucleic Acids useful in selecting a chelating group that should be present in the Res. 35, D198–D201. inhibitor of a metalloprotein, which should be the starting point in Pierce, M. M., Raman, C. S. & Nall, B. T. (1999). Methods, 19, 213–221. designing new metalloprotein inhibitors. MetLigDB is freely acces- Sang, Q. X., Jin, Y., Newcomer, R. G., Monroe, S. C., Fang, X., Hurst, D. R., sible at http://silver.sejong.ac.kr/MetLigDB. Users are invited to Lee, S., Cao, Q. & Schwartz, M. A. (2006). Curr. Top. Med. Chem. 6, 289–316. Shi, J., Ye, J.-W., Song, T.-Y., Zhang, D.-J., Ma, K.-R., Ha, J., Xu, J.-N. & Zhang, contact us by selecting ‘Contact’ on the menu bar on that site. P. (2007). Inorg. Chem. Commun. 10, 1534–1536. Suggestions for changes in data sets and in web site features are Shi, W. & Chance, M. R. (2008). Cell. Mol. Life Sci. 65, 3040–3048. welcomed. Studies using MetLigDB should cite this paper. Williams, R. J. (2003). Chem. Commun. pp. 1109–1113.

J. Appl. Cryst. (2011). 44, 878–881 Choi, Kang and Park  MetLigDB 881