Ligand-receptor Interaction between Triterpenoids and the 11β-Hydroxysteroid dehydrogenase type 2 (11βHSD2) Enzyme Predicts Their Toxic Effects against Tumorigenic r/m HM-SFME-1 Cells*

The present study deals with in silico prediction and in vitro evaluation of the selective cytotoxic effects of triterpenoids on tumorigenic human c-Ha-ras and mouse c-myc cotransfected highly metastatic serum-free mouse embryo-1 (r/m HM-SFME-1) cells. Ligand fitting of five different triterpenoids to 11β-hydroxysteroid dehydrogenase type 2 (11βHSD2) was analyzed with a molecular modeling method, and glycyrrhetinic acid (GA) was the best-fitted triterpenoid to the ligand binding site in 11βHSD2. Analysis of antiproliferative effects revealed that GA, oleanolic acid, and ursolic acid had selective toxicity against the tumor cells and that GA was the most potent triterpenoid in its selectivity. The toxic activity of the tested triterpenoids against the tumor cells showed good correlations with the partition coefficient (logP) and polar surface area values. Time-lapse microscopy, fluorescence staining, and confocal laser scanning microscopic observation revealed that GA induced morphologic changes typical of apoptosis such as cell shrinkage and blebbing and also disrupted the cytoskeletal proteins. Furthermore, GA exhibited a strong inhibitory effect on 11βHSD2 activity in the tumor cells. Our current results suggest that analysis of the ligand-receptor interaction between triterpenoids and 11βHSD2 can be utilized to predict their antitumor effects and that GA can be used as a possible chemopreventive and therapeutic antitumor agent. To the best of our knowledge, this is the first report on in silico prediction of the toxic effects of triterpenoids on tumor cells by 11βHSD2 inhibition.

