Quantum Mechanical Models of the Resting State of the Vanadium

Apr 16, 2004 - Giuseppe Zampella,†,‡ Joslyn Yudenfreund Kravitz,†,§ Charles Edwin Webster,| Piercarlo Fantucci,‡. Michael B. Hall,| Heather A...
0 downloads 0 Views 227KB Size
Inorg. Chem. 2004, 43, 4127−4136

Quantum Mechanical Models of the Resting State of the Vanadium-Dependent Haloperoxidase Giuseppe Zampella,†,‡ Joslyn Yudenfreund Kravitz,†,§ Charles Edwin Webster,| Piercarlo Fantucci,‡ Michael B. Hall,| Heather A. Carlson,*,§,⊥ Vincent L. Pecoraro,*,§ and Luca De Gioia*,‡ Department of Biotechnology and Biosciences, UniVersity of Milano-Bicocca, Piazza della Scienza 2, 20126 Milano, Italy, Department of Chemistry, UniVersity of Michigan, Ann Arbor, Michigan 48109-1055, Department of Chemistry, Texas A&M UniVersity, College Station, Texas 77843, and Department of Medicinal Chemistry, College of Pharmacy, UniVersity of Michigan, Ann Arbor, Michigan 48109-1065 Received November 17, 2003

Density functional theory has been used to investigate structural and electronic properties of complexes related to the resting form of the active site of vanadium haloperoxidase as a function of environment and protonation state. Results obtained by studying models of varying size and complexity highlight the influence of environment and protonation state on the structure and stability of the metal cofactor. The study shows that, in the trigonal bipyramidal active site, where one axial position is occupied by a key histidine, the trans position cannot contain a terminal oxo group. Further, a highly negatively charged vanadate unit is not stable. Protonation of at least one equatorial oxo ligand appears necessary to stabilize the metal cofactor. The study also indicates that, while at rest within the protein, the vanadate unit is most likely an anion with an axial hydroxide and an equatorial plane containing two oxos and a hydroxide. For the neutral, protonated state of the vanadate unit, there were two minima found. The first structure is characterized by an axial water with two oxo and one hydroxo group in the equatorial plane. The second structure contains an axial hydroxo group and an equatorial plane composed of one oxo and two hydroxo oxygen atoms. These two species are not significantly different in energy, indicating that either form may be important during the catalytic cycle. These data support the initial crystallographic assignment of an axially bound hydroxide, but an axial water is also a possibility. This study also shows that the protonation state of the vanadate ion is most likely greater than previously proposed.

the following reaction scheme.3

Introduction Evidence of the involvement of vanadium ions in the vanadium-dependent haloperoxidase (VHPO)1 enzymes has stimulated interest in biomimetic models for catalytic halogenation of organic substrates.2 In general, haloperoxidases catalyze the two electron oxidation of halide ions via * To whom correspondence should be addressed. E-mail: carlsonh@ umich.edu (H.A.C.); [email protected] (L.D.G.); [email protected] (V.L.P). † These authors contributed equally to this work. ‡ University of Milano-Bicocca. § Department of Chemistry, University of Michigan. | Texas A&M University. ⊥ Department of Medicinal Chemistry, College of Pharmacy, University of Michigan. (1) Wever, R.; Kustin, K. In AdVances in Inorganic Chemistry; Sykes, A. G., Ed.; Academic Press: New York, 1990; Vol. 35, pp 103-137. (2) de la Rosa, R. I.; Clague, M. J.; Butler, A. J. Am. Chem. Soc. 1992, 114, 760-761.

10.1021/ic0353256 CCC: $27.50 Published on Web 04/16/2004

© 2004 American Chemical Society

H2O2 + H+ + X- f H2O + “HOX”

(1)

The enzyme is named on the basis of its ability to oxidize a halide (X-). Thus, chloroperoxidases can oxidize chloride, bromide, or iodide, whereas bromoperoxidases are incapable of using chloride as a substrate. The identity of “HOX” is dependent upon both the pH of the reaction and the halide involved. Initial studies indicated that the chloroperoxidase oxidizes a chloride ion to HOCl,4 whereas the products of bromide oxidation are generally thought to be a thermodynamic distribution of Br3-, HOBr, and Br2.5 The “HOX” (3) Everett, R. R.; Butler, A. Inorg. Chem. 1989, 28, 393-395. (4) van Schijndel, J. W. P. M.; Barnett, P.; Roelse, J.; Vollenbroek, E. G. M.; Wever, R. Eur. J. Biochem. 1994, 225, 151-157. (5) Butler, A. Coord. Chem. ReV. 1999, 187, 17-35.

Inorganic Chemistry, Vol. 43, No. 14, 2004

4127

Zampella et al. anism.11

Figure 1. Atomic structure of the active site of the vanadium containing chloroperoxidase from the fungus CurVularia inaequalis as obtained from X-ray structure analysis.7

