Isolation, semi-synthesis, docking-based prediction, and bioassay-based activity of Dolichandrone spathacea iridoids: new catalpol derivatives as glucosidase inhibitors

Dolichandrone spathacea iridoids are promising anti-diabetic inhibitors towards α-glucosidase protein (PDB-3W37) and oligo-1,6-glucosidase protein (PDB-3AJ7). Five catalpol iridoids (1, 2, 10, 13, 14) were isolated from mangrove plant D. spathacea, and their derivatives (3, 4, 5, 6, 7, 8, 9, 11, 12, 15) were obtained from reduction, acetylation, O-alkylation, acetonisation, or hydrolysation starting from naturally isolated compounds. They were identified by spectral methods such as IR, MS, and 1D and 2D NMR. Their glucosidase-related (3W37 and 3AJ7) inhibitability and physiological compatibility were predicted by molecular docking simulation and prescreened based on Lipinski's rule of five. Experimental α-glucosidase inhibition of 1–15 was evaluated using enzyme assays. Compounds 3, 4, 5, 6, and 9 are new iridoid derivatives, introduced to the literature for the first time, while all fifteen compounds 1–15 are studied for molecular docking for the first time. Regarding protein 3W37, the five strongest predicted inhibitors assemble in the order 2 > 10 > 1 > 9 > 14. In respect to 3AJ7, the corresponding order is 14 > 2 > 10 > 5 > 1 = 9. Lipinski's criteria suggest 10 as the candidate with the most potential for oral administration. The in vitro bioassay revealed that compound 10 is the most effective inhibitor with a respective IC50 value of 0.05 μM, in the order 10 > 2 > 14 > 13 > 1. The computational and experimental results show good consistency. The study opens an alternative approach for diabetes treatment based on inhibitability of natural and semi-synthesised catalpol iridoid derivatives towards carbohydrate-hydrolases.