The present study deals with in silico prediction and in vitro evaluation of the selective cytotoxic effects of triterpenoids on tumorigenic human c-Ha-ras and mouse c-myc cotransfected highly metastatic serum-free mouse embryo-1 (r/m HM-SFME-1) cells. Ligand fitting of five different triterpenoids to 11␤-hydroxysteroid dehydrogenase type 2 (11␤HSD2) was analyzed with a molecular modeling method, and glycyrrhetinic acid (GA) was the best-fitted triterpenoid to the ligand binding site in 11␤HSD2. Analysis of antiproliferative effects revealed that GA, oleanolic acid, and ursolic acid had selective toxicity against the tumor cells and that GA was the most potent triterpenoid in its selectivity. The toxic activity of the tested triterpenoids against the tumor cells showed good correlations with the partition coefficient (logP) and polar surface area values. Time-lapse microscopy, fluorescence staining, and confocal laser scanning microscopic observation revealed that GA induced morphologic changes typical of apoptosis such as cell shrinkage and blebbing and also disrupted the cytoskeletal proteins. Furthermore, GA exhibited a strong inhibitory effect on 11␤HSD2 activity in the tumor cells. Our current results suggest that analysis of the ligand-receptor interaction between triterpenoids and 11␤HSD2 can be utilized to predict their antitumor effects and that GA can be used as a possible chemopreventive and therapeutic antitumor agent. To the best of our knowledge, this is the first report on in silico prediction of the toxic effects of triterpenoids on tumor cells by 11␤HSD2 inhibition.
Triterpenoids, which are biosynthesized in plants by cyclization of squalene, are widely distributed throughout the vegeta-ble kingdom, utilized in many food products and the major components of medicinal plants used in Asian countries (1). There is a growing interest in elucidating the biological and pharmacological roles of triterpenoids in analgesic, anti-inflammatory, anti-tumor, hepatoprotective, and immunomodulatory effects (1), and we have been focusing our attention on certain triterpenoids as multifunctional agents for the prevention and treatment of cancer. Recently, we found that ursolic acid (UA) 2 from apples was selectively toxic to tumorigenic human c-Ha-ras and mouse c-myc cotransfected highly metastatic serum-free mouse embryo-1 (r/m HM-SFME-1) cells (2). We further found that glycyrrhetinic acid (GA), a licorice component, was not only selectively toxic to the tumor cells but also more potent than some clinically available antitumor agents in its selectivity (3). SFME cells, which were established by Loo et al. (4), were originally derived from a 16-day-old whole Balb/c mouse embryo and are maintained in a serum-free culture medium. These cells do not undergo growth crisis, maintain their diploid karyotype for extended passages, and are non-tumorigenic in vivo. Consequently, they are non-transformed, behave as primary cultures, have a finite lifespan, and display the characteristics of the CNS progenitor cells (5,6). Another SFME-derived cell line is the r/m HM-SFME-1 cell line (7). Although SFME cells are non-tumorigenic in vivo and require EGF for their survival, growth, and proliferation (5,6), r/m HM-SFME-1 cells are tumorigenic and do not require any growth factors such as EGF (8). Analyzing the characteristics and behaviors of normal and tumorigenic SFME cells could be of great importance in the field of toxicological studies for cancer prevention and therapy because they are of the same lineage and simple comparisons of these cells may contribute to our understanding of the behavioral differences between normal and tumor cells in the CNS in their responses to anti-tumor agents.
11␤-Hydroxysteroid dehydrogenase type 2 (11␤HSD2) requires NAD ϩ , shows dehydrogenase activity for endogenous glucocorticoids such as corticosterone (CORT; 9 -11), and has been reported to be associated most notably with pituitary adenomas in the CNS (12), but also with colonic adenomas (13) and breast (14,15) and colorectal (16) cancers. The underlying explanation for the aberrant 11␤HSD2 expression is uncertain, but it has been postulated to control glucocorticoid regulation of cellular proliferation (17). Results from in vitro studies using malignant transformed cell lines demonstrated the anti-proliferative actions of glucocorticoids. Therefore, the local inactivation of glucocorticoids such as CORT by 11␤HSD2 could be an important oncogenic process promoting cellular proliferation (18). Furthermore, 11␤HSD2 inhibition by GA prevented colon cancer without triggering adverse side effects in the cardiovascular system (13). GA also had adverse effects on the proliferation of pituitary adenomas (19), and 11␤HSD2 inhibition induced apoptosis of corticotroph tumor cells (12). Taken together, these reports suggest that 11␤HSD2 inhibition can be utilized as a potential therapeutic option in controlling cancer. Thus, prediction of the toxic effects of drugs on tumors can be achieved by investigating 11␤HSD2 inhibition in silico and also in vitro in tumorigenic r/m HM-SFME-1 cells.
In the present study, in light of our previously reported selective toxicity of GA and UA against the tumor cells, in silico prediction and in vitro evaluation of the anti-tumor effects of triterpenoids were investigated. Ligand fitting of five triterpenoids to 11␤HSD2 was analyzed with a molecular modeling method to predict their cytotoxic effects. Normal SFME and tumorigenic r/m HM-SFME-1 cells were treated with the triterpenoids to investigate their efficacy as tumor cell-selective toxic agents. Subsequently, the half-maximal inhibitory concentration (IC 50 ) values of the tested triterpenoids for the normal and tumor cells were analyzed with the partition coefficient between 1-octanol and aqueous phases (logP) and the polar surface area (PSA) values to examine whether there were correlations between them. Furthermore, the tumor cells were exposed to GA, which was the most potent triterpenoid in its selective toxicity, and time-lapse microscopy, fluorescence staining, and confocal laser scanning microscopic observation were adopted to analyze the morphologic and cytoskeletal changes. Furthermore, the effects of GA on CORT levels were analyzed to assess the inhibition of 11␤HSD2 enzyme activity by GA.
In Silico Ligand-Receptor Interaction between Triterpenoids and 11␤HSD2-The binding site selection and exploration for 11␤HSD2 were carried out as reported previously (20). In brief, 11␤HSD1 (Protein Data Bank code 3HFG) was selected as a template for the structure modeling of 11␤HSD2 (NCBI reference sequence NM_008289.2) because of its good crystal structure resolution (2.3 Å) and because its information was the latest (from 2009) among the reported 11␤HSD1 models. For the construction of the 11␤HSD2 model, 100 independent models of the target protein were built using a Boltzmann-weighted randomized modeling procedure in the Molecular Operating Environment 2009.10 (Chemical Computing Group, Inc., Montreal, Canada), which was adapted from reports by Levitt (21) and Fechteler et al. (22). The intermediate models were evaluated by a residue packing quality function, which is sensitive to the degrees to which non-polar side chain groups are buried, and hydrogen bonding opportunities are satisfied. The 11␤HSD2 model with the best packing quality function and full energy minimization was selected for further analyses. The secondary structures of the 11␤HSD2 model exhibited a central six-stranded all-parallel ␤-sheet sandwich-like structure, flanked on both sides by three helices, which are in agreement with the 11␤HSD1 model. Hydrophobic or hydrophilic ␣spheres, which were created by the Site Finder module of Molecular Operating Environment 2009.10, were utilized to define potential ligand binding sites (LBS).
The analysis of the ligand-receptor interaction between the ligands (CORT, PubChem CID 5753; AA, PubChem CID 73170; BA, PubChem CID 5318287; GA, PubChem CID 10114; OA, PubChem CID 10494; and UA, PubChem CID 64945) and the 11␤HSD2 model was performed with the ␣-sphere and excluded volume-based ligand-protein docking (ASE-Dock) module of Molecular Operating Environment 2009.10 (23). In the ASE-Dock module, ligand atoms have ␣-spheres within 1 Å. Based on this property, concave models were created, and ligand atoms from a large number of conformations were generated by superimposition with these points can be evaluated and scored by the maximum overlap with the ␣-spheres and minimum overlap with the receptor atoms. The ligand fitting scores and ASE scores are then obtained (24). The scoring function used by the ASE-Dock module is based on ligand-receptor interaction energies and the score is expressed as a U total value. The ligand conformations were subjected to energy minimization using the MMF94S force field (25), and 500 conformations were generated using the default systematic search parameters. Five thousand poses per conformation were randomly placed onto the ␣-spheres located within the LBS in 11␤HSD2. From the resulting 500,000 poses, the 200 poses with the lowest U total values were selected for further optimization with the MMF94S force field. During the refinement step, the ligands were free to move within the binding pocket.
Cell Lines and Cell Culture-SFME cells were a gift from Dr. S. Shirahata (Kyushu University, Fukuoka, Japan) and r/m HM-SFME-1 cells were taken from our cell stocks (7). The basal nutrient medium was a 1:1 mixture of Dulbecco's Modified Eagle's medium and nutrient mixture F-12 Ham (Dulbecco's Modified Eagle's medium/F-12) (26,27) and supplemented with sodium bicarbonate, sodium selenite, and gentamicin sulfate. Cells were maintained in Dulbecco's Modified Eagle's medium/F-12 supplemented with insulin, transferrin and EGF in 60-mm diameter dishes precoated with bovine fibronectin (Biomedical Technologies, Cambridge, MA) in a humidified atmosphere containing 20% O 2 and 5% CO 2 at 37°C.
Measurement of Antiproliferative Activity-Cells plated at 1 ϫ 10 4 cells/well in 96-well microplates were treated with the test compounds at half confluency. After culture for another 24 h, the cell numbers were determined by the MTT assay (28).
Calculation of logP and PSA-The Spartan 06 program (Wavefunction, Inc., Irvine, CA) was adopted for calculation of logP from the Crippen model and PSA.
Western Blotting Analysis-Proteins were extracted with PBS containing 1 mM PMSF, 1 mM EDTA, 2 mM 2-mercaptoethanol, and 1% Triton X-100 at 4°C for 3.5 h. For Western blotting analysis, aliquots of proteins were separated by SDS-PAGE, blotted onto a nitrocellulose membrane, and probed with a primary antibody followed by a secondary antibody. The primary antibodies used were mouse monoclonal anti-GAPDH (Santa Cruz Biotechnology, Santa Cruz, CA) and rabbit polyclonal anti-11␤HSD2 (Cayman Chemical, Ann Arbor, MI). The secondary antibodies used were alkaline phosphatase-conjugated anti-mouse and anti-rabbit IgG 1 (Chemicon International, Temecula, CA). Visualization of the antigen-antibody complexes was performed with 33 l of 5-bromo-4-chloro-3-indolyl phosphate, 66 l of nitroblue tetrazolium, and 40 l of 1 M MgCl 2 in 10 ml of 0.1 M Tris-HCl buffer (pH 9.5). Images of the positive bands were obtained by scanning, and the densities were determined using an LAS-3000 image analyzer (FujiFilm, Tokyo, Japan).
Time-lapse Microscopy-Cells plated on 60-mm diameter dishes were cultured with or without GA at half confluency in a mini incubator under a microscope (CK30; Olympus, Tokyo, Japan) and subjected to time-lapse microscopy using 10ϫ and 20ϫ objective lenses for 48 -96 h (1 frame/10 min) and a digital imaging system equipped with a 3ϫ magnification zoom-in camera (DIGA XP22V; Panasonic, Osaka, Japan).
Fluorescence Staining and Confocal Laser Scanning Microscopic Observation-Cells were fixed with 3.7% formaldehyde in PBS and permeabilized in PBS containing 0.1% Triton X-100. For F-actin labeling, cells were incubated with rhodamine-phalloidin. For ␤III-tubulin labeling, cells were blocked in TBS containing 0.05% Tween 20 and 3% nonfat dried milk. Next, the cells were incubated with a mouse monoclonal anti-␤III-tubulin antibody (R&D Systems, Minneapolis, MN), washed with PBS, and incubated with HiLyte Fluor TM 488-conjugated goat anti-mouse IgG (AnaSpec, San Diego, CA). After washing and mounting with ProLong Gold anti-fade reagent (Molecular Probes, Eugene, OR), the cells were observed by confocal microscopy using an LSM510 META confocal laser scanning microscope equipped with Ar and He-Ne lasers (Carl Zeiss Japan, Tokyo, Japan) or a BIOREVO BZ-9000 fluorescence microscope system (Keyence, Osaka, Japan). Images were captured using 40ϫ and 63ϫ objective lenses, analyzed, and processed with the software Image Browser (Carl Zeiss Japan).
Measurement of CORT Levels-Inhibition of 11␤HSD2 enzyme activity was assessed by measuring the increase in CORT accumulation in the culture medium supernatant. The 11␤HSD2 enzyme activity inhibition was also examined in the reaction mixture reported by Mazzocchi et al. (29) except that CORT was used for the substrate and the reaction mixture contained 10 M triterpenoids as inhibitors of 11␤HSD2. The preparation of 11␤HSD2 from the cells for the reaction mixture was prosecuted by the methods of Brown et al. (30), and samples were incubated at 37°C for 1 h. The CORT levels were determined using a corticosterone enzyme immunoassay kit (Enzo Life Sciences International Inc., Plymouth Meeting, PA) in accordance with the manufacturer's protocol. The CORT levels were normalized to the cell numbers and expressed as ng/ml/ 10 6 cells.
Statistical Analysis-Experiments were performed in triplicate and repeated at least three times. The values quoted are given as means Ϯ S.D. Ordinary or repeated-measures analysis of variance followed by Student's t test or Tukey-Kramer's multiple comparison test was used to evaluate the   OCTOBER 21, 2011 • VOLUME 286 • NUMBER 42 statistical significance of differences between groups. IC 50 values were obtained using Prism 4.0 (GraphPad Software, Inc., San Diego, CA).