species can then halogenate an organic substrate if one is present. More recent work has suggested that certain organic substrates can be regioselectively halogenated,6 presumably in the active site of the enzyme. Of this enzyme family, the vanadium-dependent chloroperoxidase (VCPO) from the pathogenic fungus CurVularia inaequalis has been investigated the most thoroughly. Its X-ray crystal structure was the first of this family to be solved and has shed light on the characteristics of the active site of the enzyme.7 Subsequent crystal structures of other VHPOs have revealed a drastic difference in the overall structure of the protein.8,9 However, the entire first and second hydrogenbonding shells for the cofactor are conserved across species, with the exception of the substitution of a histidine in the bromoperoxidase for a phenylalanine in the chloroperoxidase.8 We will, therefore, refer to protein residues using the sequence and numbering in the 2.15 Å VCPO structure (PDB code 1VNI) of Messerschmidt and Wever.7 In the native state, the vanadate cofactor (VO43-) is bound to the protein through a single covalent bond from the V to the N of His496 (Figure 1). The high negative charge of the cofactor is offset by a number of protonated amino acids in the active site which donate hydrogen bonds to the oxygen atoms of the cofactor (Lys353, Arg360, His404, and Arg490). The crystal structure of this enzyme indicates that the vanadate unit is a trigonal bipyramidal structure with a hydroxide and His496 in the axial positions and three oxo moieties in the equatorial plane.7 There is no evidence that the vanadium undergoes redox cycling during catalysis, and hence, it is believed that the vanadium ion plays the role of a strong Lewis acid with respect to the primary oxidant, H2O2.10 Kinetic studies suggest that protonation of the bound peroxo group is a crucial step in the heterolytic cleavage of the O-O bond. These proposals suggest that the protonated oxygen is then transferred to the halide according to an oxo-transfer mech(6) Martinez, J. S.; Carrol, G. L.; Tschirret-Guth, R. A.; Altenhoff, G.; Little, R. D.; Butler, A. J. Am. Chem. Soc. 2001, 123, 3289-3294. (7) Messerschmidt, A.; Wever, R. Proc. Natl. Acad. Sci., U.S.A. 1996, 93, 392-396. (8) Weyand, M.; Hecht, H.-J.; Kie, M.; Liaud, M.-F.; Vilter, H.; Schomburg, D. J. Mol. Biol. 1999, 293, 595-611. (9) Isupov, M. N.; Dalby, A. R.; Brindley, A. A.; Izumi, Y.; Tanabe, T.; Murshudov, G. N.; Littlechild, J. A. J. Mol. Biol. 2000, 299, 10351049. (10) de Boer, E.; Boon, K.; Wever, R. Biochemistry 1988, 27, 1629-1635.

4128 Inorganic Chemistry, Vol. 43, No. 14, 2004

Finally, the hypohalous acid formed can act as a stereoselective halogenating agent.12 Site-directed mutagenesis experiments (alanine scanning) have elucidated the relative importance of various activesite residues in catalytic activity.13 His496 is the primary chemical anchor for vanadium, and its replacement with Ala results in complete loss of activity. Substitution of the positively charged residues Lys353, Arg360, and Arg490 revealed less severe effects, but of these mutants, the largest decrease in activity is observed for K353A. In fact, Lys353 is the only residue which interacts with the peroxo-bound moiety directly. It has been suggested that it plays a role in polarizing the peroxo bond, making it more susceptible to nucleophilic attack.14 The crystal structures indicate that Arg490 is involved in two hydrogen bonds with oxygen atoms of the cofactor, both in the native and the peroxo forms, whereas Arg360 donates only one hydrogen bond. This difference in participation is reflected in the milder reduction of catalytic efficiency shown by the R360A mutant with respect to K353A and R490A. Peroxo-vanadium complexes are good functional models of the VHPO and have been shown to epoxidize alkenes,15 hydroxylate hydrocarbons,16 and oxidize many organic17 and inorganic compounds,18 including alcohols19 and sulfides.20 However, attempts to synthesize structural models of the VHPO active site have not been as successful, and a synthetic complex reproducing the exact coordination environment of vanadium in VHPO has not yet been obtained. It has been proposed that such a complex has not been isolated because a five-coordinate imidazole complex is unstable unless sequestered in a protein active site.21 Therefore, synthetic models have thus far provided limited insight into structural and electronic features of the vanadium cofactor. In this context, quantum mechanics calculations are useful tools to investigate structural and electronic properties of models of the active site of metalloenzymes. Though computational investigations of vanadium complexes have been reported,22,23 the ligand sets of the complexes in previous studies did not reflect the coordination environment of the (11) Colpas, G. J.; Hamstra, B. J.; Kampf, J. W.; Pecoraro, V. L. J. Am. Chem. Soc. 1996, 118, 3469-3478. (12) Clague, M. J.; Keder, N. L.; Butler, A. Inorg. Chem. 1993, 32, 47544761. (13) Hemrika, W.; Renirie, R.; Macedo-Ribeiro, S.; Messerschmidt, A.; Wever, R. J. Biol. Chem. 1999, 274, 23820-23827. (14) Messerschmidt, A.; Prade, L.; Wever, R. Biol. Chem. 1997, 378, 309315. (15) Sheldon, R. A.; Kochi, J. K. Metal-catalyzed oxidations of organic compounds: mechanistic principles and synthetic methodology including biochemical processes; Academic Press: New York, 1981. (16) Mimoun, H.; Saussine, L.; Daire, E.; Postel, M.; Fischer, J.; Weiss, R. J. Am. Chem. Soc. 1983, 105, 3101-3110. (17) Conte, V.; DiFuria, F.; Moro, S. J. Phys. Org. Chem. 1996, 9, 329336. (18) Butler, A.; Clague, M. J.; Meister, G. E. Chem. ReV. 1994, 94, 625638. (19) Bartolini, O.; Conte, V.; Di Furia, F.; Modena, G. NouV. J. Chim. 1985, 9, 147-150. (20) Smith, T. S., II; Pecoraro, V. L. Inorg. Chem. 2002, 41, 6754-6760. (21) Ligtenbarg, A. G. J. In Chemistry; University of Groningen: Groningen, Holland, 2001. (22) Buhl, M. J. Comput. Chem. 1999, 20, 1254-1261. (23) Conte, V.; Bortolini, O.; Carraro, S.; Moro, S. J. Inorg. Biochem. 2000, 80, 41-49.

Models of Vanadium-Dependent Haloperoxidase

vanadium atom observed in the enzyme. In this paper, we report results obtained using density functional theory (DFT) to investigate structural properties of complexes related to the native form of VCPO as a function of environment and protonation state.

Scheme 1. Schematic Representation of the Starting Structures Used in Calculations To Investigate Native Models of the Active Site of VCPOa

