Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
MOLECULAR CANCER THERAPEUTICS | SMALL MOLECULE THERAPEUTICS
Development of a Novel Multi-Isoform ALDH Inhibitor Effective as an Antimelanoma Agent Saketh S. Dinavahi1,2,3, Raghavendra Gowda1,2,3,4, Krishne Gowda1, Christopher G. Bazewicz2,3,5, Venkat R. Chirasani1, Madhu Babu Battu6, Arthur Berg7, Nikolay V. Dokholyan1,8, Shantu Amin1, and Gavin P. Robertson1,2,3,4,5,9,10
ABSTRACT ◥ The aldehyde dehydrogenases (ALDH) are a major family of KS100 was mitigated by development of a nanoliposomal formu- detoxifying enzymes that contribute to cancer progression and lation, called NanoKS100. NanoKS100 had a loading efficiency of therapy resistance. ALDH overexpression is associated with a poor approximately 69% and was stable long-term. NanoKS100 was prognosis in many cancer types. The use of multi-ALDH isoform or 5-fold more selective for killing melanoma cells compared with isoform-specific ALDH inhibitors as anticancer agents is currently normal human fibroblasts. NanoKS100 administered intravenously hindered by the lack of viable candidates. Most multi-ALDH at a submaximal dose (3-fold lower) was effective at inhibiting isoform inhibitors lack bioavailability and are nonspecific or toxic, xenografted melanoma tumor growth by approximately 65% with- whereas most isoform-specific inhibitors are not effective as mono- out organ-related toxicity. Mechanistically, inhibition by KS100 therapy due to the overlapping functions of ALDH family members. significantly reduced total cellular ALDH activity to increase reac- The present study details the development of a novel, potent, multi- tive oxygen species generation, lipid peroxidation, and accumula- isoform ALDH inhibitor, called KS100. The rationale for drug tion of toxic aldehydes leading to apoptosis and autophagy. Col- development was that inhibition of multiple ALDH isoforms might lectively, these data suggest the successful preclinical development be more efficacious for cancer compared with isoform-specific of a nontoxic, bioavailable, nanoliposomal formulation containing a inhibition. Enzymatic IC50s of KS100 were 207, 1,410, and 240 novel multi-ALDH isoform inhibitor effective in the treatment of nmol/L toward ALDH1A1, 2, and 3A1, respectively. Toxicity of cancer.
Introduction they often leave behind therapy-resistant cancer cells with a stem cell– like phenotype, which serve as a reservoir for disease recurrence and Malignant melanoma is an aggressive neoplasm accounting for the metastasis (3). majority of skin cancer–related deaths (1). The outlook for metastatic Cancer cells with stem cell characteristics comprise a small subset of disease remains poor, with current 5-year survival rates of 20% (1). undifferentiated cells that initiate tumor formation and generate However, treatment strategies for malignant melanoma have vastly multipotent progenitors (3). They promote tumor aggressiveness, improved with the discovery of targeted therapies to BRAF and MEK repopulation after injury, and metastasis, having intrinsic resistance along with the development of immune checkpoint inhibitors (2). to radiotherapy, chemotherapy, and targeted therapies (4). A major Although current treatment strategies may kill the bulk of tumor cells, mechanism by which these cells develop resistance is through upre- gulation of the aldehyde dehydrogenases (ALDH), which has impaired the response to preoperative chemotherapy and radiotherapy in 1Department of Pharmacology, The Pennsylvania State University College of esophageal carcinoma (5), conventional chemotherapy, erlotinib and Medicine, Hershey, Pennsylvania. 2The Melanoma and Skin Cancer Center, The gefitinib in lung carcinoma (6), olaparib in breast carcinoma (7), and Pennsylvania State University College of Medicine, Hershey, Pennsylvania. 3The cyclophosphamide in a myriad of carcinomas (8, 9). Melanoma Therapeutics Program, The Pennsylvania State University College of The 19 human ALDH isozymes are broadly defined as a superfamily 4 þ Medicine, Hershey, Pennsylvania. Foreman Foundation for Melanoma Research, of NAD(P) -dependent enzymes that participate in aldehyde metab- The Pennsylvania State University College of Medicine, Hershey, Pennsylvania. 5Department of Dermatology, The Pennsylvania State University College of olism, catalyzing the oxidation of toxic aldehydes into carboxylic Medicine, Hershey, Pennsylvania. 6Laboratory of Molecular Cell Biology, Centre acids (10–13). ALDH activity within cells is generally a composite of for DNA Fingerprinting and Diagnostics, Uppal, Hyderabad, India. 7Department the activities of multiple ALDH isoforms, which have overlapping of Public Health Sciences, The Pennsylvania State University College of Medicine, substrate specificity (14, 15). The ALDHs confer a survival advantage 8 Hershey, Pennsylvania. Department of Biochemistry and Molecular Biology, to metabolically active cancer cells, by oxidizing aldehydes that The Pennsylvania State University College of Medicine, Hershey, Pennsylvania. accumulate and cause oxidative damage, into less toxic, more soluble 9Department of Pathology, The Pennsylvania State University College of Med- icine, Hershey, Pennsylvania. 10Department of Surgery, The Pennsylvania State carboxylic acids (16, 17). Accordingly, ALDH overexpression is linked University College of Medicine, Hershey, Pennsylvania. to poorer survival in gastric, breast, lung, pancreatic, and prostate carcinomas, as well as in head and neck squamous cell carcinomas Note: Supplementary data for this article are available at Molecular Cancer Therapeutics Online (http://mct.aacrjournals.org/). (HNSCC; refs. 8, 11, 14, 18, 19). The ALDH1A1, 1A2, 1A3, 3A1, and 3A2 isozymes are particularly important in cancer progression and Corresponding Author: Gavin P. Robertson, The Pennsylvania State University College of Medicine, 500 University Drive, Hershey, PA 17033. Phone: 717-531- resistance to anticancer therapies (8, 15, 20, 21). 8098; Fax: 717-531-0480; E-mail: [email protected] Current ALDH inhibitors can be categorized into multi-ALDH isoform inhibitors and isoform-specific inhibitors, which primarily Mol Cancer Ther 2020;19:447–59 inhibit one isoform (11). Limitations of multi-ALDH isoform inhi- doi: 10.1158/1535-7163.MCT-19-0360 bitors, such as N,N-diethylaminobenzaldehyde (DEAB), which targets 2019 American Association for Cancer Research. ALDH1A1, 1A2, 1A3, 1B1, 2, and 5A1, 4-dimethylamino-4-methyl–
AACRJournals.org | 447
Downloaded from mct.aacrjournals.org on September 24, 2021. © 2020 American Association for Cancer Research. Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
Dinavahi et al.
pent-2-ynthioic acid-S-methylester (DIMATE), which targets purchased from Tocris Biosciences. Isatin and the multi-ALDH ALDH1A1 and 3A1, and citral, which targets ALDH1A1, 1A3, and isoform inhibitor DEAB were purchased from Sigma Aldrich. 2, lack bioavailability or have toxicity (11). DIMATE has tumor- inhibitory efficacy when injected i.p. but will require further preclinical ALDH structure preparation evaluation (22). More recently, the ALDH inhibitors (aldis) -1, -2, -3, The structures of ALDH1A1, 2, and 3A1 bound to the inhibitors -4, and -6 have been developed, which target ALDH1A1, 2, and 3A1, CM037, psoralen, and CB7, respectively (4X4L, 5L13, and 4L2O), were and show efficacy in killing cultured cancer cells, particularly as retrieved from the protein data bank (PDB). The three-dimensional combinatorial therapy (23–25). However, these compounds have structures of the protein complexes were prepared using the protein mainly been tested in vitro and thus require further validation in preparation wizard tool (Schrodinger, LLC, 2017); water molecules preclinical models (23–25). were deleted except those in the inhibitor-binding pocket, bond orders Isoform-specific inhibitors, such as Cpd 3 and CM037 (targeting were assigned, hydrogen atoms were added, and metal ions were ALDH1A1), CVT10216 (targeting ALDH2), and CB7 and CB29 treated as described previously (35–37). Next, the orientation of the (targeting ALDH3A1), have limited efficacy in killing cultured cancer side chain structures of Gln and Asn was flipped, if necessary, to cells, particularly when used as monotherapy, and have not been tested provide the maximum degree of H-bond interactions. The charge state in animal cancer models (18, 26, 27). NCT-501, which targets of His residues was optimized. Finally, a restrained minimization of the ALDH1A1, has been shown to be effective in inhibiting HNSCC protein structure was performed using the OPLS force field with growth in animals via intratumoral injection, suggesting poor systemic backbone atoms being fixed. The minimized protein was used for the bioavailability (28). Other more bioavailable ALDH1A1-specific inhi- docking analysis. The structure was validated using Gaia webserver bitors have been developed, such as the orally bioavailable compounds (http://chiron.dokhlab.org). NCT-505 and NCT-506, and the i.p. available compounds 13 g and 13 h, but have not yet been evaluated in preclinical animal models (29, 30). ALDH grid generation and ligand preparation Therefore, ALDH inhibitors are needed that inhibit the multiple, Prepared protein structures were used to generate scoring grids for functionally overlapping ALDH isoforms, with an acceptable phar- subsequent docking calculations as described previously (35–37). To macologic profile. each protein crystal structure, a grid box of default size (20 20 The current study describes the design and development of a novel, 20A0) was centered on the corresponding active site position. Default potent, multi-isoform ALDH inhibitor, called KS100. KS100 was parameters were used, and no constraints were included during grid developed because ALDH1A1, 2, and 3A1 overexpression was generation. The ligand preparation was then performed using the observed in a cell line melanoma progression model, and targeting ligprep module in Schrodinger (35–37). these individual ALDH isoforms did not affect cultured cell growth. KS100 potently inhibited multiple ALDH isoforms with negligible Molecular docking modeling of ALDH with inhibitors toxicity when administered in a nanoliposomal form, called The starting conformations of ligands were minimized using the NanoKS100. NanoKS100 was bioavailable and inhibited melanoma OPLS 2005 force field until the energy difference between subsequent tumor growth by approximately 65% at submaximal (3-fold lower) structures was 0.001 kJ/mol-A0. The docking study was performed doses. Most importantly, KS100 significantly reduced total cellular using GLIDE 6.6 in Maestro 10.1 (35–37). The GLIDE (Grid Ligand ALDH activity compared with several ALDH inhibitors leading to Docking with Energetics) algorithm estimates a systematic search of enhanced reactive oxygen species (ROS) generation, lipid peroxida- positions, orientations, and conformations of the ligand in the tion, and accumulation of toxic aldehydes causing increased apoptosis enzyme-binding pocket via a series of hierarchical filters. The shape and autophagy. and properties of the receptor are symbolized on a grid by various dissimilar sets of fields that furnish precise scoring of the ligand pose. The potential hit compounds with good fitness score were considered Materials and Methods for docking analysis. All the hits were subjected to the extra precision Cell lines, culture conditions, and chemicals (XP) mode of GLIDE. Default values were accepted for Van der Waals Normal human fibroblasts (FF2441) were provided by Dr. Craig scaling, and input partial charges were used. During the docking Myers, Penn State College of Medicine (Hershey, PA). The human process, the GLIDE score was used to select the best conformation melanoma cell lines WM35, WM115, WM278, WM3211, 1205 Lu, for each ligand (35–37). and A375M and normal melanocytes (NHEM) were provided by Dr. Meenhard Herlyn Wistar Institute (Philadelphia, PA). The human Modeling to assess specificity of ALDH inhibitor binding melanoma cell line UACC 903 was provided by Dr. Mark Nelson, All bound crystal water molecules and ligands were stripped out of University of Arizona (Tucson, AZ). The wild-type BRAF melanoma the crystal structures of ALDH1A1, 2, and 3A1 prior to docking. cell line C8161.Cl9 was provided by Dr. Danny Welch, University of Simultaneously, the structure of KS100 was built and optimized in Kansas (Kansas City, KS), and MelJuSo was provided by Dr. Judith Marvin sketch workspace. Because ALDH1A1, 2, and 3A1 are depos- Johnson, Institute for Immunology, Germany. Cell lines were main- ited in oligomeric states in the PDB database, monomeric conforma- fi tained in a 37 C humidi ed 5% CO2 atmosphere incubator and tions of respective structures were extracted, and missing atoms or periodically monitored for phenotypic and genotypic characteristics residues were relocated through homology modeling. The structures and tumorigenic potential to validate and confirm cell line identity. were optimized using DMD software suite, and subsequently, molec- Cell lines were authenticated and tested for mycoplasma contamina- ular docking using Medusadock suite (http://medusadock.dokhlab. tions periodically and are used within 15 passages after authentication. org/) was employed, which is known for its rapid sampling efficiency The ALDH1A1- and 3A1-specific inhibitors, Cpd 3 (31, 32) and and high prediction accuracy (38). Initially, molecular docking of CB7 (18, 33) respectively, were synthesized in-house according to KS100 to the active site of ALDH1A1 alone was attempted as the previously published procedures. The ALDH1A1-specific inhibitor, conformations of ALDH1A1, 2, and 3A1 are structurally identical CM037 (34), and ALDH2-specific inhibitor, CVT10216 (34), were (Supplementary Fig. S1). From the ALDH1A1–KS100 docked
448 Mol Cancer Ther; 19(2) February 2020 MOLECULAR CANCER THERAPEUTICS
Downloaded from mct.aacrjournals.org on September 24, 2021. © 2020 American Association for Cancer Research. Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
ALDH Inhibition for Melanoma
complex, it was evident that the KS100 binding pocket in ALDH1A1 and release the drug into the solution. The precipitated lipids were was lined by the residues: Ser-121, Phe-171, Val-174, Met-175, separated via centrifugation at 10,000 rpm for 15 minutes. The Trp-178, Glu-269, Phe-290, His-293, Gly-294, Tyr-297, Cys-302, supernatant was then used to measure KS100 concentration, Cys-303, Ile-304, Tyr-457, Gly-458, Val-460, and Phe-466. calculated from a standard curve of KS100 from 0.01 to 1 mg/mL. To identify the off-target effects of KS100, the binding scaffold of A 1:1 solution of chloroform to methanol was used as the reference KS100 as a substructure was extracted and employed in Erebus (http:// blank. The percentage of KS100 incorporated into nanoliposomes erebus.dokhlab.org), a webserver that searches the entire PDB database was calculated as (incorporated KS100/total KS100) X 100 (40, 41). for a given substructural scaffold (39). Erebus identifies off-target structures from the PDB database by matching substructures with the (b) Stability. Stability of NanoKS100 stored at 4 C was assessed weekly same amino acids and atoms segregated by identical distances (within by comparing size and zeta potential using the Malvern Zetasizer fi a given tolerance) as the atoms of the query structure (39). Finally, the and measuring IC50 ef cacy for killing UACC 903 melanoma cells prediction accuracy of a match was evaluated by the root-mean-square by MTS assay and comparing these values with that of freshly deviation (RMSD) or by the normal weight with a given variance. manufactured NanoKS100 (40, 41).
Synthesis of KS100 (c) In vitro drug release kinetics. Drug release kinetics of KS100 from KS100 was synthesized as shown in Scheme 1 of Fig. 2B. Briefly, 5,7- the liposomes were measured using 1 mL of purified NanoKS100 dibromoisatin (1) (10 mmol) was dissolved in anhydrous DMF by dialysis in 1 L of 10 mmol/L reduced glutathione at room (30 mL) and cooled on ice with stirring. Solid K2CO3 (11 mmol) temperature through a molecular weight cutoff of 25 kDa mem- was added, and the dark-colored suspension was brought to room brane (Spectra Por). NanoKS100 (0.05 mL) in the dialysis bag was temperature and stirred for 1 hour. 1,4-Bis(bromomethyl)benzene removed at 0.5, 1, 2, 4, 8, 12, 24, 36, 48, and 72 hours, and the (40 mmol) was added slowly with constant stirring until the starting amount of KS100 released at each time point was estimated using material had been consumed (monitored by TLC). The reaction UV-visible spectrophotometry as detailed previously (40, 41). mixture was poured into cold water and extracted with ethyl acetate. The ethyl acetate layer was washed with water, brine, and dried (d) Hemolytic activity. The hemolytic activity assay was performed as fl over MgSO4. The solvent was removed, and the crude product was described previously (41, 42). Brie y, fresh mouse and rat blood purified by silica gel column chromatography using 80:20 hexanes/ were drawn and placed into an EDTA test tube. Erythrocytes were ethyl acetate as the eluent to yield the intermediate 5,7-dibromo-1-(4- separated from plasma by centrifugation at 1,500 rpm for 10 bromomethylbenzyl)-1H-indole-2,3-dione (2) as orange-red crystals. minutes at 4 C using PBS. Erythrocyte pellets were diluted with To the intermediate compound (1.02 mmol), thiourea (1.02 mmol) 50 mL PBS in centrifuge tubes to give a 5% v/v solution and and ethanol (25 mL) were added and heated to reflux until the then treated with 5 mmol/L KS100 in DMSO, NanoKS100 (10– starting material had been consumed (monitored by TLC). The 40 mmol/L) in PBS, empty liposome, or 1% Triton X-100 (positive solvent was removed under vacuum. The final product (2-[4-(5,7- control). Samples were incubated at 37 C for 60 minutes and then dibromo-2,3-dioxo-2,3-dihydro-indol-1-ylmethyl)benzyl]isothiourea) centrifuged at 12,000 rpm for 10 minutes. Next, supernatants were (3) was recrystallized in ethanol-ethyl acetate to afford KS100 (yield transferred to a 96-well plate and absorption measured at 540 nm. 70%). The identity of KS100 was confirmed by nuclear magnetic The amount of hemoglobin released in the presence of 1% Triton resonance as well as mass spectra analysis, and purity (>99%) was X-100 was set as 100% lysis, and % hemolysis was calculated as: quantified by high-performance liquid chromatography analysis. (absorbance of the samples at 540 nm/absorbance of the positive control) X 100. Preparation of NanoKS100 KS100 was encapsulated into a nanoliposome by first combining siRNA transfections L-a-Phosphatidylcholine and 1,2-Dipalmitoyl-sn-Glycero-3- Duplex stealth siRNA sequences for scrambled and ALDH1A1, Phosphoethanolamine-N-[Methoxy(Polyethylene glycol)-2000] 1A2, 1A3, 1L1, 2, 3A1, 5A1, 18A1 and BRAF were obtained from ammonium salt in chloroform at 80:20 mol % for a final lipid Invitrogen. Individual siRNAs were introduced into UACC 903 cells concentration of 25 mg/mL (Avanti Polar Lipids; refs. 40, 41). Note via nucleofection using an Amaxa nucleofector with solution that 5 mg of KS100 (in methanol) was then added to the lipid mixture, R/program K-17. Nucleofection efficiency was >90% with 80% to dried under nitrogen gas, and resuspended in saline at 60 C. The 90% cell viability. siRNA knockdown was confirmed either by mixture was then sonicated at 60 C for 30 minutes followed by qRT-PCR or by Western blots. extrusion through a 100-nm polycarbonate membrane using Avanti Mini-Extruder (Avanti Polar Lipids Inc.). (40, 41). qRT-PCR analysis Total RNA was extracted by Trizol (Sigma), and cDNA was Characterization of nanoparticle-based KS100 called generated by reverse transcription kit (Applied Biosystems). qRT-PCR NanoKS100 was performed using a SYBR green kit (Qiagen). Expression of each isoform in melanoma cell lines was normalized to corresponding (a) Drug encapsulation. Efficiency of encapsulation of KS100 in the expression in fibroblasts. Primers for the ALDH isoforms were used nanoliposomal formulation was estimated by UV-visible spectro- as described previously (15). photometry (SPECTRAmax M2 plate reader; Molecular devices; refs. 40, 41). Specifically, 1 mL of NanoKS100 solution was added ALDH isoform–specific enzyme assays to a 10 kDa Centricon filter tube (Millipore) and centrifuged at ALDH1A1, 2, and 3A1 enzyme assays were performed as described 3,750 rpm for 30 minutes to remove free KS100. Next, 0.5 mL of by the manufacturer (R & D systems). Isoform-specific aldehydes were purified NanoKS100 was combined with 0.5 mL of a 1:1 solution of converted to their respective carboxylic acids along with the conver- þ chloroform to methanol to destroy the nanoliposomal structure sion of NAD to NADH (absorbance at 340 nm). Specifically, 1 mg/mL
AACRJournals.org Mol Cancer Ther; 19(2) February 2020 449
Downloaded from mct.aacrjournals.org on September 24, 2021. © 2020 American Association for Cancer Research. Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
Dinavahi et al.
of ALDH1A1 was treated with various concentrations of ALDH before measuring fluorescence at 485 nm excitation and 520 nm inhibitors for 15 minutes followed by addition of substrate mixture emission. (10 mmol/L propionaldehyde; 100 mmol/L KCl; 1 mmol/L NAD; 2 mmol/L DTT; 50 mmol/L Tris, pH 8.5), and the absorbance Lipid peroxidation of NADH was measured in kinetic mode for 5 minutes. Similarly, Lipid peroxidation was measured using the thiobarbituric acid 0.5 mg/mL of ALDH2 was used with 2 mmol/L of acetaldehyde as the reactive substances kit according to the manufacturer's instructions substrate, and 0.2 mg/mL of ALDH3A1 was used with 1 mmol/L of (Cayman Chemicals). Briefly, cells were treated with 5 mmol/L of 4-nitrobenzaldehyde as the substrate following the addition of ALDH ALDH inhibitor or DMSO for 24 hours. Cell pellets were lysed in PBS inhibitors. by sonication on ice. Lipids in the lysates were hydrolyzed in the presence of acetic acid and sodium hydroxide. Free malondialdehyde Total cellular ALDH activity assay (MDA) was measured by the reaction to thiobarbituric acid colori- Briefly, total cellular ALDH activity assays were performed on metrically at 530 nm. cell lysates. Note that 20 mg of cell lysate was treated with 1 mmol/L of ALDH inhibitors or DMSO for 15 minutes followed by Apoptosis assay addition of substrate mixture (2 mmol/L acetaldehyde; 100 mmol/L The Annexin-V-PE/7-AAD kit (Millipore) was used to distinguish KCl; 1 mmol/L NAD; 2 mmol/L DTT; 50 mmol/L Tris, pH 8.5), live cells from apoptotic cells as described previously (43). Briefly, cells and the absorbance of NADH was measured in kinetic mode for were incubated with 5 mmol/L of ALDH inhibitor or DMSO for 5 minutes. 24 hours. Cells were washed with PBS and stained with Annexin- V–PE and 7-AAD solution per the manufacturer's instructions. Cells Cell viability assays were acquired by BD Fortessa flow cytometer and gated for four Cell viability assays of UACC 903 cells transfected with siRNA, and distinct regions, namely, live cells (Annexin V-7 AAD ), early apo- þ þ þ melanoma cells (UACC 903, 1205 Lu, C8161.CI9, MelJuSo), normal ptotic (Annexin V-7 AAD ), late apoptotic (Annexin V-7 AAD ), þ human fibroblasts (FF2441), and melanocytes (NHEM) treated with and necrotic (Annexin V-7 AAD ) cells. ALDH inhibitors were performed as described previously (43–45). Briefly, 5,000 cells per well were plated in a 96-well plate and incubated Western blot analysis overnight at 37 C in a 5% CO2 atmosphere. For the siRNA knockdown Cells were harvested by the addition of RIPA lysis buffer, and experiment, cells were incubated for another 72 hours. For the ALDH samples were processed as previously described (46, 47). Briefly, 1 to 6 inhibitor experiments, cells were treated with agents at various con- 2 10 cells were incubated overnight at 37 C in a 5% CO2 atmo- centrations and incubated for 72 hours. Twenty microliter of MTS sphere. For experiments with KS100, the agent was added and protein reagent was then added into each well, and formation of tetrazolium lysates collected following 24 hours of treatment. Blots were probed was measured by absorbance after 1 hour at 492 nm. IC50 values or % with antibodies according to each supplier's recommendations: anti- cells for each experimental group were measured in three independent bodies to cleaved PARP and LC3B from Cell Signaling Technology; experiments. alpha-enolase, ALDH1A1, 2, 3A1, 18A1, BRAF, and secondary anti- bodies conjugated with horseradish peroxidase from Santa Cruz Toxicity and MTD animal studies Biotechnology. Immunoblots were developed using the enhanced All the animal experiments were conducted according to the guide- chemiluminescence detection system (Thermo Fisher Scientific). lines of Penn State, Hershey Institutional Animal Care and Use Alpha-enolase served as the loading control. Committee. To determine the effective dose for in vivo efficacy studies, KS100 and NanoKS100 were injected i.p. and i.v., respectively, into Tumor efficacy and toxicity assessment Swiss Webster mice (Jackson labs) once daily for 7 days (40, 41). Efficacy and toxicity studies were performed in nude mice as Animals were monitored for changes in body weight, behavior, and described previously (40, 46, 48, 49). Briefly, 1 million cells were physical distress compared with control (DMSO for KS100, empty injected in both flanks of 4- to 6-week-old female nude Balb/c mice liposomes for NanoKS100). Dose escalation was performed to identify (Envigo). After a week, when the tumors were vascularized, animals the MTD. were treated with either NanoKS100 (at various doses) or empty liposomes. Tumor volumes, animal weight and behavior were mon- AldeRed ALDH detection assay itored continuously every other day. Animals were sacrificed after The AldeRed ALDH detection assay (Millipore) was used to tumor volumes in the empty liposome groups exceeded 2,500 mm3 and þ distinguish ALDH cells from the ALDH cells. Briefly, cells were tumors were subsequently collected. incubated with 5 mmol/L of ALDH inhibitor or DMSO for 24 hours. Cells were washed with PBS and stained with AldeRed reagent Assessment of serum biomarkers of major organ toxicity (AldeRed 588-A) for 1 hour as per the manufacturer's instructions. At the end of the UACC 903 xenograft study, serum samples were Cells were acquired by BD Fortessa flow cytometer and gated for analyzed for levels of alanine aminotransferase, alkaline phosphatase, þ ALDH cells using DEAB as a negative control. albumin, globulin, total protein, total bilirubin, blood urea nitrogen, glucose, creatinine, amylase, and calcium (40, 46, 48, 49). Serum ROS assay analysis was performed at the Centralized Biological Laboratory, Penn To quantify ROS levels, the nonfluorescent dye DCFDA (Sigma) State, University Park. Empty liposomes served as the control. was used (43). DCFDA turns to highly fluorescent 20,70-dichlorofluor- escein upon oxidation by ROS generated in cells (43). Briefly, cells were Statistical analysis treated with 5 mmol/L of ALDH inhibitor or DMSO for 24 hours. Statistical analysis was undertaken using the one-way/two-way DCFDA (10 mmol/L) was then added and incubated for 30 minutes ANOVA GraphPad PRISM Version 7.04 software. Dunnett's post hoc
450 Mol Cancer Ther; 19(2) February 2020 MOLECULAR CANCER THERAPEUTICS
Downloaded from mct.aacrjournals.org on September 24, 2021. © 2020 American Association for Cancer Research. Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
ALDH Inhibition for Melanoma
analysis was performed when there was a significant difference. Results Subsequent analysis of The Cancer Genome Atlas (TCGA) database were considered significant at a P value of < 0.05. to determine the relationship of ALDH overexpression on survival of patients with melanoma yielded variable results. Specifically, over- expression of ALDH1A1 and 2 was associated with slightly improved Results survival (Fig. 1B), whereas high ALDH3A1 expression was associated ALDH overexpression occurs in melanoma and is associated with lower survival (Fig. 1C). However, these data are measurements with disease progression of RNA expression and thus do not take into account overall protein Cancer cell expression of ALDHs often increases with disease levels or ALDH activity in cancer cells. progression, as oxidative stress secondary to high metabolic demands To functionally determine whether targeting ALDH1A1, 2, or 3A1 leads to ROS generation, lipid peroxidation, and the accumulation of in melanoma affects cell proliferation, a rapid siRNA screen was toxic aldehydes, which can inhibit cancer cells (17, 50). Western blot undertaken (Fig. 1D). siRNA for ALDH18A1, a unique ALDH analysis of ALDH1A1, 2, and 3A1 in melanoma cells revealed that isoform that promotes melanoma cell survival through proline syn- ALDH overexpression occurs in melanoma compared with control thesis (51), and V600EBRAF were used as positive controls. Knockdown fibroblast (FF2441) and melanocyte (NHEM) cells (Fig. 1A). Further, of each respective protein by its siRNA is shown in Fig. 1E. Similar to the degree of ALDH expression correlated with melanoma stage the scrambled siRNA, individual siRNA knockdown of ALDH1A1, 2, such that metastatic melanomas exhibited the highest ALDH expres- and 3A1 did not affect UACC 903 cell survival up to 72 hours where sion levels, followed by vertical growth phase and finally radial as the positive control siRNA caused an approximately 50% reduction growth phase melanomas. ALDH expression was not dependent in cell survival (Fig. 1D). Pharmacologic inhibition of ALDH1A1, 2, on BRAF mutational status, as ALDH levels were similar between and 3A1 by isoform-specific inhibitors also had no effect on cell mutant V600EBRAF and wild-type BRAF cells. qRT-PCR analysis of survival, even when 100 mmol/L concentrations were used for 72 hours the ALDH isoforms demonstrated that the isoforms overexpressed in (Fig. 1F). In contrast, DEAB, a multi-ALDH isoform inhibitor, UACC 903 and 1205 Lu melanoma cells compared with control reduced UACC 903 cell survival by 30% at a 100 mmol/L concentration fibroblasts are ALDH1A1, 1A2, 1A3, 1L1, 2, 3A1, 5A1, and 18A1 after 72 hours. These results suggested that targeting multiple ALDH (Supplementary Table S1). Expectedly, when cells were stained with isoforms with overlapping function may be more effective for mela- þ AldeRed dye to isolate the ALDH (cells with high levels of ALDH) noma therapy specifically and anticancer therapy in general. In þ cells and ALDH cells using flow sorting, the ALDH cells had a addition, that inhibiting individual ALDH isoforms may lead to high expression of these isoforms compared with ALDH cells resistance mediated through upregulation of functionally similar (Supplementary Table S1). ALDHs.
Figure 1. The ALDH family is collectively important in melanoma. Western blot showing ALDH1A1, 2, and 3A1 expression levels in normal human fibroblasts (FF2441), melanocytes (NHEM), radial growth phase (RGP), vertical growth phase (VGP), and metastatic melanoma cell lines. ALDH expression in general increased during disease progression and was not dependent on BRAF mutational status. Alpha-enolase served as the loading control (A). Data from the TCGA database showing slightly better survival with ALDH1A1 and 2 overexpression (B) and worse survival with ALDH3A1 overexpression (C) in patients with melanoma. The data are available through the UCSC Xena Cancer Browser. Individual siRNA knockdown of ALDH1A1, 2, and 3A1 did not significantly reduce the survival of UACC 903 cells after 72 hours in an MTS assay. siRNA to BRAF and ALDH18A1 served as positive controls. Scrambled siRNA served as the negative control (D). siRNA knockdown of ALDH1A1, 2, 3A1, 18A1, and BRAF in UACC 903 cells was confirmed via Western blot. Alpha-enolase served as loading control (E). Pharmacologic inhibition of ALDH1A1, 2, and 3A1 using ALDH isoform–specific inhibitors (Cpd 3, CVT10216, and CB7, respectively), and the multi-ALDH isoform inhibitor, DEAB, revealed multi-ALDH isoform inhibition was most effective in inhibiting UACC 903 cell survival (F).
AACRJournals.org Mol Cancer Ther; 19(2) February 2020 451
Downloaded from mct.aacrjournals.org on September 24, 2021. © 2020 American Association for Cancer Research. Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
Dinavahi et al.
Identification and development of the novel, potent, would bind and interact more effectively in the ligand-binding pocket multi-isoform ALDH inhibitor, called KS100 of the ALDHs, using the backbones of Isatin and Cpd 3. A series of To create a multi-ALDH isoform inhibitor, an in silico screen was compounds were tested through in silico modeling to determine initially undertaken based on the x-ray crystal structure of ALDH1A1 whether they had optimal docking in the ligand-binding pocket of using various natural products. Isatin was identified during this screen ALDH1A1, and KS100 was selected as the best candidate (Fig. 2A). It as weakly binding to ALDH1A1 compared with the ALDH1A1- was also found to fit well into the ligand-binding pockets of ALDH2 specific inhibitors Cpd 3 and CM037 (Fig. 2A). A medicinal chemistry and 3A1. KS100 had docking scores of 10.247, 8.716, and 13.851 approach was subsequently undertaken to design compounds that for ALDH1A1, 2, and 3A1, respectively (Table 1), compared with
Figure 2. Design, synthesis, and toxicity analysis of the novel, ALDH1A1, 2, and 3A1 inhibitor, called KS100. Based on the structure and binding of Isatin, Cpd 3, and CM037, a medicinal chemistry approach was undertaken to design KS100, which exhibited more effective binding to ALDH1A1, 2, and 3A1 (A). KS100 was synthesized from 5,7-dibromoisatin followed by benzylation as detailed in the Supplementary Materials and Methods (B). A 7-day repeated dose study was conducted for KS100. KS100 was administered i.p. daily, whereas animal body weight, physical and behavioral changes, and mortality were monitored. KS100 was toxic starting at 5 mg/kg (C).
452 Mol Cancer Ther; 19(2) February 2020 MOLECULAR CANCER THERAPEUTICS
Downloaded from mct.aacrjournals.org on September 24, 2021. © 2020 American Association for Cancer Research. Published OnlineFirst November 21, 2019; DOI: 10.1158/1535-7163.MCT-19-0360
ALDH Inhibition for Melanoma
Table 1. Docking scores and ALDH-inhibitory activity for KS100 and other ALDH inhibitors.
Docking scores IC50s (nmol/L) Compound ALDH1A1 ALDH2 ALDH3A1 ALDH1A1 ALDH2 ALDH3A1
Isatin 5.46 6.398 5.819 15,635 1,821 168,661 28,679 5,047 304 Cpd 3 7.686 9.839 7.695 44 12 72,136 1,640 11,866 548 CM037 11.276 7.137 8.137 98 34 2,278 250 1,774 303 CVT10216 7.892 11.809 8.924 2,427 194 53 2 2,719 608 CB7 8.159 7.846 14.576 139,016 16,934 144,409 11,470 298 29 DEAB 9.154 10.026 11.211 89 23 833 277 15,119 4,091 KS100 10.247 8.716 13.851 207 10 1,410 248 240 50
Note: Isatin, Cpd 3, CM037, CVT10216, CB7, DEAB, and KS100 were docked into the active site pockets of ALDH1A1, 2, and 3A1. Docking scores were calculated using
GLIDE 6.6 in Maestro 10.1. Enzyme inhibition studies for ALDH1A1, 2, and 3A1 were performed as described in the Materials and Methods. KS100 displayed low IC50 values for all three ALDH isoforms tested. Preexisting isoform-specific ALDH inhibitors (Cpd 3 and CM037 for ALDH1A1, CVT10216 for ALDH2, CB7 for ALDH3A1) and the multi-ALDH isoform inhibitor, DEAB, were used as positive controls.