1 Department of Microbiology and Molecular Cell Biology, Eastern Virginia Medical School, Norfolk, VA 23507 and 2 Department of Pharmaceutical Sciences, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Keywords: ganciclovir/gene therapy/thymidine/thymidine kinase/thymidylate kinase
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Because of its long history as the chemotherapeutic target for antiviral agents, the structure/activity of HSV-1 TK has been extensively studied via characterization of drug-resistant viral strains (Larder et al., 1983; Darby et al., 1986
), site-directed or random sequence mutagenesis (Liu and Summers, 1988
; Munir et al., 1992
; Michael et al., 1994
; Black and Loeb, 1996
), photoaffinity labeling (Rechtin et al., 1995
, 1996
) and X-ray crystallography (Brown et al., 1995
; Champness et al., 1998
). As a result of these studies, mutant HSV-1 TKs have been generated that display selective preferences for different HSV-1 TK substrates. For example, we had previously described mutations at Gln125 to Asn that resulted in poor utilization of thymidine as a substrate, but retained the ability to phosphorylate GCV (Drake et al., 1999
; Hinds et al., 2000
). Another approach has been the use of random oligonucleotide mutagenesis to produce drug-specific HSV-1 TK mutant enzymes (Black and Loeb, 1996
; Black et al., 1996
). Two of these HSV-1 TK mutants, termed TK30 and TK75, were highly sensitive to GCV when expressed in mammalian cells, resulting in more cell killing compared with cells expressing the wild-type HSV-1TK (Black et al., 1996
; Kokoris et al., 1999
). In this study, we hypothesized that the combination of the Asn mutations at Gln125 and the mutated regions in TK30 or TK75 would create a GCV-specific enzyme by eliminating the endogenous substrate competition of thymidine and improving the efficacy of GCV phosphorylation.
Based on this hypothesis, 16 mutant HSV-1 TK enzymes were generated using site-directed mutagenesis. These mutants and also wild-type HSV-1 TK, were expressed and partially purified from Escherichia coli, then screened for enzymatic activity with thymidine and GCV. Mutant HSV-1 TK enzymes deficient in thymidine phosphorylation activities, but active for GCV phosphorylation, were further characterized for their kinetic properties and expressed in human colon carcinoma cell lines for cell killing and drug metabolism studies. The cumulative substrate binding information generated from the enzymatic and cell culture studies were used in constructing an HSV-1 TK molecular model of active-site/substrate binding interactions with FlexiDock, a SYBYL 6.7-based program. Of the mutant enzymes evaluated, two were identified that possessed unique substrate properties consistent with the goal of generating HSV-1 TKs with minimal thymidine and maximal GCV phosphorylation activities. The correlation between the functional enzymatic and cell culture data and the modeled ligand/protein structures for these HSV-1 TK mutants is discussed and has resulted in a new predictive structural modeling approach for generating substrate specific HSV-1 TKs.
![]() |
Materials and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
All reagents and nucleotides were purchased from Sigma Chemical unless specified otherwise. [Methyl-3H]thymidine (69 Ci/mmol) and [8-3H]ganciclovir (17.4 Ci/mmol) were purchased from Moravek Biochemicals. DE81 filter-paper disks, isopropyl-ß-D-thiogalactopyranoside (IPTG), dithiothreitol (DTT), heparin resin and DE-52 resin were purchased from Fisher Scientific. Automated DNA sequencing of plasmids was done using a Model 377 DNA sequencer from Applied Biosystems. Oligonucleotide primers were synthesized by Clontech. Restriction endonucleases and T4 DNA ligase were purchased from New England Biolabs.
Site-directed mutagenesis
Site-specific mutations were introduced into the flanking regions of conserved sequences within the HSV-1 TK gene using the PCR-based, QuikChange Site-directed Mutagenesis Kit as per the manufacturers instructions (Stratagene). Four pairs of mutagenic primers were used to alter the amino acid sequence of three genes previously subcloned into a pET9a (Promega) expression construct: wild-type HSV-1 TK, Q125N TK and Q125E TK (Hinds et al., 2000). Two different sets of mutations, TK30 and TK75 (Black et al., 1996
), were introduced at the flanking regions of conserved site 3 and site 4. Each of the site 3 (30-3 and 75-3) and site 4 (30-4 and 75-4) mutants were generated independently. For mutants incorporating both site 3 and site 4 mutations, the corresponding site 3 mutant was used as the template. Specifically, in the TK30 set, three amino acid substitutions at site 3, L159I, I160L and F161A, and two mutations at site 4, A168Y and L169F, were introduced. These amino acid substitutions were incorporated into the vector templates using the following primers: CCGGCCCTCACC ATCCTGGCTGACCGCCATCCC and GGATGGCGGTCAGCCAGGATGGTGAGGGCCGG for site 3 (30-3) and CATCCCATCGCCTACTTCCTGTGCTACCCG and GTAGGGTAGCGGATGAAGCAGACGATGGGC for site 4 (30-4). In the TK75 set, two amino acid substitutions at site 3, I160L and F161L, and two substitutions as site 4, A168V and A169M, were introduced. These substitutions were incorporated into the template vectors using the following sets of mutagenic primers: CCGGCCCTCACCCTGCTGCTGGACCGCCATCCC and GGGATGGCGGTCCAGCAGCAGGGTGAGGGCCGG for site 3 (75-3) and CATCCCATCGCCGTGATGCTGTGCTACCCG and CGGGAGCACAGCATCACGGCGATGGGATG for site 4 (75-4). The mutant plasmids generated by this procedure were transformed into JM109 competent E.coli cells. Positive transformants were identified by growth on LB medium containing kanamycin (50 µg/ml). Plasmid DNA was purified from positive transformants using the Qiagen Mini Prep kit as per the manufacturers instructions. To confirm the presence of the specified mutated residues and fidelity of the coding sequence of HSV-1 TK, automated sequencing of each construct was done using a primer directed at the T7 promoter located on the pET-9a expression vector.
Enzyme expression and purification
Following sequencing analysis, the mutant plasmids were transformed into competent E.coli SY211 TK cells and expression induced with IPTG as described previously (Rechtin et al., 1995; Hinds et al., 2000
). Soluble protein lysate was loaded on to a DE 52 cellulose column coupled to a heparin column and purified as described previously (Hinds et al., 2000
). This procedure results in 0.53 mg of total protein with a minimum purity of 70%. The mutant enzymes can be stored at 80°C for at least 6 months.
Photoaffinity labeling of lutant HSV-1 TK enzymes
The TMP photoaffinity analog [32P]-5-azido-dUMP (5N3dUMP), was synthesized enzymatically using HSV-1 TK, 5-azidodeoxyuridine and [32P]ATP as described previously (Rechtin et al., 1995
). For photocrosslinking, 5 µg of the partially purified enzymes were incubated with 10 µM [32P]-5N3dUMP, 10 mM MgCl2 and 40 mM KCl for 10 s, followed by UV irradiation for 90 s with a hand-held UV lamp (254 nm UVP-11, Ultraviolet Products) at a distance of 3 cm. Samples were separated by SDSPAGE and subsequent autoradiography as described previously (Rechtin et al., 1995
).
Enzymatic assays
Partially purified HSV-1 TK mutants were screened in triplicate for thymidine and GCV phosphorylation using a standard reaction mix, 20 mM potassium phosphate, pH 7.6, 1 mM DTT, 5 mM ATP, 5 mM MgCl2, 25 mM NaF, 40 mM KCl and 0.5 mg/ml bovine serum albumin (BSA) in a total volume of 25 µl and the following linear velocity conditions established for wild-type HSV-1 TK: for thymidine, 10 ng enzyme, 2 µM [3H]thymidine, 10 min incubation at 37°C; for GCV, 250 ng, 40 µM [3H]GCV, 20 min incubation at 37°C (Hinds et al., 2000). The entire reaction mixture (25 µl) was spotted on to an ion-exchange filter, DE-81, washed and quantitated by scintillation counting as described previously (Hinds et al., 2000
). Control reactions incubated without enzyme were spotted and processed to correct for backgrounds. These mutant enzymes were also screened for thymidylate kinase (TMPK) activities. For these experiments, 20 µM [3H]TMP was substituted for deoxythymidine in the above assay mixture and [3H]TDP product was separated from [3H]TMP on thin-layer PEI-cellulose plates developed in 0.4 M LiCl (Hinds et al., 2000
). The TDP product was quantified by scintillation counting of the excised TMP and TDP spots.
Those HSV-1 TK mutant enzymes that were active and demonstrated unique pharmacological properties were characterized further by determining the kinetic parameters (Km, kcat, Vmax) for each enzyme. The protein concentrations, substrate ranges and incubation times used in determining the Km of thymidine and GCV for each mutant were done with linear velocity conditions established for each enzyme. Each Km determination was assayed in triplicate using the following conditions: for [3H]thymidine, wild-type HSV-1 TK, 10 ng, 0.255 µM, 15 min; Q30 TK, 500 ng, 260 µM, 10 min; Q75 TK, 50 ng, 0.540 µM, 10 min; Q75-3 TK, 100 ng, 0.540 µM, 5 min; N7530 TK, 400 ng, 10100 µM, 20 min; and Q7530 TK, 100 ng, 10100 µM, 5 min for [3H]GCV, wild-type, 150 ng, 20100 µM, 20 min; Q30 TK, 300 ng, 20120 µM, 20 min; Q75 TK, 100 ng, 10100 µM, 20 min; Q75-3 TK, 200 ng, 270 µM, 20 min; N7530 TK, 400 ng, 20 min , 40200 µM; and Q7530 TK, 100 ng, 10100 µM, 20 min. The kinetic values (Km, Vmax and kcat/Km) for each enzyme and substrate were determined using Microsoft Kalidegraph software version 3.6.
Pseudo-first-order rate constants for wild-type HSV-1 TK and N30-3 TK were determined using thymidine and GCV as substrates. Using the filter-binding assays described above, the experiments were carried out by preincubating a 50-fold excess of wild-type HSV-1 TK or N30-3 TK with 2 nM thymidine5 µl [3H]thymidine or 80 nM GCV5 µl [3H]GCV. Equal volumes of the enzymesubstrate mixture were added with assay buffer and incubated for a total of 10 min for thymidine or 20 min for GCV. Aliquots were taken at time zero as a baseline, then every 1530 s and spotted on a DE81 filter. The filters were washed as above and radioactivity was quantitated by scintillation counting of the washed filters. The resulting data points were plotted and fitted to a best-fit single exponential with the Kalidegraph software, from which the pseudo-first-order rate constants were obtained. Wild-type HSV-1 TK had a pseudo rate of 0.5 nM/min for thymidine and 0.8 nM/min for GCV with product yields of 65 and 70%, respectively (data not shown). The rates for N30-3 TK were 1.2 nM/min for thymidine and 0.4 nM/min for GCV; however, the product yields were lower at 3 and 16%, respectively (data not shown).
HSV-1 TK expressing cell lines
The HCT116 and SW620 human colon carcinoma cell lines were maintained in RPMI 1640 supplemented with L-glutamine and 10% fetal bovine serum (Cellgro). These two cell lines were obtained from American Type Culture Collection (Rockville, MD). To generate a Moloney murine leukemia virus-derived plasmid for each mutant, the unique BspEIMluI restriction fragment that contains the mutated flanking regions of site 3 and site 4 were excised from the respective mutant pET-TK plasmids and ligated into the same site of pLENTK (McMasters et al., 1998) as described previously (Hinds et al., 2000
). Each new pLEN-mutant TK construct was sequenced to confirm the presence of the mutated residues. Following sequence analysis, wild-type HSV-1 TK, Q30 TK, Q75 TK, Q75-3 TK, N30-3 TK, Q7530 TK and N7530 TK were transfected into the HCT-116 cell line using lipofectin reagent (Gibco-BRL), 2 µg of plasmid, 15 µl of lipid and 3x105 cells. Transfected cell lines were maintained in their respective growth medium and placed under G418 selection (200600 µg/ml) for 2 weeks. For each HSV-1 TK mutant, at least five positive clones were screened, selected and expanded based on the results of [3H]GCV metabolic labeling, GCV sensitivity (10 µM) and western blot analysis (Hinds et al., 2000
). Similar procedures were used to generate the wild-type HSV-1 TK and Q7530 TK expressing SW620 cell lines. The resulting clonal isolates from each mutant HSV-1 TK cell line were screened by western blot analysis to create a panel of representative cell lines that expressed their respective HSV-TK isoform atsimilar protein levels to one another (Hinds et al., 2000
).
Metabolic labeling assays
For metabolic labeling with [3H]GCV, the same general procedure was followed for each HSV-TK expressing cell line in the panel: cells were plated (5x105/well) in triplicate in a six-well plate. The following day, cells were incubated for 18 h in 3 ml of medium containing 2 µCi [3H]GCV (17.4 Ci/mmol) (Moravek) and 4 µM unlabeled GCV. The next day, the cells were rinsed twice with PBS, trypsinized, pelleted and the nucleotides were extracted with 0.150 ml of 70% methanol on ice for 10 min (Hinds et al., 2000). For metabolic labeling with [3H]thymidine, confluent cell monolayers in six-well plates were incubated for 1 h in 3 ml of medium containing 2 µCi [3H]thymidine (69 Ci/mmol) and 2 µM unlabeled thymidine prior to methanol extraction. If the thymidine labelings were performed on sub-confluent, dividing cultures, the metabolite numbers would reflect growth rates and therefore cellular kinase activities, not specific to HSV-1 TK activities. For both labeling experiments, an aliquot of each methanol-soluble supernatant was analyzed for radioactivity by scintillation counting.
Cell viability and clonal dilution assays
For MTT cell viability assays, cells expressing wild-type, Q30 TK, Q75 TK, Q75-3 TK, N30-3 TK and N7530 TK were plated in 0.1 ml (5000 cells/well) in a 96-well microtiter plate. The following day, GCV was added to the cells within the dose range 0.0110 µM. After 4 days, 3-(4,5-dimethylthiazol-2-yl)-2,5-dephenyltetrazolium bromide (MTT) was added to each well (50 µg/well) for 3 h followed by the addition of dimethyl sulfoxide to solubilize the precipitated crystals. Cell viability was determined by using a plate reader measuring optical density at 540 nm. For clonal dilution assays, cell lines were seeded (105 cells/well) in 24-well plates in 1 ml of medium. The following day, 0.0001, 0.001, 0.01, 0.1 or 1 µM GCV was added to the appropriate wells in triplicate. After 24 h, the cells in each well were sequentially diluted from 1:10 to 1:10 000 in 1 ml of fresh medium as described previously (Drake et al., 1999). The resulting cell colonies were counted after 7 days.
Molecular modeling of substrate interactions with mutant HSV-1 TKs
All of the molecular modeling analyses were performed using SYBYL version 6.6/6.7 (Tripos, St. Louis, MO) on Silicon Graphics R4400 and R5000 computers in the UAMS Biomedical Visualization Laboratory. The coordinates of two X-ray crystal structures, HSV-1 TK bound with thymidine and HSV-1 TK bound with GCV (Champness et al., 1998), were downloaded from the Brookhaven protein database and loaded into the FlexiDock module of the molecular modeling program SYBYL 6.7/UNITY 4.2. The FlexiDock program consists of two main components, a genetic algorithm and an energy evaluation function for scoring the resulting solutions. The genetic algorithm continually evolves possible ligandreceptor interactions selecting for the structure with the lowest energy form. The energy evaluation scoring functions in the program are based on the Tripos force field, van der Waals interactions, electrostatic interactions and torsional and constraint energy terms. For simplicity, one monomer was used in these studies. Initially, all ligands and water molecules, except for four water molecules important in thymidine binding to the active site, were deleted from the wild-type structures before the amino acid substitutions for Q75-3 TK, Q30-4 TK and Q7530 TK were introduced separately. Next, the active site, designated by a 6 Å radius around Y101, site 3 (L159H164) and site 4 (Y168P173), was fully minimized using an annealing function that relaxes the local strain introduced by the amino acid substitutions. Using the FlexiDock program, the ligandreceptor complex for wild-type thymidineTK and GCVTK were optimized according to the following parameters: initial population, generation number and flexible bonds. For thymidine, the C
Cß bond for I160, F161, D162, A168, L169 and L170, the C
Cß bond and the OH side chain for Y101, Y172 and Y132, the C
C
bond for R163 and Q125 and the CßC
bond for E83 and E225 were marked as flexible. For, GCV, the C
Cß bond for I160, F161, D162, A168, L169 and L170, the C
Cß bond and OH side chain for Y101, Y132 and Y172 and the C
C
bond for R163 and Q125 were identified as flexible. The resulting ligandreceptor geometries obtained from a FlexiDock run using the above parameters were tested against the ligandreceptor geometry seen in the crystal structure. By using these parameters and the wild-type HSV-1 TK structure, the FlexiDock program was able to match the ligandprotein geometry and hydrogen-bonding pattern seen in the crystal for both substrates with root mean square values of 0.5048 and 0.5750 for thymidine and GCV respectively. These parameters for thymidine and GCV were then used to evaluate ligand binding in the active sites of mutant HSV-1 TK enzymes. An initial population of 2000 individuals, analyzed for 50 000 generations, was chosen as the optimal conditions for analyzing thymidine and GCV binding in the active site of HSV-1 TK. Because thymidine and GCV have different binding properties, the hydrogen donoracceptor pairs and the flexible bonds varied. On average, these conditions required at least 15 h to complete one FlexiDock analysis per structure. Further details about this algorithm will be provided freely by the authors upon request.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
To introduce amino acid substitutions in the flanking regions of site 3 and site 4 in HSV-1 TK corresponding to the TK30 and TK75 enzymes (Black et al., 1996), the coding sequences of three HSV-1 TK genes were used encoding wild-type HSV-1 TK, Q125N TK and Q125E TK that had been previously cloned into a bacterial expression vector, pET-9a (Hinds et al., 2000
). Using the QuikChange site-directed mutagenesis protocol, 16 HSV-1 TK mutants were developed, expressed and purified from E.coli as described in Materials and methods (Hinds et al., 2000
). Initially, the HSV-1 TK mutants containing the full set of mutations, Q30 TK, Q75 TK, N30 TK, N75 TK, E30 TK and E75 TK were assayed for phosphorylation of thymidine (2 µM) and GCV (40 µM) using enzymatic conditions optimized for wild-type HSV-1 TK (Table I
). For comparative purposes, the results obtained for each mutant were normalized as a percentage relative to wild-type HSV-TK activities (set at 100%). As expected, wild-type HSV-1 TK and Q75 TK effectively phosphorylated both substrates and Q30 TK was more GCV-specific. However, substitution of Gln125 with Asn or Glu into Q30 TK or Q75 TK led to comparatively low phosphorylation activities for both substrates for the N30, E30, N75 and E75 HSV-1 TK mutants.
|
Following the activity screens, Km (appt) and kcat (appt) for thymidine and GCV were determined for each of the following HSV-1 TK mutant enzymes: Q30 TK, Q75 TK, Q75-3 TK, Q7530 TK, N30-3 TK and N7530 TK. Linear velocity conditions for each individual enzyme and substrate were established prior to Km (appt) determinations (see Materials and methods). The kinetic data obtained from each experiment were analyzed using a double-reciprocal plot and the Km (appt), kcat and kcat/Km (appt) values were derived from these graphs and are compiled in Table II. For each mutant evaluated, the Km (appt) values for thymidine were 7100 times higher than those for wild-type HSV-1 TK. Of note for GCV, the Q75-3 TK and Q7530 TK enzymes had lower Km (appt) values of 34 and 40 µM, respectively, combined with the higher Km (appt) values of 15 and 53 µM, respectively, for thymidine. The Q7530 enzyme actually has a lower Km (appt) for GCV than thymidine. The determined kcat values were also informative, in that for thymidine there were relatively minor changes compared with wild-type HSV-1 TK, suggesting alteration of thymidine binding as the effect of the mutations. In contrast, there was a large range (20100-fold) of decreased kcat values for GCV as a substrate relative to wild-type HSV-1 TK. Although N30-3 TK had detectable activities, linear velocity conditions could not be established. Under the conditions described in Materials and methods, N30-3 TK was evaluated in comparison with wild-type HSV-1 TK to determine a pseudo-first-order rate for thymidine and GCV. Although their respective pseudo-first-order rates were not very different, the product yields were drastically changed, 65% for thymidine and 70% for GCV with wild-type HSV-1 TK compared with 3% and 16%, respectively, for N30-3 TK (data not shown). This apparent product inhibition for N30-3 TK is similar to that reported for the Glu125 TK mutant (Drake et al., 1999
; Hinds et al., 2000
).
|
|
The cDNAs of a subset of HSV-1 TK mutants (Q30 TK, Q75 TK, Q7530TK, N7530TK and N30-3 TK) were subcloned into a retroviral plasmid, pLENTK, and used to transfect stably a human colon cancer cell line, HCT116. Additionally, wild-type HSV-1 TK and Q7530 TK retroviral plasmids were transfected into SW620, another human colon tumor line. From a panel of multiple HSV-1 TK expressing clones, a subset of HCT-116 clones expressing either wild-type HSV-1 TK, Q30 TK, Q75 TK, Q7530 TK N30-3 TK or N7530 TK were selected for equivalent levels of HSV-1 TK expression based on western blot determinations (data not shown). Despite numerous attempts, a clone expressing equivalent levels of N30-3 TK was not identified. Therefore, the N30-3 TK clone utilized in the following studies had 20% of the expression level of protein compared with the rest of the HSV-TK expressing cell panel.
The HCT116 cell lines expressing the different HSV-1 TK isoforms were evaluated for dose-dependent GCV sensitivity by using an MTT cell viability assay and a more sensitive clonal dilution assay. As shown in Figure 2, the Q30 TK, Q75 TK and Q7530 TK expressing cells were more sensitive to the lowest GCV concentration tested than wild-type HSV-1 TK expressing cells. At the 0.110 µM GCV doses, these same mutant TK expressing cell lines were similar in effectiveness to wild-type HSV-1 TK cells. The cell lines expressing N7530 TK and N30-3 TK, which were both poor GCV kinases in the activity screens, displayed only minor sensitivity to GCV at higher concentrations (110 µM). These cell lines were further evaluated using a more sensitive clonal dilution assay as described in Materials and methods. Each cell line was treated with low GCV doses (0.00011 µM) for 24 h, followed by drug removal and sequential dilution. As shown in Figure 3
, the number of surviving colonies plotted versus GCV doses indicated IC50 values of 0.01 and 0.008 µM, respectively, for wild-type HSV-1 TK and Q7530 TK. Q30TK expressing cell lines had GCV sensitivities (0.009 µM) analogous to the Q7530TK and wild-type HSV-TK cells (data not shown). Even the poor GCV metabolizing N7530 TK and N30-3 TK HCT-116 cell lines had therapeutically low IC50 values of 0.08 and 0.5 µM, respectively.
|
|
These same HCT-116 cell lines were evaluated for intracellular metabolism of [3H]thymidine and [3H]GCV. Following an 18 h incubation with 4 µM [3H]GCV, total methanol-soluble metabolites were extracted from each cell line. As shown in Table III, the Q7530 TK and Q75TK cell lines metabolized GCV at levels approaching those extracted from wild-type HSV-1 TK cells, whereas the Q30TK cells generated half of these levels. The N7530 TK and N30-3 TK cells generated GCV metabolite levels that were 15- and 200-fold lower than wild-type HSV-1 TK cells, consistent with the MTT cell killing data presented in Figure 2
. For [3H]thymidine, the same cell lines were grown to confluency and labeled with 2 µM [3H]thymidine for only 1 h to minimize cellular thymidine metabolism associated with DNA replication. Under these conditions, wild-type HSV-1 TK and Q75 TK cells generated twice the levels of thymidine metabolites compared with parent HCT-116 cells. In contrast, the Q30 TK, Q7530 TK, N30-3 TK and N7530 TK cell lines generated thymidine metabolites at levels comparable to the parent HCT-116 cells. An additional metabolic labeling comparison was done in another colon tumor cell line, SW620, expressing wild-type HSV-1 TK and Q7530 TK. As presented in the second column of Table III
, analogous results to the HCT-116 cell lines were obtained.
|
The FlexiDock component of the molecular modeling program SYBL 6.7 was used with the coordinates of two wild-type HSV-1 TK crystal structures (Champness et al., 1998) to dock thymidine or GCV in the active site of Q7530 TK and the corresponding half-site mutants, Q75-3 TK and Q30-4 TK. The FlexiDock program uses a genetic algorithm and energy scoring functions to evaluate thousands of potential ligandreceptor interactions, optimized for HSV-1 TK according to the following parameters: initial structure population, generation number, flexible bonds and hydrogen bonding interactions (see Materials and methods). When the native HSV-1 TK structures were modeled with thymidine (Figure 4A
) or GCV (Figure 4B
), there was a good correlation between the crystal and FlexiDock model. Root mean square (r.m.s.) values of 0.5594 for thymidine and 0.5750 for GCV were obtained, indicative of convergence in the thousands of structures evaluated by the genetic algorithm. As seen in Figure 5A
, modeling of the full set of mutations in Q7530 TK with thymidine predicted that it would bind backwards in the active site relative to wild-type HSV-1 TK. Modeling of GCV with Q7530 TK generated conformations analogous to that of GCV bound to wild-type HSV-1 TK (Figure 5B
). Modeling of the half-site Q75-3 TK bound with thymidine was similar to that of wild-type HSV-1 TK, which was expected since the Q75-3 TK possessed robust thymidine phosphorylation activities (Figure 6A
). In contrast, for modeling of Q30-4 TK, an enzyme that lacked catalytic activity, the most stable predicted conformation positioned thymidine with its 5-methyl group rotated 180° with respect to its position in the wild-type enzyme (Figure 6B
). Binding of GCV was not predicted to be significantly altered in the active site of Q75-3 TK or Q30-4 TK (data not shown). The cumulative modeling data indicate that the mutations near site 4 are most responsible for the observed loss of thymidine phosphorylation activities.
|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The Q7530 TK enzyme possessed the most interesting properties with regard to the development of a GCV-specific kinase. In comparison with wild-type TK, it had a decreased Km (appt) for GCV to 40 µM and a greater than 50-fold increase in the Km (appt) for thymidine to 53 µM. To our knowledge, the Q7530 TK is the first HSV-TK derived mutant to have a Km for GCV that is lower than that of thymidine. Q7530 TK also lacked TMPK activity and only trace amounts of [32P]-5-azido-dUMP were photoincorporated (Figure 1), results consistent with the inability of TMP to bind properly in the active site of this enzyme. In support of the kinetic data, Q7530 TK expressed in the HCT116 cell line generated levels of [3H]thymidine metabolites comparable to parental cells lacking HSV-1 TK, yet it was able to generate GCV metabolites comparable to those of expressed wild-type HSV-1 TK (Table III
). In cell viability assays, Q7530 TK expressing cells were equally sensitive to low GCV concentrations as cells expressing wild-type HSV-1 TK or Q30 TK. A more rapid initial rate of GCV metabolism, incorporation into DNA and accelerated S-phase cell cycle arrest have been observed in the Q7530 TK cells relative to those expressing wild-type enzyme (data not shown). It remains to be determined with in vivo studies whether the properties of Q7530 TK confer any therapeutic advantage in ablating tumor growth, as has been reported previously for Q30 TK relative to wild-type HSV-1 TK (Kokoris et al., 1999
). Another mutant HSV-1 TK possessing improved GCV metabolism, termed sr39 TK, has demonstrated improved efficacy of tumor regression in animal models and more sensitive killing of sr39TK transduced T-lymphocytes (Pantuck et al., 2002
; Qasim et al., 2002
). Hence it is feasible that Q7530 TK will function similarly in in vivo studies and experiments to evaluate these mutants in a murine prostate tumor model are in progress.
The other mutant enzyme with GCV kinase properties, N7530 TK, was not as efficient as Q7530 TK in utilizing GCV as a substrate. Despite a Km (appt) of 71 µM for GCV similar to wild-type HSV-1 TK, N7530 TK possessed an 800-fold decrease in the catalytic efficiency for GCV phosphorylation. However, the change in the Km (appt) of thymidine to 91 µM for N7530 TK was the most drastic change for any of the mutants analyzed. Like Q7530 TK, N7530 TK also lacked the ability to phosphorylate TMP (Table I, Figure 1
). When expressed in the HCT116 cells, extracted GCV metabolites relative to wild-type HSV-1 TK for N7530 TK were 40-fold lower. In cell viability assays, the N7530 TK expressing cells were insensitive to lower GCV doses, but displayed limited cytoxicity at higher doses (10 µM) (Figure 2
). We hypothesize that these GCV cytotoxicities would not have been observed had the deoxyprymidine activities of this enzyme not been decreased. This is because normal intracellular concentrations of thymidine exceed the Km for HSV-1 TK, while intracellular GCV concentrations are considerably lower than its Km. Thus, thymidine would normally out-compete GCV for binding to the active site (Markham and Faulds, 1994
; Kokoris et al., 1999
).
In order to evaluate why the different structural mutations comprised in Q7530 TK resulted in a more GCV-specific kinase, a FlexiDock molecular modeling algorithm was generated and tested. This FlexiDock program allows both the ligand and binding pocket to flex during docking, so that an induced fit can be more thoroughly evaluated among the thousands of different possible structures analyzed. As shown in Figures 46, the structural differences caused by mutation at Ala168 were the primary contributors to the different thymidine activities observed between wild-type HSV-1 TK, Q7530 TK, Q75-3 TK and Q30-4 TK. For Q7530 TK, the loss in thymidine kinase activity was directly related to the Tyr168 mutation near site 4 (Figure 5A
). In both Q7530 TK and Q30-4 TK, the bulky phenyl group of Tyr168 sterically blocks the binding pocket where the 5-methyl group of thymidine normally is bound. In a previous molecular modeling analysis of Q30 TK based on a temperature annealing algorithm (Kokoris et al., 1999
), the same conclusion of a steric interference role for Tyr168 was drawn. Also, an early structural study examining HSV-1 TK mutants isolated from BVdU-resistant viruses revealed that a single amino acid substitution at A168T conferred BVdU resistance (Larder et al., 1983
). The likelihood that thymidine actually is bound with its 5-methyl group flipped 180° relative to wild-type enzyme, seen in Q30-4 TK (Figure 6B
), or rotated 180°, seen in Q7530 TK (Figure 5A
), is unknown since the FlexiDock modeling program assumes that the substrate had entered the binding site. The most probable scenario is that thymidine has limited access to the active site of these particular enzymes, as supported by the activity, photoaffinity and cellular data. Although not shown or described, the molecular modeling algorithm could also accurately predict poor or altered thymidine binding in N7530 TK, N30-3 TK and Q30TK.
For modeling with GCV, more subtle variations in binding were detected. Although thymidine and GCV share the same binding site, GCV is positioned higher, closer to Arg176. As a result, the steric interactions between GCV and the phenyl group of Tyr168 were not significant enough to alter GCV binding in Q7530 TK (Figure 5B). For N7530 TK, the extra substitution of Asn125 may be a destabilizing factor to account for its decreased activity relative to wild-type or Q7530 TK. One inconsistency with the modeling and activity data is that of the Q30-4 TK enzyme. Even though the purified Q30-4 TK enzyme lacked apparent GCV phosphorylation activity, it was positioned in the modeled active site similar to that of the highly active Q7530 TK. One reason for this inconsistency may be that as determined in the kinetic evaluations for these mutants (Table II
), the overall changes in the Kms for GCV (less than 2-fold) were not as high as those for thymidine (1590-fold). Hence there was a greater probability that the modeling program would more readily predict greater changes for thymidine binding. Also, the modeling program does not allow the visualization of how the substrate enters the active site, as the FlexiDock program makes an assumption that the substrate can gain access to the active site. It is possible that other structural perturbations within Q30-4 TK prevent both thymidine and GCV binding, and hence activity, in the purified enzyme. Lastly, the current HSV-1 TK model is weighted for the influence of hydrogen bonding on substrate binding, interactions which have consistently proven to be important for thymidine binding (Brown et al., 1995
; Champness et al., 1998
; Hinds et al., 2000
) and is reflected in the models of the mutant enzymes. As additional HSV-1 TK mutants are generated, there are many parameters within the FlexiDock program that can be refined to improve analysis of GCV binding, especially an emphasis on hydrophobic and steric interactions.
In conclusion, from the 16 different HSV-1 TK mutants generated, purified and assayed, two mutant enzymes, Q7530 TK and N7530 TK, have been characterized as GCV-specific kinases. Originally we proposed to develop these novel GCV kinases by incorporating mutations known to minimize deoxypyrimidine substrate utilization (Hinds et al., 2000) with mutations known to improve the catalytic efficiency of GCV phosphorylation (Black et al., 1996
). The Q7530 TK essentially achieved our initial goal by possessing a Km for GCV lower than that of thymidine. Introduction of an Asn or Glu for Gln125 in the Q75 TK or Q30 TK mutant backgrounds eliminated thymidine metabolism of HSV-1 TK and, with the exception of N7530 TK, also proved to be equally detrimental for GCV phosphorylation. However, these studies do demonstrate that even though all nucleoside/nucleotide substrates bind in the same active site region, discrimination of phosphorylation of guanine analogs such as GCV can be achieved while eliminating activity for thymidine/TMP. Besides the well documented importance of Gln125 (Brown et al., 1995
, Champness et al., 1998
; Hinds et al., 2000
), Ala168 was again identified (Kokoris et al., 1999
) as being a critical participant in substrate binding in the active site. It is clear that the half-site mutations at site 3 in the 753 TK are responsible for the improved thymidine and GCV activities, while the site 4 half-site mutations (that included mutations to Ala168) were generally detrimental to substrate binding. By analyzing these new HSV-1 TK mutant enzymes with the FlexiDock program, we have refined a computer modeling algorithm that can be used to design and test substrate binding properties before generating the mutants. Our future plan is to use this modeling approach as a screen to identify new HSV-1 TK mutants, in an effort to develop enzymes that lack deoxypyrimidine binding and possess Kms for GCV that are lower than 10 µM.
![]() |
Notes |
---|
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Black,M.E., Newcomb,T.G., Wilson,H.M. and Loeb,L.A. (1996) Proc. Natl Acad. Sci. USA, 93, 35253529.
Boucher,P.D., Ruch,R.J. and Shewach,D.S. (1998) Hum. Gene Ther., 9, 801814.[ISI][Medline]
Brown,D.G., Visse,R., Sandhu,G., Davies,A., Rizkallah,P.J., Melitz,C., Summers,W.C. and Sanderson,M.R. (1995) Nature Struct. Biol., 2, 876881.[ISI][Medline]
Champness,J.N., Bennett,M.S., Wien,F., Visse,R., Summers,W.C., Herdewjin,P., de Clercq,E., Ostrowski,T., Jarvest,R.L. and Sanderson,M.R. (1998) Proteins, 32, 350361.[CrossRef][ISI][Medline]
Chen,M.S. and Prusoff,W.H. (1978) J. Biol.Chem., 253, 13251327.[Abstract]
Chen,M.S., Walker,J. and Prusoff,W.H. (1979) J. Biol. Chem., 254, 1074710753.[Medline]
Darby,G., Larder,B.A. and Inglis,M.M. (1986) J. Gen. Virol., 67, 753758.[Abstract]
Drake,R.R., Wilbert,T.N., Hinds,T.A. and Gilbert,K.M. (1999) J. Biol. Chem., 274, 3718637192.
Freeman,S.M., Whartenby,K.A., Freeman J.L., Abboud,C.N. and Marrogi,A.J. (1996) Semin. Oncol., 23, 3145.[ISI][Medline]
Halloran,P.J. and Fenton,R.G. (1998) Cancer Res., 58, 38553865.[Abstract]
Hinds,T.A., Compadre,C., Hurlburt,B.K. and Drake,R.R. (2000) Biochemistry, 39, 41054111.[CrossRef][ISI][Medline]
Klatzmann,D., Cherin P., Bensimon,G. and Boyer,O. et al. (1998) Hum. Gene Ther., 9, 25852594.[CrossRef][ISI][Medline]
Kokoris,M.S., Sabo,P., Adman,E.T. and Black,M.E. (1999) Gene Ther., 6, 14151426.[CrossRef][ISI][Medline]
Larder,B.A., Derse,D., Cheng,Y. and Darby G. (1983) J. Biol. Chem., 258, 20272033.
Liu,Q. and Summers,W.C. (1988) Virology, 163, 638642.[ISI][Medline]
Machida,H. (1986) Antimicrob. Agents Chemother., 29, 524526.[ISI][Medline]
Markham,A. and Faulds,D. (1994) Drugs, 48, 455484.[ISI][Medline]
McMasters,R., Jones,K.E., Saylors,R.L., Hendrix,E.M., Moyer,M.P. and Drake,R.R. (1998) Hum. Gene Ther., 8, 18251835.[ISI]
Michael,M., Fetzer,J. and Folkers,G. (1994) Eur. J. Biochem., 226, 219226.[Abstract]
Munir,K.M., French,D.C., Dube,D.K. and Loeb,L.A. (1992) J. Biol. Chem., 267, 65846589.
Pantuck,A.J,, Matherly,J., Zisman,A., Nguyen,D., Berger,F., Gambhir,S.S., Black,M.E., Belldegrun,A. and Wu,L. (2002) Hum. Gene Ther., 13, 777789.[CrossRef][ISI][Medline]
Qasim,W., Thrasher. A.J., Buddle,J., Kinnon,C., Black,M.E. and Gaspar,H.B. (2002) Gene Ther., 9, 824827.[CrossRef][ISI][Medline]
Rechtin,T.M., Black,M.E., Mao,F., Lewis,M. and Drake,R.R. (1995) J. Biol. Chem., 270, 70557060.
Rechtin,T.M., Black,M.E. and Drake,R.R. (1996) Anal. Biochem., 237, 135140.[CrossRef][ISI][Medline]
Rubsam,L.Z., Boucher,P.D., Murphy,P.J., KuKuruga,M. and Shewach,D.S. (1998) Cancer Res., 58, 38733882[Abstract]
Sterman,D.H., Treat,J., Litzky,L.A., Amin,K.M. et al. (1998) Hum. Gene Ther., 9, 10831092.[ISI][Medline]
Thust,R., Tomicic,M., Klocking,R., Wutzler,P. and Kaina,B. (2000a) Mutagenesis, 2, 177184.[CrossRef]
Thust,R., Tomicic,M., Klocking,R., Voutilainen,N., Wutzler,P. and Kaina,B. (2000b) Cancer Gene Ther., 7, 107117.[CrossRef][ISI][Medline]
Wei,S., Chao,Y., Lin,W., Yang,D., Shih,Y. Chang,L., Whang-Peng,J. and Yang,W.K. (1998) Exp. Cell. Res., 241, 6675.[CrossRef][ISI][Medline]
Received May 9, 2002;