Computational Methods All calculations presented in this paper were performed using Gaussian 98.24 Quantum chemical calculations were carried out using models of four sizes, based on the first shell of residues in contact with the vanadate cofactor. Bond distances and angles for all models were taken from the crystal structure of the enzyme.7 The smallest models used an ammonia in place of His496 and contained no external hydrogen-bonding molecules. The second set of models was expanded to include contributions from one hydrogen-bonding partner in the active site, represented by ammonia or ammonium, depending on the charge state desired. The third set of models was expanded once again, and small molecule models which more closely resembled the active-site residues were included. His496 was modeled with an imidazole, and Lys353 and Arg490 were modeled by methylamine and guanidine, respectively. Guanidinium and ammonium can be considered acceptable structural models for the polar heads of arginine and lysine, respectively, but they are expected to be slightly less basic than the corresponding amino acids.25 Finally, the largest scaffolds were expanded to include the six amino acids from the crystal structure that appear to be donating hydrogen bonds to the cofactor. In addition to Arg490, another guanidinium and a glycyl-serine dipeptide were added to represent Arg360, Ser402, and Gly403, respectively. The representative group for Lys353 was reduced to an ammonium cation, and the representative group for His496 was changed to be a methylimidazole molecule. Additionally, a crystallographically derived water was included to act as a hydrogen-bonding partner for the axial oxygen ligand. Due to the size of the scaffolds (6162 atoms), a lower level of theory was required. With the exception of the largest models, all calculations were carried out using DFT.26 They employed the hybrid exchangecorrelation functional, B3LYP.27 The basis set used for vanadium was the effective core potential basis set of Hay and Wadt (LanL2DZ)28 as modified by Couty and Hall (341/341/41 for V), where the two outermost p functions have been replaced by a (41) split of the optimized 4p function.29 The cc-pVDZ30 basis set was (24) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Zakrzewski, V. G.; Montgomery, J. A., Jr.; Stratmann, R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A. D.; Kudin, K. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Cioslowski, J.; Ortiz, J. V.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Gonzalez, C.; Challacombe, M.; Gill, P. M. W.; Johnson, B. G.; Chen, W.; Wong, M. W.; Andres, J. L.; Head-Gordon, M.; Replogle, E. S.; Pople, J. A. Gaussian 98, revision A.11.1; Gaussian, Inc.: Pittsburgh, PA, 1998. (25) Tantillo, D. J.; Chen, J.; Houk, K. N. Curr. Opin. Chem. Biol. 1998, 2, 743-750. (26) Parr, R. G.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: Oxford, U.K., 1989. (27) Lee, C.; Yang, W.; Parr, R. G. Phys. ReV. B 1988, 37, 785-789. (28) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 83, 270-283. (29) Couty, M.; Hall, M. B. J. Comput. Chem. 1996, 17, 1359-1370. (30) Woon, D.; Dunning, T. H. J. J. Chem. Phys. 1993, 98, 1358-1371.

a N is either an ammonia molecule or an imidazole. Complexes which minimize to biologically relevant structures are shown in bold.

used for N and O atoms, and D9531 for hydrogens. Full geometry optimizations were carried out, and minima were verified through frequency calculations. The largest models underwent a series of partial optimizations at the Hartree-Fock level of theory and using LanL2DZ for V and D95V31 for all other atoms. All hydrogens were minimized as well as the vanadate oxygens and the imidazole ring in the axial position. Vanadium and the heavy atoms of the protein and water were held fixed. Figures were drawn using ChemEdit.32

Results In Scheme 1, the figures, and tables, the models used in this study have been labeled with an alphanumeric designation defined as follows: The first part of the name refers to the charge state of the vanadate cluster (a, n, or c, to designate anionic, neutral, or cationic structures, respectively). A preceding number indicates the charge of the system if it is not one or zero (i.e., 3a for the -3 charge state). The letter is followed by a number which corresponds to the starting structure of the vanadate unit as defined in Scheme 1. For complexes with one or more external hydrogen-bonding molecules (Figures 3-5), the vanadate structure code will be followed by a superscript for each external molecule. The superscript will be either K, n, or c representing a potassium ion, a neutral hydrogen-bond donor, or positively charged hydrogen-bonding molecule, respectively. For species with (31) Dunning, T. H. J.; Hay, P. J. In Modern Theoretical Chemistry; Schaefer, H. F., III, Ed.; Plenum: New York, 1976; Vol. 3, pp 1-28. (32) Lim, C. E.; Jorgensen, W. L. In Encyclopedia of Computational Chemistry; Schleyer, P. v. R., Ed.; John Wiley and Sons: New York, 1998; Vol. 5, p 3295.

Inorganic Chemistry, Vol. 43, No. 14, 2004

4129

Zampella et al.

Figure 3. Calculated minima of anionic and neutral structures with one added hydrogen-bonding molecule. Atom labels correspond to those in Table 2. The panels in the top row show the anionic minima in order of increasing energy from left to right. The first two panels of the second row show the minimum structures for a neutral vanadate unit with an axial water. The last panel in the second row shows the minimum structure for a neutral vanadate unit with an axial hydroxo moiety. The panels of the third row show other low-energy structures of this type with energy increasing from left to right. Figure 2. Calculated minima of non-hydrogen-bonded structures. Atom labels correspond to those in Table 1. Panel 2aKK shows the minimized structure of the dianionic vanadate unit with two potassium atoms included in the calculation. Panels a1 and n1 show the structures of the global minima for the anionic and neutral units, respectively. Panels a3, n3, and n5 show the biologically relevant minima for their respective charge states. Panels c5 and c6 show the low energy cationic structures.

more than one local energy minimum, the alphanumeric designation will also contain the subscript i, ii, iii, or iv, depending on the number of minima observed. The bound amine ligand, which mimics His496, will be referred to generically as N. The identity of N (imidazole or ammonia) will be specified as appropriate. No External Hydrogen-Bonding Group (Figure 2). (VO4-3)N and (HVO4-2)N (Scheme 1: 3a, 2a1, and 2a2). Initially, we investigated the structural properties of models of the active site without the stabilizing influence of the protein environment. Optimization of unprotonated (3a) and singly protonated (2a1 and 2a2) structures led to separation of each into a tetrahedral vanadate and an amine. The results were the same when either imidazole or ammonia occupied the axial N position. The calculation was repeated for the imidazole complexes using an extended basis set (triple-ζ + polarization on all atoms) and two different functionals (B3LYP and BP8633). Similar results were obtained for all methods (data not shown). To investigate whether neutralizing the negative charge of the complex would improve its stability, two K+ ions were (33) Becke, A. D. Phys. ReV. A 1988, 38, 3098-3100. Perdew, J. P. Phys. ReV. B 1986, 33, 8822-8824; 1986, 34, 7406.

