Tran Thi Phuong Thao‡
ab,
Thanh Q. Bui‡c,
Nguyen Thi Thanh Haic,
Lam K. Huynh
dh,
Phan Tu Quye,
Nguyen Chi Baof,
Nguyen Thi Dunga,
Nguyen Linh Chiab,
Tran Van Locab,
Irina E. Smirnovag,
Anastasiya V. Petrovag,
Pham Thi Ninhab,
Tran Van Sung*ab and
Nguyen Thi Ai Nhung
*c
aInstitute of Chemistry, Vietnam Academy of Science and Technology (VAST), 18 Hoang Quoc Viet Road, Cau Giay, Ha Noi, Vietnam. E-mail: tranvansungvhh@gmail.com
bGraduate University of Science and Technology, VAST, 18 Hoang Quoc Viet Road, Cau Giay, Ha Noi, Vietnam
cDepartment of Chemistry, University of Sciences, Hue University, Hue City, Vietnam. E-mail: ntanhung@hueuni.edu.vn
dInternational University, Quarter 6, Linh Trung Ward, Thu Duc District, Ho Chi Minh City, Vietnam
eDepartment of Natural Sciences & Technology, Tay Nguyen University, Buon Ma Thuot, Vietnam
fHue University, Hue City, Vietnam
gUfa Institute of Chemistry-Subdivision of the Ufa Federal Research Centre of the Russian Academy of Sciences, Prospekt Oktyabrya 71, Ufa, Russian Federation
hVietnam National University, Ho Chi Minh City, Vietnam
First published on 4th November 2021
Dipterocarpus alatus-derived products are expected to exhibit anti-diabetes properties. Natural dipterocarpol (1) was isolated from Dipterocarpus alatus collected in Quang Nam province, Vietnam; afterwards, 20 derivatives including 13 oxime esters (2 and 3a–3m) and 7 lactones (4, 5, 6a–6e) were semi-synthesised. Their inhibitory effects towards diabetes-related proteins were investigated experimentally (α-glucosidase) and computationally (3W37, 3AJ7, and PTP1B). Except for compound 2, the other 19 compounds (3a–3m, 4, 5, and 6a–6d) are reported for the first time, which were modified at positions C-3, C-24 and C-25 of the dipterocarpol via imidation, esterification, oxidative cleavage and lactonisation reactions. A framework based on docking-QSARIS combination was proposed to predict the inhibitory behaviour of the ligand-protein complexes. Enzyme assays revealed the most effective α-glucosidase inhibitors, which follow the order 5 (IC50 of 2.73 ± 0.05 μM) > 6c (IC50 of 4.62 ± 0.12 μM) > 6e (IC50 of 7.31 ± 0.11 μM), and the computation-based analysis confirmed this, i.e., 5 (mass: 416.2 amu; polarisability: 52.4 Å3; DS: −14.9 kcal mol−1) > 6c (mass: 490.1 amu; polarisability: 48.8 Å3; DS: −13.7 kcal mol−1) > 6e (mass: 549.2 amu; polarisability: 51.6 Å3; DS: −15.2 kcal mol−1). Further theoretical justifications predicted 5 and 6c as versatile anti-diabetic inhibitors. The experimental results encourage next stages for the development of anti-diabetic drugs and the computational strategy invites more relevant work for validation.
Inhibition of glucose- and insulin-based enzymes is the main strategy for the therapeutic treatment of type 2 diabetes. The former attempts to reduce gut glucose absorption by controlling postprandial hyperglycemia by mitigating the activity of glucosidases; meanwhile, the latter involves insulin signaling regulation by compensating defects in insulin secretion and insulin action.4 Firstly, α-glucosidase, an exoenzyme found in animals, plants, bacteria, and fungi, breaks down starch and disaccharides to glucose. A study suggested that this enzyme can only yield monosaccharides by catalysing the hydrolysis of α-(1→4) and α-(1→6) bonds,5 confirming the sources of α-glucosidase from sugar beet seeds.6 Data on the α-glucosidase crystal structure can be found at the Worldwide Protein Data Bank under entry PDB-3W37 (DOI: 10.2210/pdb3W37/pdb). Secondly, oligo-1,6-glucosidase, also known as isomaltase, is a debranching endoenzyme, hydrolysing only the α-1,6 linkage in starch and glycogen to produce sugars with an α-configuration.7 Besides some bacterial species such as Bacillus cereus, oligo-1,6-glucosidase is present mainly in the animal kingdom.8 In humans, it can be found on the small intestine brush border.9 Protein structural data of isomaltase can be download directly from the Worldwide Protein Data Bank database under entry PDB-3AJ7 (DOI: 10.2210/pdb3AJ7/pdb). Thirdly, protein tyrosine phosphatase 1B (PTP1B) is a major glucose-homeostasis and energy-metabolism regulator, and thus regarded as an attractive drug target for therapeutic intervention in type 2 diabetes and obesity.10 It blocks insulin receptor substrate-1 and dephosphorylate phosphotyrosine residues, thus causing insulin insensitivity or even cut-off of intracellular insulin signaling.11 Conversely, regarding the leptin signaling pathway, this protein binds and dephosphorylates leptin receptor Janus kinase 2 (JAK2), thereby resulting in the malfunctioning of energy balance.11 The PTP1B structure has already been well determined and its information is available for public reference at UniProtKB, archived under entry ID: UniProtKB-A0A0U1XP67. Therefore, 3W37, 3AJ7, and PTP1B (Fig. 1) are considered highly promising drug targets for the effective treatment of type 2 diabetes. In principle, multiple inhibition of the protein tyrosine phosphatase 1B (PTP1B) and glycoside hydrolase proteins (3W37 and 3AJ7) is a promising strategy to suppress hyperglycemia and improve insulin sensitisation simultaneously.
![]() | ||
Fig. 1 (a) α-Glucosidase protein 3W37. (b) Oligo-1,6-glucosidase protein 3W37. (c) Protein tyrosine phosphatase 1B PTP1B. |
Natural products are considered highly valuable for drug discovery due to their diverse structures and biological activities. In particular, triterpenes belong to a large group of naturally occurring compounds derived from the precursor squalene.12 Dammarane is a tetracyclic triterpene, existing as the main components in dammar resin13,14 and ginseng.15 The basic skeleton of dammarane is comprised of a tetracyclic moiety and a side chain at C-17.16,17 Recently, the number of publications and patents on the anti-diabetic activity of dammarane triterpenoids has increased rapidly.18–21 Besides, clinical trials and animal model experiments have proven the ability of dammarane triterpenoids to prevent and treat diabetes by lowering the blood glucose, regulating lipid metabolism, and increasing insulin sensitivity.22 Thus, considering the advantage of the diversity and novelty of the dammarane structure, a series of chemical modifications of leading compounds has been achieved. Consequently, many synthetic dammarane derivatives have been reported to improve hyperglycemia, glucose tolerance and α-glucosidase inhibition.22–25 Dipterocarpol, a dammarane-type triterpenoid, is the main metabolite of the resin from Dipterocarpus alatus (Dipterocarpaceae), a tropical plant growing in countries in South East Asia such as Cambodia, Laos, Thailand, the Philippines and Vietnam26 (Fig. 2). It was reported that dipterocarpol and its derivatives exhibit antituberculosis,27 antimicrobial,28 hepatoprotective,29 and cytotoxic activity.30
![]() | ||
Fig. 2 (a) Dipterocarpus alatus tree. (b) Collection of dammar resin from Dipterocarpus alatus. (c) Structure of dipterocarpol. |
In this work, two different approaches to modify the structure of dipterocarpol are proposed based on their advantages reported in the literature. The modification of the ring at C-3 of dipterocarpol to obtain oxime and its acyl derivatives (Fig. 3a) was designed. The activating group ketone at C-3 can be easily converted to oxime, followed by O-acylation to obtain acyloxyimino compounds (–CRN–O–COOR). This type of compound (containing acyloxyimino functional groups) is known to exhibit important biological activities, e.g., antiviral, anticancer and anti-diabetes.31 According to our knowledge, several studies on the synthesis of oxime esters derived from lupane or oleane triterpenoids have been carried out;31 nevertheless, similar transformation regarding dammarane triterpenoids has not reported to date. Another modification was at the side chain, focusing on the formation of a lactone ring, instead of a long-chain linear of dipterocarpol (Fig. 3b). Previous studies showed that the location of a lactone ring at C-17 of some dammarane derivatives helps promote the antihyperglycemic effects significantly.20,25,32 Besides, a series of ester derivatives at C-3 was also introduced.
In silico research has been developed extensively in an attempt to accurately predict molecular behaviour at the biological scale. In particular, molecular dynamics gives a view of the dynamic “evolution” of atoms and molecules in a system, whose trajectories are determined by numerically solving Newton's equations of motion for a system of interacting particles. This technique is widely known as one of the most powerful computer-aided tools for pharmacophore development and drug design. However, its downside of long computation times, high cost and significant power consumption highly limit its applicability. Conversely, although successfully complementing the dynamics-related disadvantage by introducing a view of static inhibitory stability, the molecular docking technique is based on a precursor condition that the ligands are already placed within a certain radius to the targeted protein, and thus there is a lack of information on their preceding transport in or interaction with plasmatic components in biological media. Accordingly, a strategy based on the combination of the static docking technique and molecular dynamics has been proposed in an attempt to retrieve meaningful models for ligand–protein inhibition.33–35 Nevertheless, the cost of its implementation is still considered high in low-based research institutes. QSARIS,36 which was released by SciVision Biotech Inc., is a leading-edge database system for predicting the physiochemical properties of lead compounds in drug and agrochemical research. This software promises sufficient compensation for the information lost in the pre-docking presumption. Therefore, a QSARIS-docking combination can be a potential approach for the prediction of ligand–protein inhibition, in general.
Herein, both experimental and theoretical studies were conducted. The isolation of dipterocarpol and semi-synthesis of its twenty derivatives were carried out on Dipterocarpus alatus from Vietnam. Their α-glucosidase inhibition was investigated by enzyme assays. The experimental evidence collected was used to retrieve a comprehensive theoretical model framed by molecular docking simulation coupled with QSARIS-based analysis. The experiment-theory retrieval was applied to all the diabetes-related protein structures.
The synthesis of the oxime esters is summarised in Scheme 1.
The synthesis of the lactone derivatives is summarised in Scheme 2.
α-Glucosidase inhibition (%) = [Abs(control) − Abs(sample)]/Abs(control) × 100% |
To the best of our knowledge, there was only one study reported on the oxidative cleavage at the double bond (C-24 and C-25) of dipterocarpol. The authors used ozone/CH2Cl2 or ozone/MeOH to obtain hemiacetal or acetal, respectively.21 In our case, other oxidative cleavage reagents were utilised such as CrO3, SeO2, KMnO4. Consequently, the use of CrO3/CH3COOH (Jones oxidation) afforded the lactone products (compound 4) under mild conditions (room temperature, 40 min) in good yield (85.5%). Other oxidative reagents gave several unwanted products, which could not be purified after the reaction. Under Jones oxidation, the aldehyde formed at C-24 was further oxidised to carboxylic acid, which continuously underwent an intramolecular reaction with the hydroxyl group at C-20 to form lactone 4 in acidic condition during oxidation. The reduction of the ketone group at C-3 of compound 4 to hydroxyl group was carried out with NaBH4 in MeOH for 20 min. Compound 5 was obtained in 74% yield. Due to the steric hindrance in the β-site of the A-ring of 4, the BH4− anion can only approach the C-3-keto group from the α-site, and thus the 3β-hydroxyl product 5 is the major product. Compound 5 was the only product obtained after column chromatography of the reduction mixture. Further, 5 was reacted with anhydric acid/pyridine to produce 6a, with chloride acid to give 6b and 6c, and with carboxylic acid/DCC/DMAP to obtain 6d and 6e in good yields.
Instead of a ketone signal (δC 218.08) in the 13C NMR spectrum of 1, a signal at δC 167.40, which is characteristic of an oxime carbon at C-3, was observed for compound 2. The formation of oxime esters 3a–3m was confirmed by the appearance of carbonyl ester groups in the range of δC 178.02–174.48 in their 13C NMR spectra. The signals representing the aromatic carbons of compounds 3a–3i are resonated at δC 163.03–124.37 and δH 9.17–7.12 in their NMR spectra. The 1H NMR spectrum of compound 3k provides an additional methyl signal at δH 2.17 The signals at δH 2.42
and 0.92
demonstrate the presence of a pentanoyl group in compound 3l. The presence of an isobutyryl group is inferred by the signal at δH 1.23 (6H, d, J = 6.0) in the 1H NMR spectrum of compound 3m. The NMR spectra of compound 4 indicate the disappearance of the olefine bond at C-24 and C-25 and the two methyl groups at C-26, C-27, compared with compound 1. The presence of a carbonyl lactone at C-24 (δC 176.69), together with the much more down-field chemical shift of the tertiary oxygenated carbon (δC 90.02) suggests the formation of a lactone ring at the side chain of the dipterocarpol skeleton in compound 4. The structure of 4 was further confirmed by 2D NMR (COSY, HSQC, HMBC, and NOESY) spectra. The HMBC correlations between H-23 (δH 2.55–2.50), H-22 (δH 2.49–2.47)/C-24 (δC 176.69), C-20 (δC 90.02), H-21 (δH 1.30)/C-20 (δC 90.02) and C-17 (49.30) confirm the location of the lactone ring at C-17. The orientation of the methyl group (C-21) was established as α by the NOESY correlation between H-17 (δH 2.06–2.04) and H-21 (δH 1.30).
The 1H NMR spectrum of compound 5 presents a hydroxymethine group at δH 3.19 (1H, dd, J = 11.5, 5.0, H-3). The large coupling constant (J = 11.5) suggests the axial orientation of proton H-3 and the β-configuration of the hydroxyl group at C-3. An oxygenated methine carbon (δC 78.93) was observed in the 13C NMR spectrum of compound 5, instead of a ketone group (δC 217.88) in compound 4. The 13C NMR spectra of compounds 6a–6e contain additional signals of carbonyl ester in the range of δC 176.63–170.98. The signals at δC 144.29–113.09 and δH 7.83–6.43 indicate the presence of aromatic rings in compounds 6d and 6e. A methyl carbonyl ester at δH 2.04 a pentanoyl group at δH 2.29
and an isobutyryl group at δH 1.17 (3H, d, J = 7.5), 1.16 (3H, d, J = 7.5) were determined from the 1H NMR spectra of compounds 6a, 6b and 6c, respectively.
To the best of our knowledge, the twelve oxime esters and seven lactones are new dammarane derivatives that have not been reported in the literature to date. The analysis of the HR-ESIMS spectra data of all these newly synthesised compounds provided in Section 3.2 (characteristics) also further confirms the elucidation of their structures.
Number | Sample/compound | IC50 value ± SD (μM) | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
a Acarbose (IC50 = 183.44 ± 3.78 μM) and NA: no activity (IC50 > 390 μM). | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 | 1 | 264.91 ± 9.64 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2 | 2 | 36.94 ± 2.80 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3 | 3a | NA | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
4 | 3b | 33.12 ± 5.37 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
5 | 3c | 16.38 ± 1.00 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
6 | 3d | NA | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
7 | 3e | 36.60 ± 1.26 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
8 | 3f | 23.98 ± 1.29 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
9 | 3g | NA | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
10 | 3h | NA | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
11 | 3i | 46.79 ± 5.33 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
12 | 3k | 59.15 ± 2.08 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
13 | 3l | 43.51 ± 1.68 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
14 | 3m | NA | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
15 | 4 | NA | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
16 | 5 | 2.73 ± 0.05 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
17 | 6a | 220.5 ± 4.77 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
18 | 6b | 31.95 ± 0.69 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
19 | 6c | 4.62 ± 0.12 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
20 | 6d | 228.09 ± 4.62 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
21 | 6e | 7.31 ± 0.11 |
The starting material dipterocarpol (1) was also tested for comparison with its synthetic derivatives. Acarbose was used as the standard compound in this assay. The results showed that dipterocarpol possesses weak activity (IC50 264.91 μM) compared with acarbose (IC50 183.44 μM). Nevertheless, most of its oxime ester derivatives (2, 3b, 3c, 3e, 3f, 3i, 3k, and 3l) exhibited better the bio-activity than the standard compound acarbose towards α-glucosidase, as shown by their IC50 values of 16.38–59.15 μM. The benzoate derivative (3c) is the most active compound in the oxime ester series with an IC50 of 16.38 μM. It is followed by the parachlorobenzoate (3f, IC50 of 23.98 μM) and cinnamate (3b, IC50 of 33.12 μM) derivatives. The substitution of monoiodine (3h) and diodine benzoate (3d) led to the loss of their activity. Compound 3i with meta-fluorobenzoate exhibited good activity with an IC50 value of 46.79 μM. However, the installation of para-fluorobenzoate (3g) resulted in no α-glucosidase inhibition. The above-mentioned inhibitory effect suggests that the presence of the halogen substitutions in the aromatic rings increase the overall inhibition in the order of Cl > F > I. The heterocyclic oxime ester with a nicotinic ring (3e) showed good activity with an IC50 of 36.60 μM; in contrast, the pyrazine derivative was observed to be inactive. The acetate (3k) and pentanoate (3l) derivatives were found to exhibit good activity with the corresponding IC50 of 59.15 and 43.51 μM, respectively. The results suggest that the chain length of the alkoxy substitution may improve the α-glucosidase inhibitory activity. The isobutyrate (3m) oxime ester was inactive in this test. This means that the bulkiness of the alkoxy oxime ester group seems to lower the activity overall.
Compound 4 with a lactone ring located at C-17 was inactive to the α-glucosidase inhibition; however, the reduction of the ketone to hydroxyl group at C-3 led to potential inhibitory activity. Compound 5 with a free hydroxyl group was registered as the most active compound (IC50 of 2.73 μM) among the ester derivatives, as shown in Scheme 2. This evidence suggests that the hydroxyl group at C-3 plays an important role in α-glucosidase inhibition. Sufficiently increasing the chain length and bulkiness of the alkanoate moieties may enhance the activity. This was evidenced by the increase in activity, which followed the order of: 6a (IC50 220.5 μM) < 6b (IC50 31.95 μM) < 6c (IC50 4.62 μM) (Table 1). The cinnamate derivative showed significantly better activity than the salicylate one. This observation suggests that the electron density of the benzene rings may have a negative correlation with the inhibitory activity of the host molecules towards α-glucosidase. Also, the dose–response curves of the most active inhibitors, i.e., 5, 6c, and 6e, are provided in the ESI (Section 3; pages 203–205†).
In summary, it is obvious that the most significant inhibitors towards α-glucosidase, which correspond to IC50 values lower than 10 μM, follow the order of 5 > 6c > 6e. In particular, 5 demonstrates pronounced efficacy with the corresponding value of 2.73 ± 0.05 μM, accounting for the concentration required for 50% inhibition of the enzyme. This is followed by 3c > 3f > 6b > 3b > 3e > 2, whose IC50 values are in the range of 10–40 μM, and this considered highly effective. This group includes a mildly effective inhibitor, namely 3e, which is expected to be very different according to the docking-based simulation. Also, the experiments found that 3a, 3d, 3g, 3h, 3m, and 4 exhibit no inhibition activity towards α-glucosidase; nevertheless, the simulation gave the opposite prediction given 3d. The inconsistencies observed are readdressed and discussed in detail later in this report.
C | 3W37 | 3AJ7 | PTP1B | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Site 1 | Site 2 | Site 3 | Site 4 | Site 1 | Site 2 | Site 3 | Site 4 | Site 1 | Site 2 | Site 3 | Site 4 | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
E | H | E | H | V | R | E | H | E | H | E | H | V | R | E | H | E | H | E | H | E | H | V | R | E | H | E | H | E | H | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
a C: compound; E: docking score energy (kcal mol−1); H: hydrogen bonds; V: van der Waals interactions; and R: root-mean-square deviation (Å). | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 | −8.1 | 1 | −10.9 | 2 | 14 | 1.24 | −9.2 | 2 | −9.8 | 2 | −10.2 | 2 | 15 | 1.24 | −9.4 | 1 | −8.7 | 1 | −9.3 | 1 | −10.5 | 1 | 10 | 1.19 | −8.2 | 0 | −9.1 | 1 | −7.4 | 0 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2 | −11.2 | 3 | −12.9 | 4 | 5 | 1.67 | −9.8 | 3 | −10.4 | 3 | −11.9 | 2 | 18 | 1.67 | −8.7 | 1 | −8.2 | 0 | −10.6 | 1 | −13.8 | 3 | 7 | 1.17 | −11.4 | 2 | −10.8 | 1 | −9.7 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3a | −8.1 | 1 | −9.7 | 2 | 17 | 1.89 | −8.7 | 1 | −9.1 | 1 | −9.5 | 3 | 27 | 1.89 | −8.2 | 2 | −7.3 | 1 | −7.9 | 1 | −9.3 | 2 | 11 | 1.25 | −7.2 | 1 | −8.1 | 1 | −7.3 | 0 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3b | −12.1 | 2 | −13.4 | 3 | 13 | 1.78 | −11.7 | 2 | −10.9 | 1 | −12.2 | 3 | 20 | 1.78 | −9.8 | 1 | −9.4 | 1 | −11.7 | 2 | −11.5 | 2 | 10 | 1.98 | −8.0 | 0 | −7.6 | 0 | −9.3 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3c | −8.2 | 1 | −10.6 | 2 | 11 | 1.29 | −9.0 | 2 | −9.3 | 2 | −11.8 | 3 | 25 | 1.29 | −10.1 | 2 | −8.4 | 1 | −9.6 | 2 | −9.6 | 2 | 6 | 1.92 | −8.8 | 1 | −7.5 | 1 | −6.7 | 0 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3d | −10.9 | 3 | −12.3 | 4 | 13 | 1.29 | −11.3 | 2 | −11.8 | 3 | −10.9 | 3 | 20 | 1.29 | −8.2 | 1 | −9.1 | 2 | −9.5 | 2 | −13.7 | 8 | 8 | 1.89 | −11.9 | 4 | −12.3 | 3 | −12.7 | 4 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3e | −12.4 | 3 | −13.8 | 6 | 11 | 1.12 | −11.6 | 3 | −12.6 | 4 | −15.3 | 5 | 22 | 1.12 | −12.6 | 3 | −13.1 | 3 | −11.9 | 2 | −12.9 | 4 | 8 | 1.43 | −10.7 | 2 | −11.3 | 3 | −9.6 | 2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3f | −7.9 | 2 | −10.5 | 2 | 19 | 1.94 | −8.3 | 1 | −9.5 | 1 | −12.2 | 3 | 21 | 1.94 | −9.7 | 1 | −10.1 | 1 | −11.3 | 2 | −10.1 | 3 | 8 | 1.75 | −8.9 | 1 | −9.1 | 2 | −7.5 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3g | −7.6 | 1 | −9.8 | 2 | 13 | 2.05 | −8.0 | 1 | −7.8 | 1 | −9.9 | 2 | 22 | 2.05 | −7.1 | 1 | −6.8 | 0 | −8.2 | 1 | −9.0 | 3 | 8 | 1.90 | −8.4 | 2 | −7.6 | 1 | −8.1 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3h | −8.2 | 1 | −9.9 | 2 | 15 | 1.92 | −8.3 | 1 | −8.7 | 2 | −10.4 | 2 | 13 | 1.92 | −9.7 | 1 | −8.3 | 1 | −9.6 | 1 | −9.4 | 3 | 5 | 1.36 | −7.0 | 1 | −8.4 | 1 | −8.9 | 2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3i | −10.6 | 3 | −11.3 | 4 | 13 | 1.81 | −9.6 | 2 | −10.2 | 3 | −11.7 | 2 | 23 | 1.81 | −10.2 | 1 | −9.3 | 1 | −8.0 | 0 | −10.3 | 3 | 9 | 1.14 | −9.0 | 1 | −8.9 | 1 | −9.4 | 2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3k | −9.8 | 2 | −11.6 | 3 | 12 | 1.03 | −8.9 | 1 | −10.1 | 2 | −12.7 | 5 | 14 | 1.03 | −10.8 | 2 | −11.0 | 3 | −9.5 | 2 | −9.7 | 2 | 7 | 1.35 | −7.6 | 1 | −8.0 | 1 | −7.3 | 0 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3l | −7.9 | 1 | −10.9 | 2 | 15 | 1.70 | −8.2 | 1 | −9.0 | 1 | −9.6 | 1 | 14 | 1.70 | −7.2 | 0 | −6.9 | 0 | −8.1 | 1 | −9.9 | 2 | 11 | 1.60 | −8.1 | 1 | −7.9 | 1 | −8.7 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3m | −7.7 | 0 | −9.1 | 1 | 14 | 1.95 | −7.3 | 0 | −8.2 | 1 | −10.3 | 3 | 21 | 1.95 | −8.5 | 1 | −9.0 | 2 | −7.6 | 1 | −9.3 | 1 | 9 | 1.34 | −6.9 | 0 | −7.8 | 1 | −6.3 | 0 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
4 | −7.5 | 1 | −8.5 | 1 | 7 | 1.87 | −6.7 | 0 | −7.1 | 0 | −10.2 | 3 | 14 | 1.87 | −9.2 | 1 | −8.1 | 1 | −7.5 | 1 | −11.2 | 4 | 6 | 1.23 | −9.2 | 2 | −8.1 | 1 | −10.7 | 2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
5 | −10.2 | 1 | −14.9 | 3 | 9 | 1.30 | −11.3 | 2 | −11.9 | 2 | −15.9 | 4 | 13 | 1.30 | −11.9 | 2 | −12.7 | 2 | −13.5 | 3 | −13.4 | 4 | 7 | 1.37 | −10.2 | 1 | −8.9 | 1 | −9.7 | 2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
6a | −8.1 | 0 | −9.7 | 1 | 14 | 1.19 | −8.8 | 1 | −7.9 | 0 | −11.6 | 3 | 19 | 1.19 | −9.2 | 2 | −8.9 | 1 | −10.2 | 2 | −9.8 | 2 | 8 | 1.40 | −7.3 | 1 | −6.9 | 0 | −7.0 | 0 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
6b | −9.3 | 1 | −10.1 | 2 | 14 | 1.22 | −8.9 | 1 | −9.6 | 1 | −10.6 | 2 | 20 | 1.22 | −8.8 | 1 | −7.9 | 1 | −7.6 | 0 | −11.5 | 4 | 9 | 1.85 | −10.1 | 2 | −10.9 | 2 | −9.5 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
6c | −11.9 | 1 | −13.2 | 2 | 14 | 1.12 | −10.7 | 1 | −11.3 | 2 | −13.7 | 3 | 21 | 1.12 | −11.3 | 2 | −9.1 | 1 | −10.6 | 2 | −15.8 | 4 | 9 | 1.99 | −12.1 | 2 | −10.8 | 2 | −10.2 | 1 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
6d | −10.8 | 2 | −13.7 | 4 | 10 | 1.31 | −11.4 | 2 | −12.1 | 2 | −12.5 | 4 | 20 | 1.31 | −9.7 | 2 | −10.5 | 2 | −10.9 | 3 | −10.7 | 3 | 5 | 1.70 | −9.2 | 1 | −8.9 | 1 | −9.7 | 2 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
6e | −12.6 | 3 | −15.2 | 4 | 11 | 0.31 | −11.2 | 2 | −12.1 | 2 | −14.4 | 3 | 16 | 0.31 | −10.8 | 1 | −9.3 | 1 | −8.1 | 1 | −10.3 | 2 | 11 | 1.54 | −8.9 | 1 | −12.2 | 3 | −9.5 | 1 |
Regarding the investigated structure representing α-glucosidase, viz., 3W37, the computational simulation is considered to be in good agreement with the experimental results collected from the bioassays. The docking-based computation predicted that the ligands with the highest inhibitability follow the order of 6e (DS −15.2 kcal mol−1) > 5 (DS −14.9 kcal mol−1) > 3e (DS −13.8 kcal mol−1) = 6d (DS −13.7 kcal mol−1) > 3b (DS −13.4 kcal mol−1) = 6c (DS −13.2 kcal mol−1). The significance is highly justified by the values calculated using a similar computing configuration for three commercialised drugs, i.e., acarbose (DS −14.2 kcal mol−1), voglibose (DS −11.1 kcal mol−1), and miglitol (DS −12.7 kcal mol−1), which were reported in our previous study.47 The corresponding RMSD values for the derivatives are all under 2 Å, indicating a biologically rigid conformation.48 Especially, that of 6e is 0.31 Å, meaning a significantly short average distance between its (backbone) atoms, and thus it can possibly be viewed as a quasi-crystallised structure. All the compounds were determined to be effective α-glucosidase inhibitors based on the assay experiments, except for 3e, which exhibited a moderate IC50 of 36.60 ± 1.26 μM in vitro although it was predicted to have the third lowest DS energy in silico. Furthermore, if the intermolecular interaction ability is considered, this ligand can create 6 strong hydrogen bonds, storing a total Gibbs free energy of −10.3 kcal mol−1 (ESI: Section 2.1 and Table S1†), with the in-pose amino acids of its targeted protein structure and should be highly conducive to the stability of the 3e-3W37 complex. Although observed in-discrete, this experiment-theory inconsistency is still notable, thus requiring further explanation. Alternatively, the compounds showing unnoticeable inhibition towards the enzyme in the in vitro investigation, viz., 3a, 3g, 3h, 3m, and 4, are also expected to form unstable ligand-protein complexes according to the static in silico study, given their significantly low DS values (over −10 kcal mol−1). This excludes 3d, whose IC50 value is over 390 μM (i.e., no activity) but its corresponding DS value is −12.3 kcal mol−1 (i.e., considerable stability). In addition, 6d was experimentally demonstrated to be a weak α-glucosidase inhibitor (IC50 of 228.09 ± 4.62 μM) but computationally predicted to exhibit good affinity towards the 3W37 structure (DS of −13.7 kcal mol−1). Thus, these discrepancies also require appropriate explanations. Otherwise, the docking-based static models correspond well with the assay-based inhibition experiments.
However, due to the information lost from the pre-docking conditions, we performed a QSARIS-based analysis on the physicochemical properties of the studied ligands with reference to Lipinski's rule of five. The results are summarised in Table 3.
Ligand (compound) | DSaverage (kcal mol−1) | Mass (amu) | Polarizability (Å3) | Dispersion coefficients | Hydrogen bond (3W37/3AJ7/PTP1B) | ||
---|---|---|---|---|---|---|---|
Log![]() |
Log![]() |
H-Acceptor | H-Donor | ||||
1 | −10.5 | 442.8 | 33.7 | 6.2 | −5.6 | 1/1/1 | 1/1/0 |
2 | −12.9 | 458.0 | 38.7 | 4.2 | −5.4 | 2/1/3 | 2/1/0 |
3a | −9.5 | 577.9 | 29.6 | 9.0 | −9.2 | 1/1/1 | 1/2/1 |
3b | −12.4 | 589.0 | 36.8 | 4.7 | −5.3 | 2/1/0 | 0/1/1 |
3c | −10.7 | 562.3 | 39.2 | 3.3 | −5.1 | 1/1/2 | 1/1/0 |
3d | −12.3 | 827.9 | 30.1 | 11.6 | −12.1 | 3/2/5 | 0/1/0 |
3e | −14.0 | 598.4 | 36.9 | 4.1 | −5.2 | 2/1/2 | 1/1/0 |
3f | −10.9 | 597.9 | 41.0 | 3.7 | −5.0 | 1/0/3 | 0/3/0 |
3g | −9.6 | 598.7 | 32.8 | 9.1 | −9.6 | 1/1/3 | 1/1/0 |
3h | −9.9 | 687.4 | 31.7 | 9.8 | −10.2 | 1/0/2 | 1/1/1 |
3i | −11.1 | 578.3 | 32.8 | 5.3 | −6.6 | 4/0/3 | 0/2/0 |
3k | −11.3 | 499.5 | 33.5 | 4.9 | −5.4 | 2/2/1 | 1/2/1 |
3l | −10.1 | 542.8 | 34.8 | 4.6 | −5.3 | 1/0/2 | 1/0/0 |
3m | −9.6 | 528.6 | 31.1 | 8.6 | −8.4 | 1/2/0 | 0/1/0 |
4 | −10.0 | 416.8 | 29.8 | 6.4 | −6.6 | 1/2/3 | 0/0/1 |
5 | −14.7 | 416.2 | 52.4 | 3.1 | −3.8 | 3/2/4 | 0/0/0 |
6a | −10.4 | 458.9 | 33.3 | 4.5 | −5.7 | 1/1/2 | 0/0/0 |
6b | −10.7 | 501.2 | 37.2 | 3.9 | −4.8 | 2/1/3 | 0/1/1 |
6c | −14.2 | 490.1 | 48.8 | 3.6 | −3.7 | 1/3/3 | 1/0/0 |
6d | −12.3 | 537.8 | 37.3 | 5.1 | −5.7 | 2/3/2 | 1/1/1 |
6e | −13.9 | 549.2 | 51.6 | 3.8 | −3.5 | 3/1/1 | 0/0/1 |
The first criterion regarding molecular mass (<500 amu) can be used to explain the docking-assay inconsistency. The corresponding values for ligands 3d and 3e are 827.9 and 598.4 amu, respectively. Their significant mass implies that even though they can exhibit inhibition towards the 3W37 structure when in the pose, their heavy mass is more likely to deter their transport from the medium of either the in vitro experiments or real physiological environment. A similar approach can explain the mild discrepancy observed in the case of 6d due to its slightly heavy mass of 537.8 amu, cf., Lipinski's criterion. In contrast, although having a lower-bound mass weight (mass of 461.8 amu), 4 exhibited no inhibition activity (IC50 > 390 μM) in the bioassay-based experiments, which is likely because of its low in-pose stability (DS of −8.5 kcal mol−1). Otherwise, most of the assay-based inactive compounds, viz., 3a, 3d, 3g, 3h, and 3m, also possess a molecular mass over 500 amu, thus justifying this correlation. However, although 6e is predicted by the static model in the docking simulation to be the most promising inhibitor (DS of −15.2 kcal mol−1) towards α-glucosidase, the experimental observations reveal that it is only third-ranked (IC50 of 7.31 ± 0.11 μM), after 5 (IC50 of 2.73 ± 0.05 μM) and 6c (IC50 of 4.62 ± 0.12 μM). In this case, the differences can be explained based on the slightly high mass of 6e, i.e. 549.2 amu, compromising its capability to form a highly stable 6e-3W37 inhibitory system, but not prohibitively. Besides molecular mass, according to Lipinski's criteria, a good membrane-permeable molecule should satisfy the requirements of no more than 5 groups for hydrogen bonds; no more than 10 groups receiving hydrogen bonds; and a value of logP of less than +5 (log
P < 5). Hence, it is obvious that the three most promising candidates, viz. 5, 6c, and 6e, are ready for further development in pharmaceutical applications, especially for oral administration. In addition, their polarisability constants are significant, i.e. 52.4 Å3 for 5, 48.8 Å3 for 6c, and Å3 51.6 for 6e, likely inducing the formation of molecular dielectric moments.49 This property is thought to be highly conducive to both inhibitory effects towards a protein structure given that the polypeptide molecule is made of polarised amino acids and fluidity in hydrophilic environments of the body.
In summary, the molecular docking technique coupled with analysis of the QSARIS-based physicochemical properties can provide a reasonably accurate model to predict ligand-protein inhibitability in general, and dammarane-α-glucosidase inhibition in particular. The latter determines the mobility of the ligands, whereas the former provides the stability of the ligand–protein complexes. This theoretical model suggests that the Dipterocarpus alatus-derived compounds are promising drug-like inhibitors towards α-glucosidase, following the order of 5 (mass: 416.2 amu; polarisability: 52.4 Å3; DS: −14.9 kcal mol−1) > 6c (mass: 490.1 amu; polarisability: 48.8 Å3; DS: −13.7 kcal mol−1) > 6e (mass: 549.2 amu; polarisability: 51.6 Å3; DS: −15.2 kcal mol−1), which is in excellent agreement with the experiments. However, this model is still considered premature, thus requiring more solid-evidenced validations. Molecular dynamics simulation can monitor the simultaneous “behaviour” of all atoms in the systems, thus validating the mechanism of inhibition, and surface plasmon resonance can identify the correct ligand-protein complex structures.
Relatively, the inhibitability of the studied ligands towards the other diabetes-related protein structures, i.e., 3AJ7 and PTP1B, was predicted without experimental evidence. Firstly, the former is still expected to be most vulnerable at its site 1 and with the strongest interactions derived by 5 (DS: −15.9 kcal mol−1), 6c (DS: −13.7 kcal mol−1), and 6e (DS: −14.4 kcal mol−1). Although the 3e-3AJ7 inhibitory system has favourable stability, given its DS value of −15.3 kcal mol−1, number of hydrogen bonds of 5, and number of van der Waals interactions of 22, its formation its expected to be significantly hindered due to the heavy mass (598.4 amu) and low polarisability (36.9 Å3) of the ligand. Similarly, the formation of the 2-PTP1B and 3d-PTP1B complexes, especially the latter, is almost unexpected despite the favourable stability of their intermolecular interaction. Besides, the dispersion coefficient (logP) of 3d is 11.6, which is much higher that the Lipinski's threshold, suggesting its insignificant dispersity. Meanwhile, 6e does not seem to be a promising inhibitor towards the PTP1B structure given its corresponding moderate DS values, regardless of the inhibited sites. Therefore, the preliminary analysis suggests that 5 and 6c are versatile inhibitors regarding all the diabetes-related protein structures, i.e., 3W37, 3AJ7, and PTP1B. The theoretical implications can be verified by appropriate experimental investigations, e.g., enzyme inhibition assays or clinical trials.
Finally, the 3D morphology and 2D interaction map of all the ligand-protein inhibitory complexes are visually presented in Fig. 5–7. In the design of specific protein inhibitors, the size- and shape-compatibility between the ligands and the protein sites is also important. Although the Dipterocarpus alatus-derived compounds are projected to be size-fitting (from the 3D observation) and shape-complementary (by 2D continuous proximity contours) with the in-pose topographical features of the diabetes-related protein structures, the sites are unlikely to be conducive to either the entry of significantly larger candidates or simultaneous inhibition. Therefore, a size- and medical-reference for the further development of pharmaceutical medicine are also given. In addition, the configurations of the in-pose interactions were projected in detail.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra04461c |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2021 |