Introduction
Diabetes mellitus (DM) has become prevalent in the past decades, especially in the middle-income countries. According to the International Diabetes Federation, it was responsible for about 4.2 million deaths in 2019 1 and is predicted to be the seventh leading cause of death by 2030 (WHO). 2 Also, there is evidence suggesting a relationship between this disease and the causes of cardiovascular diseases, blindness, kidney failure, stroke, and limb amputation. 3 In particular, type 2 DM is a noninsulin-dependent disorder, stemmed from the ineffective use of insulin in the body and accounting for 90-95% of total diabetes cases. 4 There are many kinds of medicines for the treatment of type-2 diabetes, based on insulin sensitizers, sodium glucose transporter-2 (SGLT-2) inhibitors, glucagon like peptidase-1 (GLP-1) analogues and dipeptidyl peptidase-4 (DPP-IV) inhibitor mechanisms. However, these therapies have signicant side effects such as weight gain, cardiovascular risks, and atherosclerosis as well as limited tolerability. Therefore, the approach for type 2 DM changes into the inhibition of glucosidase enzymes, which slows down of glucose absorption through intestinal wall. There are two main types of glucosidase responsible for carbohydrate break-down processes in the body, regarding their sources: a-glucosidase and oligo-1,6glucosidase. First, a-glucosidase is an exoenzyme found in animal, plant, bacterial, or fungal organisms, 5 especially from sugar beet seeds. 6 This type of glucosidase was demonstrated only catalysing hydrolysis of a-(1/4) and a-(1/6) bonds to yield monosaccharides. 7 Protein structural data of a-glucosidase can be referenced at Worldwide Protein Data Bank under entry PDB-3W37 (DOI: 10.2210/pdb3W37/pdb). Second, oligo-1,6-glucosidase, oen called isomaltase, is a debranching endoenzyme, hydrolysing only the a-1,6 linkage in starch and glycogen to produce sugars with an a-conguration. 8 Unlike aglucosidase, they are not responsible for a-1,4 linkage hydrolysis. Oligo-1,6-glucosidase is mainly found in animals although some bacterial species, such as Bacillus cereus, can also synthesise it. 9 In humans, it is located on the small intestine brush border. 10 Information for isomaltase crystal structure can be downloaded directly from Worldwide Protein Data Bank database under entry PDB-3AJ7 (DOI: 10.2210/pdb3AJ7/pdb). Therefore, 3W37 and 3AJ7 are considered as highly promising drug targets for effective treatment of type-2 diabetes, especially medication through oral administration. In principle, multiple inhibition of these enzymes (3W37 and 3AJ7) is considered as an efficacious strategy to suppress the hyperglycemia and improve insulin sensitisation, simultaneously. Crystal structure of the proteins are shown in Fig. 1.
In silico technique, based on computational simulation and computing calculation, is currently seeing a gain of their popularity for medical science as prescreening research. It is for reducing cost and time of wet laboratory experiments by predicting the compounds with undesirable properties and the most promising candidates. The former substances are deemed to be eliminated from the next analysis or further-developed research, while the latter justies the selection. 11 Regarding ligand-protein interaction, molecular docking simulation is an effective method to investigate the potency of a ligand as an inhibitor towards its targeted protein. The method can estimate ligand-target binding energy and intermolecular interaction, thus predicting static stability of the inhibitory systems. 12 Effectively inhibited by external ligands, the enzyme is likely to undergo conformational changes, thus loss of enzymatic functionality inevitably ensuing. This might help mitigating the amount of glucose catalytically synthesised and excreted to bloodstream. According to the methodology of molecular docking simulation applied for ligand-protein inhibitory systems in general and to the evaluation from.
Molecular Operating Environment (MOE)-based algorithms in particular, an associated docking score (DS) is considered as the main indicator for inhibitory effects, of which value lower than À3.2 kcal mol À1 indicates good binding capacity. [13][14][15] In principle, the gure is yielded by the free-energy sum of all individual intermolecular interactions, of which affinity stems from hydrophilic bonding, i.e. various hydrogen-bond types, and hydrophobic binding, i.e. van der Waals forces. Besides, a root-mean-square deviation (RMSD) value over 3Å reveals that the inhibition expects failure; meanwhile, the threshold of docking success is widely acceptable if #2Å. 16 Also, visual illustration for inhibitory morphology and interaction description is oen provided. Specication to drug design for diabetes, molecular docking technique demonstrated its contributive signicance regarding a variety of ligand families, e.g. thiazolidinedione, 17 pyrazole-rhodanine, 18 or wedelolactone. 19 Catalpol and catalpol analogues, containing epoxy ring at C-7 and C-8, are widely found in many plants belonging to some families such as Plantaginaceae, Lamiaceae, Bignonicaceae. 20 Recently, there is an increasing numbers of papers reported about the anti-diabetic efficacy of catalpol derivatives in the management of diabetes and its complications. 20,21 It has been showed that catalpol-type iridoids are able to restore diabetic liver function by improving glucose and lipid metabolism, ameliorating oxidative stress, and restoring mitochondrial function. They improve glucose uptake and lipolysis by upregulating the expression of glucose transporter as well as downregulating the expressions of peroxisome proliferators in mice. 20 Therefore, catalpol derivatives have become a very promising scaffold for the development of anti-diabetic drug candidates. Our previous study on chemical constituents of the mangrove plant Dolichandrone spathacea led to the isolation of four catalpol iridoids. 22 This seems reasonable because according to folk experiences, the plant is currently reported as a natural symptom reliever for diabetes, especially in Mekong Delta, Vietnam. Nevertheless, its chemical properties and medical potentiality still remain almost untouched possibly because of their narrow distribution.
In this paper, ten derivatives from the isolated compounds were synthesized, together with natural compounds they were  studied in molecular docking experiment on sugar beet aglucosidase protein (PDB-3W37) and oligo-1,6-glucosidase protein (PDB-3AJ7). Three commercial medicines acarbose (D1), miglitol (D2), and voglibose (D3) had been reported able to lengthen the duration time of carbohydrate absorption and atten the blood-glucose concentrations, thus widely used in the treatment for patients with type 2 diabetes. 23 Therefore, they were specied as references in this study. Their structural formula is presented in Fig. 2 13 C and DEPT) and 2D (HSQC, HMBC, COSY, NOESY) NMR spectra were recorded on a Bruker Avance 500 MHz using tetramethylsilane (TMS) as the internal standard for 1 H and solvent signal for 13   To a stirred solution of 1 or 2 (0.1 mmol) dissolved in MeOH (5 mL) was added NiCl 2 $6H 2 O (0.5 mmol). The reaction mixture was stirred for 10 minutes at room temperature and cooled to 0 C. NaBH 4 (18.9 mg, 0.5 mmol) was slowly added to the mixture, which was then allowed to stir at 0 C for 30 minutes and one additional hour at room temperature. The reaction mixture was then removed the solvent to obtain a white solid. The solid was dissolved in CH 2 Cl 2 and the solution was washed with water (2 Â 30 mL), dried over anhydrous Na 2 SO 4 , and evaporated the solvent. The crude product was puried over a silica gel column (CH 2 Cl 2 /MeOH 9 : 1) to obtain 3 (46 mg, yield: 83%) or 4 (44.5 mg, yield: 80.6%).
2.1.3.3 Synthesis of 6. To a solution of 2 (55.2 mg, 0.1 mmol) in dried CHCl 3 (3 mL) and pyridine (0.5 mL) was added triphenylmethyl chloride (55.6 mg, 0.2 mmol). The reaction mixture was allowed to stir for 12 hours at room temperature. The solvent was then removed under reduced pressure to give white solid which was dissolved in CHCl 3 (30 mL) and washed with water (2 Â 30 mL). The organic layer was dried over anhydrous Na 2 SO 4 , evaporated the solvent to afford a yellow solid, which was puried by silica gel column chromatography (CH 2 Cl 2 /MeOH/pyridine 98 : 2 : 0.5) to obtain 6.
2.1.3.4 General procedure for the synthesis of 7, 8, 11, 12. To a solution of 1 or 2 or 13 or 14 (0.1 mmol) in pyridine (1 mL) was added anhydride acetic acid (1 mL). The solution was stirred for 24 h at room temperature then ice cold water (10 mL) was added and extracted with EtOAc (3 Â 30 mL). The EtOAc extract was combined, dried with anhydrous Na 2 SO 4 , evaporated the solvent to give yellow solid. The crude product was subjected to a Sephadex LH-20 column (100% MeOH) to obtain the corresponding acetylated catalpol iridoids.
2.1.3.6 Synthesis of 15. Compound 2 (110.4 mg, 0.2 mmol) was dissolved in a mixture of (H 2 O/THF/EtOH 1 : 1 : 1) and TEA (1 mL). Hydroxyl amine hydrochloride (69.5 mg, 1 mmol) was then added to the solution. The reaction mixture was heated to 100-110 C for 6 hours. The solvent was evaporated to give a residue which was then suspended in water and partitioned with EtOAc (3 Â 20 mL). The organic layers were combined, dried over anhydrous Na 2 SO 4 and evaporated the solvent to afford a white solid, this was subjected to a silica gel column chromatography (CH 2 Cl 2 /MeOH 9 : 1) to obtain 15 (65.2 mg, 90%).
2.1.4 Enzyme assay of a-glucosidase. Each assay was performed on 96 well-plates (200 mL per well). The test samples were dissolved in DMSO 10% and further diluted to the concentrations: 1024, 256, 64, 16, 4, 1 mg.mL À1 . Each concentration of the test sample (10 mL) was incubated at 37 C with phosphate buffer 100 mM (pH 6.8, 40 mL), a-glucosidase 0.4 U mL À1 (25 mL), p-nitrophenyl a-D-glucopyranoside 2.5 mM (25 mL). Aer 30 minutes, the reaction was terminated by adding 100 mL of 0.1 M Na 2 CO 3 solution. The enzyme activities were determined by measuring the absorbance at l max 410 nm. The control was prepared using the same procedure replacing the test sample by distilled water, while activity of the reference was tested by replacing the test sample with acarbose 30 Percentage of inhibition was calculated as followed: 31 a-glucosidase inhibition ð%Þ ¼ Abs control À Abs sample Abs control Â 100%