4130 Inorganic Chemistry, Vol. 43, No. 14, 2004

added to the HVO42--imidazole model (Figure 2, 2aKK). In the optimized complex, the N-V bond was not cleaved, but it was significantly longer (2.216 Å) than that reported in the X-ray structure of the enzyme (1.96 Å).14 H2VO4NH3-1 and H3VO4NH3 (Scheme 1: a1-4 and n1-5). In order to efficiently examine all possible anionic and neutral species, a minimal model was used. There were four possible monoanionic complexes with the formula H2VO4NH3-1 and five possible neutral complexes with the formula H3VO4NH3 (Scheme 1, a1-4 and n1-5). In both charge states, the lowest-energy structure was a tetrahedral vanadate with a dissociated ammonia (Figure 2, a1 and n1). None of the structures with an oxo group axial to the ammonia resulted in stable five-coordinate structures (Scheme 1, a1, a2, n1, and n2). The lowest-energy, five-coordinate species in each charge state had a hydroxo group in the axial position (Figure 2, a3 and n3). A neutral species with water in the axial position was also a stable minimum (Figure 2, n5) but was 3.2 kcal/mol higher in energy than the axial hydroxo structure. Important bond distances and angles for these complexes are shown in Table 1. The final structures for a2, a4, n2, and n4 are neither biologically relevant nor low in energy so they are not shown (their structures are available in the Supporting Information). H4VO4NH3+1 (Scheme 1: c1-6). There were six cationic isomers of the model, all of which remained five-coordinate upon minimization (Scheme 1, c1-6). However, only half maintained the trigonal bipyramidal arrangement with the ammonia in one of the axial positions. The lowest-energy

Models of Vanadium-Dependent Haloperoxidase

Figure 4. Calculated minima of cationic structures with one added hydrogen-bonding molecule. Atom labels correspond to those in Table 3. Panel c1ni /n1c shows a high-energy minimum with an axial oxo group. Panels n3ci and n5ci show the two minimum structures of the overall cationic species, which are almost equal in energy. Panel n5cii shows a higher energy minimum which is an isomer of n5ci.

Figure 5. Calculated minima of structures with two added hydrogenbonding molecules. Atom labels correspond to those in Table 4. G stands for guanidine, and MA stands for methylamine. Panel a3nn shows the minimum structure for an overall anionic species. Panel a1nc shows the minimum structure for an overall neutral species containing an anionic vanadate unit and a protonated guanidinium molecule. Panel n5nc shows the minimum structure of a cationic species containing a neutral vanadate unit and a protonated guanidinium molecule. The panels in the bottom row show three nearly energetically equivalent overall neutral structures.

cation (Figure 2, c5) had water and ammonia in the axial positions, and an oxo and two hydroxos in the equatorial plane. The second lowest-energy cation (Figure 2, c6) contained two waters, one of which was trans to the ammonia, and two oxos in the equatorial plane. Of the remaining four cations (c1-4, structures available in Supporting Information), the highest-energy structure had four

hydroxo groups (c3), and although it remained trigonal bipyramidal, it was extremely high in energy (over 18 kcal/ mol higher than c5) and is not considered to be biologically relevant. In c1, c2, and c4, the vanadate rearranged so that there was no longer a clear axis made up of one oxo ligand and the amine. The bonding parameters and relative energies of c5 and c6 are in Table 1. One External Hydrogen-Bonding Molecule (Figures 3 and 4). H2VO4NH3-1 + NH3 and H3VO4NH3 + NH3 (Scheme 1: a1-4 and n1-5). External hydrogen-bonding molecules were added to the models. In the anionic case, addition of an ammonia molecule to each unique position in the equatorial plane did not significantly change the behavior or structure of the molecules. Loss of the bound ammonia to yield a tetrahedral vanadate again led to the lowest-energy structure (a1n, not shown). Addition of a hydrogen-bond donor to a3 led to three five-coordinate structures, all within 1 kcal/mol, which differed only in the position of ammonia with respect to the equatorial ligands (Figure 3, a3ni-iii). The various a2 and a4 models rearranged to four-coordinate vanadium species which were much higher in energy than the five-coordinate minima (structures not shown). The bond angles and distances for the vanadate unit in structures a1n and a3ni-iii are virtually the same as those for a1 and a3, respectively (Table 1), so only the hydrogenbond interactions and relative energies for a1n and a3ni-iii are listed in Table 2. In the case of the neutral species, 20 starting structures were minimized (two for n1, six for n2, and four each for n3-5, each with the additional ammonia molecule hydrogen bonding to each distinct oxygen ligand in both the axial and equatorial positions). The only biologically relevant structures were found to be the same as in the previous section (Figure 3, n3ni-iv and n5ni-ii). The lowest-energy structures, n5ni and n3ni, had an axial ligand donating a hydrogen bond to the free ammonia, but several other five-coordinate conformers were found to be stable minima. For the axial water case, the other minimum had the free ammonia donating a hydrogen bond to the equatorial hydroxide, rather than the oxo (n5nii). This configuration was within 1 kcal/mol in energy of n5ni. In the axial hydroxide case, there were three minima (n3nii-iv) approximately 3 kcal/mol higher in energy than n3ni. These four structures have roughly the same geometry about the vanadium but differ in the identities of the donors and acceptors of the hydrogen bonds (Figure 3). The majority of the bond angles and distances for structures n3ni-v and n5ni-ii are virtually the same as those in Table 1, so only hydrogen-bond interactions and relative energies are given in Table 2. Like the anionic case, the vanadium would not remain fivecoordinate with an oxo group trans to the bound ammonia. However, the lowest-energy, neutral species was fivecoordinate, unlike the anionic case. H3VO4NH3 + NH4+ and H4VO4NH3+ + NH3 (Scheme 1: n1-5 and c1-6). The hydrogen bonds for overall cationic species were modeled as (1) a neutral vanadate complex with an external ammonium ion and (2) a cationic vanadate complex with an external ammonia molecule. The Inorganic Chemistry, Vol. 43, No. 14, 2004

