eISSN: 1897-4309
ISSN: 1428-2526
Contemporary Oncology/Współczesna Onkologia
Current issue Archive Manuscripts accepted About the journal Supplements Addendum Special Issues Editorial board Reviewers Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
3/2023
vol. 27
 
Share:
Share:
Original paper

Virtual docking screening and quantitative structure-activity relationship studies to explore AKT and PI3K inhibitors acting on mTOR in cancers by theoretical biology and medical modeling

Ilham Kandoussi
1, 2
,
Hanane Abbou
1, 2
,
Ghyzlane EL Haddoumi
1, 2
,
Mariam Mansouri
1, 2
,
Lahcen Belyamani
1, 2, 3
,
Azeddine Ibrahimi
1, 2

  1. Biotechnology Laboratory (MedBiotech), Bioinova Research Center, Rabat Medical & Pharmacy School, Mohammed V University in Rabat, Morocco
  2. Centre Mohammed VI de la Recherche et de l’Innovation, Rabat, Morocco
  3. Emergency Department, Military Hospital Mohammed V, Rabat, Morocco
Contemp Oncol (Pozn) 2023; 27 (3): 155–162
Online publish date: 2023/12/13
Article file
Get citation
 
PlumX metrics:
 

Introduction

The phosphatidylinositol 3-kinase (PI3K)/AKT/mechanistic target of rapamycin (mTOR) signaling pathway dominates a wide range of cellular processes including survival, proliferation, and growth, is regulated by many signaling proteins upstream and regulates many effectors downstream. This pathway is hyperactivated or altered in many types of cancer [1].

mTOR is a serine/threonine kinase that acts through two structurally and functionally distinct protein complexes, the mTOR1 complex (mTORC1) and the mTOR 2 complex (mTORC2), to detect and integrate multiple intracellular and environmental signals [2, 3].

It is excessively over-activated in more than 70% of cancers [4]. In recent years, it has been widely demonstrated in animal models and cancer patients that mTOR dysfunction contributes to tumorigenesis [5]. Thus, the components of the mTOR pathway are among the most frequently mutated genes in cancer [6].

Several types of mTOR inhibitors, such as rapamycin, its rapalogs and mTORC1/2 kinase inhibitors, have been studied in various cancer models, including breast cancer, lung cancer, colorectal cancer, and others. However, the effects of mTOR inhibitors used as monotherapy in cancer are sometimes mitigated by several mechanisms of resistance [7].

Rapamycin is an allosteric inhibitor of mTOR and has been approved as an immunosuppressant, but interest is focused on its anticancer potential. However, the performance of rapamycin and its analogs (rapalogs) has not been distinguished despite isolated success in subsets of cancers, suggesting that the full therapeutic potential of mTOR targeting has not yet been exploited. A new generation of competitive adenosine triphosphate (ATP) inhibitors that directly target the mTOR catalytic site exhibits potent and complete inhibition of mTOR and is in early clinical trials [8].

Therefore, new therapeutic strategies based on mTOR inhibition need to be further developed, including combination therapies targeting other pathway inhibitors [9]. Inhibition of the PI3K/AKT/mTOR pathway increases antitumor activity [10]. Studies have shown that the combination of compounds simultaneously targeting different molecules of the PI3K/AKT/mTOR pathway leads to synergistic activity[11].

Computational design has indeed become a powerful and essential tool in drug discovery, especially in the context of developing dual-targeted drugs for novel kinase targets. Kinases are enzymes involved in various cellular processes and have been implicated in many diseases, including cancer. Designing drugs that can target multiple kinases simultaneously can offer significant advantages in terms of therapeutic efficacy [12].

Based on these results, we provide an in silico strategy for the exploration of competitive dual PI3K/mTOR and AKT/mTOR inhibitors of ATP, after two previous explorations of AKT/PI3K and mTOR/PI3Kinhibitors, and PI3K/AKT and mTOR/AKT inhibitors [13, 14].

