Gerrit E. Benary,
Frank Kilgenstein,
Sascha Koller and
Jürgen Scherkenbeck
*
University of Wuppertal, School of Mathematics and Natural Sciences, Gaussstrasse 2042119, Wuppertal, Germany. E-mail: scherkenbeck@uni-wuppertal.de
First published on 10th January 2025
Betulinic acid and other herbal pentacyclic triterpenes have attracted interest in cancer research as these natural products induce apoptosis and suppress tumor progression. However, the molecular basis of the antitumor effect is still unknown. Here we show that monophthalates of betulinic acid and related triterpenes inhibit GDP/GTP exchange in oncogenic K-RAS4B proteins via the PI3K/AKT downstream cascade. According to a binding model based on molecular modelling, these derivatives act like a molecular glue that stabilizes an unproductive K-RAS4Ballo:SOS complex. This represents a new mode of action and could be an attractive route for targeting RAS-related cancers.
![]() | ||
Fig. 1 Structure of a K-RAS4Bcat:SOS:K-RAS4Ballo complex (a) PDF code: 7KFZ and structure of K-RAS4B (b) PDB code: 4OBE. |
Pentacyclic triterpenes constitute a large class of complex herbal natural products, that can be divided into the lupane, ursane, and oleane types (Fig. 2).6,7 Common plant families that are rich in these compounds include Asteraceae, Fabaceae, and Rosaceae.8 A notable example of a plant containing pentacyclic triterpenes is Centella asiatica, commonly known as gotu kola, which has been used in traditional Chinese medicine for its wound healing and anti-inflammatory properties.9,10 Several pentacyclic triterpenes (Fig. 2) have attracted attention in cancer research due to their cytotoxic effects.11–13 Mechanistic studies have demonstrated the ability of betulinic acid (2) to induce apoptosis, or programmed cell death, in a variety of cancer cell lines. A recent target identification-based analysis of cervical cancer cell apoptosis identified eight potential key targets for betulinic acid (2), including the long-known AKT.14 Other reported cellular effects of betulinic acid (2) encompass inhibition of angiogenesis and suppression of tumor progression as well as disruption of mitochondrial function leading to release of cytochrome c and activation of caspases.15–17 Ursolic acid (4) has been found to inhibit the proliferation of cancer cells by modulating multiple signaling pathways.18,19 Among others it suppresses the activation of NF-κB, a key regulator of inflammation and cancer.20,21 Furthermore, oleanolic acid (6) was shown to exert its antiproliferative effect on K-RAS-transformed cells by autophagy.22 In recent years, numerous derivatives of natural pentacyclic triterpenes have been prepared that exhibit improved or altered anticancer activities.14,23–25 Several literature reports indicate that in particular ursolic and betulinic acid derivatives also inhibit RAS dependent malignant cell-lines, in particular those with G12D and G13D mutations.12,16,17,26–28 In addition, an in vivo antitumor effect in K-RASG12D mice has been observed for SYK023, the p-methoxy-phenylacetic ester of betulinic acid.29 However, the molecular basis of the RAS-inhibitory effect is still unclear.
In the context of our ongoing search for new RAS inhibitors, we were interested in the more general question of whether pentacyclic triterpene acid derivatives could inhibit K-RAS4B by interfering with the guanine exchange process.
![]() | ||
Fig. 3 Pentacyclic triterpenes derivatives prepared for biological evaluation. Grey color: compounds regarded as inactive in the GDP/GTP exchange assay (IC50 > 100 μM). |
Betulin contains an equatorial secondary hydroxy group at position 3 and an axial hydroxymethyl group at position 28, whose reactivity is further reduced by the neopentyl environment. Nevertheless, the primary OH-group is still the more reactive one and can be selectively mono-acetylated to derivative 29.38
A similar chemoselectivity was observed with maleic and succinic anhydride in the preparation of compounds 30 and 31.24,39 Only a small amount (10%) of the double succinylated product 32 was isolated.40
In addition to the BA derivatives, the C-17 carboxy group of OA was converted into the bioisosteric tetrazole derivative 39 (23%).41 The reaction of OA with N-Boc-ethylendiamine and subsequent deprotection with dilute HCl yielded the ethylamine derivative 40 (78%).42 In order to further evaluate the role of the carboxy group in position 17, hydroxamic acid derivative 41 was prepared from oleanolic acid via a four step sequence, including acetylation of the 3-hydroxy group, hydroxamic acid formation with hydroxylamine, deprotection of the 3-acetyl group, followed by reaction with 4,5-difluorophthalic anhydride.
The oleanoyl 4,5-difluoro-2-(hydroxycarbamoyl)benzo-ate 42 was obtained in low yield by direct conversion of 4,5-difluorophthalic monoester 43 with hydroxylamine. The ketolactone 49 was obtained in one-step by ozonolysis of OA in 77% yield. Subsequent stereoselective reduction with NaBH4 afforded in 92% yield the hydroxylactone 50 which was finally converted into the 3-monoester derivatives 51 and 52.43,44
Cmpd | Wild-type IC50 [μM] | G12D IC50 [μM] | G12V IC50 [μM] | Solubility [μM] |
---|---|---|---|---|
a Values were calculated from eight data points. At least three independent determinations were performed. The asterisk (*) indicates a single determination. n.d.: not determined. | ||||
Derivatives of betulinic acid | ||||
14 | 37.3 ± 19.1 | 48.5 ± 12.7 | 51.1 ± 6.2 | 447 |
18 | 1.27 ± 0.69 | 1.83 ± 1.45 | n.d. | 255 |
19 | 15.3 ± 3.8 | 10.4 ± 4.7 | 16.6 ± 5.2 | 191 |
23 | 23.3 ± 2.9 | 11.8 ± 1.0 | 9.61 ± 4.82 | 434 |
24 | 15.7 ± 4.3 | 8.63 ± 3.72 | 15.4 ± 2.9 | 378 |
25 | 39.1 ± 5.4 | 5.00 ± 0.91 | 9.79 ± 3.41 | 240 |
26 | >100 | 36.1 ± 14.7 | 26.4 ± 3.9 | 158 |
27 | n.d. | 13.8 ± 4.2 | n.d. | n.d. |
28 | 2.15 ± 0.17 | 0.59 ± 0.19 | 0.61 ± 0.14 | n.d. |
![]() |
||||
Derivatives of betulin | ||||
30 | 69.2 ± 6.8 | 6.80 ± 2.32 | 3.19 ± 0.19 | n.d. |
31 | 25.6 ± 4.3 | 10.4 ± 0.8 | 12.7 ± 5.5 | n.d. |
32 | 1.01 ± 0.31 | 0.71 ± 0.18 | 1.37 ± 0.28 | n.d. |
33 | 2.09 ± 0.59 | 1.45 ± 0.93 | 2.13 ± 0.87 | 416 |
34 | 1.18 ± 0.02 | 0.66 ± 0.13 | 2.85 ± 0.14 | 410 |
![]() |
||||
Derivatives of ursolic acid | ||||
36 | 3.46 ± 1.94 | 2.64 ± 1.75 | 2.52 ± 0.99 | 424 |
37 | 3.94 ± 1.02 | 2.53 ± 0.63 | 2.15 ± 0.31 | 220 |
38 | 53.6 ± 7.3 | 55.6 ± 6.5 | n.d. | 383 |
![]() |
||||
Derivatives of oleanolic acid | ||||
39 | 12.8 ± 2.9 | 27.8 ± 2.6 | 9.18 ± 0.99 | n.d. |
41 | 4.68 (*) | 7.61 ± 0.70 | 5.47 (*) | n.d. |
43 | 1.50 ± 0.08 | 2.22 ± 0.29 | 2.71 ± 0.39 | 394 |
44 | 1.33 ± 0.24 | 2.88 ± 0.89 | 2.51 ± 0.24 | 424 |
45 | 1.52 ± 0.50 | 3.06 ± 0.79 | 5.62 ± 0.93 | 390 |
46 | 4.15 ± 1.32 | 3.67 ± 0.35 | 8.12 ± 1.24 | 421 |
47 | 2.71 ± 1.30 | 1.61 ± 0.85 | 2.90 ± 0.51 | 302 |
51 | 7.18 ± 0.68 | 26.5 ± 5.9 | 8.57 ± 1.09 | 99 |
52 | 10.6 ± 2.9 | 27.1 ± 4.2 | 8.70 ± 0.59 | 320 |
![]() |
||||
Reference compounds | ||||
53 | <0.3 | <0.03 | <0.03 | 322 |
54 | 13.4 ± 2.1 | 13.7 ± 5.1 | 10.9 ± 0.6 | 114 |
55 | 0.13 ± 0.08 | <0.03 | <0.03 | 277 |
Disappointingly, all tested natural pentacyclic triterpenes 1–7 (Fig. 2) did not show any notable inhibition of GDP/GTP exchange for the most frequent K-RAS4BG12D mutant. In particular, for betulinic and ursolic acid IC50 values of >300 μM were measured. However, a significant inhibition of GDP/GTP exchange was found for the semi-synthetic betulinic monophthalic ester 19 (IC50 10.4 ± 4.7 μM) which has already been described as a cytotoxic compound with IC50 values in the range of up to 10 μM against a series of cancer cell lines.29
To establish an initial structure activity model, we prepared several betulinic acid derivatives, all of which were tested for their inhibition of the G12D and the more challenging G12V mutant, as well as their selectivity for K-RAS4Bwt (Table 1). Acylation of the 3-OH group (8, 9) and formation of esters or amides of the 17-carboxy group (10, 21) led to inactive betulinic acid derivatives. The inactivity of these lipophilic derivatives cannot be a mere solubility effect, as shown by the glycine-derived acid 12, the polar ester 13, and the amine 17, all of which were inactive, although their solubility is comparable to that of the reference compound BI2852 (53).
The relevance of the 3-OH group is further corroborated by betulonic acid 20, the 3-keto derivative of BA, which was also inactive. The platanic acid derivatives 15 and 16, obtained by oxidation of the 19-isopropenyl group of betulinic acid, had a limited solubility and IC50 values could be measured only up to 30 μM. The well soluble monophthalic acid ester 14 (IC50 48.5 ± 12.7 μM) showed lower inhibition of the GDP/GTP exchange than the parent compound 19 (IC50 10.4 ± 4.7 μM), suggesting that an aldehyde function is not a suitable substitute for the isopropenyl group. On the other hand, oxidation of the isopropenyl group led to the corresponding carboxylic acid 27, which showed GDP/GTP exchange inhibition (13.8 ± 4.2 μM) at the same level as BA derivative 19, even with the fully protected 3-hydroxy and 17-carboxy group. Taking these initial data together, it becomes evident that the SAR of betulinic acid is quite narrow.
In a second step, we therefore focused on phthalic mono-esters. Unfortunately, even the change from phthalic (19) to 2,3-naphthoic monoester 22 led to a loss of GDP/GTP exchange inhibition. The electron-rich 4,5-dimethoxy phthalate 26 exhibited only moderate activity in the G12D mutant. On the other hand, the replacement of the phenyl ring by an electron-deficient pyrazine-2,3-dicarboxylic monoester (23) or the introduction of electron-withdrawing fluorine atoms (24) led to betulinic acid derivatives with similar activity as the parent compound 19. Both betulinic acid derivatives had similar activities for the G12D and G12V mutant.
The corresponding 4,5-difluorophthalic monoester 25 showed a somewhat improved inhibition (IC50 5.00 ± 0.91 μM) of GDP/GTP exchange for K-RAS4BG12D and a slightly lower activity for the G12V mutant but a significant selectivity by a factor of 4–6 for K-RAS4Bwt. The thiophene derivative 18 has an impressive IC50 value (G12D mutant) of 1.83 ± 1.45 μM, which corresponds to an improvement of at least a factor of five compared to compound 19, but completely at the expense of selectivity for K-RAS4Bwt. Most strikingly, the HIV inhibitor bevirimat (28) inhibited the GDP/GTP exchange at high nanomolar levels for both the G12D (IC50 0.59 ± 0.19 μM) and G12V (IC50 0.61 ± 0.14 μM) mutant, with a selectivity of 3–4 fold (IC50 (wt) 2.15 ± 0.17 μM). This finding indicates, that betulinic acid derivatives modified in position 3 have the potential for nanomolar GDP/GTP exchange inhibitors. This is especially important for the G12V mutant, which is notoriously difficult to address due to the non-functionalized valine side-chain.
Those results encouraged us to extend our studies to related pentacyclic triterpenes such as betulin, oleanolic acid, and ursolic acid. In general, the structure activity relationships (SAR) were similar to betulinic acid (Fig. 3 and Table 1). Blocking of the hydroxy- and carboxy groups by esterification or by introduction of an amine function (35, 40) always resulted in inactive derivatives. Within the betulin series GDP/GTP exchange inhibition could be brought to a level equal to the best betulinic acid derivatives. In particular, betulin derivative 30 exhibits a good activity for the G12V mutant (IC50 3.19 ± 0.19 μM) together with a remarkable selectivity by a factor of 15–20 (IC50 (wt) 69.2 ± 6.8 μM) for wildtype K-RAS4B. Introduction of a second succinic monoester or a 4,5-difluorophthalic monoester at position 3 resulted in betulin derivatives 32 and 34, both of which inhibited the GDP/GTP exchange in K-RAS4BG12D in the high nanomolar range (0.71 ± 0.18 μM and 0.66 ± 0.13 μM), unfortunately again at the expense of selectivity. Dihydroxylation of the oleanolic acid Δ12-double bond and subsequent reduction of the 3-keto group with NaBH4 resulted in the formation of the highly stable but inactive hydroxylactone 50, demonstrating again the indispensability of a free carboxylic acid in position 17. Expectedly, only a reduced biological activity was found for the 4,5-difluorophthalic acid monoester 51 (IC50 26.5 ± 5.9 μM) prepared from the hydroxylactone 50. In the oleanolic and ursolic acid series, best inhibitory activities in the low micromolar range were obtained with the same monoesters in position 3 as for betulinic acid. The oleanolic acid derivatives 43, 44, and 47 showed a promising activity against the G12V mutant. Unfortunately, all oleanolic acid derivatives exhibited no selectivity in favor of the oncogenic K-RAS4B mutants. Worthwhile to mention, the oleanolic 4,5-dimethoxy-phthalic monoester 47 retained activity while the corresponding derivatives of betulinic (26) and ursolic acid (38) were found significantly less active. The relevance of the free phthalic carboxy group is further corroborated by hydroxamic acid 42 and diester 48 which both turned out completely inactive.
For selected betulinic acid derivatives quantitative ERK and AKT phosphorylation was measured in SNU-1 (K-RAS4BG12D) and A375 cells (K-RAS4Bwt) to determine whether pentacyclic triterpenes act via the RAS-RAF-MEK-ERK or the RAS-PI3K-PDK1-AKT downstream cascade (Table 2). The betulinic acid derivatives tested, had no activity on the ERK pathway, but inhibited AKT phosphorylation with micromolar IC50 values in agreement with literature reports.50,51 All betulinic acid derivatives showed selectivity for the K-RAS4BG12D mutant cell line SNU-1. No activity was detected in A375 cells carrying a BRAF but no K-RAS mutation. Our results demonstrate, that betulinic acid derivatives not only inhibit GDP/GTP exchange in a biochemical assay but disrupt specifically downstream signalling via the RAS-PI3K-PDK1-AKT pathway in RAS mutant cells.
Cmpd | SNU-1 pERK | SNU-1 pAKT | A375 pERK | A375 pAKT |
---|---|---|---|---|
a In general, two independent determinations were performed. | ||||
19 | >300 | 19.0 | >1000 | >1000 |
23 | >1000 | 61.7 | >1000 | >1000 |
24 | 134 | 13.0 | >1000 | >1000 |
25 | >300 | 18.1 | >1000 | >1000 |
26 | >300 | 83.7 | >1000 | >1000 |
However, neither the GDP/GTP exchange assay nor the ERK and AKT phosphorylation data provide any further information on the exact mode of action of the pentacyclic triterpene derivatives. In principle, the inhibition of GDP/GTP exchange can be achieved either by binding of the ligands to K-RAS4B, to SOS or to a K-RAS4B:SOS complex. K-RAS4B and even more K-RAS4B:SOS complex inhibitors with low micromolar or nanomolar IC50 values are still considered highly interesting.
To obtain further information on the mode of action, selected GDP/GTP exchange inhibitors were studied in additional assays (Table 3). A target engagement assay, which is based on bioluminescence resonance energy transfer (BRET), was used to measure target engagement at protein complexes in live cells (Fig. 5a). The BRET effect results from an energy transfer between the target luciferase fusion protein and a cell-permeable fluorescent tracer in intact cells. We used this assay to study the displacement of fluorescently labelled BI2852 (54) from its SI/SII-binding site.52 Expectedly, the SOS-binder BI3406 (55) did not show any effect. Noteworthy, the strong SII-site ligand MRTX1133 (53) also exhibited a nanomolar KD value for the displacement of BI2852 (54) at the SI/SII-site. Apparently, binding of MRTX1133 (53) to the SII-site induces conformational shifts at the nearby SI/SII-site that release BI2852 (54).
Cmpd | K-RAS4Bwt:SOS PPIb KD [μM] | K-RAS4Bwt NanoBRET TEc KD [μM] | K-RAS4BG12D NanoBRET TEc KD [μM] |
---|---|---|---|
a In general, two independent determinations were performed. n.d.: not determined.b Protein–protein interaction assay of (GTP)K-RAS4Bwt and SOS.c Target engagement cell-assay: displacement of BI2852 from the K-Ras4B SI/II binding-site. | |||
18 | >50 | 21.8 | 20.1 |
32 | >100 | 40.8 | 29.1 |
33 | 37.4 | n.d. | 26.4 |
34 | >100 | n.d. | >100 |
36 | >50 | n.d. | >50 |
37 | >50 | n.d. | 16.6 |
47 | >50 | 39.0 | >50 |
53 | 0.09 | 0.068 | 0.042 |
55 | 0.014 | n.d. | >100 |
Noteworthy, the majority of the particularly active GDP/GTP exchange inhibitors (32, 34, 36, 47) displace BI2852 (54) from the SI/SII-site only at higher concentrations. Since all these pentacyclic triterpene derivatives have a higher solubility than BI2852 (54), this is not a plausible explanation for the reduced activity in the NanoBRET assay. Rather, it seems more likely that the monophthalate derivatives bind to a different site.
The interaction of (GTP)K-RAS4Bwt and SOS was measured with a HTRF assay using an SOS-specific antibody labeled with terbium cryptate as HTRF (homogeneous time-resolved fluorescence) donor and a K-RAS4Bwt specific antibody labeled with an HTRF acceptor (XL665).52 Upon (GTP)K-RAS4Bwt:SOS complex formation, excitation of the donor antibody triggers fluorescence resonance energy transfer (FRET) towards the acceptor antibody, which emits specifically at 665 nm. Therefore, K-RAS4B or SOS inhibitors, which cause dissociation of the K-RAS4B:SOS complex, reduce HTRF signaling (Fig. 5b). The K-RAS4B site-II ligand MRTX1133 (53) and especially the SOS ligand BI3406 (55) are potent inhibitors of K-RAS4Bwt:SOS complex formation (Table 3).
Remarkably, a significant discrepancy in activity between the GDP/GTP exchange assay and the K-RAS4Bwt:SOS interaction assay was observed for all pentacyclic triterpenes. The KD values for dissociation of the K-RAS4B:SOS complex are a factor of 20–50 higher for all potent, single-digit micromolar GDP/GTP exchange inhibitors, indicating that GDP/GTP exchange is inhibited at a concentration where the K-RAS4B:SOS complex is still intact. Furthermore, surface plasmon resonance (SPR) measurements showed that the potent GDP/GTP exchange inhibitor bevirimat (28) does not bind to either K-RAS4BG12D (KD > 300 μM) or SOS (KD > 1000 μM) alone.
These data imply that the pentacyclic triterpenes studied here do not act directly on K-RAS4B or SOS but potentially by interaction with the K-RAS4B:SOS complex. Such a mode of action would be very interesting to the field, as only a limited number of small molecules is known to interfere with the K-RAS4B:SOS complex.53–57 In order to develop a binding hypothesis, various molecular modeling studies were carried out with selected triterpene derivatives.
Prior to in silico docking with Glide (Schroedinger molecular modeling suite), a sufficient number of energetically favored conformers was generated by a conformational search (mixed torsional/low mode sampling) procedure. Unfortunately, all known interaction sites on K-RAS4B, SOS, and at the K-RAS4Bcat:SOS interface turned out to be too small for the pentacyclic triterpene scaffolds, leading to low docking scores (4–5 kcal mol−1) with questionable significance.50,51 As expected, insignificant docking scores were also found for K-RAS4Bcat in a K-RAS4Ballo:SOS:K-RAS4Bcat complex (PDB code: 7KFZ).
Mg2+ is tightly bound to the triphosphate group of GTP through an ionic interaction with the γ-phosphate group, a coordinative bond to the neighboring β-phosphate, and a water mediated hydrogen bond to the α-phosphate residue. Together with a second water molecule Mg2+ forms an octahedral complex when bound to GTP (Fig. 6). An overlay of the energy-minimized dihydrate of magnesium mono (4,5-dimethoxy-2-(methoxycarbonyl)benzoate) and GTPγS·Mg2+·2H2O reveals an excellent match with both Mg2+ cations occupying almost identical positions. In this superposition the carboxy group of the monophthalate moiety agrees well with the γ-phosphate, and the ester carbonyl group with the phosphor oxygen double bond of the β-phosphate. Noteworthy, the quality of the overlay is considerably lower for the corresponding Mg2+·GDP·4H2O complex. The monophthalate building block, which is present in almost all active triterpene derivatives, can therefore be regarded as a Mg2+-chelator with a bite-angle similar to the terminal diphosphate in GTP. This finding prompted us to examine the monophthalate triterpenes for their interaction with the (Mg2+·GTP)K-RAS4Ballo:SOS interface.
![]() | ||
Fig. 6 Superposition of GTPγS·Mg2+·2H2O (PDB code: 4DSO) and monomethyl Mg phthalate·2H2O. |
A SiteMap search identified a smaller binding site (blue color) in the guanidine pocket region and a large major binding site (orange color), encompassing switch I, switch II, and the α2-helix of K-RAS4Ballo (Fig. 7a). While the small binding-site is too far away from the Mg2+ cation, the large site consists of two sub-pockets, both of which are able to accommodate a large pentacyclic triterpene scaffold and allow the formation of magnesium complexes.
![]() | ||
Fig. 7 Interaction of oleanolic acid derivative 47 and related pentacyclic triterpene derivatives with the allosteric K-RAS4B:SOS interface (PDB code: 7KFZ). Cyan: K-RAS4Ballo, magenta sphere: Mg2+, green carbons: GTP, orange carbons: pentacyclic triterpenes. (a) With SiteMap calculated binding sites (blue and orange). (b) Overlay of GDP/GTP exchange assay hits 34, 37, 38, 43, and 47 in the high-affinity SWII:SOSRem pocket. (c) Ligand protein interaction diagram of derivative 47 with the K-RAS4B:SOSRem domain; colors of the boxes indicate polarities of the amino acids. (d) MD simulation based lowest free binding energy (ΔG = −124 kcal mol−1) snapshot of compound 47 at the SWII:SOSRem interface. |
In silico dockings at the K-RAS4B:SOSCdc25 subpocket (Fig. 7a) revealed only low docking-scores due to steric clashes with the guanosine triphosphate in close proximity. In contrast, the docking scores for the K-RAS4B:SOSRem subpocket (Fig. 7a) were good to excellent (Table 4). In particular, compounds 38 and 47 showed docking-scores comparable to GDP (PDB code 4DSU), but were inferior to the stable GTP analogue GppNHP (PDB code 7KFZ), which is known to bind significantly stronger to K-RAS4B than GDP. All GDP/GTP exchange assay hits examined (Table 4) bind to the high-affinity K-RAS4Ballo:SOSRem site in a similar manner (Fig. 7b).
To obtain additional information on the binding mode, a 500 ns molecular dynamics simulation was conducted for oleanolic acid derivative 47 bound to the high-affinity K-RAS4Ballo:SOSRem site of the crystal structure 7KFZ. A trajectory with a total of 2000 frames was created by taking snapshots every 250 ps. MM/GBSA calculations were performed for a total of 1000 frames previously extracted from the trajectory. The mean ΔG binding energy was −60.2 ± 12.4 kcal mol−1 with a maximum value of −124 kcal mol−1.
The MD simulation revealed that a stable ligand–protein complex was maintained throughout the complete simulation (Fig. 8). The lowest energy poses of compound 47 bind almost identically to a narrow channel formed by residues of K-RAS4Ballo, SWI (Val29, Glu31, Tyr32, Asp33, Thr35, Ile36, Tyr40), and SWII (Gln61, Tyr64, Ser65, Ala66, Met67) as well as the residues Gln683, Pro684, Leu687, Met726, Lys728, and Trp729 of the SOSRem domain (Fig. 7c and d). The phthalic acid moiety is accommodated in a small lipophilic pocket, formed by RAS-residues of switch I, with the two methoxy groups located close to the hydroxy groups of the GTP ribose (Fig. 7d). The phthalic acid displaces two water molecules from the GTP·Mg2+ complex and is tightly anchored to its binding site by strong ionic interactions of Mg2+ with the negatively charged carboxylate anion and the ortho-ester carbonyl group (Fig. 7c). The octahedral coordination of the Mg2+ ion, as well as the interactions with the β- and γ-phosphate groups, are not affected by the binding of compound 47, in contrast to Ser17 and especially Thr35, which are displaced from their original positions by the aromatic ring. Further stabilization of the phthalic acid is achieved through water-mediated hydrogen bonds with Asp13, Asp30 and Glu31, as well as through a lipophilic interaction with Tyr32 (Fig. 8). The hydrocarbon scaffold of the oleanolic acid derivative 47 is additionally anchored by the 17-carboxy group, which forms strong water-mediated hydrogen bonds with the SOSRem residues Lys727, and especially Lys728, as well as a direct hydrogen bond with Met726 (Fig. 8). Apparently, the monophthalates of oleanolic acid and related pentacyclic triterpenes act as a molecular glue that stabilizes the K-RAS4Ballo:SOSRem complex.
At first glance, the proposed mode of action appears counter-intuitive, as stabilization of a K-RAS4Ballo:SOSRem complex should accelerate the rate of GDP exchange and thereby increase the level of active (GTP)RAS. However, the activation of SOS mediated by K-RAS4Ballo and the subsequent exchange of GDP/GTP in K-RAS4Bcat represents a complex feedback loop that depends on several conformational changes in all three proteins involved. The first step in GDP/GTP exchange is the insertion of the helical hairpin from helices αH and αI of SOSCdc25, which leads to move-ments of the switch I and II regions that open the nucleotide binding site and facilitate the release of GDP.58
The nucleotide-free K-RAS4B:SOS complex is stabilized by the side chain of SOSLeu938, which prevents Mg2+ binding, and SOSGlu942, which overlaps with the binding site of the α-phosphate. It is also worth noting that during the entire MD simulation time of 500 ns, no significant conformational changes occurred within these two critical amino acids and all K-RAS4Bkat residues (Phe28, Val29, Asp30, Asn116, Lys117, Asp119, Ala146) that interact with guanosine and are thus important for GTP reentry.
In contrast, the activation process of SOS upon binding of K-RAS4Ballo is less well understood. For H-RAS, it has been shown that the formation of the ternary complex (GppNHp)H-RASY64A(allo):SOS:H-RASwt induces a profound conformational change of the Rem domain relative to the CDC25 domain, which allosterically enhances the interactions between the switch I region of nucleotide-free K-RAS4Bcat and SOS, thereby promoting the rate of nucleotide exchange.3 Furthermore, it was shown that a SOS protein containing the Rem-domain preceding DH (Dbl homology) and PH (Pleckstrin homology) domains is autoinhibited by the DH domain, which packs against the face of the Rem domain (Fig. 9).59,60 The SOSRem residues 724–729, which are an essential part of the triterpene binding pocket, are completely occluded by the DH domain in the autoinhibited SOS state (Fig. 9). Our results suggest that oleanolic acid derivative 47 and related pentacyclic triterpenes act in a similar manner by blocking the same residues essential for the activation of SOS, thereby stabilizing an unproductive K-RAS4Ballo:SOS complex or preventing the formation of highly active SOS states.5 The limited selectivity of several monophthalate-modified pentacyclic triterpenes can also be explained by the proposed mode of action, since oncogenic (GTP)K-RAS4B proteins activate SOS by binding to the allosteric site to a similar extent as the wild-type form.61
![]() | ||
Fig. 9 Autoinhibition of SOS-DH-PH. PDB code: 1XD4, green carbons: residues involved in compound 47 binding, PH domain: orange color, DH domain: yellow color, Rem domain: magenta color, CDC25 domain: cyan color, EPS8 linker: red color. |
Disrupting the K-RAS4Ballo:SOSRem interaction appears an attractive new strategy for combating RAS-related cancers. However, the potential of this mode of action for the development of mutant-selective K-RAS4B therapeutics needs to be further investigated as a prerequisite for effective clinical translation.
For each sample, 5 μL K-RAS working solution in assay buffer (50 mM HEPES pH 7.5, 4 mM MgCl2, 2 mM EGTA, 0.01% Brij35, 1 mM TCEP) were transferred into a suitable assay plate (e.g. Greiner #784075). The test compound was added with an echoacoustic dispenser (Beckman Coulter, Brea, CA, USA) in a concentration range from 3000 μM to 3 μM (8-point dilution). After addition of the test compound 5 μL of SOS1-GTP mix in assay buffer were added. The reaction mixture was incubated overnight at room temperature followed by the addition of 10 μL GPD detection reagent. After a second incubation period of 1 h at room temperature the luminescence signal was measured with an Envision spectrophotometer (PerkinElmer, Waltham, MA, USA). IC50 values were determined from the sigmoidal dose response curves with the software Quattro Workflow (Quattro GmbH, Munich, Germany).
To filter out compounds which interfere with the nucleotide exchange assay independent from K-RAS and SOS1, a control assay was developed. For this, GDP was titrated and detected with the GDP Glo bioluminescent GDP detection assay for glycosyltransferases (Promega, Madison, WI, USA) to generate a luminescence signal comparable to the positive control of the nucleotide exchange assay. Compounds were checked for assay interference by performing dose response curves in the same concentration range as for the nucleotide exchange assay. For each sample, 10 μL GDP (125 nM final assay concentration) in assay buffer (50 mM HEPES pH 7.5, 4 mM MgCl2, 2 mM EGTA, 0.01% Brij35, 1 mM TCEP) was transferred into a suitable assay plate (e.g. Greiner #784075). The test compound was added with an Echo acoustic dispenser (Beckman Coulter, Brea, CA, USA) in a concentration range from 3000 μM to 3 μM (8-point dilution). The reaction mixture was incubated over night at room temperature followed by the addition of 10 μL GPD detection reagent. After a second incubation period of 1 h at room temperature the luminescence signal was measured with an Envision spectrophotometer (PerkinElmer, Waltham, MA, USA). IC50 values were determined from the sigmoidal dose response curves with the software Quattro Work-flow (Quattro GmbH, Munich, Germany).
For every sample, 1.5 μL assay buffer and 3 μL of Tag2-SOS1 protein was transferred into a suitable assay plate (e.g. Corning #4513). Compound was added via Echo acoustic dispenser (BeckmanCoulter, Brea, CA, USA) in a suitable concentration range e.g. from 50 μM to 0.115 μM (8-point dilution). After an incubation step for 10 min at room temperature 3 μL Tag1-KRAS protein was added followed by 7.5 μL detection mix containing anti-Tag2 antibody labeled with Terbium cryptate and anti-Tag1 antibody labeled with XL665. After an incubation for 2 h at room temperature the HTRF signal was measured at 340 nm excitation, 665 nm and 615 nm emission with an Envision microplate reader (Revvity, Waltham, MA, USA). IC50 values were determined from the sigmoidal dose response curves with the Scigilian Analyze software (Scigilian, Monteral, Canada).
Footnote |
† Electronic supplementary information (ESI) available: Full experimental details including spectroscopic and MS data for all compounds are provided. See DOI: https://doi.org/10.1039/d4ra08503e |
This journal is © The Royal Society of Chemistry 2025 |