4131

Zampella et al. Table 1. Bonding Parameters for Important Non-Hydrogen-Bonding Vanadate Species namea

V-Oab

V-Obb

V-Ocb

V-Odb

V-Nb

Ob-V-Oac

Ob-V-Nc

Ob-V-Occ

Ob-V-Odc

Oa-V-Nc

∆Ed

2aKK

2.11 1.63 1.88 1.59 1.81 2.22 2.11 2.12

1.71 1.65 1.92 1.79 1.59 1.62 1.57 2.12

1.67 1.82 1.64 1.74 1.83 1.82 1.76 1.59

1.67 1.85 1.64 1.79 1.80 1.62 1.76 1.59

2.26 3.53 2.29 3.71 2.26 2.14 2.12 2.12

85.7 113.3 88.8 110.7 102.4 86.3 98.8 75.2

87.1 55.3 72.6 85.9 88.5 94.6 93.4 85.7

118.0 108.0 119.1 109.5 117.6 125.1 114.3 125.1

118.3 108.2 121.2 108.1 112.5 114.4 110.9 125.1

172.8 146.9 161.2 151.6 165.4 -179.8 167.8 160.9

N/A 0 +13.6 0 +5.14 +8.34 0 +6.37

a1 a3 n1 n3 n5 c5 c6

a Complexes as defined in Scheme 1. b Distances in Å. c Angles in deg. d Difference in energy, in kcal/mol, as compared to lowest-energy structure of a given charge state.

Table 2. Hydrogen-Bonding Interactions in Structures as in Figure 3 structure a1n a3ni a3nii a3niii n5ni n5nii n3ni n3nii n3niii n3niv a

O‚‚‚H-N distancea

O‚‚‚H-N angleb

1.95 2.22 (to d) 2.23 (to c) 2.36 (to b) 2.56 (to c) 2.52 (to c) 2.73 (to d) 2.08 (to b) 2.15 (to c) 2.28 (to b) 2.61 (to b)

160.1 144.3 143.7 139.5 127.6 127.2 126.1 136.9 137.1 127.9 124.5

2.30 (to d)

135.5

N‚‚‚H-Oeq distancea 2.02

Oax‚‚‚H‚‚‚N angleb

153.1

1.74 (to d) 1.72 (to c)

∆Ec 0 12.05

165.5 161.3

2.35 (to a)

134.9

12.28

2.27 (to a)

145.1

12.90

1.65 (from a) 1.66 (from a) 1.68 (from a)

162.6 158.5 164.5

2.19 (to a) 1.77 (from a)

137.0 159.4

0 0.979 0.786 2.40 2.72 3.53

Distances in Å. b Angles in deg. c Difference in energy, in kcal/mol, as compared to lowest-energy structure of a given charge state.

Table 3. Hydrogen-Bonding Interactions in Cationic Structures as in Figure 4 structure c1ni (n1c) n3ci n5ci n5cii a

Oax‚‚‚H‚‚‚N distancea

N‚‚‚H-Oeq angleb

O‚‚‚H-N distancea

O‚‚‚H-N angleb

1.79 (to b) 1.79 (to c) 1.79 (to c) 1.88 (to b) 1.73 (to d) 1.79 (to c) 1.82 (to d)

149.8 150.0 144.3 140.8 150.1 147.0 148.9

Oax‚‚‚H-N distancea

Oax‚‚‚H-N angleb

∆Ec vs n3ci +5.37

1.64 N/A

149.3 N/A

0 +0.320 +4.19

Distances in Å. b Angles in deg. c Energy difference in kcal/mol.

external molecules were placed in all unique positions in the equatorial plane to create a series of minima in both states. After minimization, it was found that the type 1 structures were significantly lower in energy than the type 2 complexes. In fact, some type 2 starting structures transferred a proton from the vanadate unit to the ammonia resulting in type 1 minima. Specifically, one such complex started as a c5n structure and minimized to a structure similar to n3c (this structure is not shown in Figure 4 because it was too high in energy to be considered biologically relevant, but both the starting and minimum coordinates are available in the Supporting Information). The type 1 structures yielded several five-coordinate minima that were similar to the overall neutral minima (Figure 4, n3ci and n5ci-ii). However, these biologically relevant structures were not the global minima. The lowestenergy, biologically relevant structure had a hydroxo in the axial position (Figure 4, n3ci), and the structure with a water trans to the ammonia (Figure 4, n5ci) was very close in energy. Related structures were significantly higher in energy (Figure 4, n5cii and n3cii, not shown).

4132 Inorganic Chemistry, Vol. 43, No. 14, 2004

Table 4. Hydrogen-Bonding Interactions in Structures with Abbreviations as In Figure 5 G-H‚‚‚O G-H‚‚‚O G‚‚‚H-O G‚‚‚H-O MA‚‚‚H‚‚‚O MA‚‚‚H‚‚‚O distancea angleb distancea angleb distancea angleb n3nn 1.82 (to a) n5nni 2.21 (to b) 2.00 (to d) n5nnii 1.93 (to b)

117.6 178.9 154.2 174.1

n5nc 1.71 (to b) 1.73 (to d) a3nn 1.88 (to a) 1.90 (to c) a1nc 1.58 (to a) 1.54 (to b)

177.4 178.8 173.1 167.8 179.2 178.6

a

Distances in Å

b

1.71 (to d)

176.4

1.80 (from c) 1.81 (from c)

170.2 175.5

1.75 (to c)

177.1

2.36 (to b) 2.34 (to b) 1.68 (from c)

135.9 137.7 176.8

2.11 (to d)

158.5

2.12 (to c)

169.7

Angles in deg.