Molecular docking simulation.
The simulation was carried out on MOE 2015.10 soware. The structural information of targeted proteins (3W37 and 3AJ7) and docked ligands (iridoid compounds and referenced drugs) were required as the precursors. As conducted, this served as inputs to simulate intermolecular interaction between the agents, thereby yielding information on the respective ligand-protein inhibitory complex. The output included inhibitory congurations, docking score (DS) energy, root-mean-square deviation (RMSD), interaction types, and distances between ligands and proteins. In a typical procedure, a molecular docking simulation follows three steps. [13][14][15]32 2.2.1.1 Pre-docking preparation. Structural data of carbohydrate-hydrolase proteins 3W37 (DOI: 10.2210/ pdb3W37/pdb) and 3AJ7 (DOI: 10.2210/pdb3AJ7/pdb) were downloaded from Worldwide Protein Data Bank. The proteins and their 3D protonation were prepared by MOE QuickPrep functionality. Determination of protein active zones based upon a radius set at 4.5Å from their amino acids and the inhibitory ligands. The preparatory protein structures obtained were saved in format *.pdb, ready for docking simulation. Independently, the ligands, including iridoid compounds and referenced drugs, were structurally optimised by Conj Grad for energy minimisation. The optimisation conguration was set: termination for energy change ¼ 0.0001 kcal mol À1 ; max interactions ¼ 1000; charge modication by Gasteiger-Huckel method.
2.2.1.2 Docking investigation. Simulation on intermolecular interaction between the investigated agents was performed on MOE 2015.10 and the obtained inhibitory structures were saved in format *.sdf. Docking simulation parameters were congured as: poses retaining for further inhibition analysis ¼ 10; maximum solutions per iteration ¼ 1000; maximum solutions per fragmentation ¼ 200.
2.2.1.3 Post-docking analysis. The inhibitability of a certain duo-system was primarily predicted by docking score (DS) energy, which represents binding affinity of ligands and their targeted proteins in the site-site distance. Intermolecular interactions formed between the ligands and in-pose amino acids of the proteins were also probed. These include hydrophilic binding and hydrophobic interaction. The former is a variety of hydrogen-bond types, e.g. electron-transferring (Hacceptor/donor), cation-arene (H-p), arene-arene (p-p), and ionic; while, van der Waals forces constitute the latter. The results provide in-bonding amino acids, bonding lengths, and their Gibbs free energy. Besides, binding stability of the complex static conformation was evaluated by the respective values of root-mean-square deviation (RMSD) calculated. This is based on the fact that RMSD represent the average between neighbouring atoms; therefore, a smaller value means a more tightly bound conformation is formed. In addition, ligand conformation and orientation in its inhibited-protein active site was visualised on 2D and 3D planes.
2.2.2 Physicochemical and pharmaceutical compatibility. The docking parameters including DS average (kcal mol À1 ), molecular mass (Da), polarizability (Å 3 ) and volume or size (Å), and dispersion coefficients (log P and log S) were achieved by Gasteiger-Marsili method using QSARIS system. 33 They were then prescreened in an attempt to evaluate their orally pharmacological compatibility. This was based on Lipinski's rule of ve, a well-known set of indicators to predict drug-likeness. 34 According to Lipinski's criteria, a well membrane-permeable molecule should satisfy the requirements: (1) molecular mass < 500 Da; (2) no more than 5 groups for hydrogen bonds; (3) no more than 10 groups receiving hydrogen bonds; (4) the value of log P is less than +5 (log P < 5). 35,36 3. Results and discussion