Ligand-Receptor Interaction between Triterpenoids and
11␤HSD2-The ASE-Dock module revealed that, in addition to CORT (the substrate of 11␤HSD2), GA, UA, OA, BA, and AA could bind to the LBS in the 11␤HSD2 model and exhibited similar binding orientations (Fig. 1A). The ligand fitting scores for GA, UA, AA, OA, and BA were Ϫ34.7, Ϫ28.6, Ϫ25.5, Ϫ24.8, and Ϫ24.7, respectively. The ASE scores for GA, OA, UA, AA, and BA were Ϫ12.0, Ϫ10.6, Ϫ6.8, 0.6, and 2.1, respectively. Furthermore, to create ligand-receptor interaction plots for each triterpenoid-11␤HSD2, the Ligand Interactions module of the Molecular Operating Environment 2009.10 was used, which provided a clearer arrangement of putative key intermolecular interactions that aid in interpretation of the three-dimensional juxtaposition of the ligands and the LBS in 11␤HSD2 (Fig. 1B). Our results revealed the presence of hydrophobic interactions between GA and Leu-229, Phe-265, Trp-276, Leu-282, and Leu-283. UA had hydrophobic interactions with Leu-229, Phe-265, and Leu-282, and OA had hydrophobic interactions with Pro-220, Leu-229, Phe-265, Leu-282, and Leu-283. BA and AA also showed some hydrophobic interactions with the residues in the LBS, but the hydrophobic and hydrophilic interactions were rather sporadic, which may result in unstable binding of the ligands.
Effects of Triterpenoids on 11␤HSD2 Enzyme Activity-The inhibition of 11␤HSD2 enzyme activity was assessed by measuring the increase in CORT in the presence of 10 M triterpenoids. GA, UA, and OA significantly inhibited 11␤HSD2 enzyme activity (Fig. 2), which indicates that the tested triterpenoids bind to the enzyme at the catalytic site and have inhibitory effects on the enzymatic function. The inhibition rates for GA, UA, and OA were 94, 70, and 40%, respectively.
Antiproliferative Effects of Triterpenoids on Normal SFME and Tumorigenic r/m HM-SFME-1 Cells-The normal and tumor cells were treated with AA, BA, GA, OA, or UA at 2.5-20 M for 24 h, and the effects on cell growth were examined. As shown in Fig. 3A, Ͼ80% of the normal cells survived the 24-h triterpenoid treatments at 2.5-10 M, whereas GA treatment at 10 M were not toxic to the cells at all. In contrast, ϳ20 and 50% of the proliferative capability of the tumor cells was inhibited by OA and UA at 10 M, respectively, and GA was more potent, inhibiting ϳ80% of the tumor cell growth at 10 M (Fig. 3B). The structures of the tested triterpenoids are also shown in Fig.  3 for reference.
Correlations between Cytotoxic Activity of Triterpenoids and Their logP and PSA Values-The mean 24-h IC 50 values of the cells treated with AA, BA, GA, OA, or UA were determined, and the logP and PSA values of each triterpenoid were calculated (Table 1). Among the triterpenoids tested, GA, OA, and UA showed lower IC 50 values for the tumor cells than for the normal cells. The difference in the IC 50 values for GA between the normal cells and tumor cells was greater than those for OA or UA. Next, the IC 50 values of the triterpenoids for the normal and tumor cells were analyzed with the logP and PSA values to examine whether there were correlations between them. The toxic activity of the tested triterpenoids against the normal cells showed no correlations with the logP and PSA values (Fig. 4, A and C). In contrast, the IC 50 values for the tumor cells exhibited good correlations with the logP and PSA values (Fig. 4, B and D).  Table 1.  Time-lapse Microscopy for GA-treated Tumor Cells-Thus far, we have analyzed the 24-h toxic effects of triterpenoids on tumor cell growth and found that GA was the most cytotoxic triterpenoid (Fig. 3B and Table 1). However, even at the very effective GA concentration of 10 M for selective toxicity against the tumor cells, Ͼ15% of the tumor cells survived the 24-h treatment (Fig. 3B). If the GA exposure is prolonged, and the surviving tumor cells can gain resistance against GA and recover their growth and proliferative abilities, GA will lose its advantage as a selective antitumor agent. Therefore, we treated the tumor cells with 10 M GA and performed time-lapse microscopy for 96 h to determine whether the tumor cells could still be affected after a long period of GA exposure. As shown in Fig. 5A and supplemental Movie S1, the tumor cells without GA treatment proliferated vigorously and reached confluency at 48 h. In contrast, the GA-treated tumor cells stopped cell growth immediately and exhibited morphologic changes typical of apoptosis, such as cell shrinkage and blebbing, and cell toxicity seemed to be induced in most of the cells within 96 h ( Fig. 5B and supplemental Movie S2).