Of the type 2 species, only two were biologically relevant (c5n and c6n, structures not shown). These two complexes were both over 10 kcal/mol higher in energy than the lowestenergy, type 1 species, and the additional ammonia molecule accepted a hydrogen bond from the bound ammonia rather than donating hydrogen bonds to the oxygen ligands. Hydrogen-bonding parameters for low-energy, biologically relevant overall cationic structures are in Table 3. The starting and minimized coordinates, as well as the final energies for all structures, can be found in the Supporting Information. Two External Hydrogen-Bonding Molecules (Figure 5). External molecules that were more representative of the protein active site were incorporated: N was an imidazole molecule and Lys353 and Arg490 were modeled in their crystallographic positions by methylamine and guanidine, respectively.

Models of Vanadium-Dependent Haloperoxidase Table 5. Bonding Parameters for a3ncc and a3nnc a V-X distancec

Ob-V-X angled

H-bond 1 distancec

Oax Ob Oc Od Nb

1.87 1.84 1.62 1.60 2.15

84.0

2.55 1.67 (K353) 1.77 (R360) 1.82 (S402)

a3ncc 160.8 168.2 164.4 168.0

2.01 (G403) 2.10 (R490) 2.00 (R490)

164.7 175.4 172.8

Oax Ob Oc Od Nb

1.90 1.81 1.64 1.58 2.12

2.25 1.83 (K353) 2.45 (R360) 1.62 (S402)

a3nnc 173.2 169.0 165.3 176.2

1.91 (to R360) 1.75 (R490) 2.09 (G403)

163.3 150.3 118.3

114.1 136.5 88.8 87.8 126.6 121.9 82.8

H-bond 1 angled

H-bond 2 distancec

H-bond 2 angled

H-bond 3 distancec

H-bond 3 angled

2.00 (R490)

141.1

a Oax and Ob are hydroxo groups. Oc and Od are oxo groups. Oax accepts a hydrogen bond from water. b Oa-V-N angle ) 171.2° for a3ncc and 166.6° for a3nnc. c Distances in Å. d Angles in deg.

HVO4N2- + Externals (Scheme 1, 2a1-2). Minimization of the HVO42--imidazole adduct resulted in V-N bond cleavage (structures available in Supporting Information). H2VO4N1- + Externals and H3VO4N + Externals, (Scheme 1: a3, n1, n3, and n5). There was one stable, fivecoordinate, anionic species (Figure 5, a3nn). In this complex, the guanidine became involved in hydrogen bonds with the axial hydroxo and one equatorial oxo groups. Hydrogenbonding parameters are listed in Table 4. Calculations with either a protonated methylamine or guanidine resulted in proton transfer to an oxo group. To verify that this result was not an artifact of the higher acidity of methylamine and guanidine with respect to lysine and arginine, models which included the full amino acid side chains were also optimized, and proton transfer to the oxo groups was observed (data not shown). Protonation of a3nn resulted in three neutral vanadate isomers (Figure 5, n3nn, n5nni, and n5nnii). The stabilities of these species were similar, with energy differences on the order of a few kcal/mol. One anionic vanadate species with a cationic guanidinium group, yielding an overall neutral structure (Figure 5, a1nc) was also investigated. In this case, it was again observed that the imidazole did not completely dissociate, but it remained at a weakly bonding distance (2.47 Å). Only one stable, overall monocationic species was found (Figure 5, n5nc). The neutral vanadate group was composed of an axial water, one hydroxo, and two oxo groups. The cationic guanidine did not protonate the vanadate during minimization. Hydrogen-bonding interactions are listed in Table 4. Hydrogen-Bonding Scaffold (Figure 6). Calculations of the largest, most biologically relevant scaffold showed that a vanadate unit carrying a highly negative charge (-3 or -2) was not stable and lost the imidazole. In the case of a dianionic vanadate (2a2ccc and 2a2ncc, structure not shown), allowing the protons to minimize while all the heavy atoms were held fixed resulted in one proton from the ammonia being almost immediately transferred to an equatorial oxo group. This transfer generated the two structures, both with an anionic vanadate unit, that were used in the remaining calculations: an overall neutral structure in which one guanidine was deprotonated (Figure 6, a3nnc), and an overall

Figure 6. Stereoview of calculated partial minima of anionic vanadate units contained in overall cationic (a3ncc) and overall neutral (a3nnc) representative hydrogen-bonding scaffolds.

cation in which both guanidinium molecules were protonated (Figure 6, a3ncc). In both final structures, the axial hydroxo accepted a hydrogen bond from the nearby crystallographic water, and each equatorial ligand accepted two hydrogen bonds from the scaffold. In the overall cationic structure (a3ncc), the equatorial oxygens minimized to positions that closely resembled the crystal structure (average displacement was only 0.29 Å). In the overall neutral structure (a3nnc), the equatorial ligands moved on average 0.93 Å from their crystallographic positions, rotating approximately 34° within the equatorial plane to allow the equatorial hydroxide to donate a hydrogen bond to the deprotonated guanidine. This resulted in one additional hydrogen bond between the cofactor and the scaffold. Bonding parameters for these models are listed in Table 5. Inorganic Chemistry, Vol. 43, No. 14, 2004

4133