Synthesis of catalpol derivatives
The synthetic derivations are schematically summarised in Fig. 3. During our investigation on the phytochemical constituents of Dolichandrone spathacea, ve iridoid compounds (1, 2, 10, 13, 14) have been isolated and structural elucidated. Four of them (1,2,13,14), with the exception of 10, due to small amount, were converted into their derivatives by acetylation, reduction, O-alkylation, acetonisation and hydrolysation to obtain 3, 4, 5, 6, 7, 8, 9, 11, 12 and 15. Initially, the reduction of the olen group in cinnamoyl and iridoid moiety in the structure of 1 and 2 was designed. The reduction was carried out by different kind of reductive reagents such as NaBH 4 /MeOH, LiAlH 4 /MeOH. However, many unwanted products were afforded under these condition. This might due to the strong base condition formed during the working up, causing a damage of iridoid ring. The trial of reduction with hydrogen catalyzed by Pd/C did not work, even though the reaction was heated to high temperature. We therefore used an additional catalyst Lewis acid NiCl 2 $6H 2 O in combination with NaBH 4 reagent. The reaction was carried out under mild condition at 0 C for 30 minutes. In this case, the products 3, 4 were obtain in good yield (83% and 80.6%, respectively). Interestingly, this reductive condition gave regioselective reduction only at olenic bond of cinnamoyl group, but not at the D-3,4-double bond of iridoid ring. The dihydro derivative 4 was continuously acetylated with acetic anhydride in pyridine, giving the corresponding pentaacetate 5 in 80% yield. In order to investigate the effect of free hydroxyl groups on hyperglycemic activity in catalpol iridoids, compound 2 was O-alkylated with triphenylmethyl chloride and compound 1 was reacted with 2,2-dimethoxypropane/ TsOH$H 2 O to provide product 6 and 9, respectively. The trityl groups were attached only to the primary oxymethylene groups at C-10 in catalpol ring and at C-6 0 in glucose unit as well. In case of 9, cyclic ketal was formed only at C-4 0 and C-6 0 of the glucose unit, due to the steric hindrance of the iridoid and cinnamoyl moiety. Iridoids 1, 2, 13 and 14 were reacted with acetic anhydride in pyridine to give the corresponding acetyl products 7, 8, 11 and 12. To nd out the role of cinnamoyl moiety, this moiety was tried to remove under different basic condition such as NaOH, LiOH, or KOH. However, these conditions gave many unwanted products. Therefore, an another milder basic condition using NH 2 OH. HCl in TEA was applied to obtain the hydrolysation derivative 15 (catalpol) from compound 2.