Triterpenoids SFME cells r/m HM-SFME-1 cells logP PSA
Effects of GA on Cytoskeletal Disruption-The GA-treated tumor cells exhibited cell shrinkage and blebbing (Fig. 5B and supplemental Movie S2), suggesting that GA could disrupt cytoskeletal proteins because they provide the cell shape and maintain the cellular structure (31). Therefore, the tumor cells were treated with GA at 10 M for 24 h, and the effects on F-actin and ␤III-tubulin in the cells were analyzed. As shown in Fig. 6 (upper left Cont panel), the tumor cells showed broad lamellipodia and F-actin extensions. In contrast, the GAtreated tumor cells exhibited disappearance of the F-actin extensions, and instead, staining was rather granulated and condensed (upper right ϩGA panel). The loss of F-actin extensions at the periphery of the cell membrane was particularly obvious, and the tumor cells were no longer capable of maintaining the characteristic broad lamellipodia. ␤III-tubulin staining in the control cells (middle left panel) was homogeneous and mainly located in the cytoplasm, although its localization also extended into the broad lamellipodia. In contrast, the GA-treated cells showed non-homogeneous ␤III-tubulin staining (middle right panel with an enlarged image), and voids or puncture-like disorganizations of ␤III-tubulin were observed in the cytoplasm.
Effects of GA on 11␤HSD2 Inhibition-Normal and tumor cells were treated with GA at 10 M for 2, 8, and 24 h, and its effects on 11␤HSD2 expression and the CORT levels were analyzed to assess the inhibition of 11␤HSD2 enzyme activity. As shown in Fig. 7, A and B, endogenous 11␤HSD2 expression was higher in the tumor cells, and this was also the case in the GAtreated tumor cells. The endogenous CORT level was also higher in the tumor cells, and GA treatment strongly elevated it further in a time-dependent manner (Fig. 7C).