Zampella et al. Discussion Although the vanadium-dependent haloperoxidases comprise a rather rare family of enzymes, their ability to stereoselectively halogenate natural organic molecules makes them very important.6 These compounds are often the leads for new pharmaceuticals or precursors to polymeric materials. A better fundamental understanding of how these enzymes work, and more specifically, what role the vanadate cofactor plays in the catalytic cycle, will be immensely beneficial to both the natural products and inorganic chemistry communities. The DFT calculations of the haloperoxidase presented in this paper allow us to look more closely at this elusive metal center and make suggestions as to its role in the catalytic process. Although the coordinates for the starting structures in this study are taken from the crystal structure of the chloroperoxidase,7 the results can be extended to the entire enzyme family since only residues which are conserved in all haloperoxidases8,9 were used in the models. It has been suggested that the structure of all VHPO active sites is dianionic and trigonal bipyramidal. It purportedly contains three oxo groups in the equatorial plane and one hydroxo ligand trans to His496, which provides the only covalent bond to the protein (Figure 1).14 However, the resolution of the highest quality crystal structure currently available is 2.0 Å, and the protons are not visible,8 so the exact ligand set and protonation state of the cofactor are not strictly known. Therefore, in our studies we explored all permutations of oxo, hydroxo, and water groups as ligands for the vanadium in various charge states (Scheme 1). In order to make this complete search tractable, we used the smallest system possible: HxVO4NH3y, where 0 e x e 4 and -3 e y e 1. The arrangement proposed in the crystal structure is x ) 1, y ) -2. Minimization of a tri- or dianionic vanadate unit resulted in dissociation of the ammonia and a tetrahedral vanadate (Scheme 1, 3a1 and 2a1-2). Replacing the ammonia with imidazole did not stabilize the structure. The addition of one, two, or seven hydrogen-bonding groups in the vicinity of the complex was also insufficient to stabilize the complex. However, addition of two potassium ions to the HVO42system resulted in a stable, overall neutral complex (Figure 2, 2aKK). These results highlight the importance of buffering the negative charge on the cofactor, supporting previous proposals suggesting this role for the protein residues forming the active site.13 Upon minimization of the four anionic, five neutral, and six cationic structures shown in Scheme 1, we found that under no circumstances was a structure with a terminal oxo trans to the ammonia stable (Figure 2). The weakening of the bond to the ammonia is most likely a result of the sizable electron donation from the axial oxo to the vanadium (the strong trans influence of the oxo). In all anionic and neutral cases initiated with that arrangement (a1-2 and n1-2), one neutral ligand (water or ammonia) dissociated from the complex and donated a hydrogen bond back to the tetrahedral vanadate. While not biologically relevant, structures a1 and n1 are the global minima for their respective charge states.

4134 Inorganic Chemistry, Vol. 43, No. 14, 2004

In the cationic cases, the ligands rearranged, without dissociation, to nonbiologically important geometries (c1-2). It would seem that dissociation of a neutral ligand to leave a cationic vanadium species is less favorable. Only in the cationic case is the biologically relevant, trigonal bipyramidal species the lowest-energy structure (Figure 2, c5). In the anionic and neutral cases, the lowest-energy, biologically relevant structures have hydroxide in the trans position (Figure 2, a3 and n3), which supports the proposed axial hydroxo group from the crystal structure study. It should be noted, however, that this is dependent upon the ligands in the equatorial plane. The presence of three equatorial oxo ligands, as proposed in the crystal structure, leads to decomposition of the vanadate unit. There is also a stable neutral complex with an axial water (Figure 2, n5) which is 3.2 kcal/mol higher in energy than n3, but there is no stable anion with water in any position. In order to determine if any of the unstable structures could be stabilized by hydrogen bonds, we introduced an external ammonia molecule to each structure in various positions. In the anionic case, there was no significant change. The tetrahedral conformations were always the global minima for the anions. In the neutral case, the vanadate did not remain five-coordinate with a terminal oxo group trans to the bound ammonia. However, the interaction with the external ammonia did make the five-coordinate structure with an axial water the global minimum for the neutral species (Figure 3, n5i). There were five such low-energy complexes that maintained the trigonal bipyramidal geometry upon minimization (n5i-ii and n3i-iii). For investigation of cationic systems, two types of structures were used: a cationic vanadate unit with hydrogen bonds to an ammonia molecule and a neutral vanadate unit with hydrogen bonds to an ammonium cation. The overall cationic system was more stable if the vanadium unit was neutral than if the vanadium unit carried the positive charge. The two lowest-energy structures, which were almost energetically identical, both contained neutral vanadate (Figure 4, n3ci and n5ci). Interestingly, the difference in energy between the two neutral vanadate minima decreases as the complexity of the system is increased. The observed trend is the following: n3 is more than 3 kcal/mol lower in energy than n5; n5ni is 0.8 kcal/mol lower in energy than n3ni, but n3ci is only 0.3 kcal/mol lower in energy than n5ci. On the basis of the above results, we wanted to determine the effects of multiple hydrogen-bond donors that simultaneously complement several of the oxygen ligands on vanadium. This led us to investigate our next level of models which incorporated an environment that was more representative of the protein active site. In place of two ammonia molecules, we used an imidazole, a methylamine, and a guanidine to mimic His496, Lys353, and Arg490, respectively. Under all previous conditions there were a limited number of stable configurations of the vanadate cofactor. It was not feasible to study all possibilities in the large system, so we based the conformations of the large systems on the stable smaller models (Figure 5).

Models of Vanadium-Dependent Haloperoxidase

