ORIGINAL RESEARCH article
Sec. Molecular Diagnostics and Therapeutics
Volume 8 - 2021 | https://doi.org/10.3389/fmolb.2021.625979
The Multi-Level Mechanism of Action of a Pan-Ras Inhibitor Explains its Antiproliferative Activity on Cetuximab-Resistant Cancer Cells
- 1Department of Biotechnology and Biosciences, University of Milan-Bicocca, Milan, Italy
- 2SYSBIO-ISBE-IT–Candidate National Node of Italy for ISBE, Research Infrastructure for Systems Biology Europe, Milan, Italy
- 3Bicocca Bioinformatics, Biostatistics and Bioimaging Centre - B4, Milano, Italy
- 4Department of Industrial Engineering and Innovation Sciences, Eindhoven University of Technology, Eindhoven, Netherlands
Ras oncoproteins play a crucial role in the onset, maintenance, and progression of the most common and deadly human cancers. Despite extensive research efforts, only a few mutant-specific Ras inhibitors have been reported. We show that cmp4–previously identified as a water-soluble Ras inhibitor– targets multiple steps in the activation and downstream signaling of different Ras mutants and isoforms. Binding of this pan-Ras inhibitor to an extended Switch II pocket on HRas and KRas proteins induces a conformational change that down-regulates intrinsic and GEF-mediated nucleotide dissociation and exchange and effector binding. A mathematical model of the Ras activation cycle predicts that the inhibitor severely reduces the proliferation of different Ras-driven cancer cells, effectively cooperating with Cetuximab to reduce proliferation even of Cetuximab-resistant cancer cell lines. Experimental data confirm the model prediction, indicating that the pan-Ras inhibitor is an appropriate candidate for medicinal chemistry efforts tailored at improving its currently unsatisfactory affinity.
Ras proteins are small guanine nucleotide-binding (G) proteins with low intrinsic GTPase activity, cycling between a GDP-bound inactive state and a GTP-bound active state. They act as molecular switches in signaling pathways regulating many cellular processes, including cell proliferation, growth, survival, adhesion, migration, energy, and redox homeostasis (Simanshu et al., 2017). Ras activity is regulated in response to specific extracellular stimuli, by the competitive action between Guanine nucleotide Exchange Factors (GEFs) promoting the nucleotide dissociation and GDP/GTP exchange, and GTPase Activating Proteins (GAPs), which provide an essential catalytic group for GTP hydrolysis (Scheffzek et al., 1997; Boriack-Sjodin et al., 1998; Bos et al., 2007). In human cells, three RAS genes encode four homologous but functionally distinct isoforms (HRas, NRas, and KRas4A and K-Ras4B) (Omerovic et al., 2007; Lu et al., 2016a). Gain-of-function missense mutations, mainly located at codons 12, 13, and 61, constitutively activate Ras proteins and can be detected in approximately one-third of all human cancers. Oncogenic Ras mutants contribute to tumor onset, maintenance, progression, and influence the efficacy of both cytotoxic and targeted therapies (Li et al., 2018). For this reason, many efforts, mostly promoted by the RAS initiative (https://www.cancer.gov/research/key-initiatives/ras), have been devoted to investigating the mechanistic role of RAS oncogenes in cancer and to explore different strategies for attenuating the aberrant Ras oncoproteins signaling, as widely reviewed (Sacco et al., 2012c; Welsch et al., 2017; Gorfe and Cho, 2021; Ni et al., 2019; Spencer-Smith and O’Bryan, 2019; Khan et al., 2020; Tisi et al., 2020).
Notably, each oncogenic mutation occurring in RAS genes induces conformational changes in the encoded protein that alter the residence time of the protein in the GTP-bound active state (Hunter et al., 2015) and make the oncoprotein surface more or less prone to the functional binding not only with modulators and effectors but also with specific pharmacophore groups or classes of molecule drugs. The RasG12V mutant presents a weak intrinsic and GAP-mediated GTP hydrolysis, and it is particularly aggressive and refractory to exchange inhibitors (Hunter et al., 2015). We first proved that the RasG13D mutant shows self-sufficiency in nucleotide dissociation (Palmioli et al., 2009b). Structural and functional studies (Smith et al., 2013; Hunter et al., 2015; Lu et al., 2016b; Johnson et al., 2019; Rabara et al., 2019) indicate that this mutant remains sensitive to the catalytic activity of GEFs and of at least one GAP, Nf1. Active and selective inhibitors for these oncogenic mutants are not yet available. On the contrary, compounds that covalently bind the highly reactive cysteine in the KRasG12C mutant selectively inhibit its function (Ostrem et al., 2013; Lito et al., 2016; Patricelli et al., 2016; Hansen et al., 2018a; Janes et al., 2018). After optimization for clinical use, they show a promising anti-tumor effect in KRASG12C-positive patients (Canon et al., 2019; Hallin et al., 2020).
We previously demonstrated that a class of small water-soluble molecules (cmp2-4), specifically binds the Switch II (β-3/α-2) region of wild type HRas-GDP. These compounds inhibit GEF-catalyzed nucleotide exchange, attenuate Ras signaling, and reduce Ras-dependent cell proliferation in mouse fibroblasts (Palmioli et al., 2009a; Sacco et al., 2011). Here we demonstrate that cmp4 binds an extended Switch II pocket on HRas and KRas proteins harboring different mutations. cmp4 decreases the intrinsic and GEF-mediated nucleotide dissociation and exchange on wild type and G13D mutated Ras proteins, interferes with Ras binding to GEFs (RasGRF1 and Sos1) and the Raf1 effector, and reduces mitogen-activated protein kinases signaling and cell viability of KRasG13D cancer cells. A mathematical model of Ras signaling (Stites et al., 2007; McFall et al., 2019), appropriately modified according to recent data (Johnson et al., 2017; Johnson et al., 2019), predicts the ability of cmp4 to inhibit the proliferation of different Ras-driven cancer cells. In keeping with the model prediction, experimental data on human cancer cell lines expressing different Ras oncoproteins confirm that cmp4 is a pan-Ras inhibitor able to cooperate with Cetuximab to inhibit proliferation of Cetuximab-resistant cell lines. Although cmp4 currently has an unsatisfactory affinity for Ras, targeted medicinal chemistry efforts could turn it into a valuable and needed clinical drug.
2 Materials and Methods
2.1 Compounds and Recombinant
cmp4 was synthesized as described (Palmioli et al., 2009a). Recombinant N-terminal His-tagged wild type and G13D mutated H-Ras proteins (residues 1-166 of the mature protein) and Sos1 catalytic domain (aa553-1024 of the mature protein) were expressed in M15 [pREP4] E. coli strain harboring a pQETM-derived plasmid (Qiagen) and purified by affinity chromatography using a Ni2+-NTA column (Qiagen), as described (Palmioli et al., 2009b; Palmioli et al., 2017; Sacco et al., 2012a). The N-terminal GST-tagged RasGRF1 catalytic domain (residues 976–1262 of the mature protein), was expressed in BL21 [pLysE] E.coli strain harbouring a pGEX2T-derived plasmid and purified by glutathione–sepharose chromatography (Amersham Bioscience) as described (Palmioli et al., 2017).
2.2 Mass Spectrometry Experiments
Mass-spectrometry measurements were performed on a hybrid quadrupole-Time-of-Flight (Q-TOF) instrument (QSTAR ELITE, Applied Biosystems, Foster City, CA, United States), equipped with a nano-ESI sample source. Metal-coated borosilicate capillaries (Proxeon, Odense, DK), with medium-length emitter tip of 1-mm internal diameter, were used to infuse the sample. The instrument was calibrated using the renine-inhibitor (1757.9 Da) (Applied Biosystems, Foster City, CA, United States) and its fragment (109.07 Da) as standards. Spectra were acquired in the 1500–3000 m/z range, with accumulation time of 1 s, ion-spray voltage of 1200–1500 V, declustering potential of 80 V, and instrument interface of 50°C. Spectra were averaged over a time period of at least 3 min. Data analysis was performed by the program Analyst QS 2.0 (Applied Biosystems, Foster City, CA, United States). The samples were prepared in 5 mM ammonium acetate pH 6.5.
2.3 NMR analysis
For the experiments with the free ligand, cmp4 was dissolved in a [D11]-Tris buffer at pH = 7.3, 5 mM MgCl2. COSY and HSQC experiments were performed by using the standard sequences. For the binding experiments, wild type or G13D mutated HRas was dissolved in 500 μL of the same [D11]-Tris buffer, containing an amount of GDP equimolar to the protein, and transferred into a 5 mm NMR tube; 50 μL of the ligand solution dissolved in the same buffer were added slowly. Final protein concentration was 50 µM, final ligand concentration was 1 mM.
STD experiments were performed without saturation of the residual HDO signal and with spin-lock to avoid the presence of protein resonances in the spectra. A train of Gaussian-shaped pulses of 50 ms each was employed, with a total saturation time of the protein envelope of 2 s. An off-resonance frequency of δ = 40 ppm and on-resonance frequency δ = −1.5 ppm (protein aliphatic signals region) were applied. Spectra were acquired with a Varian Mercury 400 MHz instrument and processed using the program Mestre-Nova 9.
2.4 Flexible Docking Algorithm
Docking analyses were performed in Maestro 10.1 suite (Schrӧdinger) (https://www.schrodinger.com/citations#Maestro). All docking calculations were performed using the Glide software (Glide, version 6.7, Schrödinger, LLC, New York, NY, 2015). The receptor-based molecular docking was carried forward after preparing ligands and proteins as suggested by the developer’s protocols. For HRas and KRas, the pockets corresponding to the residues identified by experimental data on HRas were used as the input for grid receptor definition in induced-fit docking (IFD) workflow with flexible ligand option. The protocol generates alternative cmp4 poses not considering clashes with amino acids side-chains, then optimize the structures obtained by allowing the protein to undergo sidechain or backbone movements during the process (Schrödinger Suite 2015-2 Induced Fit Docking protocol; Glide version 6.7, Schrödinger, LLC, New York, NY, 2015; Prime version 4.0, Schrödinger, LLC, New York, NY, 2015). The IFD extended sampling protocol was employed, generating up to 20 poses per ligand on each iteration. The OPLS 2005 force field (Jorgensen et al., 1996) was used for the minimisation stage, in which residues within 10 Å of each ligand pose were optimised. All other parameters were set to their default values. GLIDE molecular docking output GScore (empirical scoring function) is reported, which is calculated by calculating ligand–protein interaction energies, root mean square deviation (RMSD), hydrogen bonds, hydrophobic interactions, internal energy, π–π stacking interactions, and desolvation. GLIDE Emodel was used to choose the best pose for the ligand in each structure, while IFD Score is based on the Prime calculation of energy content of the structure, and also considers the strain in the receptor and ligands.
2.5 Dissociation and Exchange Reactions
Intrinsic and GEF-mediated dissociation and exchange of mant-guanine nucleotides (mant-GXP, GXP being GDP or GTP; Molecular Probes; Invitrogen) assays were performed essentially as described in (Lenzen et al., 1995; Sacco et al., 2012b). Briefly, for dissociation reactions HRas protein was pre-loaded with mant-GXP by incubating for 30 min 250 µM HRas with 750 µM mant-GXP in 40 mM Hepes pH 7.5, 1 mM MgCl2, 2.5 mM DTE, 20 mM EDTA. Then 30 mM MgCl2 was added and the solution was incubated for further 30 min. Free nucleotides were removed by gel filtration using PD10 desalting columns (Amersham Bioscience) equilibrated with Lenzen buffer (40 mM Hepes, pH 7.5, 5 mM DTE, 10 mM MgCl2), and HRas-mant-GXP complex was concentered using centricon 10 KDa (Merck Millipore). The exchange reactions on Ras protein were performed by adding directly in an UV-cuvette 0.25 μM HRas-GXP, and an opportune concentration of cmp4 in Lenzen buffer. After 300 seconds of incubation, a 5-fold excess of mant-GXP (1.25 μM) and a specific concentration of the exchange factor (0 or 0.0625 µM as indicated) were added. The fluorescence measurements were carried out at 25°C using a LS45 fluorescence spectrometer (Perkin-Elmer) with an excitation wavelength of 366 nm and emission wavelength of 442 nm. The reactions were monitored for at least 1500 s. The dissociation reactions were performed in a UV-cuvette by adding to 0.25 μM HRas-mant-GXP, preincubated for 300 s with the opportune concentration of cmp4, 200 μM GXP and a specific concentration of the exchange factor (0, 0.0125, 0.025, 0.0416, 0.125 −0.25 µM), as indicated. Exchange data were fitted to a nonlinear “growth-sigmoidal Hill” curve (n = 1), while dissociation data were fitted to an “Exponential decay” curve, using the OriginPro 8.0 software (OriginLab Corporation, MA United States). The initial exchange or dissociation rate for each reaction (initial slope) was determined by computing the first derivative at time zero of the corresponding fitted curves. In the graphs, the maximum value of relative fluorescence (100 on Y-axis) represents the fully loaded Ras status obtained as a start point in dissociation reaction and plateau of an exchange curve obtained in the absence of cmp4.
To measure the affinity for entering nucleotide, a plate-based GDP/GTP titration assay was adapted from the method previously described (Ostrem et al., 2013): 1 µM HRas-mant-GDP complex was added to 96-well black plates in 40 mM Hepes, pH 7.5, 5 mM DTE, 1 mM MgCl2. The fluorescence was measured on a Variant Cary Eclipse fluorescence spectrometer (Agilent), with 360 nm excitation and 440 nm emission, before and after 2 h incubation at 25°C with a 5 mM EDTA solution with different concentrations of GDP or GTP, as indicated.
Results for each nucleotide were fitted to a sigmoidal curve using the OriginPro 8.0 software.
2.6 Surface Plasmon Resonance Analysis and G-LISA
Surface Plasmon Resonance experiments were carried out by using a BIAcoreX system (BIAcore, GEHealthcare). His-tagged HRas-GDP was immobilized onto a NTA-sensor chip surface (carboxymethylated dextran matrix pre-immobilized with NTA; BIAcore, GEHealthcare), obtaining a surface density of about 4500 resonance units. No nickel solution was injected over the reference cell. The binding with GST-fused RasGRF1 was monitored in real time in the presence of increasing concentrations of cmp4 (0-500 µM). All experiments were performed in HBS-P+ buffer (BIAcore, GE Healthcare) at a flow rate of 10 µL/min. Surface regeneration was accomplished by injecting EDTA (350 mM) in the flowing buffer (30 s contact) two or three times. The evaluation of binding kinetics was performed by using the Biaevaluation software, v. 3.0 (BIAore) and by considering a 1:1 Langmuir interaction. Notably, the value of koff measured in the SPR experiments cannot correspond to the physiological dissociation constants because the absence of free nucleotide in the experiments substantially affects this parameter.
Ras G-LISA Activation assay kit (Cytoskeleton, Inc. BK131) was used to measure the levels of HRas-GTP bound to the Ras binding domain of Raf1 (RBD-Raf1) in the presence of increasing concentrations of cmp4 (range 0.08-500 μM). HRas-GTP 0.4 nM was preincubated in batch with cmp4 for 5 min at RT and then transferred in 96-well coated with RBD-Raf1. After incubation at 4°C for 30 min, the plate was washed three times with washing buffer before the addition of antigen-presenting buffer. The captured HRas-GTP was incubated with the anti-Ras antibody followed by HRP-conjugated secondary antibody. Ras activity was quantified by measuring absorbance at 490 nm.
2.7 Cell Lines and Proliferation Assay
Human breast cancer cell line MDA-MB-231, obtained from the American Type Culture Collection, was routinely grown at 37°C in a humidified atmosphere of 5% CO2 in Dulbecco’s modified Eagle's medium (D-MEM) (Sigma D6429) supplemented with 10% Newborn Calf Serum (NCS), 2 mM glutamine, 100 units/ml penicillin and 100 mg/ml streptomycin. Human colon adenocarcinoma cell line SW48 (KRASWT/WT) and the isogenic SW48 expressing heterozygous KRasG13D (KRASWT/G13D) or KRasG12V (KRASWT/G12V) were obtained from Horizon Discovery Ltd. Cells were cultured in humidified atmosphere of 5% CO2 at 37°C in RPMI 1640 (Sigma R0883) supplemented with 10% Fetal Bovine Serum, 2 mM glutamine, 100 units/ml penicillin and 100 mg/ml streptomycin. Cells were passaged using trypsin–EDTA.
For growth kinetics and RealTime-GloTM MT Cell Viability Assay (Promega, #G9713) cells were plated into respectively 6-well or 96-well flat-bottomed culture plates at the density of 3000 cells/cm2. At 18h after seeding, predetermined concentrations of cmp4 (or water) were added to the cell culture. After 24, 48, and 72 h from treatment, cells were harvested and counted by Coulter Counter to obtain growth curves or treated with 500 X NanoLuciferase and 500 X MT cell viability substrate. The luminescence at different time points after treatment was recorded by using a Victor Multilabel Plate Reader (Perkin Elmer). The viability of cells treated with increasing concentrations of cmp4 was tested relative to the viability of the same cells treated with vehicle (water). Viability results were analyzed by using OriginPro 8.0 software and a nonlinear growth/sigmoidal Hill curve (n = 1) to calculate the relative IC50 values.
2.8 MAPK Activity
Breast cancer MDA-MB231 were plated (6000 cells/cm2) in 60-mm tissue culture dishes. After 18 h different concentrations of cmp4 (or vehicle) were added to the cell culture. After 48 h from treatment, both plate-adherent and in suspension cells were harvested in lysis buffer from PathSscan Sandwich ELISA kit (Cell Signaling). The detection of endogenous levels of Phospho-p44/42 MAPK was performed according to manufacturer's instructions, and the results were normalized on total protein content measured by Bradford analysis.
2.9 Mathematical Model
The computational analysis was performed starting from the mechanistic model presented in McFall et al. (2019), where a system of Ordinary Differential Equations (ODEs) is introduced to describe the Ras signaling network. The system of ODEs corresponds to the reactions reported in Supplementary Table S1, where reactions R1-R8 follow the mass-action kinetics, while reactions R9, R10, R11 follow the Michaelis-Menten kinetics; reaction R9 describes the GAP activity, while reactions R10 and R11 describe the GEF activity. The 11 reactions can be used to simulate both the wild type and mutant proteins by assuming different values of the kinetic parameters. In particular, the kinetic parameters of RasG13D and RasG12V mutants were obtained by scaling the wild type parameters (4th column) according to the corresponding alpha factors reported in the 5th and 6th columns of Supplementary Table S2. The scaling factors of RasG13D, related to reaction R9, were modified according to the results presented in (Johnson et al., 2019).
Specifically, the computational investigation presented in this work was performed with COPASI (Hoops et al., 2006) (version 4.27), exploiting the LSODA numerical integrator (Petzold, 1983). LSODA is an efficient simulation algorithm capable of dealing with stiff systems by automatically switching between explicit (the Adams’ method) and implicit integration methods (backward differentiation formulae). The accuracy in the description of the solution of the system of ODEs is controlled by the relative tolerance, that is the maximum error allowed in the solution, and absolute tolerance, which is the maximum error allowed in case the solution approaches zero. In the simulations performed here, we considered the following setting: relative tolerance 1e-6, absolute tolerance 1e-12, maximum number of steps executed to generate the solution, at each iteration, 1e5. COPASI was also exploited to perform a parameter sweep analysis (PSA) to investigate the effect of the parameter variations on the emergent dynamics and on the steady-state values of pivotal components of the model. The simulations concerning the PSA have been run by generating a set of different initial conditions for the model, considering a fixed range of variation of the parameter under investigation, and then executing the corresponding simulations with LSODA. In particular, the PSA was performed by varying a single kinetic parameter, considering a logarithmic sampling of values within the specified range. The responsiveness of the RasG13D mutant variant to GEF activity was analyzed by performing a PSA where the Vmax of GEF-mediated exchange reactions (R10 and R11) was multiplied for a parameter gamma, which varied in the range 0-1, where the top value represents the maximal activation of the GEF and the lower value represents the loss of GEF function. The basal level of unstimulated GEF activity is set as corresponding to a gamma value of 0.1.
The effect of different concentrations of Cetuximab and cmp4 was simulated by perturbing the reference parameterization (4th column of Supplementary Table S2) of the model as reported in Supplementary Table S1. In detail, the maximal action that could be obtained by an inhibitor acting by rescuing EGFR hyperactivation was simulated by dividing KM,10 and KM,11 by 10. The effectiveness of Cetuximab-like inhibitors was analyzed by a PSA performed by multiplying the Vmax of GEF-mediated exchange reactions (R10 and R11) for a parameter gamma. This parameter was varied in the range 0-1, where the absence of EGFR stimulation is represented by a 0.1 value. cmp4 (at 100 µM) expected effect was simulated by multiplying KM,10 and KM,11 by 0.5, k2-5 by 0.5, and k6 by 0.23 (yielding a half amount of Ras-GTP-Eff complex formation).
3.1 cmp4 Binds to both GDP-Bound Wild Type and G13D Mutated HRas Proteins
By NMR analysis, we previously showed that cmp4 (Figure 1A) binds HRas-GDP (Sacco et al., 2011) and –mainly through its aromatic moiety–in a binding pocket located between the α2-(Switch II) and α3-helices (Figure 1B). Flexible docking indicates that cmp4 binds to an extended Switch II pocket (here referred to as SII-EP) of HRas and KRas (Figures 1C,D; Supplementarys S1A,B and Table S3). This pocket partially overlaps with the Switch II groove (SII-G) identified on KRas by structural analysis in Gentile et al. (2017) (Supplementary Figures S2A,D,G). Results of an STD-NMR analysis of HRas-GDP with cmp4 (Figure 1B) and additional data collected on similar compounds (Palmioli et al., 2009b; Palmioli et al., 2009a; Palmioli et al., 2017; Colombo et al., 2010; Sacco et al., 2011) support the pivotal role of the benzyl group and the pyrocatechol group for Ras binding. We used these results to filter the top 10 poses in this and other docking experiments.
FIGURE 1. Reconstruction of a binding pose for cmp4 on wild-type and G13D mutated HRas-GDP complexes. (A) Chemical structure of cmp4; (B) NMR analysis traces report including aromatic resonance (region I), sugar resonance, and aliphatic CH2 (region II): 1)1H NMR spectrum of 1 mM cmp4; 2) 1H NMR spectrum of a sample containing 1 mM cmp4 and 50 μM HRas-GDP wt; 3) STD-NMR spectrum of a sample containing 1 mM cmp4 and 50 μM HRas-GDP wt; 4) 1H NMR spectrum of a sample containing 1 mM cmp4 and 50 μM HRasG13D-GDP; 5) STD-NMR spectrum of a sample containing 1 mM cmp4 and 50 μM HRasG13D-GDP. (C) Docking pose of cmp4 on PDB structure of HRas-GDP (PDB ID: 4q21). The image shows switch I (red), switch II (blue) and displays in pink the residues of Ras that undergo significant chemical shift perturbations after binding with cmp4 (Sacco et al., 2011); (D,E) Molecular detail of the selected pose of cmp4 on: (D) HRas-GDP (PDB ID: 4q21); (E) HRasG13DGDP (PDB ID:6dzh). Ras residues that are directly involved in binding with cmp4 are indicated. The backbone in the Switch I region is colored in red while the backbone in the Switch II region is colored in blue. The GDP nucleotide (grey) and cmp4 (yellow) are drawn in sticks. Heteroatoms are in red (oxygen) and blue (nitrogen).
cmp4 is a much bulkier molecule than the compound reported in Gentile et al. (2017) and occupies a larger pocket than the one there described (Supplementary Figures S2B,E,H), protruding towards the Gly12 P-loop. The cmp4 pyrocatechol group, as obtained in all of the docking best scoring poses, is much farther from this loop than the G12C binding compounds first described to target an allosteric switch II pocket (Ostrem et al., 2013; Patricelli et al., 2016) (Supplementary Figures S2C,F,I). Notably, catechol interacts with residues not only in α2-(switch II) (Glu62, Tyr64, Arg68) and α3-helices (Tyr96, Arg102) but also with the backbone of Gly10 in the P-loop (see the ligand interactions plot in Supplementary Figure S1A).
STD analysis on HRasG13D mutant protein saturated with cmp4 shows that cmp4 also interacts with the mutant protein. Flexible docking indicates that cmp4 maintains a similar positioning within the binding pocket of HRasG13D-GDP, or in KRasG13D-GDP as well, despite the partial switch II unfolding observed in the oncoprotein (Figures 1E,D; Supplementary Figures S1C,D). The top docking scores were slightly lower than obtained on the wild type proteins (Supplementary Table S3).
Since the pathological effect of Ras hyperactivity is due to the active, GTP-bound form, and the phenol-derived compounds occupying the Switch II groove (SII-G), identified by structural analysis in Gentile et al. (Supplementary Figure S1A), were reported to target Ras active form as well (Gentile et al., 2017), we also assessed whether the SII-EP pocket in GTP-bound HRas and HRasG13D. is available to cmp4 interaction. Due to the different conformation of Switch II, this pocket seems to be less available in the GTP-bound complex (Supplementary Figure S3), leading to a maximal docking score decreased in comparison to that observed in the GDP-bound form (Supplementary Table S3), but still consistent with data previously reported for compounds binding to analogous pockets (Ostrem et al., 2013; Lito et al., 2016; Patricelli et al., 2016).
3.2 cmp4 Inhibits the Intrinsic and GEF Mediated-Nucleotide Dissociation and Exchange on Wild Type and G13D Mutated HRas in a Dose-Dependent Manner
mant-GDP is a nucleotide analog whose fluorescence increases upon Ras binding. The decrease in fluorescence following incubation of the Ras-mant-GDP complex with an excess of unlabeled GTP allows us to follow nucleotide exit (Figure 2A, left) The increase in fluorescence obtained after incubation of a Ras-GDP complex with an excess of mant-GDP directly monitors nucleotide entry (Figure 2A, right). In the normal Ras activation cycle, the entry of a new nucleotide immediately follows the nucleotide exit.
FIGURE 2. cmp4 counteracts nucleotide dissociation from Ras complex. (A) Scheme depicting the different experiments for nucleotide dissociation or exchange, using nucleotides (GDP or GTP, namely GXP) conjugated with the fluorescent moiety MANT; (B,C) inhibitory efficacy of cmp4 on GEF-mediated nucleotide dissociation (B,C) exchange on HRas (black) and HRasG13D (red). (D) Mass spectrometry analysis of HRas-GDP in presence of cmp4. The dashed peaks correspond to the expected positions of the nucleotide free-Ras and Ras complexed with cmp4 without GDP. (E,F) Inhibitory efficacy of cmp4 on intrinsic nucleotide dissociation (E) and exchange (F) on HRasG13D; the initial dissociation or exchange rate of each reaction was determined computing the first derivative at time 0 of the fitted curves reported in Supplementary Figure S4. (G,H) EDTA-mediated competition between mant-GDP loaded on H-Ras and free unlabelled GDP (G) or GTP (H).
We previously demonstrated that cmp4 interferes with the function of the exchange factor RasGRF1 on HRas (Sacco et al., 2011). Here we show that cmp4 inhibits the GEF-catalyzed nucleotide dissociation and exchange reaction on wild type and G13D mutated HRas with similar efficiency (Figures 2B,C; Supplementary Table S4). Supplementary Figures S4A–D show the actual dissociation and exchange curves. We used the initial rates of each reaction (mean of at least three independent experiments) for calculating the IC50 reported in Figures 2B,C and Supplementary Table S4. The inhibitory effect of cmp4 on both dissociation and exchange reactions is independent of the GEF hSos1 vs. RasGRF1, (Sacco et al., 2011) and of the entering nucleotide, GDP or GTP (Supplementary Figure S5).
The docking results presented in Figure 1 suggest that cmp4 may form a stable Ras-nucleotide-cmp4 ternary complex, without promoting dissociation of the Ras-bound nucleotide, similar to peptide Ras inhibitors developed in our laboratory (Sacco et al., 2012b). The deconvoluted mass spectrum of 10 µM HRas-GDP in the presence of a 10-fold excess of cmp4 (Figure 2D) shows no signal corresponding to the nucleotide-free Ras/cmp4 complex. Except for a minor fraction of Ras-GDP binding a second inhibitor molecule at a low affinity, non-specific site, the HRas-GDP-cmp4 ternary complex is the most abundant species.
It was therefore of interest to monitor whether cmp4 can inhibit intrinsic (i.e., non GEF-catalyzed) nucleotide dissociation and exchange. We first tested the effect of cmp4 on HRasG13D, whose intrinsic nucleotide exchange rate is much higher than that of wild-type HRas (Palmioli et al., 2009b; Smith et al., 2013; Hunter et al., 2015; Johnson et al., 2019). Supplementary Figures S4E–H reports the actual dissociation and exchange curves. The inhibitor efficiently reduces the abnormally fast intrinsic nucleotide dissociation and exchange reactions on HRasG13D in a dose-dependent manner (Figures 2E,F). The inhibitory effect is also appreciable on the intrinsic activities of wild type HRas, which are very slow per se (see the inserts in Supplementary Figures S4E,G).
Titration with unlabeled GDP and GTP of a HRas-mant-GDP complex in the presence of EDTA allows monitoring whether a drug alters the affinity for the entering nucleotide. Figures 2G,H indicate that cmp4 alters the entry of both nucleotides without discriminating between GDP and GTP, unlike the SII-P binding molecules described by Ostrem et al. (2013).
These results suggest that cmp4 binding to the Switch II extended pocket (SII-EP) counteracts nucleotide release, even in conditions favoring nucleotide release, as observed in HRasG13D (Johnson et al., 2019), and/or in the presence of EDTA or a GEF catalytic domain.
3.3 cmp4 Reduces the Affinity of HRas-GDP for RasGRF1 and Raf1 Ras Binding Domain in a Dose-Dependent Manner
The inhibitory efficiency of cmp4 on the nucleotide dissociation rate on both wild type and G13D mutated HRas decreases with increasing RasGRF1 concentration (Figure 3A), suggesting that the GEF could force the nucleotide dissociation even on cmp4-bound Ras, counteracting the inhibitor action. In order to bind the GEF catalytic domain with the highest affinity, HRas has to undergo a conformational change that allows nucleotide release (Boriack-Sjodin et al., 1998), as evidenced by the superposition of HRas structures respectively in GDP-bound and nucleotide-free Sos1cat-bound form (Figure 3B). Notably, the same kind of interaction is also envisioned for the catalytic domain of RasGRF1, due to homology with Sos1 (Freedman et al., 2006).
FIGURE 3. cmp4 affects HRas binding to GEF (RasGRF1) and effector (Raf1-RBD) in a dose-dependent manner. (A) Inhibition of nucleotide dissociation rate on both HRas and HRasG13D (1 μM) in the presence of 100 μM cmp4 and increasing concentrations of RasGRF1 (range 0.01–0.25 μM). (B) Best fitting pose of cmp4 on HRas-GDP (pink) was superimposed to the structure of nucleotide-free Ras (HRasnf, in grey) from the crystal structure of the hSos1 catalytic domain associated with HRas (PDB ID: 1bkd). Switch I and II regions are stained darker. GDP is in pink, cmp4 in yellow; (C) Biacore-based direct measurement of 0.5 μM GEF (GST-RasGRF1) binding to His-HRas-GTP in the presence of increasing concentrations of cmp4 (25–1000 μM). In the insert kinetics analysis of RasGRF1 binding to HRas-GDP in the presence of different concentrations of cmp4, relative to SPR curves. All points for initial association rate (von, closed symbols, voff, open symbols) were fitted respectively to a nonlinear ‘growth-sigmoidal Hill’ curve (n = 1), which is reported in the graph as a thin line; (D) Levels of HRas-GTP bound to a Ras binding domain (RBD) of Raf1 in the presence of increasing concentrations of cmp4 (range 0–500 μM), detected with the G-LISA® kit (Cytoskeleton, Inc. BK131). Data were normalized to Ras-GTP levels measured in the absence of cmp4 (control). All data are significant at 99%, as calculated by Student’s t-test in comparison to control. In the inset, the percentage of inhibition of Ras-GTP bound to RBD as a function of cmp4 concentration, relative to the G-LISA experiment. All points were fitted respectively to a nonlinear ‘growth-sigmoidal Hill’ curve (n = 1).
SPR binding experiments analyzed the interaction between Ras and GEF in the presence of increasing concentrations of cmp4 (Figure 3C). cmp4 affects GEF (RasGRF1) binding to Ras-GDP in a dose-dependent manner, with an estimated EC50 of 170 μM. In particular, cmp4 dose-dependently reduces the association rate, and so the kon of the interaction (Figure 3C, inset), suggesting that the compound reduces the formation of the Ras/GEF complex, a key intermediate in Ras activation cycle. This finding agrees with the observation that cmp4 stabilizes the nucleotide-bound HRas conformation by bridging Switch I and Switch II (Figure 1D). This stabilized connection between Switch I and II would make Ras more refractory to the formation of the high-affinity complex with the GEF and to its catalytic action (Boriack-Sjodin et al., 1998).
The aberrant mitogenic signaling in Ras-driven cancer cells largely depends on the increased recruitment of the downstream effectors Raf1, from the constitutively active Ras oncoproteins (Metcalfe et al., 1993; Warne et al., 1993). Accordingly, molecules disrupting Ras/Raf1 association block KRas downstream signaling and impair Ras-mediated tumorigenic proliferation (Waldmann et al., 2004; Athuluri-Divakar et al., 2016; Trinh et al., 2016; Liu et al., 2017; McGee et al., 2018; Wiechmann et al., 2019). To assess the ability of cmp4 to affect Ras-GTP/Raf1 binding, we performed an ELISA assay with increasing concentrations of cmp4 (from 0 to 500 µM). Figure 3D shows that cmp4 reduces in a dose-response manner the amount of Ras-GTP complex bound to the effector Ras binding domain of Raf1 (RBD-Raf1) with an EC50 value of about 0.45 μM (IC50 about 250 μM, Figure 3D, insert).
3.4 cmp4 Reduces Cell Proliferation and MAPK Activation in KRasG13D Expressing Cancer Cells
KRas–the predominantly Ras isoform mutated in cancer–presents a different amino acid in front of the binding pocket (glutamine instead of histidine in position 95) and a more disordered Switch II region even in the active conformation (Johnson et al., 2019) when compared to HRas. Docking poses and their scores (Supplementarys Figures S1B,D; Table S3) suggest that the pocket in KRas and KRasG13D is equally available for cmp4 binding, consistently with the inhibitory effect exerted by cmp4 on KRas-transformed mouse fibroblasts (Sacco et al., 2011).
Here we evaluated the effect of cmp4 on MDA-MB-231, human breast cancer cells expressing KRasG13D. cmp4 reduces the proliferation of MDA-MB-231 cells in a dose-dependent manner (Figure 4A, IC50 of about 125 µM at 72 h), causing a significant cell detachment (see microscopy images in Figure 4B). MTT assays (Figure 4C) show that cmp4 significantly affects the viability of MDA-MB-231 cells already after 24 h-treatment. The reduced proliferative potential of cells treated with cmp4 correlates with a dose-dependent decrease of the level of activated/phosphorylated mitogen-activated protein kinases (MAPKs), as revealed by an ELISA assay performed on cell lysates collected after 24-h treatment with cmp4 (Figure 4D). Since high doses of cmp4 were administered to the cells, due to its low Ras affinity, we cannot exclude that the inhibition of proliferation is ascribable to off-target effects. However, the correlation between MAPK and cellular proliferation is consistent with a predominant specific effect on Ras activity.
FIGURE 4. Effect of cmp4 on cell viability and Ras signaling of human breast cancer MDA-MB-231 cell line expressing KRasG13D. (A) Growth curves of MDA-MB-231 treated with increasing concentration of cmp4 or vehicle (deionized water) supplemented in the growth medium. After 24, 48, and 72 h of treatment cells were trypsinized and counted with a Burker chamber. (B) Microscopy analysis of MDA-MB-231 treated for 48 h with different concentrations of cmp4. (C) Cell viability of MDA-MB-231 cells treated with increasing concentration of cmp4, or vehicle (deionized water) for 24, 48, and 72 h as measured by MTT assay; data were normalized on cells treated with vehicle imposed as equal to 1. (D) Phosphorylated MAPK level in cell lysates from MDA-MB-231 cells no treated or 24 h-treated with cmp4 or vehicle. Data were normalized on the phospho-MAPK level in MDA-MB-231 treated with vehicle imposed as equal to 1. Data shown are mean and standard deviation of two independent experiments, each performed in triplicate. Single and double asterisk above histograms indicates a statistical significance of 95% and 99% respectively, calculated by Student’s t-test in comparison to cells treated with vehicle.
3.5 Validation of the Mechanism of Action of cmp4 in Isogenic Cell Lines Expressing Different KRas Oncoproteins
Different Ras mutants produce a spectrum of distinct phenotypic effects and may display a significant difference in their ability to respond to therapies targeting the Ras pathway (Johnson et al., 2019). A recent computational model of the Ras activation cycle allows to explain and reproduce some of these different phenotypic traits, such as the peculiar sensitivity of KRAS mutants to Cetuximab, a drug targeting EGFR hyperactivation (McFall et al., 2019). Since results presented above and literature data (Sacco et al., 2011) indicate that cmp4 may interfere with multiple steps of the Ras activation cycle, we decided to use this model together with experiments on isogenic cell lines expressing different Ras mutant proteins to validate the mechanism of action of cmp4.
The model of the Ras activation cycle (Figure 5A) consists of 11 reactions (Supplementary Table S1). The first 8 reactions follow the mass-action kinetics, with a single kinetic parameter, while reactions R9, R10, R11 follow the Michaelis-Menten kinetics and require two different parameters. Parameter values can be changed to tailor the model to different cell systems. Supplementary Table S2 reports parameters used in this paper, that have been partially modified compared to McFall et al. (2019), by taking into account recent literature (Johnson et al., 2019; Rabara et al., 2019) and our own data. GEF activation induced by the interaction of a Growth Factor with its cognate receptor (reaction not included in the model) is simulated by an abrupt increase (up to 10-fold) of the Vmax of the GEF-catalyzed reactions, i.e., Vmax,10 and Vmax,11 (grey arrow pointing to GEF in Figure 5A). Figure 5B (left panel) reports the results of a simulation of virtual cells in the absence of growth factor stimulation. Starting from nucleotide-free Ras, a rapid association of Ras with the available nucleotides is observed (guided by the fast reactions R4 and R5), then the level of Ras-GTP (grey line) and of the Ras-GTP-effector complex (dotted line) reach a steady state over the course of the simulation, characterized by a low level for both the species. When Vmax,10 and Vmax,11 are increased (simulating growth factor stimulation, Figure 5B right panel), both Ras-GTP and the Ras-GTP-Effector complex reach a steady-state level that is higher than the basal level observed in the absence of growth factor. In the following, we will use the steady-state level of the Ras-GTP-Effector complex to estimate the proliferation state of the simulated cell lines and to compare simulated and experimental data.
FIGURE 5. In silico modeling of Ras signalling network. (A) Scheme depicting the species and the reactions constituting the mathematical model. The parameters affected by cmp4 are indicated by red arrows. (B) Example of the output obtained simulating growth factor unstimulated (left) and stimulated (right) conditions of a wild type system. The level of effector activity, that is the level of the complex Ras-GTP-eff, is evaluated as the steady-state value reached during the dynamic simulation.
The small red arrows in Figure 5A indicate the steps within the Ras activation cycle affected by cmp4. They include the reactions describing the intrinsic association to, and dissociation from, the nucleotide (R2-R5), reactions describing association to the effector (R6), and GEF-mediated reactions allowing nucleotide exchange (R10 and R11). To study the effect of cmp4 on the Ras activation cycle we instantiated three different models representing a cell line endowed with a constitutively active EGFR mutant (EGFRG719S). This mutant receptor constitutively recruits GEFs to the plasma membrane causing an aberrant Ras activation. We simulated this mutation by imposing the maximal value for Vmax,10 and Vmax,11. The wild type cell line carries two wild type KRAS alleles, while two mutant cell lines express KRasG13D and KRasG12V in heterozygosis. Simulation of these virtual cell lines shows that the KRASWT/G12V heterozygous mutant is the most aggressive based on the level of total KRas-effector complex, followed by the wild type and by the KRASWT/G13D (Figures 6A–C). Although surprising at first sight, this result likely reflects the lower affinity of the KRasG13D mutant protein for Raf1 (Johnson et al., 2019).
FIGURE 6. Simulation of the effect of Cetuximab and cmp4 on Ras signaling in human colorectal cancer SW48 isogenic cellular models. (A–I)In silico modeling of Ras signaling network in SW48 isogenic cellular models, expressing a hyper-activated EGF receptor mutant in combination with different Ras variants: either wild type Ras (A,D,G; SW48 KRASWT/WT) or KRasG13D (B,E,H; SW48 KRASWT/G13D) or KRasG12V (C,F,I; SW48 KRASWT/G12V). The different cellular systems were simulated under untreated condition (A–C; CTRL), or treated with the following drugs: Cetuximab (D–F), used at an ideal concentration completely blocking EGFR activity, which represent the maximal effect obtainable with the single mechanism of action based on GEF-mediated nucleotide exchange inhibition; cmp4 (G–I), used at the concentration of 100 µM, which is around IC50 for this compound on multilevel mechanisms of action (Supplementary Table S1). For each panel, the dimension and colour of the characters are indicative of the level of the components or their activity in the simulation. The resulting effector activity is illustrated as a histogram on the right of each panel, and its fold change normalized on wild type unstimulated cells is reported on top. For panels (D–I), the inhibition efficacy is calculated with respect to the untreated corresponding model.
As confirmed by our results (Figure 7A), the presence of the GAP-insensitive KRasG12V mutant confers resistance to the treatment with Cetuximab (Burgess et al., 2003; Seshacharyulu et al., 2012). Computational results predict that the theoretical maximal effect exerted by Cetuximab (i.e., a complete reversion of GEFs activation) leads to a reduction of virtual proliferation (i.e., a reduction in the level of the Ras-GTP-effector complex) of 95% in the SW48 KRASWT/WT model, of 87% in the SW48 KRASWT/G13D model and only of 20% in the SW48 KRASWT/G12V model (Figures 6D–F). These simulation results are consistent with RasG13D being responsive to GEFs action (Palmioli et al., 2009b; Smith et al., 2013; Johnson et al., 2019) and Supplementary Figure S4, whereas KRasG12V is fully active even if GEFs are not activated (Supplementary Figure S6).
FIGURE 7. Effect of 72 h-treatment with cmp4 and/or Cetuximab (CTX) on cell viability of human colorectal cancer SW48 isogenic cell lines. (A–C) Relative cell viability of SW48 KRASWT/WT, SW48 KRASWT/G12V, and SW48 KRASWT/G13D cells treated for 72 h with different concentrations of CTX (A) or cmp4 (B). (C) Results from the simulations of the SW48 KRASWT/WT, SW48 KRASWT/G12V and SW48 KRASWT/G13D mathematical models either untreated or treated with the following drugs: 0.5 nM CTX (corresponding to an inhibition of nearly 70% of GEF activity); 100 µM cmp4, which is around IC50 for this compound on multilevel mechanisms of action; a combination of both. (D) Relative cell viability of SW48 KRASWT/WT, SW48 KRASWT/G12V, and SW48 KRASWT/G13D cells treated for 72 h either with 0.5 nM CTX, 100 μM cmp4, or a combination of both. Data were normalized on cells treated with vehicle taken equal to 1. Single, double, and triple asterisk above histograms in (A,B,D) indicates a statistical significance of 95%, 99%, and 99.9% respectively, calculated by Student’s t-test in comparison to cells treated with vehicle.
The potential inhibitory effect of cmp4 was tested on all the models, in the hypothesis that it could behave as a panRas inhibitor. The appropriate constants (Figure 5A) were modified with respect to the untreated case, by considering the biochemical effect induced by treatment with 100 µM cmp4 in the appropriate in vitro assay (see Supplementary Table S1 for actual values used in simulation experiments). Both experimental cell viability assays and simulation results indicate that all three virtual cell lines are sensitive to cmp4 (Figures 6G–I, 7B), the SW48 KRASWT/G12V cell line being the less sensitive (Supplementary Table S5).
These results prompted us to test whether the combined use of both drugs could improve the pharmacological treatment of the G12V mutant. Simulation results indicate that the combined treatment is additive or nearly additive in the three cell line models SW48. The effect is striking in the KRASWT/WT and KRASWT/G13D models (Figure 7C, black and grey bars, respectively), but nevertheless noticeable also in the KRASWT/G12V model, where complete inhibition of the EGFR cascade (i.e., leaving Vmax,10 and Vmax,11 at their basal level) has only a 10% effect on the level of the Ras-GTP-effector complex (Figure 7C, white bars). We fully confirmed these simulation results by measuring the inhibition in cell proliferation of the three cell lines treated with a combination of the two drugs (Supplementary Figure S7) and in particular with 0.5 nM Cetuximab (CTX), 100 μM, cmp4, or a combination of the two drugs (Figure 7D), validating the multi-level mechanism of action of cmp4 suggested by the molecular docking and biochemical assays described above.
Reported success in the direct targeting of Ras proteins, long postulated as undruggable, has paved the way to the possible pharmacological inhibition of Ras in anti-cancer therapy. Best results, so far, were obtained with mutation-specific inhibitors, such as irreversible inhibitors binding mutant RasG12C proteins (Ostrem et al., 2013; Lito et al., 2016; Patricelli et al., 2016; Hansen et al., 2018b; Janes et al., 2018). These molecules target the SII-P allosteric cavity in the GDP-bound form and prevent the GEF-mediated nucleotide exchange and, indirectly, effector engagement. Notably, the in vivo efficacy of these inhibitors depends on the fact that RasG12C does not permanently remain in a GTP-bound form, likely because of relevant retained intrinsic GTPase activity (Hunter et al., 2015). Ras oncoproteins with impairment of both intrinsic and GAP-mediated GTP hydrolysis, such as RasG12R, RasG12V, and RasQL61 (Hunter et al., 2015), would be refractory to this inhibitory action mechanism. Other promising compounds targeting a cryptic phenol-capturing groove near SII-P in both GDP- and GTP-bound forms of non-G12C Ras mutants were identified (Gentile et al., 2017). They are reversible inhibitors preventing the GEF-mediated nucleotide exchange and PI3K engagement, but they do not affect the binding to Raf1. These inhibitors seem particularly interesting for targeting the HRas isoform, which is a more potent activator of PI3K than KRas isoform (Yun et al., 1998). Although new powerful approaches for inhibiting Ras signaling in cancer have been recently developed (Gilardi et al., 2020; Shin et al., 2020) the challenge for the identification of inhibitors effective on the non-G12C pathological Ras variants is still open.
Here we show that cmp4 is a water-soluble pan-Ras inhibitor with a complex, multi-level mechanism of action. cmp4 is the product of rational design from a lead compound in which a 3,4-dihydroxyphenyl group (catechol) and a benzyloxy group, are interconnected by a linear linker (Palmioli et al., 2009b). cmp4 targets the SII-G pocket, as the compounds identified by Gentile et al. (2017). Since cmp4 is bulkier, it occupies a more extended region protruding towards the G12(P-)loop, here named SII-EP (Supplementary Figures S2E,H). The catechol group of cmp4 can undergo several different interactions with residues not only in α2-(switch II) (such as Glu62, Tyr64, and Arg68) and α3-helices (such as Tyr96 and Arg102), but also with the backbone of Gly10 in the P-loop (see ligand interactions plot in Supplementary Figure S1). Natural compounds containing a pyrocatechol group also target this pocket: 5-O-caffeoylquinic acid (5-CQA) takes contact with HRas-GDP through its aromatic caffeic acid moiety but is less efficient in inhibiting RasGRF1 binding (Palmioli et al., 2017). Although the residues in Switch II are the most affected upon cmp4 binding according to NMR analysis, the residues revealing a change in their chemical environment are more widespread along Ras protein (Sacco et al., 2011) suggesting that the binding of the compound could induce a deeper conformational rearrangement that cannot be reproduced by any docking protocol, in agreement with the effects observed for other compounds binding to this area (Ostrem et al., 2013; Gentile et al., 2017).
Treatment with cmp4 prevents intrinsic and GEF-mediated nucleotide exchange, both in wild type and in the G13D-mutated Ras protein, which is self-sufficient in nucleotide dissociation although remaining sensitive to GEF catalytic activity. In addition, cmp4 reduces Ras/Raf1 binding. This effect suggests that cmp4 is able to accommodate in the Switch II pocket of GTP-bound Ras proteins, either interfering with the hydrogen bonds network involved in stabilizing the State 2 Switch I conformation, required for Raf1 binding (Buhrman et al., 2007), or at least destabilizing the ordered Switch II conformation (R state) which allows high-affinity binding to Raf1. A more disordered T State is indeed adopted whenever α3 helix is shifted towards Switch II (Johnson et al., 2019). The presence of cmp4 in the SII-EP site could counteract the shift to the R state, which is characterized by a narrower pocket (Buhrman et al., 2010).
Simulation of the multi-level action mechanism of cmp4 in a computational model describing the Ras activation cycle in conditions designed to represent cells with a constitutively active EGF Receptor, suggests that the compound can work on different Ras oncoproteins, including KRasG13D and KRasG12V. cmp4 effectively cooperates with compounds blocking the Ras signaling cascade at the level of the EGF Receptor, such as Cetuximab. A near additive effect is observed even in the presence of the RasG12V mutant that makes virtual and real cells insensitive to the inhibition of GEF activity resulting from treatments with Cetuximab. In vitro growth inhibition induced by cmp4 and Cetuximab (administered individually or in combination) on isogenic SW48-derived cell lines expressing different Ras mutant proteins fully confirm the simulation results.
With its multi-level mechanism of action that is only minimally superimposed with that of Cetuximab, cmp4 is a good candidate for medicinal chemistry efforts tailored at improving its currently unsatisfactory affinity for Ras proteins.
As a pan-Ras inhibitor, cmp4 is able to inhibit not only Ras oncoproteins but also the wild type variant when activated in a stimulus-dependent way. This would allow cmp4-based drugs to be effectively used in combination therapies with Cetuximab to reduce the proliferation of tumor cells expressing the constitutively activated RTK receptor, but also suggests certain cytotoxicity on proliferating cells in general, given that proliferation in mammalian cells is essentially promoted by Ras signaling. A low affinity is desired when dealing with treatments affecting wild-type Ras, in order to avoid general toxicity to non-proliferating cells, although the affinity of cmp4 still needs some improvement for this aim. It is noteworthy that the development of drugs specific for a pathogenic Ras variant could be achieved by adding chemical groups that can efficiently interact with the variant molecular features, such as the glutamate residue present in G12D or G13D KRas mutants, gaining in specificity and affinity for the targeted oncoprotein and allowing the administration of lower doses, ineffective on wild-type Ras proteins.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
RT, SA, and LG planned, performed, and discussed molecular docking experiments; MS, EM, and ES performed and discussed laboratory experiments (biochemistry and cell biology); PC, DB, MN, and RT planned, performed, and discussed simulation experiments, PC wrote the first draft of simulation methods and results; AP, CA, and FP planned, performed and discussed organic chemistry and NMR experiments; RG planned, performed and discussed MS experiments; MV and ES conceived the experimental plan, discussed all results and figures. RT, MV, and ES wrote the paper.
This work was supported by funds from the Italian Ministry of University and Research (MIUR) through grants “Research facilitation fund (Fondo per le Agevolazioni alla Ricerca–FAR)” to ES and MV, and “Dipartimenti di Eccellenza 2017” to University of Milano Bicocca.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank Annalisa D’Urzo to have implemented a preliminary version of the Ras-cycle mathematical model.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.625979/full#supplementary-material.
Athuluri-Divakar, S. K., Vasquez-Del Carpio, R., Dutta, K., Baker, S. J., Cosenza, S. C., Basu, I., et al. (2016). Divakar 2016 cell small molecule RAS‐mimetic disrupts RAS association with effector proteins to block signaling. Cell 165, 643–655. doi:10.1016/j.cell.2016.03.045 |
Boriack-Sjodin, P. A., Margarit, S. M., Bar-Sagi, D., and Kuriyan, J. (1998). The structural basis of the activation of Ras by Sos. Nature 394, 337–343. doi:10.1038/28548 |
Bos, J. L., Rehmann, H., and Wittinghofer, A. (2007). GEFs and GAPs: critical elements in the control of small G proteins. Cell 129, 865–877. doi:10.1016/j.cell.2007.05.018 |
Buhrman, G., Holzapfel, G., Fetics, S., and Mattos, C. (2010). Allosteric modulation of Ras positions Q61 for a direct role in catalysis. Proc. Natl. Acad. Sci. U.S.A 107, 4931–4936. doi:10.1073/pnas.0912226107 |
Buhrman, G., Wink, G., and Mattos, C. (2007). Article transformation efficiency of RasQ61 mutants linked to structural features of the switch regions in the presence of Raf. Structure. 15, 1618–1629. doi:10.1016/j.str.2007.10.011 |
Burgess, A. W., Cho, H. S., Eigenbrot, C., Ferguson, K. M., Garrett, T. P. J., Leahy, D. J., et al. (2003). An open-and-shut case? recent insights into the activation of EGF/ErbB receptors. Mol. Cell 12, 541–552. doi:10.1016/S1097-2765(03)00350-2 |
Canon, J., Rex, K., Saiki, A. Y., Mohr, C., Cooke, K., Bagal, D., et al. (2019). The clinical KRAS(G12C) inhibitor AMG 510 drives anti-tumour immunity. Nature 575, 217–223. doi:10.1038/s41586-019-1694-1 |
Colombo, S., Palmioli, A., Airoldi, C., Tisi, R., Fantinato, S., Olivieri, S., et al. (2010). Structure-activity studies on arylamides and arysulfonamides Ras inhibitors. Curr. Cancer Drug Targets 10, 192–199. doi:10.2174/156800910791054185 |
Freedman, T. S., Sondermann, H., Friedland, G. D., Kortemme, T., Bar-Sagi, D., Marqusee, S., et al. (2006). A Ras-induced conformational switch in the Ras activator Son of sevenless. Proc. Natl. Acad. Sci. U.S.A. 103, 16692–16697. doi:10.1073/pnas.0608127103 |
Gentile, D. R., Rathinaswamy, M. K., Jenkins, M. L., Moss, S. M., Siempelkamp, B. D., Renslo, A. R., et al. (2017). Ras binder induces a modified switch-II pocket in GTP and GDP states. Cell Chem. Biol. 24, 1455–1466.e14. doi:10.1016/j.chembiol.2017.08.025 |
Gilardi, M., Wang, Z., Proietto, M., Chilla, A., Calleja-Valera, J., Goto, Y., et al. (2020). Tipifarnib as a precision therapy for HRAS-mutant head and neck squamous cell carcinomas. Mol. Cancer Ther. 19, 1784–1796. doi:10.1158/1535-7163.MCT-19-0958 |
Gorfe, A. A., and Cho, K. J. (2021). Approaches to inhibiting oncogenic K-Ras. Small GTPases, 1–10. doi:10.1080/21541248.2019.1655883
Hallin, J., Engstrom, L. D., Hargi, L., Calinisan, A., Aranda, R., Briere, D. M., et al. (2020). The KRASG12C inhibitor MRTX849 provides insight toward therapeutic susceptibility of KRAS-mutant cancers in mouse models and patients. Cancer Discov. 10, 54–71. doi:10.1158/2159-8290.CD-19-1167 |
Hansen, R., Peters, U., Babbar, A., Chen, Y., Feng, J., Janes, M. R., et al. (2018a). The reactivity-driven biochemical mechanism of covalent KRASG12C inhibitors. Nat. Struct. Mol. Biol. 25, 454–462. doi:10.1038/s41594-018-0061-5 |
Hansen, R., Peters, U., Babbar, A., Chen, Y., Feng, J., Janes, M. R., et al. (2018b). The reactivity-driven biochemical mechanism of covalent KRASG12C inhibitors. Nat. Struct. Mol. Biol 25, 454–462. doi:10.1038/s41594-018-0061-5 |
Hoops, S., Sahle, S., Gauges, R., Lee, C., Pahle, J., Simus, N., et al. (2006). COPASI–a complex pathway simulator. Bioinformatics 22, 3067–3074. doi:10.1093/bioinformatics/btl485 |
Hunter, J. C., Manandhar, A., Carrasco, M. A., Gurbani, D., Gondi, S., and Westover, K. D. (2015). Biochemical and structural analysis of common cancer-associated KRAS mutations. Mol. Cancer Res. 13, 1325–1335. doi:10.1158/1541-7786.MCR-15-0203 |
Janes, M. R., Zhang, J., Li, L. S., Hansen, R., Peters, U., Guo, X., et al. (2018). Targeting KRAS mutant cancers with a covalent G12C-specific inhibitor. Cell 172, 578–589.e17. doi:10.1016/j.cell.2018.01.006 |
Johnson, C. W., Reid, D., Parker, J. A., Salter, S., Knihtila, R., Kuzmic, P., et al. (2017). The small GTPases K-Ras, N-Ras, and H-Ras have distinct biochemical properties determined by allosteric effects. J. Biol. Chem. 292, 12981–12993. doi:10.1074/jbc.M117.778886 |
Johnson, C. W., Lin, Y.-J., Reid, D., Parker, J., Pavlopoulos, S., Dischinger, P., et al. (2019). Isoform-specific destabilization of the active site reveals a molecular mechanism of intrinsic activation of KRas G13D. Cell Rep. 28, 1538–1550.e7. doi:10.1016/j.celrep.2019.07.026 |
Jorgensen, W. L., Maxwell, D. S., and Tirado-rives, J. (1996). Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 7863, 11225–11236. doi:10.1021/ja9621760
Khan, I., Rhett, J. M., and O’Bryan, J. P. (2020). Therapeutic targeting of RAS: new hope for drugging the “undruggable. Biochim. Biophys. Acta Mol. Cell Res. 1867, 118570. doi:10.1016/j.bbamcr.2019.118570 |
Lenzen, C., Cool, R., and Wittinghofer, A. (1995). Analysis of intrinsic and CDC25-stimulated guanine nucleotide exchange of p21ras-nucleotide complexes by fluorescence measurements. Methods Enzym. 255, 95–109. doi:10.1016/s0076-6879(95)55012-7
Li, S., Balmain, A., and Counter, C. M. (2018). A model for RAS mutation patterns in cancers: finding the sweet spot. Nat. Rev. Cancer 18, 767–777. doi:10.1038/s41568-018-0076-6 |
Lito, P., Solomon, M., Li, L., Hansen, R., and Rosen, N. (2016). Allele-specific inhibitors inactivate mutant KRAS G12C by a trapping mechanism. Science 351, 604–608. doi:10.1126/science.aad6204 |
Liu, L. J., Wang, W., Huang, S. Y., Hong, Y., Li, G., Lin, S., et al. (2017). Inhibition of the Ras/Raf interaction and repression of renal cancer xenografts in vivo by an enantiomeric iridium(III) metal-based compound. Chem. Sci. 8, 4756–4763. doi:10.1039/c7sc00311k |
Lu, S., Jang, H., Muratcioglu, S., Gursoy, A., Keskin, O., Nussinov, R., et al. (2016a). Ras conformational ensembles, allostery, and signaling. Chem. Rev. 116, 6607–6665. doi:10.1021/acs.chemrev.5b00542 |
Lu, S., Jang, H., Nussinov, R., and Zhang, J. (2016b). The structural basis of oncogenic mutations G12, G13, and Q61 in small GTPase K-Ras4B. Sci. Rep. 6, 1–15. doi:10.1038/srep21949
McFall, T., Diedrich, J. K., Mengistu, M., Littlechild, S. L., Paskvan, K. V., Sisk-Hackworth, L., et al. (2019). A systems mechanism for KRAS mutant allele–specific responses to targeted therapy. Sci. Signal. 12, eaaw8288. doi:10.1126/scisignal.aaw8288 |
McGee, J. H., Shim, S. Y., Lee, S. J., Swanson, P. K., Jiang, S. Y., Durney, M. A., et al. (2018). Exceptionally high-affinity Ras binders that remodel its effector domain. J. Biol. Chem 293, 3265–3280. doi:10.1074/jbc.M117.816348 |
Metcalfe, C., Biochem, J., Owens, G. K., Bjorkerud, S., Lyons, R. M., Gentry, L. E., et al. (1993). Complexes of Ras. GTP with Raf-1 and mitogen-activated protein kinase kinase. Science 260, 1658–61. doi:10.1126/science.8503013 |
Ni, D., Li, X., He, X., Zhang, H., Zhang, J., and Lu, S. (2019). Drugging K-RasG12C through covalent inhibitors: mission possible? Pharmacol. Ther. 202, 1–17. doi:10.1016/j.pharmthera.2019.06.007 |
Omerovic, J., Laude, A. J., and Prior, I. A. (2007). Ras proteins: paradigms for compartmentalised and isoform-specific signalling. Cell. Mol. Life Sci. 64, 2575–2589. doi:10.1007/s00018-007-7133-8 |
Ostrem, J. M., Peters, U., Sos, M. L., Wells, J. A., and Shokat, K. M. (2013). K-Ras(G12C) inhibitors allosterically control GTP affinity and effector interactions. Nature 503, 548–551. doi:10.1038/nature12796 |
Palmioli, A., Ciaramelli, C., Tisi, R., Spinelli, M., De Sanctis, G., Sacco, E., et al. (2017). Natural compounds in cancer prevention: effects of coffee extracts and their main polyphenolic component, 5-O-Caffeoylquinic acid, on oncogenic ras proteins. Chem. Asian J 12, 2457–2466. doi:10.1002/asia.201700844 |
Palmioli, A., Sacco, E., Abraham, S., Thomas, C. J., Domizio, A. D, Gioia, L. De, et al. (2009a). First experimental identification of Ras-inhibitor binding interface using a water-soluble Ras ligand. Bioorganic Med. Chem. Lett. 19, 4217–4222. doi:10.1016/j.bmcl.2009.05.107
Palmioli, A., Sacco, E., Airoldi, C., Di Nicolantonio, F., D’Urzo, A., Shirasawa, S., et al. (2009b). Selective cytotoxicity of a bicyclic Ras inhibitor in cancer cells expressing K-RasG13D. Biochem. Biophys. Res. Commun. 386, 593–597. doi:10.1016/j.bbrc.2009.06.069 |
Patricelli, M. P., Janes, M. R., Li, L. S., Hansen, R., Peters, U., Kessler, L. V., et al. (2016). Selective inhibition of oncogenic KRAS output with small molecules targeting the inactive state. Cancer Discov. 6, 316–329. doi:10.1158/2159-8290.CD-15-1105 |
Petzold, L. (1983). Automatic selection of methods for solving stiff and nonstiff systems of ordinary differential equations. SIAM J. Sci. Stat. Comput. 4, 136–148. doi:10.1137/0904010
Rabara, D., Tran, T. H., Dharmaiah, S., Stephens, R. M., McCormick, F., Simanshu, D. K., et al. (2019). KRAS G13D sensitivity to neurofibromin-mediated GTP hydrolysis. Proc. Natl. Acad. Sci. U.S.A. 116, 22122–22131. doi:10.1073/pnas.1908353116 |
Sacco, E., Abraham, S. J., Palmioli, A., Damore, G., Bargna, A., Mazzoleni, E., et al. (2011). Binding properties and biological characterization of new sugar-derived Ras ligands. Medchemcomm 2, 396–401. doi:10.1039/c0md00264j
Sacco, E., Farina, M., Greco, C., Lamperti, S., Busti, S., DeGioia, L., et al. (2012a). Regulation of hSos1 activity is a system-level property generated by its multi-domain structure. Biotechnol. Adv. 30, 154–168. doi:10.1016/j.biotechadv.2011.07.017 |
Sacco, E., Metalli, D., Spinelli, M., Manzoni, R., Samalikova, M., Grandori, R., et al. (2012b). Novel RasGRF1-derived Tat-fused peptides inhibiting Ras-dependent proliferation and migration in mouse and human cancer cells. Biotechnol. Adv. 30, 233–243. doi:10.1016/j.biotechadv.2011.05.011 |
Sacco, E., Spinelli, M., and Vanoni, M. (2012c). Approaches to Ras signaling modulation and treatment of Ras-dependent disorders: a patent review (2007 present). Expert Opin. Ther. Pat. 22, 1263), 1287. doi:10.1517/13543776.2012.728586 |
Scheffzek, K., Ahmadian, M. R., Kabsch, W., Wiesmüller, L., Lautwein, A., Schmitz, F., et al. (1997). The Ras-RasGAP complex: Structural basis for GTPase activation and its loss in oncogenic ras mutants. Science 277, 333–338. doi:10.1126/science.277.5324.333 |
Seshacharyulu, P., Ponnusamy, M. P., Haridas, D., Jain, M., Ganti, A. K., and Batra, S. K. (2012). Targeting the EGFR signaling pathway in cancer therapy. Expert Opin. Ther. Targets 16, 15–31. doi:10.1517/14728222.2011.648617 |
Shin, S. M., Kim, J. S., Park, S. W., Jun, S. Y., Kweon, H. J., Choi, D. K., et al. (2020). Direct targeting of oncogenic RAS mutants with a tumor-specific cytosol-penetrating antibody inhibits RAS mutant–driven tumor growth. Sci. Adv 6, eaay2174. doi:10.1126/sciadv.aay2174 |
Simanshu, D. K., Nissley, D. V., and McCormick, F. (2017). RAS proteins and their regulators in human disease. Cell 170, 17–33. doi:10.1016/j.cell.2017.06.009 |
Smith, M. J., Neel, B. G., and Ikura, M. (2013). NMR-based functional profiling of RASopathies and oncogenic RAS mutations. Proc. Natl. Acad. Sci. U.S.A. 110, 4574–4579. doi:10.1073/pnas.1218173110 |
Spencer-Smith, R., and O’Bryan, J. P. (2019). Direct inhibition of RAS: quest for the holy grail?. Semin. Cancer Biol. 54, 138–148. doi:10.1016/j.semcancer.2017.12.005 |
Stites, E. C., Trampont, P. C., Ma, Z., and Ravichandran, K. S. (2007). Network analysis of oncogenic Ras activation in cancer. Science 318, 463–467. doi:10.1126/science.1144642 |
Tisi, R., Gaponenko, V., Vanoni, M., and Sacco, E. (2020). Natural products attenuating biosynthesis, processing, and activity of ras oncoproteins: state of the art and future perspectives. Biomolecules 10, 1535. doi:10.3390/biom10111535
Trinh, T. B., Upadhyaya, P., Qian, Z., and Pei, D. (2016). Discovery of a direct ras inhibitor by screening a combinatorial library of cell-permeable bicyclic peptides. ACS Comb. Sci. 18, 75–85. doi:10.1021/acscombsci.5b00164 |
Waldmann, H., Karaguni, I. M., Carpintero, M., Gourzoulidou, E., Herrmann, C., Brockmann, C., et al. (2004). Sulindac-derived ras pathway inhibitors target the ras-raf interaction and downstream effectors in the ras pathway. Angew. Chemie Int. Ed. 43, 454–458. doi:10.1002/anie.200353089
Warne, P. H., Viciana, P. R., and Downward, J. (1993). The amino-terminal region of Raf-1 in vitro. Nature 364, 352–355 |
Welsch, M. E., Kaplan, A., Chambers, J. M., Stokes, M. E., Bos, P. H., Zask, A., et al. (2017). Multivalent small-molecule pan-RAS inhibitors. Cell 168, 878–889.e29. doi:10.1016/j.cell.2017.02.006 |
Wiechmann, S., Maisonneuve, P., Grebbin, B. M., Hoffmeister, M., Kaulich, M., Clevers, H., et al. (2019). Competitive inhibitors of Ras effector binding. bioRxiv:548750.
Yun, C., Yan, I., and Greene, L. A. (1998). Prevention of PC12 cell death by N -acetylcysteine requires activation of the Ras pathway. J Neurosci. 18, 4042–4049. doi:10.1523/JNEUROSCI.18-11-04042.1998 |
Keywords: RasG13D, RasG12V, anti-cancer agent, exchange factor, intrinsic nucleotide dissociation and exchange, Raf1 binding, mathematical modeling & simulation, cetuximab
Citation: Tisi R, Spinelli M, Palmioli A, Airoldi C, Cazzaniga P, Besozzi D, Nobile MS, Mazzoleni E, Arnhold S, De Gioia L, Grandori R, Peri F, Vanoni M and Sacco E (2021) The Multi-Level Mechanism of Action of a Pan-Ras Inhibitor Explains its Antiproliferative Activity on Cetuximab-Resistant Cancer Cells. Front. Mol. Biosci. 8:625979. doi: 10.3389/fmolb.2021.625979
Received: 04 November 2020; Accepted: 11 January 2021;
Published: 17 February 2021.
Edited by:Piero Crespo, Consejo Superior De Investigaciones Científicas, Spain
Reviewed by:Sayan Chakraborty, Institute of Molecular and Cell Biology (A∗STAR), Singapore
Amancio Carnero, Sevilla University, Spain
Copyright © 2021 Tisi, Spinelli, Palmioli, Airoldi, Cazzaniga, Besozzi, Nobile, Mazzoleni, Arnhold, De Gioia, Grandori, Peri, Vanoni and Sacco. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elena Sacco, firstname.lastname@example.org
†These authors have contributed equally to this work