Structural elucidation
Elucidations of the structures for compounds 1, 2 and 13, 14 using spectral methods as FT-IR, MS, HR-MS, NMR (1D and 2D) were reported in ref. 24-29. In-detail data for the characterisations of each compound is provided in ESI † for in-depth reference and cross-checking.
Structural elucidation of compound 10: The 1 H and 13 C NMR spectra of 10 exhibit typical signals for a catalpol moiety quite similar to that of compounds 1, 2, 13 and 14. The side chain seems to be a monoterpenoid with ten carbon signals in the 13 Fig. 4. To the best of our knowledge, the structures elucidated regarding compounds 3, 4, 5, 6, and 9 are in respect of new semi-synthesised iridoid-derivatives that found having unreported in the literature, given our reachable referencing. The HR-ESIMS spectra of all the new synthetic compounds in Section 3.1 (characterisation) also further conrmed the elucidation of their structure.

Simulation of enzyme inhibitability
Static interaction between the studied inhibitory agents (iridoid compounds and referenced drugs) and the targeted proteins (3W37 and 3AJ7) was simulated by molecular docking technique. List of descriptive symbols are given in ESI † (Section 2.3).
3.4.1 Molecular docking simulation for iridoid-3W37. Regarding 1, 2, 3, 4, 5, 9, 10, 14, and 15, the docking results related to their inhibitory complexes with protein 3W37 are summarised in Table 1 (in-detail parameters included in ESI † -Section 2.4) and the complex structures are virtually represented in Fig. 5. The former provides inhibition parameters (DS value and RMSD), and in-depth bonding details building up the host complexes; while, the latter visualises the optimised orientation of each ligand in an active site of protein 3W37. DS value for 2-3W37 registers À15.1 kcal mol À1 , the lowest value calculated. This means the complex is at the lowest free-energy level, thereby most stable. In detail, compound 2 seems to exhibit great affinity with a specic amino acid in 3W37 structure, aka. arginine, by creating seven ionic and four H-acceptor bonds with four different arginine molecule, i.e. Arg 814, Arg 670, Arg 676, and Arg 699, and one ionic bond with Asp 666. They altogether compose À39.2 kcal mol À1 in respect of hydrophilic attraction, of which bonding lengths are in 2.78-3.95Å. This stability is followed by those of 10-3W37 and 1-3W37 with the DS values À14.7 and À14.1 kcal mol À1 , respectively, both composed by seven hydrogen bonds to a variety of amino acid types. However, their hydrophilic energy distribution sees different patterns. While each 10-3W37 ligand-amino acid hydrogen bond contributes rather evenly to the total hydrophilic free energy of À11.9 kcal mol À1 , C-Agr 670 H-acceptor bonding holds À8.7 kcal mol À1 of free energy accounting for over a ha of 1-3W37 total gure, À17.2 kcal mol À1 . The former contains the binding performances in distance 2.71-3.30Å, while the corresponding gures for the latter are slightly elevated, i.e. 2.75-3.80Å. This also justies the slight outweighing of 10-3W37 stability, previously evaluated by DS values. On the other hand, hydrophobic binding contributes fourteen van der Waals interactions to each 2-3W37, 10-3W37 or 1-3W37 inhibitory structure. Otherwise, 3-3W37, 4-3W37, 5-3W37, 9-3W37, 14-3W37, and 15-3W37 seem less stable given their inferior DS values, i.e. À12.9, À12.3, À11.6, À13.3, À13.0, and À12.3 kcal mol À1 . Therefore, the inhibitability of the regarded compounds accords with an order: 2 > 10 > 1 > 9 > 14 > 3 > 4 ¼ 15 > 5. In addition, all RMSD values registering under 2Å indicates that the associated complexes are tolerable as biological rigid bodies. In particular, although exhibiting a moderate DS value, the average distance between neighbouring atoms in 15-3W37 structure is signicantly short given its RMSD value of 0.67Å. This suggests a ligand-bound conformation if the inhibitory structure is formed. The reason might relate to ligand-pose spatial complementarity. Also, illustrative description for the intermolecular interactions detected in each active site is projected on a 2D diagram for each complex. It is also clear that all iridoid ligands are in high exposure to their targeted protein shown by densely blurredviolet zones, while long and continuous proximity contour bounding large in-pose areas likely indicate unintermittent interactability between the docking agents. Finally, visual observations reinforce the inhibitability since the docked sites are still spacious aer the inhibition of 1, 2, 3, 4, 5, 9, 10, 14, or 15.
In respect of 7, 8, 11, 12, and 13, DS values associated with their ligand-3W37 inhibitory complexes register over À10.5 kcal mol À1 , thus deemed as less competent inhibitors towards the protein. The respective calculated RMSD varies between 1.41 and 1.84Å, acceptable values for a rm conformation. Furthermore, in-depth results of the inhibitory complexes are also summarised in Table 1 while their corresponding 3D-and 2D-virtual renderings are included in ESI † (Section 2.1). In exceptional, 6 was found not able to perform inhibition towards protein 3W37.
3.4.2 Molecular docking simulation for iridoid-3AJ7. Regarding 1, 2, 3, 4, 5, 9, 10, 14, and 15, their inhibitory complexes with protein 3AJ7 are virtually presented in Fig. 6 and their respective docking results are summarised in Table 2 (indetail parameters included in ESI † -Section 2.4). The former predicts the ligand conformation and orientation in an active site of the protein; while, the latter provides in-depth docking parameters of the complexes, including DS value, RMSD, and interaction detail. The highest DS value obtained is À14.5 kcal mol À1 , responsible for the stability of 14-3AJ7 by representing the lowest free-energy level. However, there are only ve hydrophilic interactions detected as the complex forms, including four ionic and one H-donor bonds. The targeted amino acids are aspartic acid (Asp 242, Asp 352) and glutamic acid (Glu 277), which create hydrogen-bond binding with 14 carbon atoms in distance 2.93-3.82Å. The total hydrophilic potential calculated is À9.9 kcal mol À1 , insignicantly. Hence, 14-3AJ7 stability seems accumulated from van der Waals forces arisen between its two components as they include twenty different 3AJ7 amino acids in involvement. This is followed by 2-3AJ7 of which stability is expressed via À13.5 kcal mol À1 of DS value. The value is likely cultivated from both hydrophilic and hydrophobic contribution. There are nine hydrogen bonds formed between molecule 2 and 3AJ7 amino acids in distance 2.58-4.0Å, accounting for À28.8 kcal mol À1 of Table 2 Molecular docking simulation results for inhibitory complexes between the compounds and the protein 3AJ7 with amino acids: 1-3AJ7, 2-3AJ7, 3-3AJ7, 4-3AJ7, 5-3AJ7, 6-3AJ7, 7-3AJ7, 8-3AJ7, 9-3AJ7, 10-3AJ7, 11-3AJ7, 12-3AJ7, 13-3AJ7, 14-3AJ7, and 15-3AJ7 a the former free energy; while, the latter is from van der Waals attraction constituted by een different in-pose amino acids. In contrast, 10-3AJ7 binding affinity is likely based more on hydrophobic interactions than on hydrophilic counterparts as the former includes twenty-one amino acids while the latter results in a total energy just À7.5 kcal mol À1 . Its DS value registers À12.4 kcal mol À1 , the third highest gure. Otherwise, DS values of 1-3AJ7, 3-3AJ7, 4-3AJ7, 5-3AJ7, 9-3AJ7, and 15-3AJ7 are À11.4, À10.6, À10.2, À11.7, À11.5, and À11.0 kcal mol À1 , respectively. Therefore, the ligand inhibitability towards protein 3AJ7 follows an order: 14 > 2 > 10 > 5 > 1 ¼ 9 > 15 > 3 > 4. In addition, by possessing RMSD value under 2Å, all the inhibitory complexes are accepted performing biological rigidity. Moreover, visual illustration of interacting activity between the ligand and in-pose amino acids clarify the predominance of hydrophobic interactions in the establishment of ligand-3AJ7 inhibitory structures. This is justied by a myriad of amino acids presenting around the docked compounds. Resembling to observations on ligand-3W37, spaciousness inside the active sites are highly conducive to motion of 1, 2, 3, 4, 5, 9, 10, 14, or 15 when inhibition processes are on, thus their inhibitability towards protein 3AJ7 affirmed.
In respect of 7, 8, 11, 12, and 13, DS values associated with their ligand-3AJ7 inhibitory complexes are also included in Table 2. The respective gures are À9.3, À8.8, À6.4, and À7.5 kcal mol À1 , clearly expressing inferiority of overall stability. Their RMSD values are all under 1.99Å. Also, their corresponding 3D morphology and 2D interaction map are given in ESI † (Section 2.2). Similar to investigation on ligand-3W37, failure was also recorded in an attempt to establish a 6-3AJ7 inhibitory structure. No possible intermolecular interaction is given, thus no inhibitability of molecule 6 towards protein 3AJ7 ensuing. The reason for both could be referred to spatial restraint due to 6 structural bulkiness and branching, which possibly prohibit its entry into the host active site. However, more in-depth investigations, either theoretical or experimental, are needed to reach a more accurate conclusion.
3.4.3 Molecular docking simulation for drug-protein complexes. The inhibitability of three commercial medicines acarbose (D1), miglitol (D2), and voglibose (D3) towards the carbohydrate-hydrolase proteins (3W37 and 3AJ7) were investigated for referencing purposes. The inhibitions are virtually rendered in Fig. 7, while in-detail docking results are summarised in Table 3 (in-detail parameters included in ESI † -Section 2.4). All drugs seem fully tting with their inhibiting structural topography given their continuous proximity contours. This indicates high degree of complementarity. Even quasi-macromolecule D1 still fulls this topographical tting. However, different duo-systems see different patterns of intermolecular interaction between the adducts. In particular, D1 is able to create a myriad of interactions with its in-pose amino acids regardless of either the host protein or types of binding, thus conducive to its prominent inhibitability in comparison to the others. DS values of D1-3W37 and D1-3AJ7 are À14.2 and À15.3 kcal mol À1 , respectively. These were followed by the corresponding gures for D3-protein structures, i.e. À12.7 and À13.6 kcal mol À1 , respectively. Calculation on D2-protein inhibition results in relatively inferior DS values, i.e. 11.1 kcal mol À1 for D2-3W37 and À10.7 kcal mol À1 for D2-3AJ7.
Nevertheless, these complexes are considered possessing ligand-bound conformations given their short average distance between in-active atoms represented by RMSD of 0.68 and 1.03 A. Other RMSD values register under 2Å and most binding length are predicted under 3Å. The inhibitability of the medicines towards carbohydrate-hydrolase proteins (3W37 and 3AJ7) can be briefed as: acarbose (D1) > voglibose (D3) > miglitol (D2). In referencing, there are no signicant differentials of the inhibitory indicators between these commercial medicines and the selective iridoid compounds. This implies promising applicability of these substances as carbohydrate-hydrolase inhibitors in particular and in diabetes treatment in general. 1, 2, 10 regarding 3W37 inhibition and 2, 10, 14 regarding 3AJ7 inhibition are even highly competent.