Results obtained investigating this model indicate that hydrogen bonds between methylamine or guanidine and an HVO42- cofactor (Scheme 1) were not sufficient to withdraw enough electron density from the metal center to stabilize the five-coordinate complex. Remarkably, optimization of any dianionic vanadium complex initially complemented by charged hydrogen-bonding groups resulted in proton transfer to an oxo group to give a monoanionic vanadate structure (Figure 5, a3nn). This was similar to the anions minimized with only one added ammonia. In fact, if structures a3ni and a3nii are overlaid, the resultant three hydrogen bonds are the same as those seen in a3nn. In the neutral case, we do not see quite the same correspondence between models. In the large model, neither n5nni nor n5nnii, shown in Figure 5, has any hydrogen bond to the axial water, while that is the only hydrogen-bonding interaction in n5ni and n5nii. Further, in n3nn (Figure 5), the axial hydroxo group is positioned over the equatorial hydroxos and accepts a hydrogen bond, but in the smaller model (Figure 3, n3n), the axial hydroxo is rotated so that it is above the equatorial oxo and is donating a hydrogen bond. We also calculated the minimum structures and their energies for models with one external ammonia molecule that more closely resembled the large one. These calculations included models where the axial hydroxo either accepted or donated a hydrogen bond, but both cases were significantly higher in energy than the original models (data given in Supporting Information). Comparison of structures in Figure 5 to the depiction of the crystal structure in Figure 15 shows that in n3nn there was a large shift in the position of the guanidine to make the axial hydrogen bond. This shift did not occur in the n5nn cases. A similar rearrangement of R490 is observed in the H404A mutant form of the protein,34 suggesting that this large movement is connected to the lack of the imidazole ring interacting with the axially bound hydroxide. This shift probably occurred in n3nn in order to make the axial position more waterlike, which was not necessary for n5nn. The above behavior would seem to indicate that an axial water is preferred over an axial hydroxide, which agrees with results obtained using smaller models. Careful examination of Figure 4 shows that the axial hydroxo in n3ci accepts a very tight hydrogen bond from the external ammonium (2.61 Å), making that ligand quite waterlike and almost identical in energy to n5ci (in which the ammonium did not interact with the axial water at all). Conversely, n5ni was 0.8 kcal/ mol lower in energy than n3ni whose axial hydroxo was not waterlike at all. We have also investigated one anionic vanadate species with a cationic guanidinium, yielding an overall neutral structure. This complex (Figure 5, a1nc) was related to the anionic structure that rearranged to be tetrahedral in previous cases (Figure 2, a1). Hence, we again saw that the imidazole was lost, but like n1nii, the amine remained trans to the oxo but at a nonbonding distance. Remarkably, a1nc was less stable than n3nn by only 0.66 kcal/mol. (34) Macedo-Ribeiro, S.; Hemrika, W.; Renirie, R.; Wever, R.; Messerschmidt, A. J. Biol. Inorg. Chem. 1999, 4, 545-559.

Figure 7. Structure of active site, as in Figure 1, with proposed resting protonation state.

A stable minimum which contained a protonated guanidinium was found for one overall monocationic species (Figure 5, n5nc) where its proton was not transferred to the neutral metal cofactor. This is in agreement with the smaller models where protonation of the external hydrogen-bonding group led to lower-energy complexes than the corresponding cationic vanadate units. The structure of n5nc resembled n5nni, differing by longer V-O bond distances and shorter hydrogen bonds due to the strong Coulombic interaction of the oxo groups with the positively charged guanidinium residue. As in n5nni, we did not see a hydrogen bond to the axial water in n5n. The partial optimization of the largest scaffolds (Figure 6) resulted in very similar structures that are consistent with the smaller models. Even with the plethora of hydrogen bonds available to the vanadate, neither a highly charged unit nor an axial oxo was stable. Additionally, the axial hydroxo accepted a hydrogen bond from the nearby crystallographic water in both cases. The main difference between the two structures is that minimizing the overall neutral structure (Figure 6, a3nnc) resulted in a significant rotation of the oxygens within the equatorial plane which allowed the hydroxide to donate a hydrogen bond to the neutral guanidine. The six hydrogen bonds from the scaffold were redistributed among the cofactor ligands so that none of the interactions were lost. Thus, the net number of hydrogen bonds for the equatorial oxygens was seven, rather than six as seen in the crystal and overall cationic scaffold (Figure 6, a3ncc) structures. Since the minimized, cationic scaffold most closely resembles the crystal structure, it is likely that the active site is fully protonated. The discussion above may shed some light on to the resting state of the VCPO. These studies indicate that the vanadate cofactor is most likely anionic with several charges from the environment complementing the cofactor and creating a net neutral or cationic state for the active site. If the vanadate cofactor were neutral in the resting state, it would probably oscillate between the n3 and n5 type structures. Calculations indicate that addition of a proton to either state would likely result in protonation of a side chain rather than the cofactor. Further, these calculations suggest three things about the anionic cofactor. First, the axial position is not occupied by a terminal oxo. Second, in a stable state, there is most likely a hydroxo trans to His496 and another hydroxo in the equatorial plane (Figure 7) as in a3. Inorganic Chemistry, Vol. 43, No. 14, 2004

4135

Zampella et al. Last, the ideal position at which this cofactor is protonated would seem to be dependent upon the kinetics of the reaction. Protonation at the equatorial hydroxo would lead to an n4 type structure which is unstable and would lose the water molecule quite quickly forming a high-energy, four-coordinate vanadate. This would be the preferred protonation site if peroxide were already in the active site when the reaction begins. However, protonation at either the axial position or at an equatorial oxo yields n3 or n5 species, respectively, which appear close enough in energy in the protein that they probably interconvert. This species could be stable long enough for a peroxide molecule to enter the active site and displace the water to begin the catalytic cycle. In fact, protonation of the axially coordinated hydroxide group in the native state of the enzyme (n5nni-ii and n5nc) has been proposed to be a crucial step of the VCPO catalytic cycle.35 Conclusion The goal of this study was to further our understanding of the early intermediates of the vanadium dependent chloroperoxidase catalytic cycle. On the basis of crystallographic data, we generated computational models of the (35) Hamstra, B. J.; Colpas, G. J.; Pecoraro, V. L. Inorg. Chem. 1998, 37, 949-955.

4136 Inorganic Chemistry, Vol. 43, No. 14, 2004

active-site vanadate unit with increasing contributions from the surrounding protein scaffold. Investigation of various charge states and ligand arrangements suggest that an anionic vanadate cofactor, consisting of an axial hydroxo group with a hydroxo group and two oxo moieties in the equatorial plane, is most stable when complemented by the hydrogen bonds of the protein active site in its resting state. The initial protonation of the catalytic cycle may occur at the axial hydroxo or an equatorial oxo leading to reasonable intermediate structures of varying reactivity. Acknowledgment. The authors wish to thank the following funding agencies for support of this work: NIH GM42703 (V.L.P.), NIH GM65372 (H.A.C.), Beckman Young Investigator Program (H.A.C.), NIH GM008270 Michigan Molecular Biophysics Training Grant (J.Y.K.), NSF CHE 98-00184 (M.B.H.), and the Welch Foundation, Grant A-648 (M.B.H.). Supporting Information Available: Starting and final coordinates, as well as final energies, for all minimized complexes. Sample input file and figure showing the minimized structures for the nonbiologically relevant structures. This material is available free of charge via the Internet at http://pubs.acs.org. IC0353256