DISCUSSION
For cancer prevention and therapy, selectivity is an important issue. An ideal antitumor agent should be toxic toward malignant cells, with minimum toxicity toward normal cells. However, there are currently only limited numbers of such agents available for clinical use (32). An example of such an agent is GLEEVEC, which targets the oncogenic breakpoint cluster region-abelson tyrosine kinase responsible for chronic myeloid leukemia (33). However, mutations and overexpression of the target molecules often lead to drug resistance, owing to multiple genetic and epigenetic alterations in tumor cells (34,35). Moreover, tumor cells in advanced disease stages usually exhibit genetic instability and metabolic malfunction and are often resistant to conventional anticancer drugs (32). Therefore, with the intensive need for the development of more effective and safer agents for chemoprevention of cancer, natural products from plants, and their synthetic derivatives have been expected to be utilized in creating new and better chemopreventive and therapeutic agents (1). We previously reported that a naturally occurring triterpenoid, UA from apples, scarcely affected the viability of normal SFME cells, but markedly suppressed the growth of tumorigenic r/m HM-SFME-1 cells (2). Subsequently, we reported that another triterpenoid, GA from licorice, was not only selectively toxic to the tumor cells but also more potent than some clinically available antitumor agents in its selectivity (3). These abilities of UA and GA prompted us to investigate whether there were other potent triterpenoids that can inhibit tumor cell growth without hindering normal cell growth.
In the present study, the ASE-Dock module revealed that CORT (the substrate of 11␤HSD2), GA, UA, OA, BA, and AA exhibited similar binding orientations in the LBS in the 11␤HSD2 model, suggesting that triterpenoids can act as competitive inhibitors of CORT binding. The ligand fitting scores for the triterpenoid-11␤HSD2 docking model revealed that the inhibitory effects of the triterpenoids on 11␤HSD2 can be presumed in the order of GA Ͼ UA Ͼ AA Ϸ OA Ϸ BA. The ASE scores revealed that the inhibition against 11␤HSD2 can be surmised in the order of GA Ͼ OA Ͼ UA Ͼ AA Ϸ BA. Furthermore, the two-dimensional diagrams for the ligand-receptor interaction between the triterpenoids and 11␤HSD2 showed that GA, UA, and OA would be more stably placed in the LBS than AA and BA. These in silico results suggest that 11␤HSD2 can be inhibited by the triterpenoids in the order of GA Ͼ UA Ͼ OA Ͼ AA Ϸ BA. The results for the inhibitory effects of the triterpenoids on 11␤HSD2 enzyme activity also showed that 11␤HSD2 can be inhibited by the triterpenoids in the order of GA Ͼ UA Ͼ OA Ͼ AA Ϸ BA. In fact, our current in vitro results showed that GA, UA, and OA were selectively toxic against the tumor cells and that the selectivity of GA was more potent than any of the other triterpenoids tested. The selectivity of GA was such that, at 10 M, it inhibited ϳ80% of the tumor cell growth within 24 h without affecting the normal cells and induced apoptosis in most of the tumor cells within 96 h. By analyzing the lipophilicity and toxic effects of some compounds from plants, Jiao et al. (36) reported that the logP values measured for some antitumor compounds were consistent with those predicted on the basis of their chemical structures and that logP values for measurement of the relative lipophilicity could be applied. Noshita et al. (37) reported good correlations between the PSA and the biological activities of an antitumor compound from wasabi (Wasabia japonica) and its synthetic derivatives. In the present study, we found that the IC 50 values of the tested triterpenoids against the tumor cells showed a good direct and an inverse correlation with the logP and PSA values, respectively, suggesting that a certain degree of lower lipophilicity or higher PSA value is necessary for the toxic effects against the tumor cells. These findings also indicate that a specific intracellular delivery system for triterpenoids may exist in the tumor cells. Fig. 8 shows the proposed mechanism underlying the 11␤HSD2 inhibition by GA and the tumor cell-selective toxicity. Our present results showed that the tumor cells exhibited higher 11␤HSD2 expression, which suggests that the metabolic cascades through 11␤HSD2 could be relatively easily triggered and that this may make the tumor cells sensitive to 11␤HSD2 inhibition. It has been reported that patients with specific mutations in the 11␤HSD2 gene developed severe CORT-dependent hypertension and other features of apparent mineralocorticoid excess (38). The other enzyme that could increase CORT is 11␤HSD1. Increased 11␤HSD1 and decreased 11␤HSD2 expressions up-regulated CORT (39). 11␤HSD1 acts as an oxoreductase converting 11-dehydrocorticosterone to CORT, whereas 11␤HSD2 is a dehydrogenase that inactivates CORT (40), suggesting that increased activity of 11␤HSD1 and/or decreased activity of 11␤HSD2 up-regulate CORT. GA is a well known specific inhibitor of 11␤HSD1 and 11␤HSD2 (41). Furthermore, our current in silico and in vitro results revealed that 11␤HSD2 could be inhibited and CORT was upregulated by GA. Therefore, these reports and our present findings indicate that up-regulation of CORT is strongly associated with 11␤HSD2 inhibition by GA. CORT is a known down-regulator of GSH (42), and inhibition of 11␤HSD2 activity by GA up-regulated CORT and prevented tumor growth and metastasis (13), indicating that 11␤HSD2 inhibition by GA induces down-regulation of GSH and consequently leads to suppression of tumor cell growth. This was exactly the case in our present and previous (20) studies, and the previous study also clearly revealed that the reactive oxygen species production levels were much higher, and the GSH levels were much lower in the tumor cells than in the normal cells. In the tumor cells, FIGURE 7. 11␤HSD2 expression is higher in tumor cells than in normal cells, and GA up-regulates the CORT level in tumor cells. Normal SFME and tumorigenic r/m HM-SFME-1 cells were treated with GA at 10 M for 2, 8, and 24 h, and its effects on 11␤HSD2 expression were analyzed. The CORT level was also analyzed to assess the inhibition of 11␤HSD2 enzyme activity by GA. A, images of a Western blotting analysis. GAPDH was used as a loading control. M, Size markers. B, images of the Western blotting analysis were captured, analyzed, and processed with an image analyzer. As shown in the rectangle (A) and the closed bars (B), endogenous 11␤HSD2 is higher in the tumor cells, and this is also the case in the GA-treated tumor cells. C, the endogenous CORT level is higher in the tumor cells, and the GA treatment elevates it further in a time-dependent manner. Each point is the mean Ϯ S.D. of at least three experiments. ***, p Ͻ 0.001 by Student's t test for untreated normal cells versus untreated or treated tumor cells. ϩϩϩ, p Ͻ 0.001 by Student's t test for untreated tumor cells versus treated tumor cells. FIGURE 8. Proposed mechanism underlying the 11␤HSD2 inhibition by GA and tumor cell-selective toxicity. GA is moved into the cells, possibly by a specific intracellular delivery system (Fig. 4, B and D). 11␤HSD2 is overexpressed in various tumor cells (12)(13)(14)(15)(16), including the tumorigenic r/m HM-SFME-1 cells (Fig. 7, A and B). GA inhibits 11␤HSD2 (Figs. 1, A and B, 2, and 7C; 13, 41), up-regulates CORT ( Fig. 7C; 13), and down-regulates GSH (20). GA also down-regulates Ras (3). Decreased GSH and Ras make the cells vulnerable to oxidative damage (20) and cytoskeletal disruption ( Fig. 6; 44, 45), which consequently lead to apoptosis of the tumor cells ( Fig. 5B and supplemental Movie S2).
which produced higher reactive oxygen species levels and were already under oxidative stress, addition of GA, which would further elevate the production of reactive oxygen species, could lead to apoptosis of the tumor cells. Furthermore, the decrease in the GSH level by GA, which creates a redox imbalance, would lower the antioxidant capacity of the cells, which could contribute to death of the tumor cells. To the best of our knowledge, there have been no reports directly associating disruption of cytoskeletal proteins with 11␤HSD2 inhibition by GA. Our current results could neither be interpreted as they have some association. Therefore, it remains unknown whether or not the disruption of the cytoskeletal proteins was a direct action of the 11␤HSD2 inhibition. However, our previous study revealed the potential involvement of Ras in the selective toxicity by GA against the tumor cells (3). Ras pathways have been targeted for the development of chemotherapeutic interventions against tumor cells by interfering with cytoskeletal proteins such as actin (31) and tubulin (43). Furthermore, in our previous (44,45) and current studies, we found that GA disrupted F-actin and ␤III-tubulin, which could lead to loss of functionality for cell growth or proliferation and result in apoptosis of the cells. These findings indicate that disturbances in cytoskeletal proteins may be the factors involved in the toxic effects of GA against the tumor cells. It has been reported that interference with actin or microtubule functions associated with the integrity of the cytoskeleton could be utilized as a strategy for developing novel antitumor treatments (31). It has also been reported that some agents bind directly to cytoskeletal proteins, disrupt the integrity of tumor cells, and inhibit cell growth and proliferation (43,46). Although further studies are required to elucidate whether the GA-mediated interference with the cytoskeletal proteins was a direct action or an indirect effect through some signaling pathways, our present study clearly reveals that disruption of cytoskeletal proteins could be one of the important factors that led to the induction of the selective tumor cell toxicity.
In conclusion, we utilized in silico analyses of the ligandreceptor interaction between triterpenoids and 11␤HSD2 to predict the antitumor toxicity of some triterpenoids. The results of the in silico analyses were in good agreement with the in vitro experimental evidence. To the best of our knowledge, this is the first report on in silico prediction of the toxic effects of triterpenoids on tumor cells by 11␤HSD2 inhibition. 11␤HSD2 has been reported to be associated with various cancers and identified as a possible target for the prevention and therapy of cancer (12)(13)(14)(15)(16). In recent years, molecular modeling has gained much importance in the field of drug discovery and development (47)(48)(49), and further in silico studies of 11␤HSD2 and its interactions with possible ligands are expected for the successful development of antitumor drugs targeting 11␤HSD2.