Screening of physicochemical and pharmaceutical compatibility
Several properties of studied ligands are summarised in Table 4 in order to screen their physicochemical and pharmaceutical compatibility based on Lipinski's rule of ve. The most promising candidates, aka. 1, 2, and 10, weigh ca. 500 amu and possess negative dispersion coefficients (log P and log S). Only 10 creates less than both ten hydrogen bond acceptors and ve hydrogen bond donors regardless of the targeted protein, either 3W37 or 3AJ7. This preliminarily suggests 10 as the most potential compound for oral administration in attempt to tackle diabetes. However, other promising candidates should be in great consideration. First, 507.1 amu 14 might still be suitable since the commercial product Acarbose weigh 645.6 amu, over Lipinski's rst criterion. The highly efficacy-demonstrated medicine also forms seven hydrogen bond donor with protein  3AJ7. This justies the applicability of 1 and 2. In addition, it is noticeable that their polarisability is over 10Å 3 , indicating high polarisation. The property is of signicance because it is highly conducive to protein inhibition as the polypeptide molecule is made of polarised amino acids. Therefore, the selective iridoid compounds, especially 10, are likely compatible with pharmaceutical applications in physiological medium. However, in vitro and in vivo experiments are required in order to conrm the computational predictions. In summary, the results obtained altogether are likely to justify the applicability of the novel-synthesised iridoid compounds in cytological in general, and for inhibition of carbohydrate-hydrolase proteins (3W37 and 3AJ7) in particular.

