![]()
|
|
||||||||
J. Biol. Chem., Vol. 278, Issue 47, 47110-47118, November 21, 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




¶
From the
Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853 and the
Southern Research Institute, Birmingham, Alabama 35205
Received for publication, May 2, 2003 , and in revised form, August 15, 2003.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
|
-D-[2-deoxyribofuranosyl]-6-methylpurine (MeP-dR), 9-
-D-arabinofuranosyl-2-fluoroadenine (F-araA), and 2-fluoro-2'-deoxyadenosine (F-dAdo), which are not substrates for trimeric PNPs. The toxic purine analogs (6-methylpurine (MeP) and 2-fluoroadenine (F-Ade)) generated from these prodrugs readily diffuse across cell membranes, have high bystander activity not requiring cell-to-cell contact, and can kill the complete cell population even with gene expression in only 0.1-1% of the cells (15, 20, 21). In addition, MeP and F-Ade are toxic to both proliferating and nonproliferating tumor cells (22). These characteristics distinguish PNP anticancer gene therapy from the well studied herpes simplex virus thymidine kinase anti-cancer gene therapy strategy (23-29) and suggest that the selective generation of MeP or F-Ade could elicit considerable activity against solid tumors that have a low growth fraction. To optimize PNP anticancer gene therapy, the prodrug should be completely inert with respect to human enzymes and efficiently cleaved only by tumor cells expressing E. coli PNP. But whereas excellent in vivo antitumor activity has been demonstrated with MeP-dR, F-araA, and F-dAdo, these compounds all have dose-limiting toxicities. A better understanding of the structural basis for substrate specificity of the enzyme would allow modification of the prodrugs coupled with enzyme redesign to maintain activity levels while decreasing toxicity. Therefore, we have used crystallography to study the binding of MeP-dR, F-dAdo, and other nucleoside analogs to E. coli PNP.
In our previous work (10), the purine binding site was identified in complexes with R 1-P and 6-iodopurine, whereas Koellner et al. determined the structure of complexes with the inhibitors formycin B (FMB; 8-aza-7-deazainosine) (30) and formycin A (8-aza-7-deazaadenosine) (31). We report here the structure of E. coli PNP in association with 10 purine nucleosides (Fig. 1), which vary in both the base (hypoxanthine, 8-aza-7-deazahypoxanthine, adenine, 7-deazaadenine, 2-fluoroadenine, 6-methylpurine, and 6-methylthiopurine) and sugar (ribose, 2'-deoxyribose, arabinose, and xylose) portion of the molecule. We have also performed kinetic studies on the substrates and carried out computer modeling to examine the role of sugar modifications on reactivity. This information, together with knowledge about the structure of mammalian PNPs (32, 33), should aid in the design of more effective prodrugs that can be readily cleaved by E. coli PNP. Furthermore, an understanding of substrate analog binding should allow us to rationally change the enzyme to cleave novel prodrugs that are not cleaved by either human or wild type bacterial PNPs.
|
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
X-ray Intensity MeasurementsLigands were soaked into the crystals for 24-36 h, with concentrations ranging from 5 to 50 mM depending on solubility. Crystals were flash frozen at 100 K with liquid nitrogen, using 30% glycerol as a cryoprotectant, and typically diffracted to about 2.2-Å resolution. Data were collected at CHESS stations A1, F1, and F2. Various CCD x-ray detectors were used, including Area Detector Systems Quantum 1 single module and Quantum 4 mosaic CCD detectors and a Princeton Scientific Instruments 2k CCD detector. Exposure times ranged from 30 to 50 s per frame with a frame typically consisting of 1° of rotation. A total of 50-120° of data were collected depending on beam time constraints and the type of detector used. The data were processed using the programs DENZO (35) and SCALEPACK (36). Table I summarizes the data collection statistics.
|
calc) map based on SIGMAA (42) from the CCP4 package (43). Some loops could be built at the 0.8
level, but others have missing residues in the final model. Temperature factors in this region were much higher than for the rest of the molecule. Models of inosine (Ino), FMB, 7-deazaadenosine (tubercidin (TBN)), 2-fluoroadenosine (F-Ado), 9-
-D-ribofuranosyl-6-methylthiopurine (MTP-R), MeP-dR, and F-dAdo were directly constructed into the corresponding difference electron density. The ligand densities for adenosine (Ado), 9-
-D-arabinofuranosyladenine (AraA), and 9-
-D-xylofuranosyladenine (XylA) were less clear. These maps were improved by 3-fold averaging using the program RAVE (44). Water molecules were included after the ligand model building. The final models contain three protein chains, three ligand molecules, three ions modeled as phosphates, and several hundred water molecules. B-factors for the ligands sometimes varied considerably among monomers of the same complex. The lowest average ligand B-factors within a complex ranged from 18 (F-Ado monomer C) to 76 (AraA monomer C). In 2Fo - Fc maps based on the final model, density connectivity between sugar and base was observed for all monomers of all complexes except AraA and Ado. Refinement statistics are given in Table II.
|
Modeling of Intermediate FormationTo examine the role of sugar pucker in reaching the proposed oxocarbenium intermediate reaction state, computer modeling was carried out on the Ado, AraA, Ino, and XylA ligands using the GB/SA solvation model (45) and AMBER* force field implemented in MacroModel version 7.2 (46). Because the bottom of the energy well of the crystal structure appeared to be wide, a weak restraint of 1.0 kJ/mol-Å2 was applied to the substrate to prevent drifting away from the crystallographic position, whereas the protein was restrained at 200 kJ/mol-Å2, and the phosphate was left unrestrained due to software limitations. To mimic the formation of the oxocarbenium intermediate, the force field parameter for the glycosidic bond was adjusted to a bond length 50% greater than normal, and the O-4'-C-1'-N-9 and C-2'-C-1'-N-9 angle parameters were set to 90°. The structures were energy-minimized with these modified parameters to obtain an approximate model for how the ligand moves during the reaction.
Modeling of Ribose 1-PhosphateBecause only density for the adenine base was observed in an attempted soak including R 1-P, the DockVision program (47) was utilized to position a R 1-P molecule into the E. coli PNP active site. Both the Monte Carlo simulation option and genetic algorithm within DockVision were used, but only the Monte Carlo simulation gave reasonable results. These calculations were performed using the Cornell Theory Center work station cluster. Four Monte Carlo simulations and two genetic algorithm runs were performed with a distance cut-off of 8.0 Å. For the four Monte Carlo simulations, the numbers of trials were set at 10,000, 20,000, 10,000, and 20,000, and dielectric constants were set at 1, 1, 2, and 4, respectively. For genetic algorithm runs, the number of generations was set at 10,000, and the dielectric constant was set at 1 or 2. The results of the genetic algorithm runs were inconsistent with no convergence. The Monte Carlo simulation gave much better results, and the run with the highest dielectric constant resulted in a cluster of 26 R 1-P structures located near the expected binding site. Of these 26 models, the one with the lowest energy had the best hydrogen bonding interactions and is the one described in this paper.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
-sheet core flanked by eight
-helices. The core can be divided into a large eight-stranded
-sheet and a smaller five-stranded
-sheet that pack together to form a distorted
-barrel, an arrangement also found in human PNP (2), bovine PNP (32, 33), and human 5'-deoxy-5'-methylthioadenosine phosphorylase (48), although there is little sequence homology between the E. coli enzyme and the trimeric class mammalian enzymes (10).
Active Site of E. coli PNPThe previously determined structures of E. coli PNP (10, 30, 31) identified the residues involved in substrate binding and catalysis (Fig. 3). The purine binding site consists of Ala156, Phe159, Val178, Met180, Ile206, and Asp204. The first four of these residues form a hydrophobic pocket around the purine base. Phe159 is located between the purine base and the hydrophobic face of the sugar and makes an angle of
60° with the plane of the purine ring. Met180 also lies between the purine base and the hydrophobic face of the ribosyl group. It has been proposed that in E. coli PNP, a protonated Asp204 stabilizes the transition state, in which electron density from the weakening bond is transferred to the purine ring, by donating a hydrogen bond to the purine N-7 atom (10). Evidence for protonation of Asp204 comes from fluorescence studies on the inhibitor FMB, which exists in tautomeric equilibrium between N-7-H and N-8-H forms (IUPAC conventions for ring numbering systems vary among bases; for consistency, we use only the purine ring numbering as shown for inosine in Fig. 1). Although the N-8-H form (which could hydrogen-bond to a protonated Asp204) is the minor (
20%) tautomer in solution (30, 49), fluorescence studies have detected an equilibrium shift in favor of this tautomer upon binding to E. coli PNP (49). Also consistent with the proposed role for Asp204, mutating this residue to Ala results in a
100-fold reduction in activity on MeP-dR and F-araA.2
|
The ribose binding site consists primarily of interactions with Glu181 and His4# (residues from an adjacent monomer are designated by # throughout). The Glu181 side chain accepts hydrogen bonds from the 2'- and 3'-hydroxyl groups, and His4# accepts a hydrogen bond from the 5'-hydroxyl group. The phosphate group hydrogen bonds with the 3'-hydroxyl and, in some structures, the 2'-hydroxyl.
Structures of E. coli PNP ComplexesUnlike mammalian PNPs, which are highly specific for 6-oxopurine nucleosides (or 2'-deoxynucleosides) (50, 51), E. coli PNP accepts 6-aminopurine nucleosides as substrates and shows some activity for purine nucleosides with ribosyl modifications (52, 53). The structures of the nucleosides used here (Fig. 1) include a variety of modifications to the purine bases and ribosyl group, but all are
-D-furanosides. The ligands bind to the active site in the high syn conformation about the glycosidic bond and generally adopt a C-4'-endo sugar pucker (Fig. 4A). Strong electron density shows the binding of each ligand, except AraA and Ado, where only weak density with poor sugar/base connectivity is visible. The AraA complex also shows extremely high ligand B-factors (>75 Å2 in each monomer), suggesting low ligand occupancy. Table II shows the kinetic constants for seven substrates (no cleavage of XylA by E. coli PNP was detected). The Ki values for FMB and TBN were previously determined to be 5 µM (54) and 120 µM (55), respectively. The collection of structures along with the kinetic data allows us to examine the roles of substituents at positions 1, 2, 6, and 7 of the purine base and the effects of changing stereochemistry at positions 2' and 3' of the sugar.
|
Ribose normally makes two hydrogen bonds to Glu181, one from the O-2'-hydroxyl group to one of the carboxylate oxygen atoms and one from the O-3'-hydroxyl group to the other carboxylate oxygen atom. The 3'-hydroxyl group of 2'-deoxy ligands can interact with both carboxylate oxygen atoms, as in MeP-dR (at 3.1 Å), or maintain only the usual single (but stronger, 2.6-Å) hydrogen bond, as in F-dAdo. Similar B-factors were observed for both compounds, and the Km for F-Ado was only marginally lower than for F-dAdo, suggesting that loss of the O-2'-hydroxyl group may not significantly affect substrate binding. When the glycosidic bond lengthens during the reaction, flattening of the C-4'-O-4'-C-1'-C-2' torsion requires that either the C-1' atom or the base or both must move. We investigated this motion in ribose substrates by energy-minimizing the crystal structures of Ado and Ino using a modified AMBER force field with parameters forcing the nucleosides to adopt a conformation mimicking the transition state between the substrate and oxocarbenium intermediate. For Ado, the crystallographically determined torsion angles in the three monomers range from 20 to 23°. Restrained energy minimization with parameters requiring a 50% increase in the glycosidic bond length and 90° O-4'-C-1'-N-9 and C-2'-C-1'-N-9 angles caused the C-1' atom to move toward the plane of the other atoms, giving a final torsion of 7°, even without explicitly changing this torsion parameter to require flattening. With Ino, the torsion value decreased from 14-16° to 8°, confirming that flattening of the torsion is a natural response to lengthening the glycosidic bond.
In the case of XylA, the crystallographically determined sugar pucker is C-3'-exo (Figs. 4C and 5B), which allows XylA to maintain hydrogen bonding between Glu181 and the C-2'-and C-3'-hydroxyl groups (Fig. 6E). The C-3'-exo pucker means that C-4', O-4', C-1', and C-2' are nearly coplanar before the reaction begins (with a torsion angle of 2°), making glycosidic bond lengthening incompatible with the required planar geometry. Upon AMBER energy minimization, the torsion value increased to 22°, suggesting that XylA is not a substrate because its pucker makes the oxocarbenium intermediate energetically inaccessible.
|
|
Role of Asp204 in Binding and CatalysisThe purine binding site is largely hydrophobic, but N-1, N-7, and O-6 or N-6 typically form hydrogen bonding interactions with protein side chains or water molecules. The proposed catalytic mechanism suggests that a hydrogen bond between Asp204 and N-7 should be observed in complexes with substrates. This hydrogen bond was not observed in the previous study with the inhibitor FMB (30) (but in our FMB complex, one subunit has a 2.9-Å hydrogen bond), but it was present in a complex with a formycin A analog (31). These structures left open the question of how closely purine base binding matches the binding of formycins. A comparison of the 10 present complexes shows that the overall ligand binding geometry is similar for both 6-oxo and 6-amino bases (Fig. 5), including formycin B. Although there is clearly only one major mode of base binding, the minor variations are enough to affect some contacts made by the base.
Most significantly, good substrates do not necessarily show an Asp204/N-7 hydrogen bond, despite the evidence from mutation studies that Asp204 is critical for catalysis. Although the natural substrates Ado and Ino show a hydrogen bond in one subunit each, no hydrogen bond is observed in any subunit for MeP-dR and MTP-R, both of which have similar catalytic efficiencies with PNP (Table II). In the C monomer, MeP-dR is displaced deeper into the active site pocket, farther from the active site opening and phosphate group. The 4.9-Å distance between the phosphate and sugar C-1' suggests that this conformation would not lead to catalysis. In the A monomer and B monomer (Fig. 6A), MeP-dR binds with the general conformation observed in all monomers of all of the other complexes.
A significant difference relative to other complexes is observed in all three monomers of the MTP-R complex, where the Asp204 side chain rotates toward the active site opening and cannot interact with the purine N7 (Fig. 6B). The change in position is probably due to steric conflicts between the normal Asp204 rotamer and the MTP-R sulfur atom. A similar displacement of Asp204 may be envisioned for binding 7-methyl analogs, which have previously been confirmed as substrates (57). However, these analogs already carry a positive charge at N-7, whereas MTP-R would presumably still require a nearby proton donor for catalysis. Although a water molecule is observed hydrogen-bonding with N-7 and the side chain of Ser203, Asp204 is still required, since mutations to Asn and Ala both prevent cleavage of MTP-R.2 Taken together, the structures and mutation data suggest that Asp204 is not critical for initial ligand binding but must move in concert with the substrate during catalysis to donate a proton to N-7 as the glycosidic bond breaks. Movement of Asp204 during the reaction is consistent with the work of Koellner, who proposed that Asp204 also changes conformation after protonation of N-7 (31).
Structurally Conserved Waters at N-1In 6-oxo-purines, the N-1 position is protonated, whereas in 6-aminopurines the N-1 position is unprotonated. In the PNP complexes, the N-1 position interacts with a series of water molecules, the first of which can be either a donor or acceptor for the purine N-1 position and for the second water molecule (Fig. 7). In complexes with substrates, the first water molecule must be a hydrogen bond donor to the backbone carbonyl group of Phe159, whereas the first water molecule in some of the inhibitor complex subunits is too far from this oxygen atom to form a hydrogen bond. The second water molecule must donate a hydrogen bond to the backbone carbonyl of Leu158 but can be either a donor or acceptor to a third water molecule, which is observed
50% of the time. In some cases, an additional water molecule bridges to the side chain of Asp112. The residues lining the water channel are highly conserved.
|
Effect of C-2 SubstituentsSubstitutions at C-2 can be made without affecting the water channel, as observed in both the F-Ado (Fig. 6C) and F-dAdo (Fig. 6D) complexes. E. coli PNP also accepts 2-amino substrates such as guanosine. The environment for the C-2 substituent is hydrophobic, consisting of Ala156, Val178, Met180, and Phe159. With 2-fluoro compounds, the fluorine atom packs against the Met180
-carbon and sulfur, providing more favorable van der Waals interactions than the hydrogen atom of Ado. The improved steric fit may account for the stronger electron density and significantly lower B-factors observed for F-Ado and F-dAdo relative to Ado.
Role of Arg24Because the crystals were grown from high concentrations of ammonium sulfate, and because sulfate is a known inhibitor of PNP, a sulfate ion probably occupies the phosphate binding site in the E. coli PNP complexes. However, because PNP has two binding affinities for phosphate, one of which has a Kd much lower than that for sulfate (58), the identity of the bound ion cannot be conclusively determined. The phosphate binding site consists mainly of two arginine residues (87 and 43#), backbone interactions with Gly20, and side chain interactions with Ser90. Arg24 sometimes participates as well. This positively charged cavity is buried within the active site. Upon ligand binding, the phosphate site is further shielded by the nucleoside sugar.
Arg24 may provide an explanation for the existence of two difference binding states for phosphate. Koellner observed small differences around Arg24 in open and closed protein conformations, with higher side chain B-factors and more diffuse density in the open form and stabilizing interactions between Arg24 and a moved helix in the closed form, and suggested that these differences accounted for the two different phosphate affinities (31). We consistently observe a more significant variation; in all of the complexes except AraA (where electron density indicates very poor nucleoside binding), the Arg24 guanidinium group interacts with the phosphate in only one of the three independent active sites and is far away (up to 8 Å, with clear density indicating the different conformation) in the other two. Although movement of Arg24 is the most promising structural explanation for observed phosphate affinities, note that sulfate (present at high concentrations for all E. coli PNP complexes published to date) shows only a single binding affinity (58).
Modeling Studies of Ribose 1-PhosphateIn earlier studies of E. coli PNP using monoclinic crystals (10), attempts to observe bound R 1-P were unsuccessful because of low occupancy. In the current work, high sulfate concentrations prevent R 1-P binding. Therefore, we have used modeling studies to better understand R 1-P binding.
The docking studies predict that R 1-P adopts a C-2'-exo conformation. This is the most common ribose pucker and was observed in the bovine PNP/hypoxanthine/R 1-P crystal structure (33). The C-3'-C-4'-C-5'-O-5' torsion angle is 174.3°, the C4'-O4'-C1'-O torsion angle is 102.8°, and the O-4'-C-1'-O-P torsion angle is 63.9°. Comparison with the Ino complex shows that the ribose ring conformation flips after cleavage such that the C-1' atom moves away from the purine base and about 0.9 Å toward the phosphate binding site (Fig. 8). The phosphate group moves 1.0 Å toward the sugar ring. The purine base moves 1.1 Å along the line of the glycosidic bond but retains the same general interactions. R 1-P maintains hydrogen bonds with Gly20, Arg87, Ser90, Arg43#, and Glu181, but the O-5' contact with His4# is lost. The major movement of the ribosyl group is a rotation about a line drawn from O-4' to the midpoint of the C-2'-C-3' bond. In contrast to trimeric bovine PNP, in which almost all movement is in the sugar (33), the modeling predicts that structural changes necessary to break the bond in the E. coli enzyme are shared roughly equally by the sugar, base, and phosphate. Based on the present structures, the movement of the base would be required for some substrates (such as MTP-R) to accept a proton from Asp204.
|
Because the enzyme active site allows little room for modification of the sugar, we have also begun to turn our attention to protein engineering as a way of improving catalytic efficiency. By matching a given substrate to an engineered PNP, it may be possible to optimize the conversion of the prodrug. This series of 10 complexes suggests possible strategies for achieving this goal. Residues near the sugar C-5' position (Met64, Met180, and Ile71#) are believed to be important primarily for hydrophobic packing rather than catalysis and might be mutated to smaller residues and tested with modified versions of MeP-dR, F-dAdo, or F-AraA with bulky substituents added to C-3' or C-5'. In addition, the structure of PNP complexed with nonstandard sugars suggests residues that could be modified to improve the cleavage of arbinonucleosides. Finally, combinations of the above approaches may be considered. Modeling and crystallographic studies are currently under way to explore possible PNP/prodrug pairs based on these ideas.
| FOOTNOTES |
|---|
* This work was supported by National Institutes of Health (NIH) National Cooperative Drug Discovery Grant CA-67763. This work is based upon data measured at the Cornell High Energy Synchrotron Source (CHESS), which is supported by National Science Foundation Grant DMR-9311772, using the Macromolecular Diffraction at CHESS (MacCHESS) facility, which is supported by NIH Grant RR-01646. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
¶ To whom correspondence should be addressed: Dept. of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853. Tel.: 607-255-7961; Fax: 607-255-1227; E-mail: see3{at}cornell.edu.
1 The abbreviations used are: PNP, purine nucleoside phosphorylase; Ado, adenosine; AraA, 9-
-D-arabinofuranosyladenine; F-Ado, 2 fluoroadenosine; F-Ade, 2-fluoroadenine; F-dAdo, 2-fluoro-2'-deoxyadenosine; F-araA, 9-
-D-arabinofuranosyl-2-fluoroadenine; FMB, formycin B (8-aza-7-deazainosine); Ino, inosine; MeP, 6-methylpurine; MeP-dR, 9-
-D-[2-deoxyribofuranosyl]-6-methylpurine; MTP-R, 9-
-D-ribofuranosyl-6-methylthiopurine; R 1-P, ribose 1-
-D-phosphate; TBN, tubercidin (7-deazaadenosine); XylA, 9-
-D-xylofuranosyladenine. ![]()
2 E. M. Bennett, C. Li, P. W. Allan, W. B. Parker, and S. E. Ealick, unpublished data. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. A. McCarthy and J. G. McCarthy The Structure of Two N-Methyltransferases from the Caffeine Biosynthetic Pathway Plant Physiology, June 1, 2007; 144(2): 879 - 889. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. L. Lorenzi, C. P. Landowski, A. Brancale, X. Song, L. B. Townsend, J. C. Drach, and G. L. Amidon N-METHYLPURINE DNA GLYCOSYLASE AND 8-OXOGUANINE DNA GLYCOSYLASE METABOLIZE THE ANTIVIRAL NUCLEOSIDE 2-BROMO-5,6-DICHLORO-1-({beta}-D-RIBOFURANOSYL)BENZIMIDAZOLE Drug Metab. Dispos., June 1, 2006; 34(6): 1070 - 1077. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Zang, W.-H. Wang, S.-W. Wu, S. E. Ealick, and C. C. Wang Identification of a Subversive Substrate of Trichomonas vaginalis Purine Nucleoside Phosphorylase and the Crystal Structure of the Enzyme-Substrate Complex J. Biol. Chem., June 10, 2005; 280(23): 22318 - 22325. [Abstract] [Full Text] [PDF] |
||||