The marked interest in the development of new PI3K/AKT and PI3K/mTOR inhibitors as potential agents for cancer treatment has prompted us to explore the possibility of developing these inhibitors on the basis of quantitative structure-activity relationship (QSAR) models to predict bioactivity of PI3K and AKT inhibitors towards mTOR and the interaction of the best will be evaluated by docking analysis.

Material and methods

Dataset generation

Mechanistic target of rapamycin inhibitors were taken from Binding DataBase (https://www.bindingdb.org) and their IC50 (concentration of molecules causing 50% inhibition) was converted to a logarithmic scale, pIC50. A total of 170 chemically different compounds with high activity with pIC50 greater than 8 were recruited to construct the mTOR-QSAR model. AKT and PI3K inhibitors with a pIC50 greater than 8 were also selected with the aim of predicting their activity towards mTOR using the developed QSAR model and exploring their dual activity. Compounds of significant activity were chosen to predict future effective dual inhibitors.

Quantitative structure-activity relationship model generation

184 2D descriptors available on the MOE 2008.10 (obtained from Chemical Computing Group (CCP); Montreal, QC, Canada) [15] were calculated for the 170 compounds. Invariant and insignificant descriptors were initially eliminated; then the QSAR contingency descriptor selection and intercorrelation matrices between descriptor pairs were used to extract the 52 most relevant molecular descriptors, which were employed for the distance calculation of each database entry.

All 170 selected compounds were distributed randomly to a training set with 128 compounds (75% of the data) and a test set consisting of 42 compounds (25% of the data).

Partial least squares (PLS) analysis based on the leave-one-out (LOO) method was used to correlate molecular descriptors with pIC50 values.

Quantitative structure-activity relationship model validation

The internal validation procedure evaluates the relative predictive performance of the QSAR model, firstly via the correlation coefficient (R2), which is used to measure the correlation between the experimental pIC50 and the predicted interest values with the purpose of observing the variability between the variables in the set test, and secondly by the root mean square error (RMSE), which is used to evaluate the relative error of the QSAR model.

The model is also tested by cross-validation using the LOO method and the computation of the correlation coefficient (R2) and mean squared error (RMSE), while the detection of outliers was conducted using Z-Score, $ Z-SCORE (absolute difference between the value of the model and the activity field, divided by the square root of the mean square error of the data set) and $ XZ-SCORE (absolute difference between the value of the model under a leave-one-out cross validation scheme and the activity field, divided by the square root of the mean square error of the data set).

External validation evaluates the activities of the predictions and the calculation of the numerical parameters using the model.

Activity prediction

The constructed and validated QSAR-mTOR model was used to predict the activity of two groups of AKT and PI3K inhibitors against mTOR, a first one being AKT inhibitors and second one being PI3K inhibitors. These inhibitors, retrieved from the Binding database (https://www.bindingdb.org), have a pIC50 greater than 8 with 578 inhibitors for AKT and 477 inhibitors for PI3K.

After calculating the predicted activity, the 40 inhibitors with the best predictions in each group were selected for docking into mTor in order to explore their dual activity.

Molecular docking

The 3D coordinates of the PI3K inhibitors as well as AKT that showed the best predicted activity in the QSAR-mTOR model were generated from 2D by MarvinView 5.4.1.1. For the docking analysis, 4JT6 (PDB ID), the mTOR crystallized structure recovered from the PDB database with a resolution of 3.6 Å was used. For docking studies mutual graph learning (MGL) tools 1.5.6 with AutoGrid4 and AutoDock vina (Scripps) [16] were used. The mTOR structure was hydrogenated using MGL tools and PyMol was used to visualize the results [17].

In this study we followed the same docking strategy used by the authors in previously published studies [18]. The Grid boxes were generated, using MGL tools 1.5.6, around the active site of the two 3D structures of the mTOR kinase protein. The Grid boxes were set to have 16 to 20 Å of edge with coordinates x = 49.037, y = –0.839, z = –45.349. These coordinates were determined using the potential substrate binding residues as centroids (in the hinge region and the activation loop) [19].

AZD8055 and XL388 are known mTOR inhibitors and were also docked to mTOR to serve as a control to evaluate the docking results. Their 3D structures were obtained from PubChem.

Results

Quantitative structure-activity relationship analysis

The mTOR-QSAR model was built based on 52 molecular descriptors; the validation procedure evaluates the predicted activities and the residuals for the molecules in the training set. In cross-validation, the model is tested using a LOO method. An mTOR-SAR model with a correlation coefficient (R2) of 0.80813 and an RMSE of 0.17756 was obtained after the QSAR regression analysis of 170 molecules. It was validated and evaluated by the cross-validation method LOO. The predictive performance of this model was represented by cross-validated RMSE with 0.54236 and R2 with cross-validation of 0.65460.

Figure 1 presents the plot of both the experimental and the predicted inhibitory potency inpIC50 values of the training and test set compounds, showing a comparable and similar distribution between the two groups. No outliers were detected in the test set data, and all compounds were well-predicted with a residual value less than one log unit.

Fig. 1

Relationship between observed and predicted data from QSAR-mTor model. The compounds of the training set are in red and those of the test set are in blue

/f/fulltexts/WO/52046/WO-27-52046-g001_min.jpg

Virtual screening

The mTOR-QSAR model developed and validated was used to calculate the predicted activity of AKT and PI3K inhibitors against mTOR. Then, the 40 inhibitors with the best predictions in each group were chosen to perform their docking into mTOR. The best predicted AKT and PI3K inhibitors present pIC50 activities in the range 9.36–9.95 and 9.23–9.87, respectively.

Molecular docking studies

The docking scores (affinity) obtained from the docking of AKT inhibitors and PI3K inhibitors in the catalytic site of mTor were between –7.6 and –9.6 kcal/mol for the AKT inhibitor group and between –3.0 and –10.6 kcal/mol for the PI3K inhibitor group. The mTOR reference inhibitors AZD8055 and XL388 presented score of –7.2 and –8.6 kcal/mol, respectively (Tables 1, 2).

Table 1

Docked interaction analysis of AKT inhibitors screened into mTOR

BindingDBReactan t_set_idpIC50WeightPredicted mTOR pIC50Number of H bondsActive site residues and bond length in ÅAffinity [–kcal/mol]
449888.397559.7019.74943.4–3.4 Y2226/3.1 W2239/3.2 K21719.8
3298008.167529.5869.53323.3 E2190/3.5 D23578.4
3298138.102503.5479.42923.0 S2165/3.2 Q21668.3
3298268.397547.6019.54423.0 S2165/3.1 Q21668.2
3298318.318517.5759.63733.3 S2165/2.9 Q2166/3.6 E21908.7
3298388.161555.5839.40753.3 R2348/3.3–3.4 D2252/3.4–3.2 R22518.9
3298428.387533.5739.41023.4 D2357/3.1 T22458.2
3298438.180587.5449.70023.5 D2357/3.3 N23438.1
3298668.657551.4149.71823.3 Q2167/3.4 V22408.4
3298678.552587.8759.94123.3 V2240/3.2 R22518.0
3298798.180477.5099.44133.3 S2342/3.4 R2251/3.3 V22408.3
3299068.823525.3779.85323.4 V2240/3.3 Q21668.3
3299118.677474.5359.40723.4 G2238/3.0 R22517.6
3299128.853480.9259.41823.4 V2240/3.3 Q21678.2
3299138.443506.9639.39123.4 V2240/3.3 Q21678.2
3299488.119531.6019.70023.1 S2165/3.3 Q21668.3
3299498.657537.9929.71823.2 S2165/3.2 Q21668.4
3299508.522536.0029.74232.8 V2240/3.3 Q2167/3.1 S21659.0
3299589523.9669.56823.2 Q2167/3.2 S21658.6
3299598.638583.6099.47123.3 Q2165/3.1 S21659.6
3299608.376515.5589.38523.4 S2342/3.4 T22459.1
3299668.769532.4379.53123.1 Q2167/3.2 S21658.7
3299678.537514.4469.41323.2 Q2167/3.1 S21658.3
3299688.229514.4469.41323.1 Q2167/3.3 S21658.5
3299708.387517.5759.55323.4 S2165/3.4 D22438.4
3300028.309524.0889.51623.2 S2165/3.1 Q21678.1
3300078.376564.4539.95723.4 S2165/3.3 D22438.4
3300088.744548.8919.69723.2 S2165/3.1 Q21678.6
3300208.136543.6129.83413.2 R22518.7
3300528.638542.5849.49123.0–3.0 T22459.3
3301219.096542.5849.43823.3 R2252/3.3 R23499.5
501521639461.4289.83533.5–3.1 R2251/3.2 D22527.9
501521708.522523.3779.44623.2 R2251/3.4 E21908.0
501976658.301466.5409.44012.9 Y22259.5
502126798.221466.5409.44013.1 S21658.5
504146918.337466.5409.44013.0 Y22259.5
504268068.920466.5409.44023.6 G2236/3.4 T22458.8
504537679461.4289.83523.3 D2252/3.5 R22517.8
504537758.698523.3779.36023.5 R2251/3.4 D22528.0
510614968.283429.4529.46653.1 Q2167/3.1 S2165/2.8 E2190/3.2 D2357/2.9 V22409.7
AZD80559.096465.5******33.1 V2240/2.8 S2165/2.9 Q21677.8
XL3888.004455.5******43.0 Q2167/3.3 D2357/3.2 Y2225/2.5 D21958.6

[i] Docking (affinity) scores of AKT inhibitors in the catalytic site of mTOR were set between –7.6 and –9.8 kcal/mol. AZ8055 and XL388 (mTOR reference inhibitors) showed scores of –7.8 and –8.6 kcal/mol.

Table 2

Docked interaction analysis of PI3K inhibitors screened into mTOR

BindingDBReactan t_set_idpIC50WeightPredicted mTOR pIC50Number of H bondsActive site residues and bond length in ÅAffinity [–kcal/mol]
453128.508478.3279.59543.0 K2166/2.8 S2165/3.3 D2167/3.2 Y22258.2
453189.522454.3049.42933.0 T2245/3.0 T2245/3.4 V22403.4
453478.744464.2999.63833.3 Y2225/3.1–3.1 T22458.8
2283398.698444.569.39443.5 D2357/3.4 Y2225/3.2 E2190/3.4 D22438.7
2283408.301430.5329.23443.0 E2190/3.4 Q2167/3.4 D2357/3.4C22438.6
2283418.301414.5349.43033.4 C2243/3.5 D2357/3.1 E21908.1
2349839.318360.4209.24332.9–3.4 V2240/3.4 219010.1
2350168.795390.4469.34432.9 S2160/3.0 Q2167/3.2 N23438.8
2350248.795413.4539.38333.0 V2240/2.9 E2190/3.2 D23579.6
2350908.698352.3979.31953.4 Q2167/3.0 S2165/3.3
D2357/2.9–2.8 Y2225
7.9
2350918.508348.4499.39123.4–3.1 V22419.7
2351038.638382.4239.41133.3 V2240/3.2 Y2225/3.3 D23579.7
2351208.274426.4769.41853.2 Q2165/3.1–3.1 S2165/3.1 T2245/3.2 R22518.0
2351218.585443.4799.23332.9 T2245/2.9 S2342/3.0 Q21678.7
2794608.301467.5159.54423,0 Q2167/2.9 S21658.4
501966288.318462.4619.78633.1 V2240/3.3 Q2243/3.2 S23427.5
501966408.318432.4359.73233.5 D2357/3.3 E2190/3.0 V22403.0
501966598.698476.4879.44423.2 T2245/3.3 S23427.3
501966758.568461.4729.50713.0 V22408.8
501966778.568447.4459.44733.0 Q2167/3.1 S2342/3.4 V22408.7
502911998.443458.4499.81853.3 D2357/3.3 E2190/3.2
Q2167/3.5–2.9 V2240
9.4
502912118.161462.9419.64952.9–3.4 V2240/3.5
Q2167/3.4 E2190/3.4 D2357
9.1
502912128.107423.4579.29643.1 S2165/3.2–3.1 Y2225/3.3 D21957.5
503423009.638515.9519.30543.5 E2190/3.3 Y2225/3.4 D2357/3.4 V224010.6
5034235210.154520.5209.56133.1 Q2167/2.9 S2165/3.1 G22389.0
505631738.958504.6399.50873.0–2.9–3.4 S2165/3.4 K2166/3.2–3.1
D2167/3.3 V2240/3.3 D2244
8.1
505709238.221650.7059.40633.4 Y2225/3.0 R2347/3.3 D22527.9
505713739.000435.5329.83043.3 V2240/3.0–3.1 Q2167/2.9 S21657.4
506235158.221366.3899.37743.4 D2357/3.4 N2343/3.3 V2240/3.1 T22457.6
506983278.537433.4189.48343.1 V2240/G2238 3.2/3.6 T2245/3.2 Q21679.2
507039368.698490.6139.53173.0–3.2 D2243/3.1 V2240/3.1 S23427.4
507173238.522513.6469.26143.5–2.9 V2240/3.1 S2342/3.1 H22479.0
509193488.522537.5479.87943.0 Y2225/3.5 D2357/3.1 D2167/2.9 S21659.5
509274408.657468.4429.49323.1 V2240/3.2 T22448.9
510282038.283508.0069.65552.8 V2240/3.2 D2357/3.0 D2195/3.5
Y2224/2.9 R2251
8.9
510459418.221417.4199.28613.1 Q21678.9
510653158.537445.4059.76143.1–3.5 V2240/3.2 T2245/3.1 R22517.3
510923938.468457.4619.59652.9–3.5 V2240/3.4 D2357/3.3
Q2167/3.2 E2190
9.4
510924128.619615.6579.55453.3 V2240/3.3 D2243/3.0 R2348/3.0
S2165/3.2 D2167
3.2
AZD80559.096465.5*****33.1 V2240/2.8 S2165/2.9 Q21677.8
XL3888.004455.5*****43.0 Q2167/3.3 D2357/3.2 Y2225/2.5 D21958.6

[i] Docking (affinity) scores of PI3K inhibitors in the catalytic site of mTOR were set between –3.0 and –10.6 kcal. AZ8055 and XL388 (mTOR reference inhibitors) showed scores of –7.8 and –8.6 kcal/mol.

Discussion

The PI3K/AKT/mTOR pathway is stimulated by numerous growth factors and their receptors, and it regulates metabolism, growth, survival, cell proliferation, and angiogenesis.

This pathway is one of the most frequently mutated pathways in cancer, leading to cancer progression and induces resistance to existing treatments. Combination inhibition of PI3K/AKT/mTOR signaling pathways is an approach to overcome the acquisition of resistance to cancer treatment.

In order to explore new common inhibitors of AKT/mTOR and PI3K/mTOR, an mTOR2D-QSAR model was constructed to predict the activity of one group of AKT inhibitors and another one of PI3K, then docking of those that had the best prediction (the top 40 in each group), whose results were compared to those of AZD8055 and XL388, was performed.

The mTOR inhibitors used to construct the QSAR model and those of AKT and PI3K whose activity towards mTOR was predicted were extracted from Binding Database and selected according to their IC50 bioactivity. The model was built by the PLS method using MOE software.

The AKT and PI3K inhibitor docking score compared to the AZD8055 and XL388 score predicts that the AKT inhibitors have an affinity between –7.6 and –9.8 kcal/mol and those of PI3K have an affinity between –3.0 and –10.6 kcal/mol, while AZD8055 and XL388 have affinities of only –7.8 and –8.6 kcal/mol, respectively. Therefore 39 AKT inhibitors and 30 PI3K inhibitors had an affinity higher than that of AZD8055, and 17 AKT inhibitors and 22 PI3K inhibitors had a higher affinity than that of XL388. Thus, the affinity and interaction of XL388 would be more effective than that of AZD8055.

Visualization of these interactions showed that the different compounds adapt to the ATP binding site forming hydrogen bonds, in particular with the residue V2240, which is considered as a hinge residue for mTOR. The best inhibitor of AKT and that of PI3K establish respectively 5 and 7 hydrogen bonds with mTOR, compared to AZ8055 and XL388, which have only 3 and 4 connections, respectively. The significant affinity score of some inhibitors despite the low number of hydrogen bonds can be explained by the hydrophobic interactions and that of van der Waals. Thus, the affinity and interaction of XL388 and most of the inhibitors studied are more efficient than those of AZ8055, a potent, selective and oral mTOR kinase inhibitor, which inhibits the proliferation of tumor cells in vitro and in vivo [20], but induces apoptosis [21].

A double comparison of the inhibitors studied with respect to XL388, in consideration of the affinity and the number of mTOR-mediated hydrogen bonds, predicts that 3 AKT inhibitors (Table 3) and 10 PI3K inhibitors (Table 4) have better results than XL388.

Table 3

Comparison of AKT inhibitors studied to XL388

BindingDBReactant_set_idWeightPredicted mTOR pIC50HBAFFINITY [–kcal/mol]b_1rotNa_acca_donSlogPTPSA
44988559.7019.74949.86716.001111.53
329838555.5839.40758.97635.137104.78
51061496429.4529.46659.75622.074123.69

[i] AKT inhibitors with better affinity and number of hydrogen bonds than XL388

Table 4

Comparison of PI3Kinhibitors studied to XL388

BindingDBReactant_set_idWeightPredicted mTOR pIC50HBAFFINITY [–kcal/mol]b_rotNa_acca_donSlogPTPSA
228339444.569.39448.77621.711123.33
228340430.5329.23448.66631.056134.33
50291199458.4499.81859.45623.065132.98
50291211462.9419.64959.15623.501132.98
50342300515.9519.305410.65615.54797.730
50698327433.4189.48349.24711.622101.66
50717323513.6469.26149.05822.414107.55
50919348537.5479.87949.57744.970132.22
51028203508.0069.65558.95810.213135.07
51092393457.4619.59659.45523.670120.09

[i] PI3K inhibitors with better affinity and number of hydrogen bonds than XL388

Taking into account the modified Lipinski rules (no more than 5 hydrogen bond donors, no more than 10 hydrogen bond acceptors, an octanol-water partition coefficient log P not greater than 5.5) and the Veber rule determining good bioavailability (polar surface area ≤ 140 Å (absolute polarity measurement), rotatable bonds ≤ 10 (flexibility measurement)), we found that 2 AKT inhibitors (Table 5), of which 51061496 is the most relevant (5 hydrogen bonds and an affinity of –9.7 kcal/mol), and 9 PI3K inhibitors (Table 6), with 50291199 being the most relevant (5 hydrogen bonds and an affinity of –9.4 kcal/mol), met all the criteria and could be future drugs (Fig. 2). XL388 inhibits the survival and proliferation of certain cell lines and blocks activation of the mTORC1 and mTORC2 complex [22].

Table 5

AKT inhibitors meeting Lipinski and Veber criteria

BindingDBReactant_set_idName
3298386-[4-[1-[2-(azetidin-1-yl)ethyl]-4-[4-fluoro-3-(trifluoromethyl)
phenyl]imidazol-2-yl]piperidin-1-yl]-5-(1H-pyrazol-4-yl)pyrimidin-4-amine
510614962-amino-6-(5-fluoro-6-methoxypyridin-3-yl)-8-[4-(hydroxymethoxy)
cyclohexyl]-4-methylpyrido[2,3-d]pyrimidin-7-one

[i] AKT inhibitors with better affinity and interaction with mTOR than XL388 and meeting Lipinski and Veber criteria

Table 6

PI3K inhibitors meeting Lipinski and Veber criteria

BindingDBReactant_set_idName
228339N-[5-[2-[2-amino-4-ethyl-6-(4-methoxypiperidin-1-yl)
pyrimidin-5-yl]ethynyl]-2-methylpyridin-3-yl]methanesulfonamide
228340N-[5-[2-[2-amino-4-ethyl-6-(4-hydroxypiperidin-1-yl)
pyrimidin-5-yl]ethynyl]-2-methylpyridin-3-yl]methanesulfonamide
50291199N-[5-(2-amino-4-methylpyrido[2,3-d]
pyrimidin-6-yl)-2-methoxypyridin-3-yl]-2,4-difluorobenzenesulfonamide
50291211N-[5-(2-amino-4-methylpyrido[2,3-d]pyrimidin-6-yl)-2-methoxypyridin-3-yl]
-5-chlorothiophene-2-sulfonamide
506983272-(difluoromethyl)-1-(4,6-dimorpholin-4-yl-1,3,5-triazin-2-yl)benzimidazol-4-ol
507173234-[2-(1H-indazol-4-yl)-6-[(4-methylsulfonylpiperazin-1-yl)
methyl]thieno[2,3-d]pyrimidin-4-yl]morpholine
50919348N-[2-amino-5-[3-(3,4-dimethoxyphenyl)-1H-pyrrolo[2,3-b]pyridin-5-yl]
pyridin-3-yl]-2,4-difluorobenzenesulfonamide
510282035-[5-chloro-2-[(4-methylsulfonylpiperazin-1-yl)methyl]-8-morpholin-4-ylimidazo[1,2-a]
pyrazin-6-yl]pyrimidin-2-amine
51092393N-[5-(2-amino-4-methylquinazolin-6-yl)-2-methoxypyridin-3-yl]-2,4-difluorobenzenesulfonamide

[i] PI3K inhibitors with better affinity and interaction with mTOR than XL388 and meeting Lipinski and Veber criteria

Fig. 2

Visualization of the different interactions between the compounds and the active site of mTOR via the hydrogen bonds. AZD8055 (A) and XL388 (B) respectively have 3 and 4 linkages with the catalytic site of mTOR whereas 51061496 inhibitor of AKT (C) and 50291199 inhibitor of PI3K (D) both have 5 hydrogen bonds

Red zones – oxygen atom, blue zones – nitrogen atom, green zones – other. Hydrogen bonds are represented by a dashed yellow line. Numbers represent the size of the hydrogen bonds established between the ligand and the receptor.

/f/fulltexts/WO/52046/WO-27-52046-g002_min.jpg

Conclusions

The hyperactivation of the PI3K/AKT/mTOR pathway in cancer, in association with the crucial role of mTOR signaling in tumorigenesis, has led to considerable efforts to generate inhibitors to target this pathway. Using a QSAR model, we were able to predict the activity of the powerful AKT and PI3K inhibitors towards mTOR, select the best ones and analyze their interaction with mTOR by docking while comparing them to a reference mTOR inhibitor. Promising results were obtained, and the dual inhibitors AKT/mTOR and PI3K/mTOR were predicted, pending the testing of these compounds directly on cancer cell lines in our next study.

Acknowledgments

This work was carried out with national funding from the Moroccan Ministry of Higher Education and Scientific Research (PPR program) to AI.

Notes

[7] Conflicts of interest The authors declare no conflict of interest.

References

1 

Ersahin T, Tuncbag N, Cetin-Atalay R. The PI3K/AKT/mTOR interactive pathway. Mol Biosyst 2015; 11: 1946-1954.

2 

Laplante M, Sabatini DM. mTOR signaling in growth control and disease. Cell 2012; 149: 274-293.

3 

Saxton RA, Sabatini DM. mTOR signaling in growth, metabolism, and disease. Cell 2017; 168: 960-976.

4 

Forbes SA, Bindal N, Bamford S, et al. COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res 2011; 39: D945-D950.

5 

Ciuffreda L, di Sanza C, Incani UC, Milella M. The mTOR pathway: a new target in cancer therapy. Curr Cancer Drug Targets 2010; 10: 484-495.

6 

Mayer IA, Arteaga CL. The PI3K/AKT pathway as a target for cancer treatment. Ann Rev Med 2016; 67: 11-28.

7 

Faes S, Demartines N, Dormond O. Resistance to mTORC1 inhibitors in cancer therapy: from kinase mutations to intratumoral heterogeneity of kinase activity. Oxid Med Cell Longev 2017:2017:1726078.

8 

Benjamin D, Colombi M, Moroni C, et al. Rapamycin passes the torch: a new generation of mTOR inhibitors. Nat Rev Drug Discov 2011; 10: 868-880.

9 

Tian T, Li X, Zhang J. mTOR signaling in cancer and mTOR inhibitors in solid tumor targeting therapy. Int J Mol Sci 2019; 20: 755.

10 

Vivanco I, Sawyers CL. The phosphatidylinositol 3-kinase-AKT pathway in human cancer. Nat Rev Cancer 2002; 2: 489.

11 

Polivka J, Janku F. Molecular targets for cancer therapy in the PI3K/AKT/mTOR pathway. Pharmacol Ther 2014; 142: 164-175.

12 

Sun D, Zhao Y, Zhang S, et al. Dual-target kinase drug design: current strategies and future directions in cancer therapy. Eur J Med Chem 2020; 188: 112025.

13 

Kandoussi I, Benherrif O, Lakhlili W, et al. Virtual docking screening and QSAR studies to explore AKT and mTOR inhibitors acting on PI3K in cancers. Contemp Oncol (Poznan) 2020; 24: 5-12.

14 

Kandoussi I, Benherrif O, Lakhlili W, et al. Virtual docking screening and QSAR studies to explore PI3K and mTOR inhibitors acting on AKT in cancers. J Cancer Sci Clin Ther 2021; 5: 35-48.

15 

Vilar S, Cozza G, Moro S. Medicinal chemistry and the molecular operating environment (MOE): application of QSAR and molecular docking to drug discovery. Curr Top Med Chem 2008; 8: 1555-1572.

16 

Trott O, Olson AJ. AutoDock Vina: improving the speed and accu-racy of docking with a new scoring function, efficient optimization and multithreading. J Comput Chem 2010; 31: 455-461.

17 

Seeliger D, de Groot BL. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J Comput Aided Mol Des 2010; 24: 417-422.

18 

Lakhlili W, Chevé G, Yasri A, et al. Determination and validation of mTOR kinase-domain 3D structure by homology modeling. OncoTargets Ther 2015; 8: 1923-1930.

19 

Lakhlili W, Yasri A, Ibrahimi A. Structure-activity relationships study of mTOR kinase inhibition using QSAR and structure-based drug design approaches. OncoTargets Ther 2016; 9: 7345-7353.

20 

Chresta CM, Davies BR, Hickson I, et al. AZD8055 is a potent, selective, and orally bioavailable ATP-competitive mammalian target of rapamycin kinase inhibitor with in vitro and in vivo antitumor activity. Cancer Res 2010; 70: 288-298.

21 

Cirstea D, Santo L, Hideshima T, et al. Delineating the mTOR kinase pathway using a dual TORC1/2 inhibitor, AZD8055, in multiple myeloma. Mol Cancer Ther 2014; 13: 2489-2500.

22 

Xiong Z, Zang Y, Zhong S, et al. The preclinical assessment of XL388, a mTOR kinase inhibitor, as a promising anti-renal cell carcinoma agent. Oncotarget 2017; 8: 30151-30161.

Copyright: © 2023 Termedia Sp. z o. o. This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
 
Quick links
© 2024 Termedia Sp. z o.o.
Developed by Bentus.