| Literature DB >> 32331219 |
Łukasz Tomorowicz1, Jarosław Sławiński1, Beata Żołnowska1, Krzysztof Szafrański1, Anna Kawiak2.
Abstract
A series of novel 2-[(4-amino-6-R2-1,3,5-triazin-2-yl)methylthio]-4-chloro-5-methyl-N-(5-R1-1H-benzo[d]imidazol-2(3H)-ylidene)Entities:
Keywords: 1,3,5-triazines; QSAR; benzenesulfonamide; cytotoxicity; molecular docking; synthesis
Year: 2020 PMID: 32331219 PMCID: PMC7215599 DOI: 10.3390/ijms21082924
Source DB: PubMed Journal: Int J Mol Sci ISSN: 1422-0067 Impact factor: 5.923
Figure 1Structures of Nutlin-3a and a molecular hybrid including the 2-mercaptobenzenosulfonamide fragment and imidazoline ring. The green dotted frame indicates the 2-mercaptobenzenosulfonamide fragment and blue dotted frame shows the imidazoline ring.
Scheme 1Synthesis of 5-substituted ethyl 2-{5-R1-2-[N-(5-chloro-1H-benzo[d]imidazol-2(3H)-ylidene)sulfamoyl]-4-methylphenylthio}acetate 3–5 and [(4-amino-6-R2-1,3,5-triazin-2-yl)methylthio]-4-chloro-5-methyl-N-(5-R1-1H-benzo[d]imidazol-2(3H)-ylidene)benzenesulfonamide 6–49. Reagents and conditions: (a) ethyl chloroacetate (1.1 molar equiv.), dry THF, anhydrous K2CO3, reflux 5 h; (b) appropriate 4-R1-benzene-1,2-diamine (1.05 molar equiv.), glacial acetic acid, reflux 7 h.
Cytotoxicity of compounds 6‒49 toward human cancer cell lines and non-cancerous lines HaCaT.
| Compd | IC50 [µM] | |||
|---|---|---|---|---|
| HCT-116 | MCF-7 | HeLa | HaCaT | |
|
| 49 ± 1 | 52 ± 1 | 50 ± 2 | 57 ± 3 |
|
| 75 ± 3 | 99 ± 3 | 86 ± 3 | 102 ± 4 |
|
| 23 ± 1 | 34 ± 2 | 31 ± 2 | 43 ± 2 |
|
| 25 ± 1 | 20 ± 1 | 23 ± 1 | 42 ± 2 |
|
| 32 ± 3 | 25 ± 1 | 30 ± 1 | 45 ± 2 |
|
| 37 ± 2 | 41 ± 1 | 39 ± 1 | 58 ± 2 |
|
| 36 ± 1 | 35 ± 2 | 36 ± 2 | 43 ± 2 |
|
| 32 ± 1 | 24 ± 1 | 35 ± 1 | 41 ± 2 |
|
| 54 ± 2 | 41 ± 1 | 56 ± 2 | 65 ± 3 |
|
| 22 ± 1 | 27 ± 1 | 92 ± 2 | NT * |
|
| 33 ± 1 | 32 ± 1 | 33 ± 1 | 45 ± 2 |
|
| 29 ± 1 | 31 ± 1 | 30 ± 1 | 43 ± 2 |
|
| 42 ± 2 | 82 ± 3 | 61 ± 2 | 84 ± 3 |
|
| 15 ± 1 | 19 ± 1 | 18 ± 1 | 31 ± 1 |
|
| 16 ± 0.5 | 17 ± 0.5 | 17 ± 1 | 35 ± 1 |
|
| 15 ± 1 | 24 ± 1 | 17 ± 0.5 | 32 ± 1 |
|
| 11 ± 0.5 | 24 ± 1 | 11 ± 0.5 | 34 ± 1 |
|
| 16 ± 0.5 | 17 ± 0.5 | 17 ± 1 | 35 ± 1 |
|
| 17 ± 1 | 20 ± 0.6 | 27 ± 0.5 | 42 ± 2 |
|
| 16 ± 1 | 46 ± 2 | 28 ± 1 | 52 ± 2 |
|
| 17 ± 0.5 | 20 ± 1 | 21 ± 0.4 | 35 ± 1 |
|
| 35 ± 1 | 110 ± 3 | 74 ± 1 | 98 ± 2 |
|
| 16 ± 1 | 24 ± 1 | 24 ± 1 | 35 ± 1 |
|
| 195 ± 4 | 134 ± 9 | 56 ± 2 | NT * |
|
| 17 ± 0.5 | 24 ± 1 | 18 ± 1 | 36 ± 2 |
|
| 18 ± 1 | 46 ± 2 | 39 ± 1 | 54 ± 2 |
|
| 280 ± 11 | 137 ± 3 | 250 ± 13 | NT * |
|
| 19 ± 1 | 35 ± 2 | 30 ± 2 | 48±2 |
|
| 17 ± 1 | 25 ± 1 | 30 ± 1 | 41 ± 1 |
|
| 37 ± 1 | 40 ± 1 | 39 ± 1 | 52 ± 2 |
|
| 18 ± 0.5 | 37 ± 2 | 29 ± 1 | 49 ± 1 |
|
| 17 ± 1 | 27 ± 2 | 20 ± 1 | 39 ± 1 |
|
| 14±1 | 24 ± 1 | 18 ± 1 | 38 ± 2 |
|
| 35±2 | 74 ± 2 | 28 ± 2 | 76 ± 2 |
|
| 16 ± 1 | 43 ± 1 | 18 ± 1 | 46 ± 1 |
|
| 25 ± 1 | 38 ± 2 | 30 ± 1 | 49 ± 1 |
|
| 18 ± 1 | 23 ± 1 | 24 ± 1 | 37 ± 1 |
|
| 43 ± 1 | 44 ± 2 | 45 ± 2 | 55 ± 2 |
|
| 48 ± 1 | 47 ± 2 | 50 ± 2 | 61 ± 2 |
|
| 15 ± 1 | 26 ± 1 | 28 ± 1 | 39 ± 1 |
|
| 7 ± 0.1 | 15 ± 1 | 18 ± 1 | 28 ± 1 |
|
| 21 ± 1 | 24 ± 1 | 22 ± 1 | 38 ± 1 |
|
| 34 ± 2 | 32 ± 1 | 34 ± 1 | 40 ± 2 |
|
| 30 ± 1 | 35 ± 1 | 34 ± 1 | 96 ± 1 |
|
| 2.2 ± 0.1 | 3.0 ± 0.1 | 3.8 ± 0.2 | |
NT*—not tested.
Figure 2The graph presents a comparison of the most active compounds 22 and 46 with the non-cancerous HaCaT cell line.
Summary of the QSAR equations.
|
|
| IC50 = 29.955511 a_nO + 0.002053 pmi3 + 8.132036 E_oop − 1030.177341 GCUT_SLOGP_ 1 + |
| Cell Line: MCF-7 |
| IC50 = −0.429856 a_IC + 33.052386 b_max1len + 137.464827 GCUT_SLOGP_ 2 + 1.000886 |
| Cell Line: HeLa |
| IC50 = 22.605987 ast_violation − 3.368272 a_nF − 6.575727 b_1rotN + 22.206001 h_pstrain + 0.001996 pmi + 0.547602 SlogP_VSA5 − 40.085275 |
R2—squared correlation coefficient for training set; Q2—squared correlation coefficient for test.set; F—Fisher’test; p—p-value for Fisher’s test for the whole equation; N(train)—training set; N(test)—test set. Molecular descriptors used in the models: a_nO—number of oxygen atoms (The atom count and bond count descriptors); pmi3—third diagonal element of diagonalized moment of inertia tensor (Surface Area, Volume and Shape Descriptors); E_oop—Out-of-plane potential energy (The energy descriptors); GCUT_SLOGP_1 and GCUT_SLOGP_2—descriptors using atomic contribution to logP (using the Wildman and Crippen SlogP method) instead of partial charge (Adjacency and Distance Matrix Descriptors); b_max1len—length of the longest single bond chain (The atom count and bond count descriptors); vsurf_IW6—hydrophilic integy moment (Surface Area, Volume and Shape Descriptors); a_IC—atom information content (total) (The atom count and bond count descriptors); PEOE_VSA+1—sum of vi where qi is in the range [0.05,0.10] (Partial Charge Descriptors); SMR_VSA0—adjacency and distance matrix descriptor (The Subdivided Surface Areas); std_dim3—standard dimension 3: the square root of the third largest eigenvalue of the covariance matrix of the atomic coordinates. A standard dimension is equivalent to the standard deviation along a principal component axis. (Surface Area, Volume and Shape Descriptors); ast_violation—number of Astex fragment-likeness violations (The atom count and bond count descriptors); a_nF—number of fluorine atoms (The atom count and bond count descriptors); b_1rotN—number of rotatable single bonds (The atom count and bond count descriptors); h_pstrain—the strain energy (kcal/mol) needed to convert all protonation states into the input protonation state: (kT ln 10) ( pC1 + log sum {10-pCi} ) (The Hueckel Theory descriptors); pmi—principal moment of inertia (Surface Area, Volume and Shape Descriptors); SlogP_VSA5—represent different aspects of van der Waals surface area’s contribution to compound lipophilicity (The Subdivided Surface Areas).
Figure 3Scatter plot of experimental IC50 values versus predicted IC50 values of the training set (blue) and test set (red).
Figure 43D and 2D representations of the protein–ligand (Nutlin-3a) interaction inside the active site of MDM2. Left side of figure: carbon—green, nitrogen—blue, oxygen—red, hydrogen—grey.
Docking results of the most active compounds 22, 46 with MDM2 protein (PDB ID:5C5A).
| Compound | S (kcal/mol) | Amino | Interacting Groups | Type of | Length |
|---|---|---|---|---|---|
|
| −10.1334 | Leu54 | C (piperazine) | H-donor | 3.65 |
| Met50 | Benzimidazole ring | π-H | 4.45 | ||
| Lys51 | Ph (1,2,4,5-tetrasubstituted) | π-H | 4.83 | ||
| Leu54 | 1,3,5-triazine ring | π-H | 4.06 | ||
| Leu54 | 1,3,5-triazine ring | π-H | 3.61 | ||
|
| −9.7475 | Leu54 | C (piperazine) | H-donor | 3.61 |
| Met62 | C (4-F-C6H4) | H-donor | 3.72 | ||
| Lys51 | O (S=O sulfonamide) | H-acceptor | 3.34 | ||
| Leu54 | 1,3,5-triazine ring | π-H | 4.15 | ||
| Leu54 | 1,3,5-triazine ring | π-H | 3.57 | ||
| Tyr100 | Benzimidazole ring | π-π stacking | 4.22 | ||
| Nutlin-3a | −10.8029 | Met62 | C (piperazine) | H-donor | 3.95 |
| Met62 | C ( | H-donor | 3.82 | ||
| His96 | Cl (4-Cl-C6H4) | H-donor | 3.84 | ||
| Gly58 | O [=N-(C=O)-N=] | H-acceptor | 3.29 | ||
| Gln59 | O (piperazin-2-one) | H-acceptor | 2.87 | ||
| Leu54 | Ph (4-Cl-C6H4) | π-H | 4.48 | ||
| Val93 | Ph (1,2,4-trisubstituted) | π-H | 3.75 |
Figure 53D and 2D representations of the protein–ligand interaction of compounds (A) 22, and (B) 46 inside the active site of MDM2. Left side of figure: carbon and oxygen—red, nitrogen—blue, hydrogen—grey, sulfur—yellow, chlorine and fluorine—green.
Figure 6Superimposed representation of the compounds 22, 46 and Nutlin-3a (green) in the active site of MDM2. Atoms other than carbon: sulfur—yellow, nitrogen—violet, oxygen—red, hydrogen—grey.