a-Glucosidase inhibition
The enzyme inhibitory activity of fourteen compounds (excluding 15) were assayed towards a-glucosidase. Acarbose was used as a reference whose IC 50 is determined registering 285.00 mM in the resembling assay conditions. The obtained results are shown in Table 5, which are relatively consistent with those predicted by computational simulation. In detail, 10 > 2 > 14 > 13 > 1 assemble the most effective inhibitors towards aglucosidase whose attained IC 50 values register lower than 25 mM. Compound 10 shows its marked inhibition towards the enzyme with an IC 50 value of 0.05 mM, higher than that of acarbose by 5700 times. This was followed by 2 and 14, both exhibiting signicant a-glucosidase inhibition with IC 50 value 7.24 and 8.65 mM, respectively. The corresponding gure for 13 is 20.84 mM and for 1 is 24.89 mM. The predominance of 10 inhibitability seems signifying the role of (2E,6E)-8-hydroxy-2,6dimethyl-2,6-octadienoate moiety as it is the only compound containing this functional group. Its exibility and polarity acquired from an octadienoate conguration might be advantageous to structural transformation, which in turn conducive to either protein-site entrance, polar-polar interaction, or in-pose topographical tting. However, in-depth investigations are still required in order to reach more accurate view. Although predicted exhibiting insignicant inhibitability towards aglucosidase molecule (3W37) by molecular docking simulation, experimental results reect otherwise. On the other hand, 9 also sees noticeable simulation-experiment inconsistency, in which the compound is expected to perform pronounced a-glucosidase inhibition.
Compounds 3, 4, and 5, without double bond at C-7 00 and C-8 00 , are inactive towards a-glucosidase observed from enzyme assays. This is relatively consistent with the computational prediction as these compounds are not expected to be effective 3W37 inhibitors by the simulation. Therefore, double bond at C-7 00 and C-8 00 could be subjected for further investigations on aglucosidase inhibition. Regarding acetylation, the corresponding derivatives are either inactive including 8 and 11, or unfavoured such as 7 (IC 50 167.76 mM) and 12 (IC 50 71.56 mM) as a-glucosidase inhibitors. The corresponding computing also predicts this unfavourableness. The lessening of polarity, caused by total acetyl-hydroxyl substitution, might reason this, thus justifying the important role of catalpol hydroxyl groups to the inhibition. In addition, computationally expected possessing no inhibitability towards 3W37, 6 performs assaying inhibition towards a-glucosidase with a notably high IC 50 value of 224.78 mM. This reinforces the implication that two bulky triphenyl methyl groups at C-10 and C-6 0 seem to limit the hostmolecule entrance into the protein sites.
From the obtained results, several brief views on structural activity relationship of iridoid catalpol derivatives on a-glucosidase inhibition can be speculated. Acetylation of the hydroxyl groups at C-10 and glucose unit leads to reduced activity. Reduction of the double bond in the cinnamoyl moiety as well as locking the hydroxyl groups at C-6 0 and C-4 0 of the glucose unit by a cyclic acetal group totally deteriorate the activity. The presence of a hydroxyl group at C-4 00 of cinnamoyl moiety is likely to play an important role of activity, much better than the corresponding property for methoxy or acetyl group. The attachment of two bulky triphenylmethyl groups at C-10 and C-6 0 strongly reduces the activity, 32-fold decrease of inhibition. Finally, the location of a long, linear unsaturated alkyl ester chain at C-6 (octadienonate) instead of cinnamoyl moiety signicantly rises the a-glucosidase inhibitory activity. Nemoroside (10) is the most active compound with an IC 50 value of 0.05 mM, 5.700 times better than acarbose. This could be explained by the exibility of its bonding structure, thus conducive to its topological tness while inside the inhibitory holes.

Conclusions
This study suggests a promising approach to the use of iridoid chemotype from Dolichandrone spathaceae as a-glucosidase (protein 3W37) and oligo-1,6-glucosidase (protein 3AJ7) inhibitors. Five isolated natural iridoids and ten synthetic derivatives were successfully identied by structural elucidation. These include ve new iridoid-derived structures not found having reported before in the literature, i.e. 3, 4, 5, 6, and 9. Molecular