Learn more: PMC Disclaimer | PMC Copyright Notice
Allosteric Inhibitors of the Eya2 Phosphatase Are Selective and Inhibit Eya2-mediated Cell Migration*
Background: The phosphatase activity of Eya is important for transformation, invasion, migration, and metastasis of breast cancer cells.
Results: A class of N-arylidenebenzohydrazide compounds specifically inhibits the phosphatase activity of Eya2 but not Eya3.
Conclusion: This class of compounds likely acts through an allosteric mechanism.
Significance: These inhibitors may be developed into chemical probes or anti-cancer drugs.
Abstract
Eya proteins are essential co-activators of the Six family of transcription factors and contain a unique tyrosine phosphatase domain belonging to the haloacid dehalogenase family of phosphatases. The phosphatase activity of Eya is important for the transcription of a subset of Six1-target genes, and also directs cells to the repair rather than apoptosis pathway upon DNA damage. Furthermore, Eya phosphatase activity has been shown to mediate transformation, invasion, migration, and metastasis of breast cancer cells, making it a potential new drug target for breast cancer. We have previously identified a class of N-arylidenebenzohydrazide compounds that specifically inhibit the Eya2 phosphatase. Herein, we demonstrate that these compounds are reversible inhibitors that selectively inhibit the phosphatase activity of Eya2, but not Eya3. Our mutagenesis results suggest that this class of compounds does not bind to the active site and the binding does not require the coordination with Mg2+. Moreover, these compounds likely bind within a site on the opposite face of the active site, and function as allosteric inhibitors. We also demonstrate that this class of compounds inhibits Eya2 phosphatase-mediated cell migration, setting the foundation for these molecules to be developed into chemical probes for understanding the specific function of the Eya2 phosphatase and to serve as a prototype for the development of Eya2 phosphatase specific anti-cancer drugs.
Introduction
The Eya 1–4 (eyes absent family of proteins) were first discovered as essential coactivators (1,–4) of the homeobox transcription factor Six1 (2, 5,–11). Both Six1 and Eya proteins are necessary for cellular proliferation in a number of different cell types (2, 8, 12). In support of a cooperative interaction between Six1 and Eya, the Eya1 knock-out (KO) mouse phenocopies the Six1 KO mouse, and organ defects in both KO mice are due to a decrease in proliferation and an increase in apoptosis (13, 14). Both Six1 and Eya1 are linked to branchio-oto-renal syndrome, which is characterized by branchial defects, hearing loss, and renal abnormalities (9, 10). These data suggest that Six1 and Eya cooperate in normal development of human and mouse tissues. Importantly, the Six1-Eya complex is overexpressed in many different tumor types, including Wilms' tumor, ovarian cancer, acute leukemia, malignant peripheral nerve-sheath tumors, glioblastomas, and breast cancer (8, 15,–20). We have recently demonstrated that Eya2 knockdown inhibits the ability of Six1 to induce TGF-β signaling, epithelial-to-mesenchymal transition, and tumor-initiating cell characteristics in breast cancer, properties that are associated with Six1-induced breast tumorigenesis and metastasis (17). We further demonstrated that a direct interaction between Six1 and Eya2 is required for epithelial-to-mesenchymal transition and metastasis (21). Furthermore, analyses of microarray data from 295 breast tumors demonstrated that coordinated overexpression of Six1 and Eya together significantly predict decreased time to relapse and metastasis as well as shortened survival, whereas the overexpression of individual genes does not (17). These clinical data strongly suggest that the Six1-Eya complex mediates breast tumorigenesis, and confers poor prognosis in breast cancer.
Eya proteins also contain a unique protein tyrosine phosphatase activity in their highly conserved Eya domain (ED).4 This phosphatase activity has been demonstrated to switch the Six1 transcriptional complex from a repressor to an activator for some Six1-induced genes (2). In Drosophila, although the phosphatase activity of Eya is not globally required for the ability of Six1 to induce transcription, it is required to induce transcription of a subset of genes (1). The Eya proteins therefore represent the first transcriptional co-activators with intrinsic phosphatase activity to modulate transcription (2, 3, 22). Eya proteins are the founding members of a class of non-thiol-based protein phosphatases belonging to the haloacid dehalogenase (HAD) family of enzymes. The HAD family of phosphatases use an aspartic acid (Asp) as the catalytic active-site residue instead of the more commonly used cysteine (Cys) used by most cellular phosphatases (23). Some HAD members, such as Scp1 and chronophin, dephosphorylate protein substrates, whereas substrates of many HAD phosphatases are small molecules (3, 23). All other known HAD protein phosphatases are Ser/Thr phosphatases, leaving Eya as the sole member targeting phosphorylated Tyr (23, 24). Structurally, HAD members share a highly conserved catalytic core but contain a variable cap structure. Eya members have a unique helical bundle as the cap that is not observed in any other HAD protein (23, 25).
Recently, the tyrosine phosphatase activity of Eya has been shown to be critical for transformation, migration, invasion, and metastasis of breast cancer cells (26). Although it has been argued that the phosphatase activity of Eya is not essential for its ability to induce proliferation (26), Pestell and colleagues (27) recently demonstrated that the Eya1 phosphatase activity is in fact necessary for Eya-induced proliferation of breast cancer cells. In addition, the phosphatase activity of Eya plays a role in the DNA damage response via its ability to dephosphorylate histone variant H2AX upon DNA damage, leading cells to the repair instead of apoptosis pathway (24, 28). These data suggest that inhibiting Eya phosphatases may not only inhibit tumor onset and progression but may also sensitize cells to commonly used chemotherapeutics that act through inducing DNA damage. Importantly, drugs that target the Eya phosphatases may confer limited side effects as most Eya family members are highly expressed in developing tissues but are not expressed in most normal adult tissues (29).
We have previously identified a series of structurally related small molecule compounds that specifically inhibit the Eya2 phosphatase using high throughput screening (HTS) (30). In this work, we describe characterization of the mechanism of action of these inhibitors, which suggests that this class of compounds functions as allosteric inhibitors of the Eya2 phosphatase activity. We also demonstrate that this class of compounds inhibits Eya2-mediated migration of breast cells. These studies set the foundation for further developing these inhibitors into chemical probes or potential therapeutic agents targeting the Eya phosphatase.
EXPERIMENTAL PROCEDURES
Protein Expression and Purification
Eya2 ED (residues 253–538) was subcloned into a pGEX-6p1 vector (GE Healthcare) and the protein was expressed as a GST-fusion protein in Escherichia coli. The protein was first purified on a glutathione column. GST was then cleaved off using PreScission protease and eluted Eya2 ED was further purified on a Superdex 200 (S200) column (GE Healthcare).
Eya Phosphatase Assays
Eya activity was measured in 50-μl reactions using black, 96-well, half-volume microtiter plates (Greiner Bio-one) with the substrate OMFP (3-O-methylfluorescein phosphate, Sigma-Aldrich), which is converted to a fluorescent product OMF upon dephosphorylation. The assay conditions used were 50 mm MES, pH 6.5, 50 mm NaCl, 1.25 mm MgCl2, 0.05% BSA, and 1 mm DTT, and the reaction contained 50 nm Eya2 ED and 100 μm OMFP. Reactions were started by the addition of OMFP, preceded for 1 h at room temperature followed by terminating the reaction with the addition of 75 mm EDTA. Fluorescence intensity was measured at 485/515 nm excitation/emission on a Fluoromax-3 plate reader (Horiba Jobin Yvon).
Isothermal Titration Calorimetry
Eya proteins were purified via S200 in isothermal titration calorimetry (ITC) buffer (50 mm NaPO4 buffer, pH 6.5, 50 mm NaCl). Inhibitor solutions were prepared by diluting from a DMSO stock into ITC buffer. DMSO was added to the protein solution to the same concentration as the inhibitor solution (0.2%), and all solutions were degassed. ITC experiments were performed with a VP-ITC microcalorimeter (MicroCal, Inc., GE Healthcare) at 25 °C. Due to low solubility of inhibitors in aqueous buffer, we injected Eya protein into the sample cell containing inhibitor, and the reference cell was filled with ITC buffer. After a 2-μl initial injection, 10-μl aliquots of 250 μm protein was added stepwise at 5-min intervals into the sample cell containing 20 μm inhibitor until saturation. Origin (version 7.0) was used to analyze data with embedded calorimetric fitting programs to generate the binding curves.
Enzymatic Kinetics
All kinetic experiments were performed using 384-well black medium binding plates (Greiner Bio-One). Fluorescein diphosphate (FDP, Anaspec) was prepared as a 100 mm stock solution in Eya2 ED assay buffer (25 mm HEPES, 50 mm NaCl, 5 mm MgCl2, pH 6.5) and stored in single-use aliquots at −80 °C and protected from light. Both Eya2 ED and FDP stock solutions in assay buffer were diluted to working solutions in Eya2 ED buffer with freshly added 5 mm DTT and 0.33% BSA prior to the assay. To calculate Km values and generate a Lineweaver-Burke plot for Eya2 ED phosphatase activity with FDP as substrate, we performed kinetic measurements at various concentrations of FDP and compound. Briefly, 50 nl of compound (NCGC00249327) at 10, 3.3, 1.1, 0.37, 0.12, 0.04, 0.014, or 0 mm in DMSO were added by ATS Acoustic Liquid Dispenser (EDC Biosystems) to empty wells of a plate. Immediately after, 5 μl of 0.5 μm Eya2 ED were added and the solutions incubated for 10 min at room temperature. Next, 5 μl of FDP substrate solutions at 8, 4, 2, 1, 0.5, 0.25, 0.125, or 0.063 mm were added to the wells using a CyBi®-well 384-channel simultaneous pipettor (CyBio, U.S., Inc.). The plate was briefly spun down, and fluorescence intensity was measured every 5 min for a 1-h time course using a ViewLux Imager (PerkinElmer, Inc.) with the following settings: excitation wavelength of 485 nm, emission wavelength of 525 nm, energy light of 750, and exposure time of 0.5 s. Kinetic analysis was performed using GraphPad Prism (version 4, GraphPad Software).
Reversibility Assay
The enzyme at a concentration that is 100-fold (100 nm) more than what is required for the activity assay (1 nm) is incubated with a concentration of inhibitor (40 μm) equivalent to 10-fold of the IC50. The enzyme-inhibitor complex is then diluted 100-fold and substrate is added to initiate the enzymatic reaction. The enzymatic activity at different time points (1–3 h) is compared with that of a similar sample of enzyme incubated and diluted in the absence of inhibitor. If the inhibition is reversible, the enzyme activity will recover to roughly the uninhibited level (the inhibition could be slowly reversible where the enzyme will take longer to reach full activity after dilution). If the inhibition is irreversible, the enzymatic activity will remain very low after dilution because the compounds have irreversibly inactivated the enzyme.
Cell Lines
Stable integration of full-length human Eya2 or phosphatase-dead Eya2 (D274N) in MCF10A cells was achieved through retroviral transduction. Eya2 or D274N was cloned into pMSCV-IRES-YFP backbone, and BOSC cells were used to package viral particles. YFP-positive cells were sorted 1 week after infection.
Motility Assay
Motility was measured using a gap closure assay, where a silicone μ-well insert (Ibidi, Verona, WI) in a 24-well plate was used to create an ∼500-μm gap between 40,000 cells/chamber that were plated overnight. Photos were taken of the gap immediately after removing the insert and adding 10 μm compound (or vehicle control) containing medium and then again 6 h later on a CKX41 microscope (Olympus, Tokyo, Japan). Distance migrated was determined by subtracting the size of the gap at the end time point from the size of the gap at the initial time point, using DP2-BSW software (version 2.2; Olympus). Statistics were calculated with Prism (version 5.0, GraphPad, San Diego, CA).
Thermal Shift Experiments
To provide evidence for the binding between hydrazides and Eya2 ED, thermal shift experiments were performed to evaluate whether Eya2 ED melting temperature (Tm) could be shifted in the presence of different concentrations of hydrazides. Compounds were serially diluted in DMSO in 96-well polypropylene plates (Thermo Fisher Scientific). A mixture of 49 μl of Eya2 ED and SYPRO Orange (5000× stock solution, Invitrogen) were transferred to a 96-well skirted thin-wall PCR plate (Bio-Rad) with a final concentration of 1 μm and 5×, respectively. Compounds prepared in the aforementioned mother plate were subsequently added to the protein-dye mixture, with final concentrations ranging from 0.34 nm to 200 μm (final volume, 50 μl). DMSO alone was also transferred within each compound dilution series as a control sample. DMSO concentration was consistent across each dilution series and was maintained at 2%. The plate was sealed with Optical-Quality Sealing Tape (Bio-Rad) and was centrifuged at 1000 rpm for 20 s before heating commenced. Plate heating was performed on an iQ5 real-time PCR detection system (Bio-Rad) from 20 to 95 °C with an increment of 1 °C and a ramping rate of 0.1 °C/s. A CCD camera was used to measure SYPRO Orange fluorescence through a filter set corresponding to excitation and emission wavelengths of 490 and 575 nm, respectively. Protein melting temperature was obtained from an EXCEL-based worksheet (provided by Structure Genomics Consortium) and GraphPad Prism (version 4).
NMR Spectroscopy
15N-labeled Eya2 ED was expressed in E. coli in the presence of 15NH4Cl in minimal medium and purified similarly as the unlabeled Eya2 ED. Optimal NMR buffer conditions were determined to be 50 mm Bicine, pH 7.5, 50 mm NaCl, 0.5% glycerol. Maximum Eya2 ED concentration used was 150 μm due to aggregation at higher concentrations. HSQC experiments were collected at 25 °C on a Varian 900 MHz at a concentration of 150 μm Eya2 ED. Compound was added to saturate Eya2 ED while keeping DMSO concentration below 0.1% DMSO.
UV-visible Spectra Analysis of Selected Hydrazides
Compounds were dissolved in acetonitrile at a final concentration of 25 μm with final Mg2+ concentration at 0, 0.5, 1, 5, 10, 25, 50, 100, and 200 mm. UV-visible spectra were obtained using an Agilent 8453 UV-visible Spectrophotometer (Agilent Technologies) in 0.1-cm quartz cuvette at room temperature.
Molecular Docking
The AutoDock program (version 4.0) was used to dock the compounds to the active site of the crystal structure of Eya2 ED (Protein Data Bank code 3HB1). The active site of the protein was defined by a grid of 70 × 70 × 70 points with a grid spacing of 0.375 Å centered at the Mg2+ ion. The Lamarckian genetic algorithm was applied with 100 runs, and the maximum number of energy evaluations was set to 2 × 106. Results differing by less than 1.5 Å in positional root mean square deviation of substrate were clustered, and the final binding conformations were represented by the one with the most favorable free energy of binding. The optimal binding complexes were subjected to a stepwise energy minimization and molecular dynamics simulations.
Docking at the allosteric site was performed with the flexible docking protocol (31) of Accelrys Discovery Studio (version 3.1, Accelrys Software) using the crystal structure of Eya2 ED (Protein Data Bank code 3GEB). Predicted binding-induced changes in protein structure were determined by the conjugate gradient minimization protocol using a CHARMM force field and the Generalized Born implicit solvent model (32). All minimizations converged to a root mean square gradient of <0.001 kcal/mol.
RESULTS
HTS Identified a New Chemical Series That Specifically Inhibits the Phosphatase Activity of Eya2
We have identified a series of structurally related compounds containing an N-arylidenebenzohydrazide core that inhibits the phosphatase activity of Eya2 with varying potency using an OMFP-based phosphatase assay and the Eya2 ED (structures of selected members of this class are shown in Fig. 1A) (30). Initial enzyme kinetics studies indicated that MLS000544460 exhibits a mixed mode behavior between competitive and non-competitive inhibition (30). Additional kinetic studies indicate that the competitive nature of these molecules depended on the concentration of Mg2+ in the assay. At a high concentration of Mg2+ (5 mm), the compounds display a non-competitive behavior (Fig. 1B), whereas at a low concentration (5 μm), they behave in a mixed mode manner (30). All of our enzymatic assays reported in this work were performed at high Mg2+ concentrations where the compounds are non-competitive.
![An external file that holds a picture, illustration, etc.
Object name is zbc0281486820001.jpg An external file that holds a picture, illustration, etc.
Object name is zbc0281486820001.jpg](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4047403/bin/zbc0281486820001.jpg)
Eya2 ED enzyme kinetic experiments indicate that the N-arylidenebenzohydrazide-containing compounds are non-competitive inhibitors at high Mg2+ concentrations (5 mm). A, chemical structures of representative compounds identified directly in HTS (MLS000544460 and MLS000585814), a commercially available compound that has slightly higher potency (NCGC00249327), and a structurally related but inactive compound (NCGC00241224). B, saturation curve for Eya2 ED showing initial rate with FDP concentration at each concentration of inhibitor NCGC00249327 (top). A Lineweaver-Burk plot shows non-competitive inhibition by NCGC00249327 given that the plots have the same x intercept as uninhibited enzyme and differing slopes and y intercepts. RFU, relative fluorescence units.
We have previously shown that these compounds do not inhibit three other cellular phosphatases (Scp1, PPM1A, and PTP1B), demonstrating specificity toward the Eya2 ED (30). Because there are multiple Eya family members with high homology in the ED, we tested the ability of this series of compounds to inhibit other Eya family members. Using the OMFP-based phosphatase assay, we demonstrated that this class of compounds has no discernible inhibition against the Eya3 ED (we have not been able to express and purify sufficient quantities of human Eya1 or Eya4 ED to test against these other Eya family members) (Fig. 2A). To further validate this selectivity, we measured the binding affinity of compound MLS000544460 to Eya2 ED and Eya3 ED using ITC. MLS000544460 binds to the Eya2 ED with a KD of 2.0 ± 0.3 μm, but it does not bind Eya3 ED, indicating that indeed this class is highly selective toward Eya2 (Fig. 2B).
MLS000544460 Inhibits Eya2-mediated Cell Migration
We next generated stable human Eya2, Eya2 phosphatase-dead (D274N), and YFP (as a negative control) expressing MCF10A cells. We analyzed Eya2-dependent motility using a gap closure assay, in which each of the three cell lines were treated with vehicle control (DMSO), an inactive compound NCGC00241224, and the active compound MLS000544460 at 10 μm. We found that Eya2 induces an increase in migration over the YFP control, whereas phosphatase-dead Eya2 does not, indicating that the phosphatase activity of Eya2 is important for the migration of MCF10A cells, consistent with previous reports (33) (Fig. 3, the DMSO columns). Compound MLS000544460 inhibits Eya2-mediated cell migration, whereas the inactive compound NCGC00241224 does not (Fig. 3). The above data indicate that the Eya2 phosphatase is required for cell migration and that this migration can be inhibited by our lead compound targeting the phosphatase activity of Eya2.
![An external file that holds a picture, illustration, etc.
Object name is zbc0281486820003.jpg An external file that holds a picture, illustration, etc.
Object name is zbc0281486820003.jpg](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4047403/bin/zbc0281486820003.jpg)
MLS000544460 inhibits Eya2-dependent cell migration. Compound MLS000544460 (10 μm) inhibits Eya2 phosphatase-dependent migration of MCF10A cells, in a gap closure migration assay for 6 h, whereas the structurally similar but inactive compound NCGC00241224 has no effect. *, p < 0.05; ***, p < 0.001, two-way analysis of variance.
MLS000544460 Does Not Require Mg2+ for Binding to Eya2 ED
We next explored the mechanism of action of this class of compounds. One possible mechanism of inhibition for these compounds is to serve as Mg2+ chelators because these compounds contain a number of heteroatoms that could coordinate a Mg2+ ion (Fig. 4A). We evaluated UV spectra of compound MLS00054460 in the presence of increasing concentrations of Mg2+ to examine this possibility. There is no change in the UV spectra of MLS00054460 in aqueous buffers (including the phosphatase assay buffer) with increasing concentrations of Mg2+, indicating that this compound has very low, if any, affinity to Mg2+ in aqueous solution. To evaluate the capacity of these molecules to complex Mg2+ in a more lipophilic environment, as a mimetic of the enzymatic active site, we further evaluated the UV spectra of 25 μm MLS000544460 in 100% acetonitrile (34) with increasing concentrations of Mg2+. We observed a dose-dependent shift in the maximum absorbance wavelength (λmax) from 322 (0 mm Mg2+) to 334 nm (100 mm Mg2+) (Fig. 4, B and C). Similar studies using the monovalent cation Na+ produced no discernible shift (Fig. 4D), indicating that these compounds have intrinsic affinity for Mg2+ in a more lipophilic environment similar to acetonitrile but not for other cations such as Na+. We subsequently examined the UV spectrum of a low activity inhibitor, MLS000585814 (IC50 = 73 μm) upon increasing Mg2+ concentrations. A similar, dose-dependent shift in λmax was observed from 316 nm (0 mm Mg2+) to 327 nm (50 mm Mg2+) (Fig. 4, E and F). Although this compound has low activity in the inhibition assay, it retains the ability to bind Mg2+ with similar affinity as the more potent compound MLS000544460 (Fig. 4, C compared with F). These data suggest that the binding affinity of this series of compounds to Eya2 ED is the dominating factor for inhibitory activity and not their intrinsic binding affinity to Mg2+.
The interaction between Eya2 ED and the N-arylidenebenzohydrazide-containing compound does not require Mg2+. A, structures of two active compounds with heteroatoms that can potentially coordinate Mg2+ shown in red. B, UV absorption spectrum of compound MLS000544460 in the presence of Mg2+. λmax shifted from 322 nm to 334 nm with increasing concentrations of Mg2+. C, the λmax shift observed in B is quantified as a function of Mg2+ concentration. D, the UV absorption spectrum of MLS000544460 in the presence of Na+. λmax remained at 323 nm when titrated with Na+. E, UV absorption spectrum of the low activity analog MLS000585814 in the presence of Mg2+. λmax shifted from 316 to 327 nm with increasing Mg2+ concentrations. F, the λmax shift observed in E is quantified as a function of Mg2+ concentration. G, the Tm increase caused by the binding of compounds in the presence of 5 mm EDTA correlates with their IC50 values, suggesting that these compounds stabilize Eya2 ED when Mg2+ is removed. H, MLS000544460 binds with lower affinity in the presence of 5 mm Mg2+ (KD = 6.0 ± 1.2 μm). I, MLS000544460 has higher affinity toward Eya2 ED when Mg2+ is removed with 10 mm EDTA (KD, 0.80 ± 0.04 μm). J, the structurally similar but inactive compound NCGC00241224 does not bind Eya2 ED. Abs, absorbance; a.u., absorbance units.
We also carried out protein thermal shift experiments to examine the binding interaction between this class of compounds and Eya2 ED in the presence or absence of Mg2+. These compounds stabilize Eya2 ED (reflected by an increase of Tm) when exogenous Mg2+ is removed by the addition of 5 mm EDTA (Fig. 4G), suggesting that the binding between Eya2 ED and these compounds does not require Mg2+. We observed a good correlation between the increase in Tm and the IC50 values of the compounds in our series (Fig. 4G), indicating that Eya2 binding is the dominant factor for inhibiting the enzyme. Because the increase in Tm values induced by a compound in a chemical series correlates with its affinity toward the protein (35), these data suggest that the compounds with greater affinity toward Eya2 ED have an increased inhibitory capacity (lower IC50).
We next preformed ITC experiments to further examine binding affinities between the inhibitor and Eya2 ED in the presence and absence of Mg2+. We demonstrated that MLS000544460 binds to the Eya2 ED (containing only the Mg2+ carried through purification) with a KD value of 2.0 ± 0.3 μm (Fig. 1C). When 5 mm exogenous Mg2+ is added, the KD shifts to 6.1 ± 1.2 μm, implying that Mg2+ can potentially interfere with compound binding (Fig. 4H). When Eya2 ED is dialyzed into 10 mm EDTA to remove Mg2+, it binds MLS000544460 with a higher affinity (KD = 0.80 ± 0.04 μm), indicating that this class of compounds does not require Mg2+ for its interaction with the Eya2 ED (Fig. 4I). Furthermore, we demonstrated that the structurally related but inactive compound NCGC00241224 does not bind to Eya2 ED in ITC experiments, indicating that the binding of this chemical series is driven by specific interactions of chemical moieties in the active compounds rather than nonspecific binding caused by the common chemotype within this class (Fig. 4J).
MLS000544460 Is a Reversible Inhibitor
Another possible mechanism of action for these compounds is that the hydrazide functional group provides a reactive chemical moiety that irreversibly binds to the enzyme. To evaluate this possibility, we conducted a traditional irreversibility experiment. (Details of the methods are described under “Experimental Procedures.”) We first examined the inhibition capacity of 0.4, 4.0, and 40 μm MLS000544460 after incubation with 1 nm Eya2 ED for 1 h. We observed a clear dose-dependent inhibition with an IC50 ∼ 4 μm. Next, 100 nm Eya2 ED was incubated with 40 μm compound for 10 min. The sample was diluted 100-fold to 1 nm Eya2 ED and incubated for 1 to 3 h (Fig. 5). Substrate was added, and enzyme activity was measured after 20 min. Eya2 ED essentially regains full activity after dilution compared with control (100× (40 μm) versus 100× (no compound)), indicating that these inhibitors are reversible. However, the enzyme slowly regains its activity after the dilution and requires ∼3 h for almost complete functional recovery (compare the 1- and 3-h incubation panel), indicating a slow dissociation rate and a long occupation time on ED. The reduced activity (∼80%) for the 100× (no compound) sample is likely caused by partial precipitation and potential denaturation of the enzyme, observed at high compound concentrations.
![An external file that holds a picture, illustration, etc.
Object name is zbc0281486820005.jpg An external file that holds a picture, illustration, etc.
Object name is zbc0281486820005.jpg](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4047403/bin/zbc0281486820005.jpg)
MLS000544460 is a reversible inhibitor. 1 nm Eya2 ED incubated with various concentrations of MLS000544460 demonstrates that the IC50 of MLS000544460 is ∼4 μm under this condition (first three columns in each panel). The enzymatic activity of Eya2 ED at 1 or 3 h following a 100-fold dilution of 100 nm Eya2 ED incubation with or without MLS000544460 indicate that MLS000544460 is a reversible inhibitor with a slow off rate. Cpd, compound.
Identification of Compound-binding Sites Using Crystallography and NMR
Extensive efforts trying to co-crystallize Eya2 ED with MLS000544460 and to soak MLS00544460 into crystals of Eya2 ED (25) were unsuccessful. We next carried out NMR experiments in an attempt to determine the binding site of these compounds. When the active compound MLS000544460 is added to the Eya2 ED at 200 μm concentration, significant changes are observed in the HSQC spectrum of Eya2 ED, including peak shifts and the appearance and disappearance of peaks (Fig. 6). The inactive compound NCGC00241224 does not produce significant changes in the NMR spectrum compared with vehicle control (0.1% DMSO). Unfortunately, the quality of the HSQC spectrum is insufficient (likely due to Eya2 ED aggregates at high concentrations) for us to assign the HSQC spectrum to identify specific residues in Eya2 ED that interact with the compound.
Mutagenesis Studies Suggest That Compounds Are Not Binding in the Active Site of Eya2
Because crystallography and NMR were unsuccessful in identifying the binding site of these compounds, we used molecular modeling and site-directed mutagenesis to identify potential compound binding sites. As we first embarked on these docking studies, we only had the earlier kinetic data indicating that these compounds are mixed mode inhibitors. In addition, several Eya inhibitors identified by other groups have been proposed to bind in the active site (36) or to a pocket immediately adjacent to the active site (33). We therefore first considered whether our compounds also bind to the active site. Computational docking of molecules of this series in the active site allows their accommodation in several potential binding modes. Fig. 7A shows one example of these docking models, where the compound interacts with Mg2+, forming a hydrogen bond with Arg-293 and an aromatic stacking interaction with His-283. To evaluate the validity of these models, we generated mutants that either abolish predicted specific interactions (H283A, T289A, R293A, T448A, T449A, Q450A, E505A, E506A) or sterically block inhibitor binding (T448F and S521F). In addition, three double-mutants (H283A/R293A, T289A/E505A, and T448A/E505A) were generated in case single amino acid mutations were insufficient to affect compound binding. Active site mutants were expressed and purified in a manner similar to wild type (WT) Eya2 ED and were tested for their phosphatase activity (Fig. 7B) as well as their ability to be inhibited by MLS000544460. Of these mutants, E506A and T449A were enzymatically inactive. Some of the remaining mutants displayed decreased activity, but all could still be inhibited by MLS000544460 with similar IC50 values (representative dose response curves of several mutants are shown in Fig. 7C). These results in combination with our observation that Mg2+ is not necessary for the binding of our compounds to Eya2, suggest that the site of interaction for these molecules is not the active site.
Eya2 inhibitors likely do not bind in the active site. A, a representative docking model of MLS000544460 in the active site. B, phosphatase activity of representative active site mutants of Eya2 ED in the presence of DMSO (vehicle control) or EDTA (fully inhibited). C, representative dose response curves of active site mutations demonstrate no effect on the ability of MLS000544460 to inhibit the mutant Eya2 ED. The phosphatase activity of Eya2 ED in the presence of DMSO is normalized to 1, and the activity in the presence of EDTA is normalized to 0 in these plots. RFU, relative fluorescence units.
Mutation Analyses Suggest That MLS000544460 Binds in an Allosteric Site
Because the compounds selectively inhibit Eya2 over Eya3, we next examined the amino acid sequence differences between Eya2 and Eya3 as a way to determine where the compound binds. Eya2 ED and Eya3 ED have a 68% sequence identity (Fig. 8). We mapped the amino acid differences to the Eya2 ED crystal structure to identify surface regions that potentially harbor the compound binding sites (Fig. 9A). Four sets of mutants were generated that change Eya2 residues to Eya3 on the Eya2 surface (T298P/T300V/S301V/R303G/I304S, T404P/P405Q/W410L/L411Q, L423S/T426G/H427T, and V504R/Q508I) and two additional mutants (Δ359-368 and Δ356-368) were created that deleted various portions of the unstructured region in the crystal structure that is least homologous between Eya2 and Eya3 ED (Fig. 9B). The phosphatase activity of these mutants was tested (Fig. 9C), and dose response curves with compound MLS000544460 were generated. All but one set of mutants (L423S/T426G/H427T) were inhibited by MLS000544460 with similar IC50 values as the WT. Mutant L423S/T426G/H427T significantly increased the IC50 value to 56 ± 12 μm, indicating that compound binding site is greatly affected by the changes of these residues (Fig. 9D). ITC experiments confirm that compound MLS000544460 has dramatically reduced binding to the Eya2 ED carrying mutant L423S/T426G/H427T (Fig. 9E). These results suggest that the compound is binding near the Leu-423/Thr-426/His-427 region.
![An external file that holds a picture, illustration, etc.
Object name is zbc0281486820008.jpg An external file that holds a picture, illustration, etc.
Object name is zbc0281486820008.jpg](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4047403/bin/zbc0281486820008.jpg)
Sequence alignment of Eya2 ED and Eya3 ED. Residues colored gray are not in the ED (which begins at the bar) but were in the Eya construct used for phosphatase, ITC, and NMR experiments. Symbols beneath sequences represent identical residues (*), strongly similar properties (colon), weakly similar properties (semicolon), or no similarity (no symbol). Residues highlighted in red are surface amino acids that were mutated in Fig. 9.
Mutagenesis based on differences between Eya2 and Eya3 ED identifies a potential compound binding site. A, a surface representation of the Eya2 ED structure with surface residues that differ between Eya2 and Eya3 highlighted in different colors. B, Eya2 ED mutants generated to change corresponding Eya2 residues to those of Eya3. C, phosphatase activity of representative mutants of Eya2 ED in the presence of DMSO (vehicle control) or EDTA (fully inhibited). D, dose response curves of Eya2 mutants. One mutant, L423S/T426G/H427T, demonstrated significantly reduced inhibition by MLS000544460, indicating a potential binding site. The phosphatase activity of Eya2 ED in the presence of DMSO is normalized to 1, and the activity in the presence of EDTA is normalized to 0 in these plots. E, ITC experiments demonstrate dramatically reduced affinity between the L423S/T426G/H427T and MLS000544460. RFU, relative fluorescence units.
We next modeled compound MLS000544460 in the pocket near the L423S/T426G/H427T mutant on the opposite side of the active site (Fig. 10A). Energy minimizations of the binding mode indicate that the hydrogen bond interactions involving Glu-418, Thr-426, and Thr-421 at the helix-turn-helix motif between helices 7 and 8 can be disrupted to form a pocket that allows insertion of the pyridine ring of MLS000544460. The compound sits on a largely hydrophobic patch on Eya that involves residues Ile-279, Ile-280, Leu-425, Leu-429, Leu-463, Val-466, and Phe-467. There is also a hydrogen bond interaction between the hydrazone moiety and the main chain of Gly-462. To test this binding model, we designed multiple Eya2 mutants, including L425N, which would be expected to disrupt the hydrophobic surface, and E418W and L417W, which would physically block access to the pocket. These mutants retained enzymatic activity, although at a generally lower activity level than the WT (2–4-fold lower) (Fig. 10B). All of these mutants abolished the inhibition of Eya2 ED by MLS000544460 (Fig. 10C). ITC experiments using L425N and L417W (Fig. 10D) demonstrate no apparent binding, implicating this pocket as the potential allosteric binding site. Mutants R414E and R414A do not have an effect on the inhibition of Eya2 by MLS000544460 (Fig. 10C), suggesting that the side chain of Arg-414 is flexible and not directly involved in compound binding. Substitutions of Arg-414 with bulky residues (R414Y and R414W) caused an increase in the IC50, likely through steric interference (Fig. 10C).
MLS000544460 likely binds to an allosteric site. A, the proposed allosteric compound binding site is on the opposite face of the active site. Compound is shown as a green ball-and-stick model, and the Mg2+ ion in the active site is represented by a green sphere. B, a docking model demonstrating residues in the allosteric binding pocket surrounding the compound. C, phosphatase activity of representative allosteric site mutants of Eya2 ED in the presence of DMSO (vehicle control) or EDTA (fully inhibited). Various mutants were characterized at different times so they are presented separately with their corresponding WT controls in the same experiment. D, effect of mutations in the allosteric binding pocket on the ability of MLS000544460 to inhibit the mutant Eya2 ED. The phosphatase activity of Eya2 ED in the presence of DMSO is normalized to 1, and the activity in the presence of EDTA is normalized to 0 in these plots. E, ITC experiments show that MLS000544460 does not bind Eya2 L425N and L417W.
DISCUSSION
Phosphatases play important roles in multiple disease processes, including cancer, obesity, and type II diabetes (37), yet they have typically been difficult to target, and no phosphatase inhibitors have yet been approved for therapeutic use (38). Difficulty in obtaining small molecules that inhibit phosphatases in a selective manner is a major issue. Because the active sites of phosphatases are generally highly conserved and positively charged, HTS approaches often lead to compounds lacking specificity or that are highly negatively charged, which limits their bioavailability and cell permeability (39).
Eya2 is a member of the unique HAD family of protein Tyr phosphatases that utilize an aspartic acid as the active site residue instead of the cysteine utilized by most cellular phosphatases. We previously exploited this mechanistic difference and carried out a HTS of >330,000 compounds. This screen led to the identification of a series of compounds that demonstrate specificity to Eya2 over other cellular phosphatases, including other HAD family members. This class of compounds displays unexpected high selectivity and does not inhibit or interact with Eya3. This result was surprising as the EDs of Eya2 and Eya3 have 68% sequence identity and 81% similarity, and all of the active site residues are identical or highly conserved. Mutations of residues within the active site of the Eya2 ED revealed no effect on the inhibition of Eya2 phosphatase activity by the compound, suggesting that the compound does not inhibit the phosphatase activity through binding to the active site. This is consistent with the observed compound specificity against Eya2 but not Eya3 as we would expect the compound to inhibit both Eyas if it binds to the active site. The binding between the compound and Eya2 ED also does not require Mg2+ (Fig. 4, G and I). In the absence of Mg2+ (with 5 mm EDTA), the increase in binding between Eya2 ED and compounds in the series (reflected by the increase in Tm upon compound binding) correlates with better inhibition (Fig. 4G). These observations are all consistent with the notion that the compounds do not bind to the active site and coordinate Mg2+.
We identified a potential binding site for the compound on Eya2 ED that is on the opposite face of the active site, by mutating surface residues differing between the Eya2 and Eya3 ED and evaluating the effect of mutations on the inhibitory activity of the compound. Therefore, our compound likely acts through a novel allosteric mechanism that has not been seen with other HAD inhibitors. Two other groups have identified benzbromarone (33) and metal-chelating compounds (36) as Eya2 phosphatase inhibitors. These Eya inhibitors were proposed to use an Mg2+-chelation mechanism in the active site (36) or to bind in a pocket adjacent to the Mg2+ (33). Rabeprazole, the HAD phosphatase inhibitor to Scp1, also binds in a hydrophobic pocket adjacent to the Mg2+ in the active site (40). The allosteric mechanism of our compounds is consistent with the kinetic studies under high Mg2+ concentrations, which suggest that this series of compounds are non-competitive inhibitors (Fig. 1B). At low Mg2+ concentrations (e.g. 5 μm), these compounds may also inhibit Eya2 phosphatase by chelating Mg2+. Because Eya2 requires Mg2+ for enzymatic activity, and because Mg2+ is located in the active site, the compound can also act as a competitive inhibitor under low Mg2+ conditions. This competition disappears at high concentrations of Mg2+, with the compound acting purely as non-competitive inhibitor.
Although we cannot completely rule out the possibility that mutations around the proposed allosteric site indirectly affect a yet to be identified pocket that is the true compound binding site, multiple lines of evidence suggest that the allosteric site is a direct compound binding site. First, we have made many mutations across the Eya2 ED surface, including the active site, which had no effect on the inhibition of Eya2 ED by our compound, serving as negative controls for our mutagenesis experiments in general. In particular, the T404P/P405Q/W410L/L411Q cluster of mutants contain residues Trp-410 and Leu-411, which are located on the same α7 helix as Arg-414, Leu-417, and Glu-418 whose mutations affect compound inhibition. However, the T404P/P405Q/W410L/L411Q cluster of mutants has no effect on compound inhibition. Second, we have generated four mutations at residue Arg-414. Two of these mutants changed Arg-414 to a dramatically different residue (R414A and R414E) but had no effect on compound inhibition. We only observed interference with compound inhibition when Arg-414 was mutated to bulky residues such as Tyr and Trp. Third, the mutations that do affect compound binding cluster around the binding pocket (instead of scattering all over the Eya2 ED surface) and are located on different secondary structures. It is worth noting that the compounds have a slow off rate and long residence time, as observed by reversibility experiments (Fig. 5), which is highly desirable for a drug and indicates a potential conformational change of the binding site induced by the compound. Our current docking model may not have completely captured the actual binding interaction between the enzyme and the compound after the potential induced fit. We continue to attempt to obtain a co-crystal structure of the compound and Eya2 ED, which would provide the most definitive proof and the most detailed description of the compound-binding mode.
As a member of the HAD family of enzymes, Eya contains a Rossmann-like catalytic core common to all HAD family enzymes (Fig. 11A, orange) and a unique cap structure consisting of a helix-bundle-motif (HBM) only found in Eya (Fig. 11A, cyan) (23). Structural homology searches against other HAD enzymes using the DALI server (41) demonstrate that the catalytic core contains a highly conserved structure (Z = 9.7–12.3), whereas the HBM cap shows no structural homology with other structures (Z < 4.0) (25). Furthermore, a structural homology search of the Eya2 region consisting of the proposed binding site (helix 7, helix 8, and residues forming the base of the pocket (Leu-463, Gly-462, Tyr-461)) displays no homology (Z < 3.5) with any other protein structure. The proposed allosteric compound binding site is located at the interface between the catalytic core and the HBM cap (Fig. 11, A and B). This could provide the selectively against other HAD family of phosphatases, supported by the inactivity of these compounds against Scp1. This proposed binding surface is involved in crystal contacts and is not accessible in the Eya2 ED structures (25), which would explain why we have not been able to obtain compound-containing Eya2 ED crystals in soaking experiments. In the crystal structure of the Six1-Eya2 ED complex (21), the proposed allosteric binding site is solvent-accessible; however, these crystals disintegrate upon soaking with inhibitors, which could be due to conformational changes in Eya2 ED induced by compound binding.
![An external file that holds a picture, illustration, etc.
Object name is zbc0281486820011.jpg An external file that holds a picture, illustration, etc.
Object name is zbc0281486820011.jpg](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4047403/bin/zbc0281486820011.jpg)
Possible mechanisms of allosteric inhibition by MLS000544460. A, Eya2 is a HAD family member consisting of a catalytic core (orange) and a HBM cap (cyan). The allosteric compound binding site is located at the interface between these two motifs (α7-α8 loop) (compound is shown in a green ball-and-stick model). B, surface representation of the catalytic core (orange) and HBM (cyan) showing MLS000544460 at the interface. C, a ribbon diagram demonstrating that the compound (green) binds directly behind the catalytic residues centered around the Mg2+ ion (green sphere). The compound binding at the allosteric site can potentially induce rigidity or conformational changes of the active site pocket.
Binding at the interface between the catalytic core and the cap suggests a possible mechanism of inhibition by these compounds. Some HAD members utilize their cap motif for active site solvent occlusion as part of their catalytic mechanism (23). The cap of HAD members is able to move from an “open” conformation where the active site is exposed to a “closed” conformation in which active site access is restricted after substrate binding (23). Because a portion of the active site of Eya2 (on the opposite side of the allosteric compound-binding site) is flanked by helices 1 and 2 of the HBM cap (Fig. 11C) and most HAD caps are flexible domains that require movement for catalysis (23), it is possible that these allosteric Eya2 inhibitors lock Eya2 into an open conformation in which active site solvent occlusion cannot occur when substrate binds, consequently inhibiting Eya2 activity. A similar mechanism has been identified for LTV-1, an inhibitor that targets the lymphoid tyrosine phosphatase Lyp, a canonical thiol-based, non-HAD phosphatase involved in T cell antigen receptor signaling. LTV-1 selectively prevents the capping of the “WPD” loop, and the inhibitor can only be modeled into an open conformation of Lyp (42).
Another possibility of how these compounds act as allosteric inhibitors is that these compounds may disturb the active site pocket or induce rigidity in the loops that either contain or are close to residues required for catalysis and coordination of the Mg2+ (Asp-274, Asp-276, and Asp-502). Analysis of the crystal structure of Eya2 conjugated to BeF3 demonstrates that the α4-α5 loop, β2-α9 loop, and β3-α10 loop move laterally compared with the native structure, suggesting intrinsic mobility may be required for catalysis (25). In our model, our inhibitor binds into a pocket directly opposite the active site (Fig. 11C). It is possible that the positioning of the compound could either induce rigidity in the catalytic loops and/or interfere with the geometry of the active site thereby inhibiting catalysis.
The chemical series we identified specifically inhibit Eya2 but not Eya3, although we do not yet know whether these compounds inhibit Eya1 and Eya4 because we have not been able to express and purify human Eya1 and Eya4 from E. coli. Multiple amino acids in the allosteric binding pocket of Eya2 whose mutation affected compound inhibition are conserved among all the Eyas (Arg-414, Leu-417, Glu-418, Leu-425) with the exception of Leu-423, Thr-426, and His-427 (Fig. 8). In Eya1 and -4, the Thr-426 residue is conserved but not the Leu-423 and His-427 residues. It is at this point difficult to predict whether this class of compounds will inhibit Eya1 and -4, and further experimental evidence with purified Eya1 and 4 is likely required to address this specificity question.
In summary, we have identified a series of N-arylidenebenzohydrazide compounds that specifically inhibits Eya2 but not Eya3 or other cellular phosphatases. These compounds also inhibit Eya2 phosphatase mediated cell migration, suggesting that these compounds could potentially inhibit metastasis. Our previous data suggest that Eya2 is a critical player in breast and ovarian cancers (16, 17). Eya2 is amplified in 14.8% of ovarian cancers (43), where its overexpression is associated with shortened overall survival in advanced ovarian cancer cases (43). In addition, Eya2 cooperates with Six1 to mediate poor prognosis in breast cancers. (Eya1 also shows this cooperation, but not enough breast tumors overexpress Eya3 or Eya4 to enable meaningful analysis of their cooperation with Six1 (17).) Thus, inhibiting the phosphatase of Eya2 is likely to have a significant therapeutic effect in these cancers. We have demonstrated that the N-arylidenebenzohydrazide compounds likely bind to an allosteric site on the opposite side of the active site. These molecules are promising lead compounds that can potentially be developed into chemical probes for understanding the function of Eya2 phosphatase or therapeutic agents targeting Eya2 overexpressing tumors.
Acknowledgments
The x-ray and NMR facilities were supported in part by the Structural Biology Shared Resource of the University of Colorado Cancer Center (Grant P30CA046934). We thank Drs. David Jones, Geoffrey Armstrong, Natalie Ahn, and the Ahn laboratory for advice and technical help.
*This work was supported by the Molecular Libraries Initiative of the National Institutes of Health Roadmap for Medical Research, National Institutes of Health Grant R03DA030559 (to R. Z. and H. L. F.). This work was also supported by Department of Defense Synergistic Idea Award BC084105 (to H. L. F. and R. Z.), an American Association for Cancer Research and Breast Cancer Research Foundation grant (to H. L. F. and R. Z.), a State of Colorado Proof of Concept grant (to R. Z. and H. L. F.), National Institutes of Health Small Business Technology Transfer Grant R41CA180347 (to H. L. F. and R. Z.), and a Cancer League of Colorado grant (to P. R., H. L. F., and R. Z.).
4The abbreviations used are:
- ED
- Eya domain
- HAD
- haloacid dehalogenase(s)
- HTS
- high throughput screening
- OMFP
- 3-O-methylfluorescein phosphate
- ITC
- isothermal calorimetry
- HBM
- helix-bundle-motif
- DMSO
- dimethyl sulfoxide
- FDP
- fluorescein diphosphate
- HSQC
- heteronuclear single quantum coherence.