|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 8, 6976-6985, February 20, 2004
Crystal Structure of the E2 Transactivation Domain of Human Papillomavirus Type 11 Bound to a Protein Interaction Inhibitor*![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ||
From the
Received for publication, October 16, 2003 , and in revised form, November 21, 2003.
Interaction between the E2 protein and E1 helicase of human papillomaviruses (HPVs) is essential for the initiation of viral DNA replication. We recently described a series of small molecules that bind to the N-terminal transactivation domain (TAD) of HPV type 11 E2 and inhibits its interaction with E1 in vitro and in cellular assays. Here we report the crystal structures of both the HPV11 TAD and of a complex between this domain and an inhibitor, at 2.5- and 2.4-Å resolution, respectively. The HPV11 TAD structure is very similar to that of the analogous domain of HPV16. Inhibitor binding caused no significant alteration of the protein backbone, but movements of several amino acid side chains at the binding site, in particular those of Tyr-19, His-32, Leu-94, and Glu-100, resulted in the formation of a deep hydrophobic pocket that accommodates the indandione moiety of the inhibitor. Mutational analysis provides functional evidence for specific interactions between Tyr-19 and E1 and between His-32 and the inhibitor. A second inhibitor molecule is also present at the binding pocket. Although evidence is presented that this second molecule makes only weak interactions with the protein and is likely an artifact of crystallization, its presence defines additional regions of the binding pocket that could be exploited to design more potent inhibitors.
Human papillomaviruses (HPVs)1 are the etiological agents of malignant and benign lesions of the differentiating squamous or mucosal epithelium, notably of cervical cancer. Approximately 25 HPV types replicate in mucosal tissues of the anogenital tract. HPV16, -18, and -31 are the most prevalent "high-risk" types found in pre-cancerous or malignant lesions of the cervix. HPV6 and -11 are the most common "low-risk" types, which cause benign genital warts (condyloma acuminata), a less serious condition but one of the most common sexually transmitted diseases (1). Currently, no specific antivirals are available for the treatment of HPV infections.
The small circular double-stranded DNA genome of papillomavirus is actively maintained as a multicopy episome in the nucleus of infected epithelial cells. This process is dependent on replication of the viral genome by the viral E1 and E2 proteins, in conjunction with the host DNA replication machinery. E2 is a sequence-specific DNA-binding protein that has a number of functions in the viral lifecycle. In addition to its role in the initiation of viral DNA replication, E2 is involved in regulating the transcription of viral genes (2-7), and in the segregation of the viral genome during cell division (8, 9). As a replication initiation factor, E2 binds with high affinity to specific sites located within the viral origin (ori) to help recruit it to the E1 helicase (10-13). Formation of a ternary complex between E1, E2, and the origin relies not only on the interaction of E1 and E2 with specific DNA sequences at the origin but is also critically dependent on a direct interaction between these two proteins (14-18).
The 40-kDa E2 protein exists in solution as a dimer and consists of an N-terminal transactivation domain (TAD) of
Crystal structures have been reported for the TAD of HPV16 (amino acids 1-201) (20) and for a truncated TAD of HPV18 (amino acids 66-215) lacking two N-terminal
We recently described a series of small molecules that antagonize the E1-E2 protein-protein interaction by binding reversibly to the E2 TAD (22, 23). This class of inhibitors, termed "indandiones" because they feature an indandione system spirofused onto a substituted tetrahydrofuran ring, blocks assembly of the HPV11 E1-E2-ori ternary complex at submicromolar concentrations in vitro and abrogates viral DNA replication in a transient transfection assay. To our knowledge, these are the only small molecule inhibitors of HPV DNA replication with cellular activity. Here, we report the three-dimensional structure of the HPV11 TAD and that of a TAD-indandione inhibitor complex. The HPV11 TAD structure is very similar to that of the analogous domains of HPV16 and -18. Inhibitor binding does not significantly affect the protein backbone but results in the formation of a deep pocket that makes several well defined interactions with the inhibitor. Residues that have previously been implicated in binding E1 or the indandione inhibitors, including Gln-12 (24), Glu-39 (25-27), and Glu-100 (23) are in or near this pocket.
ProteinsWild type (WT) HPV11 E2 TAD with an N-terminal polyhistidine tag (His-TAD, sequence MGHHHHHH followed by E2 amino acids 2-201) was purified from bacteria as described (23). To incorporate selenomethionine, the protein was expressed in Escherichia coli strain B834 (28), a methionine auxotroph, grown in DL30 medium (29) lacking methionine and supplemented with 2 µg/ml biotin and thiamin and 50 µg/ml seleno-D,L-methionine (Se-Met). Purification was performed as described for the WT protein, except that buffers were sparged with helium and all purification steps were completed within 1 day. Addition of a 4-lysine tail to the C terminus of His-TAD was performed by PCR using a primer encoding four lysines, followed by ligation of the PCR product into the plasmid pET15b (Novagen). The encoded His-Lys-TAD protein was purified as described for His-TAD. Electrospray mass spectrometry was used to verify the molecular weights of these three proteins and to determine that they lacked the N-terminal methionine. Complete incorporation of Se-Met was also verified by amino acid analysis. Purified proteins were >95% pure by SDS-PAGE. Protein concentrations were determined by absorbance at 280 nm using a molar extinction coefficient of 49,060 M-1 cm-1. Site-directed mutagenesis of HPV11, HPV6, and HPV31 E2 or of His-Lys-TAD was performed using the QuikChange kit (Stratagene) according to instructions supplied by the manufacturer. E2 proteins were synthesized in vitro by coupled transcription/translation using the TNT reticulocyte lysate system (Promega). InhibitorsInhibitors were synthesized by the same procedures as described for other compounds in the same series (22). Single enantiomers were isolated by chiral chromatography. Radiolabeled inhibitor 3 was synthesized from an enantiomerically pure Boc-protected amine precursor. The acetate moiety, tritiated at the methyl group, was incorporated using [3H]N-acetoxyphthalimide (American Radiolabeled Chemicals Inc.). Specific activity of the product was 9.2 Ci/mmol. Size-exclusion ChromatographySize-exclusion chromatography was performed using the SMART System (Pharmacia), on a Superdex-75 (PC 3.2/30) column in buffer composed of 25 mM Tris-HCl (pH 8.0), 500 mM NaCl, 5 mM TCEP, 0.1 mM EDTA at a flow rate of 40 µl/m. (All buffers were adjusted at room temperature to the stated pH.)
Analytical UltracentrifugationSedimentation velocity and equilibrium experiments were performed at 20 °C using an XL-A analytical ultracentrifuge (Beckman-Spinco) with an An-60 Ti rotor, in 1.2-cm charcoal-filled Epon centerpieces (double sector for velocity and six-sector for equilibrium) and quartz windows. TAD proteins were exchanged into a buffer composed of 20 mM Tris (pH 8.0), 500 mM NaCl, 2 mM TCEP, 0.1 mM EDTA. Velocity data were collected at 60,000 rpm at 0.005-cm intervals, with 5 replicate readings per point. Protein concentrations were Isothermal Titration CalorimetryTitration of inhibitor 1 into a solution of HPV11 His-Lys-TAD was performed using a VP-isothermal titration calorimetry instrument (MicroCal Inc.) in a buffer composed of 20 mM Tris (pH 7.5), 100 mM NaCl, 0.5 mM TCEP, 0.1 mM EDTA, and 1% (v/v) Me2SO, at 20 °C, and was analyzed as described previously (23). Data from a blank run, using a solution of compound 1 injected into buffer without His-Lys-TAD, were subtracted from the inhibitor-protein run prior to fitting the data to a one-site binding model. Crystallization and Data CollectionFinal crystals of apo-HPV11 His-TAD containing either Met or Se-Met were obtained at 4 °C using micro-seeding in a hanging drop vapor diffusion plate. The precipitant solution was composed of 0.1 M sodium succinate (pH 5.0), 18% PEG 5000 monomethyl ether, and 0.2 M (NH4)2SO4. Diffraction data on the Se-Met E2 TAD were collected on an ADSC Q4 CCD at the National Synchrotron Light Source beamline X4a (Brookhaven National Laboratory, NY). Four datasets were collected from a single crystal, cooled to 100 K, at four different wavelengths near the selenium absorption edge (0.9790, 0.9794, 0.9743, and 0.9879 Å). The first dataset was selected for structure refinement; data collection statistics for this dataset are listed in Table I.
The complex of His-Lys-TAD with inhibitor 1 was formed by adding 1 µl of a 60 mM solution of 1 (BILH 434) in Me2SO to 74 µl of His-Lys-TAD at 10 mg/ml in the purification buffer, and incubating the mixture on ice for 2-3 h. Crystals were then obtained by vapor diffusion, using a precipitant solution composed of 0.1 M sodium citrate (pH 5.5) and 35% methyl-2,4-pentanediol. Diffraction data were collected on the National Synchrotron Light Source beamline X25 (Brookhaven National Laboratory, NY). Data from a single crystal cooled at 100 K were recorded on a Brandeis B4 detector (Brandeis University) mounted on a -axis goniometer (Enraf-Nonius, The Netherlands). Data collection statistics are listed in Table I. Structure Determination and RefinementAll data were processed with the HKL software Denzo and Scalepack (Table I). Initial phasing of the Se-Met HPV11 His-TAD Multiwavelength Anomalous Dispersion data was done using the software SHARP (33) with selenium sites located by SHELX (34). Model building into initial electron density maps was carried out using the software O (Alwyn Jones, Uppsala University, Sweden) and model refinement was done with the software CNX (Accelrys). At a later stage, model building was also helped by molecular replacement phasing (CNX) using the published HPV16 coordinates (Protein Data Bank code 1DTO [PDB] (20)). 10% of the data was used as a test set for the Rfree calculation. The final model includes residues 3-121 and 128-195 (Table I), but 15 residues were modeled as alanines because of inadequate electron density for their side chains (Lys-6, Asp-9, Glu-17, Gln-48, Glu-66, Lys-107, Lys-120, Val-128, Lys-173, Ser-180, Thr-181, Asn-182, His-183, Ser-190, and Ser-195). The structure of the His-Lys-TAD·1 complex was solved by molecular replacement using the refined HPV11 His-TAD model. This final model includes residues 2-196 (all side chains were modeled as the correct residue), two molecules of inhibitor 1, 27 water molecules, one molecule of Me2SO, and one partially defined buffer molecule (statistics in Table I). The refined atomic coordinates and structure factors have been deposited in the Protein Data Bank (accession codes 1R6K and 1R6N for the TAD and TAD:inhibitor 1 structures, respectively). Computational StudiesThe atoms making up the binding pocket were identified by applying the Site Finder application within the MOE software package (Chemical Computing Group Inc., Montreal, Canada) to the crystal structure of the complex after both inhibitors, as well as solvent atoms, were removed. The probe radii were 1.4 and 1.8 Å. The isolated donor acceptor radius was 3.0 Å and the connection distance was 2.0 Å. Inhibitor-protein binding affinities were evaluated using the empirical scoring functions PMF (35), PLP (36), and LUDI (37). Each inhibitor-protein interaction was evaluated after removal of the other inhibitor as well as all solvent atoms (except one water molecule between His-29 and His-32 that appears to be hydrogen-bonded to both residues). For these calculations, hydrogen atoms were added to all residues and molecules.
E1-E2-Ori Complex Formation AssayThis assay was performed using in vitro translated E2 proteins as described previously (23). Briefly, binding reactions (80 µl) consisted of HPV11 E1-containing insect cell nuclear extract (titrated to give Circular Dichroism (CD) SpectroscopyFar UV CD spectra were acquired between 260 and 187 nm using a J-715 spectropolarimeter (Jasco) in a 0.05-cm cell at room temperature (approximately 23 °C). Proteins were diluted to 5 µM in a buffer composed of 20 mM sodium phosphate (pH 7.6), 50 mM Na2SO4, and 0.5 mM TCEP. A total of 5 accumulations were taken for each sample at a speed of 20 nm/m, a resolution of 0.2 nm, a response of 2 s, and a bandwidth of 1.0 nm. Differential Scanning CalorimetryMelting temperatures of TAD proteins were determined using a VP-DSC (MicroCal Inc.). TAD proteins were exchanged into a buffer composed of 50 mM Tris (pH 7.6), 150 mM NaCl, and 0.5 mM TCEP, adjusted to 5 µM, filtered, and degassed. Temperature was scanned between 10 and 65 °C at 1.5° min-1. Unfolding under these conditions was irreversible for all TAD proteins, but relative denaturation temperatures were estimated, after buffer subtraction, by fitting Cp versus temperature curves to a modified two-state model (38), using the Origin software supplied with the instrument. It is acceptable to use the equations of equilibrium thermodynamics to analyze unfolding events that are in practice irreversible, because irreversibility typically results from secondary, slower events, such as aggregation (39).
Binding of Tritiated 3 to TAD ProteinsTitrations were performed in 96-well nickel chelate Flashplates (PerkinElmer Life Sciences) using 50 nM His-Lys-TAD protein and 4-500 nM tritiated inhibitor 3, in a buffer composed of 25 mM MOPS (pH 7.0), 100 mM NaCl, 0.5 mM TCEP, 0.0025% (v/v) Tween 20, and 2% Me2SO. After a 3-h incubation to allow for equilibration of interactions between the His-Lys-TAD and both the nickel plate and inhibitor 3, signals from bound 3 were detected by scintillation counting using a TopCount NXT (PerkinElmer Life Sciences). Blank titrations of compound 3 into wells lacking His-Lys-TAD were performed in parallel and signals from nonspecific binding of the compound was subtracted prior to data analysis. Signals ranged from 30 to 4000 cpm for negative controls and 120 to 7000 for WT TAD (prior to blank subtraction). Binding data were analyzed by non-linear regression (GraFit 3.0, Erithicus Software Ltd.) using the quadratic equation,
where Max denotes the counts/min observed at saturation. ATPase AssaysATPase activity of purified recombinant E1 (40) was measured using a scintillation proximity assay as described previously (40, 41). The concentration of HPV11 E1 was 3 nM and ATP and magnesium acetate were used at 20 and 500 µM, respectively. After subtraction of background measured in reactions lacking E1, inhibition data were fit to a three-parameter logistic using SAS (release 6.12, SAS Institute Inc., Cary, NC). TAD mutants Y19A and E39A inhibited only weakly, such that their IC50 and Imax parameters could only be estimated.
Expression, Purification, and Characterization of the HPV11 E2 TADThe HPV11 E2 TAD (amino acids 1-201, Fig. 1) was expressed in E. coli as an N-terminal polyhistidine-tagged protein. The purified protein eluted in a sharp peak when analyzed by analytical size-exclusion chromatography (Fig. 2A) and comparison of its elution time to those of molecular mass standards gave an estimated mass of 24.7 kDa, close to that expected for the monomer based on its sequence (23.8 kDa).
The HPV11 TAD also appeared to be monomeric by analytical ultracentrifugation. In an equilibrium sedimentation experiment, no significant trends in apparent mass were observed at increasing rotor speed or loading concentration (Table II). Such a trend would be observed if the protein formed a mixture of monomers and dimers under experimental conditions (42). Data from a centrifugation velocity experiment were best fit to an estimated average molecular mass of 22.5 kDa for a single species, based on individual parameters of S = 1.960 ± 0.003 Svedberg and D = 8.26 ± 0.03 Frick. Collectively, these results indicate that the HPV11 E2 TAD is fully monomeric in solution over the concentration range studied, up to 20 µM.
An even more soluble version of the HPV11 TAD containing 4 lysines at its C terminus was also purified and was also found to be monomeric by size-exclusion chromatography and analytical ultracentrifugation (data not shown). Binding of inhibitor 1 (Fig. 3) to this Lys-tailed TAD was analyzed by isothermal titration calorimetry. Formation of a 1:1 complex was observed with a Kd value of 40 nM (Fig. 2C).
Crystal Structure of the Apo-HPV11 E2 TADThe structure of the HPV11 TAD was determined from Multiwavelength Anomalous Dispersion diffraction data obtained using crystals of the protein labeled with Se-Met (see "Experimental Procedures"). Overall the structure of the HPV11 TAD shown in Fig. 4 is very similar to that previously reported for the analogous domain of HPV16. The all-C superposition of the HPV11 and -16 structures (187 residues) yielded an r.m.s. deviation of 1.54 Å, and all secondary structural elements are conserved. The major differences are a slight twist (<1°) in the relative orientation of the -helical and -sheet subdomains and to variations in the position of loops, in particular those comprised of residues 140-143 and 160-165, which are involved in crystal packing in the HPV16 and HPV11 structures, respectively. The N-terminal (amino acids 3-94) and C-terminal (amino acids 95-195) subdomains can be individually superimposed with r.m.s. deviations of 0.85 and 1.14 Å, respectively. Interestingly, the HPV11 protein crystal lattice was distinct from that reported for HPV16, such that a different monomer-monomer interface was observed in the HPV11 crystal structure.
Crystal Structure of an HPV11 E2 TAD-Inhibitor ComplexWe obtained crystals of the more soluble Lys-tailed TAD in complex with compound 1, albeit using different conditions than those used to crystallize the apo-His-TAD. The space group of these crystals was also different from any of the previous TAD structures (Table I). The conformation of the HPV11 TAD within the complex (Fig. 5A) is very similar to that observed for the apoprotein; all backbone C could be superimposed with an r.m.s. deviation of 0.88 Å (187 residues). The two inhibitor 1 molecules are both in the free carboxylate form (as drawn in Fig. 3), one of two distinct forms that this class of inhibitors adopts in solution (43). The presence in the crystal structure of the free carboxylate form confirms a prediction that this form is the one with inhibitory activity (43).
Because we had demonstrated that the TAD and inhibitor 1 form a 1:1 complex in solution, we were surprised to find two molecules of inhibitor bound per protein monomer (Fig. 5A). We believe that the molecule referred to below as inhibitor A is the one that binds tightly in biological assays (see below). In contrast, inhibitor B appears to interact only weakly with a secondary binding pocket, formed not only by the TAD protein alone, but also in part by one surface of inhibitor A, as well as by residues from an adjacent protein monomer in the crystal lattice. In fact, inhibitor B has approximately equal surface area in contact with each of the two TAD molecules, as well as with inhibitor A (data not shown). The close contact between the two inhibitor molecules may distort the observed interaction between the TAD and the dichlorophenyl ring of inhibitor A (see below). The TAD-TAD interface that incorporates inhibitor B is distinct from that in the apo structure, and in fact does not result in formation of a dimer, but rather an infinite array (along the 41 screw axis of the crystal), which is unlikely to exist in solution except at very high protein and inhibitor concentrations.
To further demonstrate that the region of the TAD near inhibitor A constitutes a better binding pocket than does the one near B, a computational analysis of the protein surface was performed using the Site Finder application within the MOE software package (see "Experimental Procedures"). This program uses an In a separate computational study, binding energies for each inhibitor molecule to the TAD (in the absence of the other) were predicted using three different empirical protein-binding scoring functions (Table III). As expected, results significantly favored binding of inhibitor A relative to B, with predicted binding energies 2-fold higher for inhibitor A. Taken together, these different methods suggest that inhibitor A is solely responsible for inhibition at the low concentrations used in biological assays and that the presence of inhibitor B might be an artifact of crystallization. This suggestion is further supported by the mutational analysis of the inhibitor binding pocket, presented below. Although an artifact, it is possible that the fortuitous presence of this second molecule aided in crystallization of the complex.
Inhibitor A lies within 4.5 Å of 12 residues (marked by asterisks in Fig. 1), and the backbone superposition for these in the apo and complex structures gives an r.m.s. deviation of only 0.28 Å. However, binding of this inhibitor induces significant movement of several amino acid side chains. The indandione moiety of the inhibitor is buried in a pocket that is not present in the apo structure, but results from rotations of the side chains of Tyr-19 ( 1, 85°; 2, 40°) and His-32 ( 2, 90°) (Fig. 5C). In the apo structure, the aromatic rings of His-32 and Trp-33 form a T-shaped edge-to-face structure, with the hydroxyl group of Tyr-19 close enough to form a hydrogen bond with the -nitrogen of His-32. In the complex structure, this interaction is lost and the rings of His-32 and Trp-33 become parallel to each other. Two other side chains, Leu-94 and Glu-100, also adjust their conformations to accommodate the inhibitor (Fig. 5C). Together, these changes in side chain conformation result in the formation of the observed binding pocket (Fig. 6).
The new positions of Tyr-19 and His-32 not only open up the pocket, but also allow the indandione moiety of the inhibitor to become a middle layer of a three-tier -stacking interaction between His-32 and Trp-33. Additional Van der Waals contacts from the Leu-94 side chain result in a highly complementary shape to the deepest part of the pocket. Two other important interactions are made with the backbone nitrogen atoms of Tyr-99 and Glu-100; the distance (2.73 and 2.77 Å) and angle between these atoms and the inhibitor carboxylate oxygens suggests the formation of strong hydrogen bonds (Fig. 5C). Remaining contacts of inhibitor A with the pocket are less well defined hydrophobic interactions. The dichlorophenyl ring lies in a wider but shallower portion of the pocket without making obvious specific interactions. Although it has been found in previous studies (22) that both chlorine atoms contribute significantly to binding, neither atom is very close to the protein; the closest one being 4 Å from the Tyr-19 side chain. The observed position of the dichlorophenyl ring may be due in part to distortions induced by crystal packing or the presence of the second inhibitor molecule, as discussed above. The thiadiazolylphenyl group lies along the protein surface outside of the pocket. Notably, the phenyl-thiadiazole bond is close to the side chain of Met-101. This finding is consistent with recent photoaffinity labeling experiments2 in which the carbon atom para to the amide bond of an analogue was found to be covalently linked to the terminal carbon of Met-101 (see Fig. 6).
Mutational Analysis of the Inhibitor Binding PocketSeveral residues in or near the inhibitor binding pocket were mutated to examine the roles of specific side chains in binding the indandione, as well as gain insight into the possible overlap between the E1- and inhibitor-binding sites. Most residues were mutated to Ala, but a few residues were replaced by amino acids found in other HPV types: His-32 was changed to Tyr as found in HPV16 E2, whereas Ile-30, Met-31, and Ile-36 were changed to Val, Leu, and Met, respectively, as found in HPV6 E2. These last three residues are the only ones in the vicinity of inhibitor A that are different between HPV11 and -6 E2. Also included in this study was the amino acid substitution I73A that is located on the opposite face of the protein and was previously shown to specifically affect the transactivation activity of E2 rather than its replication or E1 binding activity (25). Each mutant protein was produced by in vitro translation and tested in our previously described E1-E2-ori complex formation assay using radiolabeled origin DNA (Table IV). The ability of each protein to interact with E1 was scored qualitatively, because we did not control for variations in the level of protein produced in the in vitro transcription-translation reaction and, in any case, could not assess the proper folding or stability of the proteins. Most mutant proteins retained at least some ability to productively interact with E1 at the origin.
We also determined the sensitivity of these mutant proteins to inhibitor 2 (Fig. 3), an analogue of compound 1 that we have described previously (23). A number of substitutions had little if any effect on interaction with E1 or on inhibitor binding. These included V64A and K68A, which interact with inhibitor molecule B in the inhibitor-TAD complex, thereby providing further evidence that binding of this second molecule is not relevant to inhibition. Substitutions I30V/M31L and I36M, which make the HPV11 TAD more similar to that of HPV6 E2, also had little effect on inhibition. These results suggest that none of these three residues are responsible for the 30-fold weaker affinity of HPV6 E2 for compound 2 (23). In further support that these residues do not play a major role in the type selectivity of the indandione inhibitors, we found that the corresponding mutations (V30I/L31M and M36I) in HPV6 E2 also had little effect on inhibition (data not shown). Other substitutions (L15A, Y19A, E39A, and S98A) produced inactive or weakly active proteins. In principle this could indicate either that these side chains are directly involved in interaction with E1, or that the mutant proteins are not properly folded. The low signals obtained with the Y19A and E39A E2 proteins in the E1-E2-ori complex formation assay were insufficient to assess inhibition. Substitution of Gln-12, His-32, Leu-43, Ser-65, and Ile-73 had a more deleterious effect on compound binding (IC50 values) than on interaction with E1, suggesting that these residues may have a specific role in binding the inhibitor. However, three of these substitutions are at the periphery of the inhibitor-binding site (Gln-12, Leu-43, and Ser-65), and Ile-73 is on the opposite face of the protein. A possible rationale for these observations is discussed below. In contrast, His-32 is in direct contact with the biologically relevant inhibitor A. The presence of a tyrosine instead of a histidine at position 32 of the E2 proteins of HPV16, -18, and -31, (Fig. 1) may contribute to the insensitivity of these proteins to the indandione class of inhibitors. Consistent with this proposal, we found in other studies3 that an HPV31 E2 carrying the double amino acid substitution D31M/Y32H was significantly more sensitive to compound 2 (IC50 = 6 µM) than its WT counterpart (IC50 = 70 µM), when tested in an assay format similar to that used to generate the results reported in Table IV. Finally, a particularly interesting result, which we have noted previously (23), was the improved affinity of E100A E2 for this class of inhibitors.
Biophysical and Functional Characterization of Purified Mutant E2 TADsBecause the roles of Tyr-19 and Glu-39 in the binding of E1 or inhibitors could not be assessed using the E1-E2-ori complex formation assay, we produced the corresponding mutant TAD proteins in E. coli to carry out additional assays and biophysical studies. We also expressed the L15A TAD, for which a moderately reduced affinity for E1 and inhibitor 2 were observed, as well as the I73A and E100A proteins, to better understand the significant effects of these two mutations on inhibitor binding. Analysis of all these proteins by circular dichroism indicated no significant change in secondary structure compared with the WT (data not shown). The thermal stability of each mutant protein was determined by differential scanning calorimetry and found to be comparable with that of the WT TAD, with the L15A protein being reproducibly the least stable (Table V).
Next, we measured the affinity of each TAD protein for compound 3 (Fig. 3) using an equilibrium binding assay. The Kd values for binding of 3 to WT and Y19A TADs were similar (Fig. 7A and Table V). Because the side chain of Tyr-19 undergoes a substantial change in position upon inhibitor binding, the lack of effect of the Y19A mutation on the affinity for inhibitor 3 may actually be due to a combination of two offsetting effects: a loss of Tyr-19-inhibitor interactions partially compensated for by a reduction in the conformational change (entropic cost) necessary to form the inhibitor binding pocket. Further mutagenesis of this residue would be necessary to address this possibility. The Kd value for 3 binding to E100A was roughly 10-fold lower, but could not be precisely determined in part because the value of Kd is well below the protein concentration used in the experiment (5 versus 50 nM). I73A, L15A, and especially E39A had significantly reduced affinities for the compound. Where a direct comparison can be made, results in Table V are in qualitative agreement with those from the E1-E2-ori complex formation assay shown in Table IV.
We assessed the relative affinity of each of these TAD proteins for E1 by taking advantage of our earlier finding that the ATPase activity of E1 is inhibited by addition of the transactivation domain. TAD binding significantly increases the apparent Km of E1 for ATP. However, we observed that TAD-saturated E1 retained some ATPase activity. We also reported a corresponding weakening of E2 binding in the presence of high ATP concentrations, which might play a role in the dissociation of E2 and assembly of E1 hexamers at the origin (40). Although not a direct binding assay, monitoring ATPase activity in the presence of E2 (or the TAD) provides a DNA-independent, solution-phase method of assessing the E1-E2 interaction. Titration of WT TAD into E1 ATPase assays inhibited activity to a maximum of 86% (Imax = 86%), with an IC50 of 4 nM (Fig. 7B and Table V). Similar results were obtained for I73A and E100A confirming that these mutations do not affect E1 binding substantially. It should be noted, however, that the IC50 values reported for I73A, E100A, and WT TAD are close to the E1 concentration used in these assays (3 nM), and hence that these IC50 values are likely to be an underestimate of the affinity of these proteins for E1. The remaining three mutant E2 proteins (L15A, Y19A, and E39A) differed from their WT counterpart in two ways. First, the maximal degree of E1-ATPase inhibition obtained at saturating concentrations of TAD was lower for these mutants than for WT E2. This was more prominent for Y19A and E39A E2, which achieved a maximum inhibition of only 32 and 38%, respectively. Although we do not completely understand why these mutant E2 TADs achieve a lower maximum level of inhibition it suggests that the geometry of the TAD-E1 complexes that do form for these two proteins is such that the complexes retain greater residual ATPase activity. The second way in which these three mutant proteins differed from WT E2 is in their lower affinity for E1, as measured by the IC50 values required to reach their respective maximal level of ATPase inhibition. By this measure, L15A was found to reduce only slightly, if at all, the affinity of E2 for E1, whereas Y19A and E39A had a larger effect. Qualitatively, these results are in good agreement with those obtained in the E1-E2-ori complex formation assay (Table IV) in that the E1-E2 interaction is little affected by I73A and E100A, moderately affected by L15A, but more significantly reduced by Y19A and, as anticipated, by E39A. Notably, all mutant proteins retained some ability to interact with E1.
Structure of HPV11 E2 TAD The tertiary structure of the HPV11 E2 TAD is nearly identical to that of the analogous domain from HPV16 (20). In particular, the orientations of the two subdomains are similar for both TADs, suggesting that the relative position of these subdomains is a fundamental characteristic of this -fold, and that the overall structure is fairly rigid even though it is far from globular. This idea is further supported by the fact that neither the binding of the indandione inhibitor nor the differing crystal packing arrangements of the three structures significantly altered the observed tertiary structure.
Based on the structure of the HPV16 E2 TAD, Antson et al. (20) proposed the existence of an additional E2 dimerization interface within the TAD, comprised of the second and third
Nature of the Inhibitor Binding Site
Conclusions Based on Studies with Mutant E2 Proteins The Side Chain of Glu-100 May Impede Inhibitor Binding
Amino Acids Outside the Pocket Can Reduce Inhibitor Binding
Importance of His-32 for Inhibitor Binding
A Role for Tyrosine 19 in Binding to E1
Analysis of the two crystal structures presented here has provided a detailed understanding of the mechanism by which the indandione inhibitors bind to E2, to disrupt the E1-E2 protein-protein interaction. Comparison of these two structures revealed the induced fit mechanism by which the inhibitor binding pocket is formed and further illustrates that it is very difficult to accurately identify binding pockets or model protein-ligand interactions using only the structure of an unliganded protein (48). The binding pocket includes both a deep, narrow cavity filled by the indandione moiety, as well as a wider shallow region that makes only a few contacts with the inhibitor. The finding that a second inhibitor molecule can bind to this shallow region of the pocket raises the possibility that molecules could be found that would retain specific interactions with the deep cavity but gain additional binding energy by interaction with the larger shallow region. The characterization of a well defined small molecule binding pocket on the E2 TAD represents an important step toward the design of HPV antivirals.
The atomic coordinates and structure factors (code 1R6K and 1R6N) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).
* 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 Biological Sciences, Boehringer Ingelheim (Canada) Ltd., 2100 Cunard St., Laval, Quebec H7S 2G5, Canada. Tel.: 450-682-4640 (ext. 4269); Fax: 450-682-4642; E-mail: pwhite{at}lav.boehringer-ingelheim.com.
1 The abbreviations used are: HPV, human papillomavirus; His-TAD, TAD with an N-terminal polyhistidine tag; His-Lys-TAD, His-TAD with a C-terminal four-lysine tail; Me2SO, dimethyl sulfoxide; MOPS, 3-(N-morpholino)propanesulfonic acid; ori, origin of replication; r.m.s., root mean squared; Se-Met, selenomethionine; TAD, E2 transactivation domain; TCEP, tris(2-carboxyethyl)phosphine; WT, wild type.
2 Davidson, W., McGibbon, G. A., White, P. W., Yoakim, C., Hopkins, J. L., Guse, I., Hambly, D. M., Frego, L., Ogilivie, W. W., Lavallée, P., and Archambault, J. (2004) Anal. Chem. 76, in press.
3 P. W. White, R. Coulombe, and M. J. Massariol, unpublished data.
4 E. Abbate and M. Botchan, personal communication.
We thank Dr. Michael Cordingley for stimulating discussions and Drs. Nathalie Goudreau, Christopher Lemke, and Steven Mason (Boehringer Ingelheim (Canada) Ltd.), as well as Prof. Liang Tong (Columbia University) for critical reading of the manuscript. We also thank Eric Abbate and Prof. Michael Botchan for communicating results prior to publication.
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||