Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, December 2005, p. 4834-4842, Vol. 49, No. 12
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.12.4834-4842.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Biological Sciences,1 Department of Chemistry, Boehringer Ingelheim (Canada) Ltd., Laval, Quebec, Canada H7S 2G52
Received 18 July 2005/ Returned for modification 22 September 2005/ Accepted 4 October 2005
|
|
|---|
|
|
|---|
All papillomaviruses have a small circular double-stranded DNA genome which encodes for only eight well-characterized proteins (for a recent review, see reference 21). The most highly conserved protein, and the only one with enzymatic activity, is the E1 helicase (reviewed in references 33 and 40). E1 is an attractive target for the development of anti-HPV drugs because it is essential for viral replication and pathogenesis (24, 34). Indeed, it has been shown in the cottontail rabbit papillomavirus (CRPV) infection model that frameshift mutations in the E1 open reading frame abrogate the ability of the CRPV genome to induce papillomas when inoculated into the skin of domestic rabbits (41).
E1 is the replicative helicase of papillomaviruses. It binds cooperatively to the origin of replication in conjunction with the E2 protein (12, 22, 27, 30). Formation of the E1-E2-origin complex involves not only the binding of both proteins to specific DNA elements in the origin but also a protein-protein interaction between the N-terminal transactivation domain of E2 and the helicase/ATPase domain of E1 (1, 3, 4, 36). We recently reported a class of small-molecule inhibitors of HPV DNA replication that bind to the transactivation domain of E2 and prevent its interaction with E1 (37, 39). Assembly of the E1-E2-ori complex facilitates the recruitment of additional E1 molecules to the origin, which assemble into hexamers in a reaction that is stimulated by ATP binding (11, 26, 35). These hexamers are the catalytically active form of E1 capable of melting the origin and unwinding the viral DNA ahead of the replication fork (28). As is the case for most helicases, the DNA-unwinding activity of E1 is powered by the hydrolysis of ATP. We previously characterized the enzymatic activities of highly purified recombinant HPV6 and -11 E1 proteins, produced with a baculovirus expression system (38). These studies revealed that HPV6 and -11 E1 proteins have similar Km values for ATP (12 and 6 µM, respectively) and that the ATPase and unwinding activities are contained within the C-terminal half of E1 (amino acids 353 to 649), the same region that binds to E2 (38). Interestingly, we found that E2 hinders the ATPase activity of E1 by raising its Km for ATP approximately sevenfold. Conversely, we observed that ATP impairs the cooperative binding of E1 and E2 to the origin, most likely by weakening the E1-E2 protein-protein interaction (38). These results and others led to the proposal that ATP not only powers the helicase activity of E1 but also acts as a molecular trigger during the initiation of viral DNA replication to help sever the interaction of E1 with E2 and promote its hexamerization at the origin (26, 38).
We recently reported on the discovery, by high-throughput screening of our compound collection, of a biphenylsulfonacetic acid inhibitor of HPV6 E1 (compound 1; Table 1) (10). This compound inhibits the ATPase activity of HPV6 E1 with a 50% inhibitory concentration (IC50) of 3 µM at an ATP concentration of 1 µM. More potent derivatives of this original hit were identified which inhibited activity with IC50 values as low as 4 nM. Here we show that these more potent inhibitors of HPV6 E1 are less active against the related HPV11 E1 protein and that a single amino acid difference between the two enzymes accounts for most of this differential inhibition. Furthermore, we describe the kinetic mechanism of action for these compounds and the activity of more potent compounds against the helicase activity of E1. Together, these results define a pocket in E1 that can be exploited for small-molecule modulation of E1 activity.
|
View this table: [in a new window] |
TABLE 1. Inhibitor potency against HPV6, 11, and 18 E1 ATPasesa
|
|
|
|---|
E1 ATPase assays in solution.
Most solution assays were performed using the previously reported scintillation proximity method (19). Briefly, 45-µl reactions were carried out in assay buffer (20 mM HEPES, pH 7.5, 1 mM dithiothreitol, 0.05 mM EDTA, 0.005% IGEPAL-CA630 [vol/vol]; Sigma [equivalent to Nonidet P-40]) supplemented with 2% dimethyl sulfoxide. Reaction mixtures contained 1 µM ATP, unless stated otherwise, and a concentration of Mg(OAc)2 500 µM greater than that of ATP. [
-33P]ATP was included either at 30 nCi/well in reactions performed at 1 µM of total ATP or at 100 µCi per µmol of unlabeled ATP in initial velocity studies. Detection was performed by adding to reaction mixtures 40 µl of a suspension of scintillation proximity beads (10 mg/ml; Amersham Biosciences) in ammonium molybdate (0.67%) and 0.8 M HCl, followed by a solution of 100 mM citric acid and 7 M cesium chloride (80 µl). Scintillation signals were quantified using a TopCount NXT microplate scintillation counter (Packard Instruments). To verify inhibitor specificity, reactions were carried out with either purified herpes simplex virus helicase primase (2 nM) (9) or SV40 T antigen (30 nM) (32) under the same conditions used for E1. Approximately 20% conversion of the ATP substrate was obtained at these concentrations of enzymes and under these reaction conditions.
Kinetic parameters were determined by initial velocity experiments. Reactions were performed in total volumes of 160 µl rather than 45 µl. Twenty-microliter samples were removed at multiple time points and quenched with 40 µl of ammonium molybdate-scintillation proximity bead suspension followed by 80 µl of citric acid-cesium chloride. The concentration of product formed at each time point was estimated by determining, for each experiment, a conversion factor relating nM of phosphate product (determined from the ratio of ATP to ADP measured by thin-layer chromatography) to cpm detected by scintillation proximity. Because of the multiple steps involved in this procedure and the relatively high variability in quantification by thin-layer chromatography, estimates of the absolute concentrations of phosphate were found to vary by a factor of 2, although the relative concentrations obtained for different time points or reaction wells were determined much more accurately. This procedure has been described in detail previously (19). Kinetic data were analyzed using the equations for linear or hyperbolic competitive inhibition (29) with the program GraFit 3.0 (Erithacus Software Ltd.).
E1 helicase assays. Helicase assays were carried out using the same buffer used for ATPase assays and the method described previously (38). Briefly, reaction mixtures contained HPV6 E1 (10 nM), a partial duplex substrate consisting of a 33P-labeled 24-base oligonucleotide covalently linked to biotin and annealed to M13 single-stranded DNA, ATP at 300 µM unless stated otherwise, and Mg(OAc)2 at a concentration 500 µM greater than the ATP concentration. Reaction mixtures were incubated at 37°C for 30 min, and products were detected by addition of streptavidin-coated scintillation proximity beads.
In the experiment used to measure inhibition by ATP-
-S, helicase assays were carried out essentially as described above, except that ATP and Mg(OAc)2 concentrations were 300 µM and 1 mM, respectively. Helicase and ATPase activities were detected by removing two 20-µl aliquots from 60-µl reaction mixtures. One aliquot was used to quantify helicase activity by measuring the amount of unwound oligonucleotide product following separation from the substrate by electrophoresis on 20% acrylamide gels, as described previously (38). The other aliquot was used to quantify ATP hydrolysis by a colorimetric method using ammonium molybdate and malachite green (18).
ATPase assays using immobilized E1. We previously showed that the helicase domain of HPV11 E1 (amino acids 353 to 649) made by in vitro transcription/translation is oligomeric (35) and has ATPase activity comparable to that of the full-length protein (38). HPV6 and -11 E1(353-649) were immobilized on plates through an affinity-purified rabbit antiserum directed against the C-terminal 14 amino acids of HPV11 E1. The purified antibody (1 µg/well) was immobilized in white Maxi-Sorp 96-well plates (Nunc) by incubation for 1 h at room temperature (or overnight at 4°C) in 100 µl antibody coating buffer (100 mM potassium phosphate, pH 6.5, 0.02% sodium azide, 0.001% glutaraldehyde). After removal of the coating solution, the plates were blocked for 1 h at room temperature (or 4°C overnight) with 200 µl/well of antibody blocking buffer (50 mM phosphate-buffered saline, 1% sucrose, 0.05% casein, 0.05% sodium azide). Blocking buffer was removed by aspiration, and coated plates were stored at 4°C until needed, up to 3 months. For reactions with FLAG-tagged E1 proteins, the monoclonal antibody M2 (Kodak) was immobilized by the same procedure. On the day of the experiment, plates were washed once with 200 µl of E1 coating buffer (20 mM Tris, pH 7.6, 500 mM NaCl, 0.02% NP-40 [vol/vol], 5 mg/ml bovine serum albumin, 1 mM dithiothreitol). In vitro-translated helicase domains (2 µl/well of the translation reaction for the HPV6 and -11 proteins or 8 µl/well for the others) were immobilized by incubation for 1 h at room temperature in E1 coating buffer (100-µl total volume). Negative control wells contained 2 µl or 8 µl per well of an in vitro translation reaction mixture programmed with pTM1 plasmid. Wells were then washed three to five times in wash buffer (same as ATPase buffer described above but supplemented with 100 mM NaCl and 0.1% [vol/vol] Triton X-100). Thorough washing was needed to remove proteins from the in vitro translation mixture, other than E1, that have ATPase activity. Efficient removal of these contaminating ATPases was verified by comparing the signals in negative control wells to those in wells in which no in vitro translation mixture was added. To determine the activity of the immobilized E1, 45-µl ATPase reactions were carried out and products were detected exactly as in the solution-reaction protocol. The activities of the HPV6 and HPV11 proteins, and of the chimeras and mutants derivatives thereof, varied by a factor of 3, but we did not determine whether this variation was caused by differences in catalytic efficiency of the enzymes or differences in the level of protein obtained by in vitro translation. Activity for the HPV16, HPV18, and CRPV enzymes was lower. In control experiments, we determined that immobilization of purified HPV6 and HPV11 E1 did not affect their activities or their levels of inhibition by the ATPase inhibitors (data not shown).
|
|
|---|
Preliminary experiments indicated that inhibitor 1 and other compounds in this series were less active at higher ATP concentrations, suggesting that they might be competitive inhibitors. Data from initial velocity kinetic experiments using a range of ATP and inhibitor concentrations were collected and analyzed by nonlinear regression as well as by visual examination of Lineweaver-Burke double-reciprocal plots. Systematic deviations of the residuals were observed after fitting to a linear competitive mechanism (Fig. 1A, dotted lines). One hallmark of the standard linear competitive mechanism is that a replot of the slopes for each inhibitor concentration from the Lineweaver-Burke plot should yield a straight line versus inhibitor concentration. This was not observed for compound 7 and other inhibitors in this series. Rather, a replot of their slopes yielded a hyperbolic curve (Fig. 1B), which is the defining characteristic of hyperbolic inhibition. As expected, a much better fit was observed for the hyperbolic competitive mechanism (Fig. 1A, solid lines), which allows formation of both the binary enzyme-inhibitor and ternary enzyme-inhibitor-substrate complexes but assumes that both complexes turn over substrate at the same rate (Fig. 1C) (29). Fitting to a more complex mechanism, in which the substrate turnover of these two complexes is different (hyperbolic noncompetitive inhibition (29), did not improve the fit. Values for the parameter
, the degree to which saturation with inhibitor raises the apparent value of Km for substrate or conversely the degree to which saturation by substrate raises the apparent Ki of the inhibitor, ranged from 10 to greater than 100 for different inhibitors in this series (Table 2). One consequence of this mechanism is that at high concentrations of substrate (i.e., near or above
Km), the enzyme is active even at saturating inhibitor concentrations. In practice, inhibitors giving high values of
may appear to act by the simpler linear competitive mechanism. An important conclusion from these results is that compounds in this series do not bind directly at the ATP-binding site, but rather affect substrate binding through an allosteric mechanism. Other examples of this kinetic mechanism have been reported for both synthetic and natural allosteric inhibitors (see Discussion).
![]() View larger version (14K): [in a new window] |
FIG. 1. Kinetic mechanism of compound 7. (A) Lineweaver-Burke double reciprocal plots for initial velocity data using 10 nM HPV6 E1 and 12.5 to 100 µM ATP. Data points correspond to observed initial velocities for inhibitor concentrations of 480 nM (black triangles), 240 nM (red squares), 120 nM (green circles), 60 nM (blue squares), 30 nM (brown asterisks), and 15 nM (black circles), as well as those observed (obs) without inhibitor (magenta diamonds). Corresponding lines were determined by nonlinear regression using the full data set to the equations for linear competitive or hyperbolic competitive inhibition (dotted or solid lines, respectively). The kinetic parameters obtained for compound 7 from the hyperbolic competitive fit are given in Table 2 (first of the two experiments). (B) Slope replot. For each inhibitor concentration in panel A, the slope of the 1/v versus 1/[S] plot was determined by linear regression, and these slopes are plotted versus inhibitor concentration (diamonds). Lines are plotted for the predicted slopes after fitting to the linear or hyperbolic competitive mechanisms (dotted or solid lines in both panels A and B). The solid line corresponds closely to the experimentally determined slopes. (C) Scheme for the hyperbolic competitive mechanism (29). The parameter indicates the degree to which inhibitor binding affects substrate binding, or vice versa. Values of greater than 1 indicate that binding of one ligand weakens binding of the other. As approaches infinity, the lower line becomes insignificant and the mechanism becomes indistinguishable from linear competitive inhibition. For a relatively slow enzyme such as E1, it is expected but not rigorously proven that the value of Km is approximately equal to the substrate dissociation constant Ks and the value of kcat is approximately equal to the rate constant for conversion of bound substrate to product, kp.
|
|
View this table: [in a new window] |
TABLE 2. Kinetic parameters for inhibition of HPV6 E1 ATPase activitya
|
-S, was 3- to 15-fold lower than that predicted solely on the basis of the higher ATP concentration. For example, we obtained an IC50 for compound 8 of 1.4 µM in a helicase assay performed at 100 µM ATP (Table 3) but a lower value of 100 nM for inhibition of ATPase activity at the same ATP concentration. Hence, other factors related to differences in assay kinetics and/or conditions (e.g., temperature or the presence of nucleic acids) probably also influence inhibitor potency. Further experiments will be required to pinpoint these additional contributing factors. Nevertheless, it remains that the main conclusion from the above results is that the biphenylsulfonacetic acid class of inhibitors is capable of abrogating the helicase activity of E1. We also tested whether these inhibitors could affect other activities of E1 that are modulated by ATP, specifically its ability to oligomerize into hexamers and its cooperative binding with E2 to the origin (see introduction). In neither case, however, did we observe any significant inhibition (data not shown), perhaps because high ATP concentrations are needed in both of these assays. |
View this table: [in a new window] |
TABLE 3. Inhibitor potency against HPV6 E1 helicase activitya
|
![]() View larger version (12K): [in a new window] |
FIG. 2. Inhibition of HPV6 E1 helicase activity. Effect of compounds 2 (closed circles) and 8 (open circles) on HPV6 E1 helicase activity. An ATP concentration of 100 µM and serial 2.5-fold dilutions of the inhibitors (160 to 0.7 µM) were used. Data points are averages of duplicate values and are connected by lines to aid in visualization. The IC50 values obtained from these curves are given in Table 3.
|
To facilitate these studies, we adapted our standard scintillation proximity method to measure the ATPase activity of the isolated helicase domain of E1 (amino acids 353 to 649), made by in vitro translation in a reticulocyte lysate. We found that the ATPase activity of the helicase domains of HPV6 and -11 E1 proteins could be measured following immobilization of the proteins in 96-well plates with an antibody specific for the C terminus of the protein and removal, by several washes, of the contaminating ATPase activities present in the reticulocyte lysate.
This methodology was used to quickly narrow down a region in the helicase domain of HPV11 E1 responsible for its lower sensitivity to the biphenylsulfonacetic acid inhibitors. Specifically, genes encoding the four HPV6-HPV11 E1 chimeric proteins shown in Fig. 3A were constructed, and the corresponding in vitro-translated proteins were tested for inhibition by compound 3 (Fig. 3B). Whereas relatively weak inhibition was observed for the HPV11 protein and for chimera A, this compound inhibited the activity of chimeras B, C, and D to a similar extent as that of the HPV6 helicase domain. Similar results were obtained with other inhibitors of the same class (data not shown). These results indicated that the region of E1 between amino acids 438 and 493 plays a substantial role in inhibitor sensitivity. Within this region, only three residues, at positions 461, 486, and 493, are not conserved between HPV6 and HPV11 E1. To test if one or more of these residues were responsible for the observed potency differences, they were individually mutated in the context of the isolated helicase domain of HPV11 E1 to the corresponding residue in HPV6 E1. The substitution of Cys-486 for Tyr substantially increased the sensitivity of HPV11 E1 to inhibition by the more potent inhibitor, 9 (Fig. 3C), as well as by other biphenylsulfonacetic acids with various substituents at C-3' or C-4' (data not shown). The IC50 value determined for HPV11 E1 C486Y was 24 nM, only fourfold higher than that measured for the HPV6 protein (5.6 nM). Note also that the IC50 of compound 9 against the helicase domain of HPV6 E1 made in vitro is very similar to the value of 4.3 nM obtained with purified full-length HPV6 E1 (Table 1), further supporting the notion that the N-terminal domain of E1 plays no significant role in ATP hydrolysis. From these results, we conclude that a tyrosine at position 486 of the E1 helicase domain is important for inhibitor binding.
![]() View larger version (11K): [in a new window] |
FIG. 3. Amino acid 486 of HPV6 and HPV11 E1 is a key determinant of inhibitor potency. (A, top) Schematic representation of the E1 helicase domain, with the four regions (labeled A to D) (5) of similarity to SV40 T antigen indicated by black boxes. Arrows indicate the positions of the highly conserved Walker A and B motifs, involved in ATP binding and hydrolysis. Positions of residues differing in the HPV6 and -11 sequences are indicated by black bars. (Bottom) Schematic representations of chimeras A to D, with portions derived from HPV11 in black and those derived from HPV6 in white. Chimeras A to D have residues corresponding to HPV6 up to amino acids 437, 493, 531, and 582, respectively. Amino acids 461, 486, and 493, the three amino acid differences between chimeras A and B, are indicated by asterisks. (B) Inhibition observed in ATPase assays using 2 µM compound 3 for HPV11 and -6 E1 as well as the four chimeras. Values are averages of duplicate data points with the standard deviations indicated. Reactions without inhibitors gave signals from 2,000 to 6,000 cpm depending on the protein. Signals from blank reactions (360 cpm) were subtracted to determine percent inhibition. (C) Inhibition observed in ATPase assays using 10 nM of inhibitor 9 for HPV11 and -6 E1 as well as HPV11 mutants L461F, C486L, C486Y, and K493S. (HPV16 E1 encodes a leucine at position 486 but the same residue as HPV11 at the other two positions.) Values are averages of duplicate or quadruplicate data points with the standard deviations indicated. In this experiment, reactions without inhibitor had 3,000 to 9,000 cpm, whereas blanks, subtracted before calculating percent inhibition, had 300 cpm. IC50 values obtained for two of the proteins in a separate experiment are also given.
|
ATPase activity observed with the high-risk E1 or CRPV protein was lower than that observed for HPV11 and -6, and significant activity could only be observed using larger amounts of protein (data not shown). The enzymatic properties of E1 from these other HPV types have not been well characterized in the literature, and further experiments will be required to determine whether this lower activity was due to lower protein expression levels, a lower degree of oligomerization, or differences in kinetic parameters. Regardless, compound 9 at 80 µM clearly inhibited the HPV18 enzyme, but had only a modest effect on CRPV E1 and no discernible activity against the HPV16 enzyme (Fig. 4). Several inhibitors had similar or even greater activity against HPV18 E1 than we had observed against the HPV11 protein (Table 1). Interestingly, for both HPV18 and CRPV, the residue corresponding to Tyr-486 is also a tyrosine, whereas it is a leucine in HPV16. These observations are consistent with the need for a tyrosine at position 486 for optimal compound binding (Fig. 3C).
![]() View larger version (13K): [in a new window] |
FIG. 4. Determination of ATPase activity and compound inhibition for HPV16, HPV18, and CRPV E1 proteins. Scintillation proximity ATPase assays were performed using in vitro-translated E1 helicase domains immobilized via the anti-FLAG antibody. Raw cpm are shown for negative controls (in vitro translation reactions programmed with empty pTM1 vector alone; white bars), and E1 reactions in the absence of inhibitor (black) or in the presence of 80 µM compound 5 (dark gray) or 9 (light gray). Error bars correspond to standard deviations for quadruplicate values. Values for percent inhibition were calculated by comparing blank-subtracted reactions in the presence and absence of inhibitor and are given above the corresponding gray bars.
|
|
|
|---|
The crystal structure of the HPV18 E1 helicase domain has recently been reported (1), and it is interesting to use this structure to speculate on the location of the biphenylsulfonacetic acid binding site. The residue identified in our studies, Tyr-486 in HPV6 E1 (corresponding to Tyr-492 in HPV18; Fig. 5A), is close to the ATP binding site, on the same
-helix as the highly conserved Lys-484 (Lys-490 in HPV18), which interacts with ATP and is essential for catalysis (15). This crystallized domain of E1 is monomeric, but based on the structure of hexameric SV40 T antigen (14, 20), which is homologous to E1, Tyr-486 should be near to but not at the monomer-monomer interface. Although we have not provided evidence that the biphenylsulfonacetic acids interact directly with Tyr-486, this is a likely possibility. The tyrosine side chain extends from the
-helix in the opposite direction as does the lysine (Fig. 5B). Inhibitors interacting with this tyrosine could cause a conformational change that would affect binding of ATP, consistent with the allosteric mechanism observed for this series. Using the Site Finder application of the MOE software package (Chemical Computing Group Inc., Montreal, Canada), we found no potential small-molecule binding sites on the surface of HPV18 E1 except for the ATP binding site, so it is conceivable that these compounds bind through an induced-fit mechanism.
![]() View larger version (47K): [in a new window] |
FIG. 5. Structure of HPV18 showing the possible binding region of biphenylsulfonacetic acid inhibitors. (A) Sequence alignment of part of the helicase domains of the E1 proteins from two low-risk (HPV6 and -11) and two high-risk (HPV16, and -18) HPV types as well as those of bovine papillomavirus 1 (BPV1) and CRPV. Fully conserved residues are red, highly conserved residues are blue, and others are black. The positions of Tyr-486 in HPV6a E1 and of the corresponding amino acids in other E1 proteins are highlighted. For the entire helicase domain, the two low-risk sequences are 93% identical, whereas the two high-risk sequences are 65% identical. For other possible pair-wise combinations, the sequence identity ranges from 50% to 60%. (B) Ribbon diagram for the HPV18 E1 helicase domain (1). The Walker A and B motifs involved in ATP binding are blue and orange, respectively. The side chains of the catalytic Lys-490 and Tyr-492 are shown extending in opposite directions from the -helix (equivalent to Lys-484 and Tyr-486, respectively, in HPV6 E1). This illustration was created using Viewer Lite version 5.0 (Accelrys).
|
The work presented here has highlighted the presence of a small-molecule binding pocket on at least some HPV E1 proteins. Guided by structural information, it may now be possible to develop compounds active against not only HPV6 and -11 E1 but also HPV18 and possibly other high-risk proteins. Although the novel ATPase inhibitors presented in this work currently lack the necessary physicochemical properties to be good drug candidates (10), they clearly show the potential for compounds targeting this site to inhibit E1 ATPase and helicase activities and thus potentially to serve as antivirals for HPV.
Present address: Institut de Recherches Cliniques de Montreal, 110 Pine Ave. West, Montreal, Quebec, Canada H2W 1R7. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»