Skip to main content

New solid phase methodology for the synthesis of biscoumarin derivatives: experimental and in silico approaches


The simple and greener one-pot approach for the synthesis of biscoumarin derivatives using catalytic amounts of nano-MoO3 catalyst under mortar-pestle grinding was described. The use of non-toxic and mild catalyst, cost-effectiveness, ordinary grinding, and good to the excellent yield of the final product makes this procedure a more attractive pathway for the synthesis of biologically remarkable pharmacophores. Accordingly, biscoumarin derivatives were successfully extended in the developed protocols. Next, a computational investigation was performed to identify the potential biological targets of this set of compounds. In this case, first, a similarity search on different virtual libraries was performed to find an ideal biological target for these derivatives. Results showed that the synthesized derivatives can be α-glucosidase inhibitors. In another step, molecular docking studies were carried out against human lysosomal acid-alpha-glucosidase (PDB ID: 5NN8) to determine the detailed binding modes and critical interactions with the proposed target. In silico assessments showed the gold score value in the range of 17.56 to 29.49. Additionally, molecular dynamic simulations and the MM-GBSA method of the most active derivative against α-glucosidase were conducted to study the behavior of selected compounds in the biological system. Ligand 1 stabilized after around 30 ns and participated in various interactions with Trp481, Asp518, Asp616, His674, Phe649, and Leu677 residues.

Peer Review reports


Coumarin compounds exhibited brilliant and remarkable pharmaceutical activities [1]. Geiparvarin (1) [2], daphsafnin (2) [3] as well as daphjamilin (3) [4] known as natural coumarin compounds with anti-cancer and anti-monoamine oxidase activities (Fig. 1).

Fig. 1
figure 1

Natural coumarin compounds

Fig. 2
figure 2

coumarin derivatives drugs

Biscoumarins derivatives comprise a diverse and interesting group of heterocyclic drugs which are extremely important for their biological activities. Some coumarin derivatives, in general, and biscoumarins, in particular, are well known for their biological activity [5]. Some approved drugs with various properties are shown in Fig. 2. Dicoumarol (4) and warfarin (5) are naturally anticoagulant that depletes stores of vitamin K. Another anticoagulant drug is acenocoumarol (6) which inhibited the reduction of vitamin K by vitamin K reductase [6]. Aminocoumarin (7) as a coumarin derivative was an antibiotic for the treatment of staphylococci infections.

Several methods have been reported for the synthesis of biscoumarins, recently Kurt et al. synthesized a series of novel bis-coumarin derivatives containing triazole moiety [7]. Bavandi et al. demonstrated the synthesis of bis-4-hydroxy coumarins derivatives in the presence of Porcine pancreas lipase as a green catalyst [8]. Zeynizadeh et al. reported the novel magnetic graphene oxide/Zn–Ni–Fe layered double hydroxide nanocomposite for the green preparation of biscoumarins [9]. Hagiwara et al. showed the use of Et2AlCl as a Lewis acid for the condensation of 4-hydroxycoumarin and aldehydes in acetonitrile or dichloromethane at room temperature [10]. In continue, other researchers reported a similar reaction using piperidine, molecular iodine [11], tetrabutylammonium bromide (TBAB) [12], heteropolyacids [13], phospho- tungstic acid [14] and sodium dodecyl sulfate (SDS)[15] as catalysts. Some of these procedures required refluxing for hours in organic solvents, use of expensive catalysts, and tedious work-up. With the increasing public concern over environmental degradation, water is commonly considered as a benign solvent because of its non-toxicity, lower cost, and abundant natural occurrence [16].

Recently, mechanochemical synthesis has received lots of attention among synthetic chemists as a brilliant standard technology [17]. Many chemical reactions such as Grignard reactions [18] reductions [19] click reactions [20] and Knoevenagel condensations [21] have been carried out using the green and mechanochemical techniques. The importance of this type of protocol is free from organic solvents.

Metal oxide nanoparticles possess huge surface areas as a useful heterogeneous catalyst traditionally for catalyzing organic reactions. They are defined as friendly, and environmentally materials with attractive physical and chemical properties [22]. Metal oxide nanoparticles are effective catalysts that can replace nonselective, unstable, or toxic catalysts [23].

MoO3-nanoparticle and its derivatives are extensively attractive due to its unique structure, it was used as an inexpensive, eco-friendly, and highly reactive catalyst. Moreover, this non-toxic catalyst can be easily applied for diverse organic synthesis, generating corresponding products in excellent yields [24]. MoO3 showed significant Lewis acid property for diverse organic transformations in the liquid phase [25]. It demonstrated quite active over a wide range of temperatures and resistant to thermal excursions [26].

Similarity search, a subcategory of ligand-based virtual screening, has emerged as a reliable, fast, and inexpensive method that finds compounds with high similarity in some ways, especially molecular features, to known bioactive molecules. The assumption that molecules that are globally similar in structure could exhibit similar biological activity is generally valid [27, 28]. Once a target has been identified, molecular docking and molecular dynamic simulation, known as structure-based virtual screening, can be applied [29,30,31]. In docking assessments databases of available molecules are docked into the region of interest of protein in silico and scored based on predicted interactions with the site. In these cases, compounds with a higher affinity toward the target can be discovered [32]. Molecular dynamic (MD) simulation is an effective technique to measure the behavior of fluids at the molecular level. In the MD simulation, the movement of system particles is determined in a certain period and the system evolution is investigated. MD simulation is a process similar to real experiments. Thus, it is an ideal tool to study the different interactions and conformation of ligands in the biological system at the nanometer scale [33, 34].

As a result, with inspiration from the remarkable biological activities of biscoumarin derivatives [35,36,37] hereby, the solid-phase pathway for Knoevenagel–Michael reactions of biscoumarins (3) by the reaction of aldehydes (2) and 4- hydroxycoumarin (1) in the presence of 10 mol% nano-MoO3 were reported. The products were achieved by grinding the reactants in a mortar with a pestle under solvent-free conditions. Similarity search analysis and molecular docking studies were performed to find an ideal biological target for this set of compounds. Finally, MD simulations were combined with the generalized-Born surface area method (MM/GBSA) to evaluate the solvation free energy of the protein and the ligand in the dynamic situations.

Experimental section


Solvents were purified by standard procedures. Reactions were followed by TLC using. Melting points were determined with an Electrothermal 9100 apparatus in open capillary tubes and are uncorrected. The chemicals used in this work were purchased from Fluka or Merck companies and were used without further purification. FT-IR spectra were recorded on Thermo Nicolet Nexus 670 and 1H NMR were determined by Bruker Avance 300 MHz spectrometers.

General procedure for the synthesis of biscoumarin derivatives using mortar–pestle grinding

An aromatic aldehyde (1 mmol), 4-hydroxycoumarin (2 mmol) and MoO3 (10 mol%) were subjected to mortar and pestle grinding for specified time. After completion of the reaction (monitored by TLC), 3 mL of distilled water was added to the reaction mixture and the product was extracted and washed with water (10 ml). Recrystallization from hot ethanol afforded pure biscumarin derivatives in good to excellent yield. The structure of all products was confirmed by appropriate spectroscopic and physical methods (Melting point or IR, 1HNMR) with those reported or with authentic samples prepared by the conventional method.

The procedure for the synthesis of 3, 3'-((4-hydroxyphenyl) methylene) bis (4-hydroxy-2H-chromen-2-one) using reflux:

4-hydroxybenzaldehyde (1 mmol) and 4-hydroxycoumarin (2 mmol) were taken in 25 mL round bottom flask followed by the addition of a different solvent (15 mL) and MoO3 (10 mol%) was stirred and reflux for the appropriate times. The progress of the reaction was monitored by TLC, the reaction mixture was washed with H2O (10 mL) and EtOAc (10 mL) to almost afford the pure product. The crude products was recrystallized from hot ethanol to obtain the pure compound.

3,3'-((4-chlorophenyl)methylene)bis(4-hydroxy-2H-chromen-2-one) (3h): FT-IR (KBr) (KBr, ν cm−1) 3072, 2923,2856, 1667, 1610, 1562, 1349, 1095, 764; 1H NMR (300 MHz, CDCl3) δ (ppm): 11.5 (s, 1H, OH), 11.3 (s, 1H, OH), 8.0 (s, 2H, arom), 7.6 (s, 2H, arom),7.1–7.4 (m, 8H, arom), 6.04 ( s, 1H, CH).

3,3'-((4-methoxyphenyl)methylene)bis(4-hydroxy-2H-chromen-2-one) (3i): FT-IR (KBr, cm−1) 3065, 2989, 2842, 2727, 1665, 1613, 1563, 1349, 1042, 769; 1H NMR (300 MHz, CDCl3) δ (ppm): 11.52 (s, 1H, OH), 11.31 (s, 1H, OH), 8.04 (s, 2H, arom), 7.65–7.63 (m, 2H, arom),7.42–7.39 (m, 4H, arom),7.27 (s,2H, arom) 7.14- 7.12 (m, 2H, arom), 6.04 ( s, 1H, CH), 3.89 (s, 3H, CH3).

3,3'-((4-nitrophenyl)methylene)bis(4-hydroxy-2H-chromen-2-one) (3j): FT-IR (KBr, ν cm−1) 3064, 2922, 2857, 2720, 1654, 1609, 1564, 1339, 1096, 764 1H NMR (300 MHz, CDCl3) δ (ppm): 11.56 (s, 2H, OH), 8.15–8.02 (m, 4H, arom), 7.65–7.26 (m, 8H, arom), 6.11 ( s, 1H, CH).

3,3'-((3-nitrophenyl)methylene)bis(4-hydroxy-2H-chromen-2-one) (3k): FT-IR (KBr, ν cm−1) 3071, 2925, 2859, 2772, 1656, 1609, 1526, 1346, 1010, 763; 1H NMR (300 MHz, CDCl3) δ (ppm): 11.58 (s, 1H, OH), 11.38 (s, 1H, OH), 8.14–7.99 (m, 4H, arom), 7.67–7.43 (m, 8H, arom),6.14 (s, 1H, CH).

3,3'-((2,4-dichlorophenyl)methylene)bis(4-hydroxy-2H-chromen-2-one)(3l): FT-IR (KBr, ν cm−1) 3071, 2925, 2861, 1653, 1311, 1096, 756; 1H NMR (300 MHz, CDCl3) δ (ppm): 11.60 (s, 1H, OH), 11.59 (s, 1H, OH), 7.98–7.34 (m, 11H, arom), 6.07–6.01 (s, 1H, CH).

3,3'-(phenylmethylene)bis(4-hydroxy-2H-chromen-2-one) (3m): FT-IR (KBr, ν cm−1) 3065,3000, 2927, 2737, 1659, 1609, 1336, 1090, 755; 1H NMR (300 MHz, CDCl3) δ (ppm): 11.43 (s, 2H, OH), 8.04 (s, 2H, arom), 7.39–7.31 (m, 9H, arom), 6.10 (s, 1H, CH).

4-(bis(4-hydroxy-2-oxo-2H-chromen-3-yl)methyl)benzaldehyde (3n): FT-IR (KBr, ν cm−1) 3068, 2923, 2806, 2730, 1701, 1607, 1097, 764, 1H NMR (300 MHz, CDCl3) δ (ppm): 11.54 (s, 1H, OH), 11.34 (s, 1H, OH), 9.93 (s, 1H, arom), 7.97–7.38 (m, 12H, arom), 6.10 (s, 1H, CH).

Insights into the biological activities of biscoumarin derivatives

In order to propose biological activities of the synthesized compounds, three steps computational process was applied, using ligand-based similarity search, molecular docking, and molecular dynamics.

Similarity-based analog searching

To find an ideal biological target for this set of compounds, ligand-based similarity search on several libraries including SwissTargetPrediction (, PubChem similarity search (, SEA Search Server (, MolTarPred (, SuperPred was performed. In these databases, the structure of compounds was uploaded and different searching approaches including fingerprint and shape-based similarity pharmacophores were applied automatically to find the most similar bioactive agents compared to the synthesized compounds.

Procedure for docking studies

In this study, the Gold docking program was used to carry out the docking calculations between the biscoumarin analogs and the binding site of human lysosomal acid-alpha-glucosidase (PDB ID: 5NN8, [38]. The protein structure was prepared using the Discovery Studio Client so that ligands and waters were removed from 5NN8 and all hydrogens were added. The box for docking calculations was built taking into account the center of the co-crystallized ligand for the enzyme with a 10 Å radius sphere around the co-crystallized ligand.

First, GOLD docking program with different functions was used for docking analyses via re-dock acarbose inside the 5NN8. All other options were set as default. The best accuracy with the lowest RMSD was seen in the ChemScore fitness function. The 3D structures of ligands were first generated by the Hyperchem and then energy minimization and optimization were performed to generate the initial confirmation via molecular mechanics (Amber) followed by molecular dynamic (AM1) approaches. These approaches are widely applied in ligand-docking and MD simulations which offer acceptable criteria to determine the proper and favorable molecular arrangement and potential energy of a molecule [39,40,41]. The derivatives were docked into the active site of protein using default parameters 10 runs for each ligand; Genetic algorithm with 100% efficacy, min ops 10,000; max ops 125,000 were chosen, All other options were set as default. The top ChemScore value was used for further analysis. The higher value confirms better interaction with the active site. Finally, protein–ligand interactions were analyzed with Discovery Studio Visualizer [42,43,44].

Molecular dynamics simulations

Starting model was obtained by imposing the best ChemScore to acid-α-glucosidase (PDB ID: 5NN8). MD simulations were conducted using the desmond operator of Schrodingers suit maestro [45,46,47]. To build the system for MD simulation, the protein–ligand complexes were solvated with SPC explicit 22,530 water molecules and placed in the center of an orthorhombic box in the periodic boundary condition. The system’s charge was neutralized by adding 84 atoms of Na and 63 atoms of Cl to simulate the real cellular ionic concentrations, respectively. The MD simulations protocol involved minimization, pre-production, and finally production MD simulation steps. In the minimization procedure, the entire system was allowed to relax for 2500 steps by the steepest descent approach. Then the temperature of the system was raised from 0 to 300 K with a small force constant on the enzyme to restrict any drastic changes. MD simulations were performed via NPT (constant number of atoms; constant pressure, i.e., 1.01325 bar; and constant temperature, i.e., 300 K) ensemble. The Nose–Hoover chain method was used as the default thermostat with 1.0 ps interval and Martyna-Tobias-Klein as the default barostat with 2.0 ps interval by applying an isotropic coupling style. Long-range electrostatic forces were calculated based on the particle-mesh-based Ewald approach with the cutoff radius for Columbia forces set to 9.0 Å. Finally, the system was subjected to produce MD simulations for 100 ns for each protein–ligand complex. During the simulation, every 1000 ps of the actual frame was stored. The dynamic behavior and structural changes of the systems were analyzed by the calculation of the RMSD and RMSF. Subsequently, the energy-minimized structure calculated from the equilibrated trajectory system was evaluated for the investigation of each ligand–protein complex interaction.

Results and discussion

Chemical synthesis

MoO3-nanoparticles were previously synthesized [26, 48]. For identification of nanostructure powder X-ray diffraction (PXRD) and scanning electron microscopy (SEM) analysis were carried out. Figure 3, displayed powder X-ray diffraction (PXRD), and the surface morphology and the diameter of the nanostructure were studied by SEM in Fig. 4. Debye–Scherrer Eq. (1) was used for calculating crystal size structure, according to this equation: D is the mean size of crystalline, k is a constant (= 0.9 assuming that the particles are spherical), λ is the X-ray wavelength, β is the line width (obtained after correction for the instrumental broadening) and θ is the angle of diffraction (Bragg angle). The average particle size obtained from XRD data is approximately ~ 50 nm [26]

$${\text{D }} = {\text{ k}}\lambda /\beta {\text{cos}}\theta$$
Fig. 3
figure 3

PXRD patterns of the synthesized MoO3 nano-crystals

Fig. 4
figure 4

Scanning electron micrographs for prepared nano-MoO3

The synthetic route to access compounds 3a- 3n are achieved by the solid phase general pathway illustrated in Fig. 5. The condensation reaction between 4-hydroxycoumarin and various aldehydes is summarized. This solvent-free reaction needed only a few minutes of reaction time. This kind of reaction is expected to be the most economical method since neither solvent is used.

Fig. 5
figure 5

General pathway for synthesis of biscoumarin derivatives

Effect of different catalyst

The amount of catalyst has shown an excellent effect on the rate and yield of the reaction, To optimize the reaction conditions, the amount of various catalysts were applied in the solid phase synthesis of 3,3'-((4-hydroxphenyl) methylene)bis(4-hydroxy-2H-chromen-2-one) 3a from the condensation of 4-hydroxycoumarin with 4-hydroxybenzaldehyde and the best result at room temperature was obtained (Table 1).

Table 1 Optimization of the solvent-free reaction

According to the result in Table 1, when the reaction was performed in the absence of catalyst, longer reaction time was required (60 min) and very low yield of product was achieved (≤ 5%) even if the reaction time was prolonged (Table 1, entry 1). To obtain the satisfactory results of (3a), the reaction using various homogeneous and heterogeneous Bronsted and/or Lewis acids was performed (Table 1, entries 2–7). Accordingly, the nano-MoO3 was the best catalyst and applied for all reactions. Using 10 mol % of nano-MoO3, 3a was isolated in 93% yield after the progress of the reaction for a few minutes (Table 1, entry7). By increasing the amount of the catalyst in the model reaction, no change was observed in the time and yield of the reaction (Table 1, entries 8–9). The decline in the use of the catalyst amount less than 10% resulted in low yields (Table 1, entries 10–11).

Under the optimized reaction conditions, a series of biscoumarins derivatives (3a–n) were synthesized, the structure of all products was confirmed by appropriate spectroscopic and physical methods (Table 2).

Table 2 Synthesis of biscoumarins derivatives (4a–n) using nano-MoO3 under mortar and pestle grinding

Effect of different solvents and temperature in reflux condition

Our next approach was to study various solvents in different temperatures on the model reaction of 4-hydroxycoumarin and 4-hydroxybenzaldehyde in the presence of (10 mol %) of MoO3 for preparation of 3a (Table 3).

Table 3 Optimization of the solvent and temperature in room temperature or reflux condition

The model reaction was performed in the different solvents and various temperatures (Table 3 entries 1–12), the results confirmed that carrying out the reaction in EtOH/ H2O (1:1) in reflux condition gave the highest yield (Table 3, entries 10). Eventually, the reaction using mortar–pestle grinding method is low cost, more efficient, and simple reaction without the usage of organic solvent that gives the desired product.

Effect of different aldehydes

After reaction optimization, different aldehydes were chosen and reacted under optimum conditions. This multi-component and solid phase approach can be used for both aromatic aldehydes with electron-withdrawing and electron-donating groups. Furthermore, a wide range of aromatic aldehydes was successfully used in this reaction with excellent results.

A proposed mechanistic route for the condensation of aldehydes and 4-hydroxycoumarin that rationalizes the formation of the products is exhibited in (Fig. 6). First, nucleophilic attack of 4-hydroxycoumarin to the activated aldehyde (by MoO3 coordination), followed by H2O elimination provides intermediate “A” that was further activated by MoO3. This, in turn, undergoes a second nucleophilic attack by another 4-hydroxycoumarin to provide the final product.

Fig. 6
figure 6

A plausible mechanism for synthesis of biscoumarins derivatives using MoO3

Similarity search and docking studies

The similarity metrics analysis on several libraries indicated that the synthetic biscoumarin derivatives can be ideal α-glucosidase inhibitors. Regarding the similarity of reported α-glucosidase inhibitors with the designed structure, molecular docking evaluations were performed to study the binding mode of all derivatives with the α-glucosidase active site. α-glucosidase enzyme catalyzes the hydrolysis of starch to simple sugars which leads to an increase in blood glucose levels. α-glucosidase inhibitors can be ideal and effective anti-diabetic agents [51,52,53]. Docking studies of the mentioned compounds were carried out with gold docking software using different score fitness functions including chem score, gold score, ChemPLP, and the best accuracy with the lowest RMSD was seen in the ChemScore fitness function. Validation of the molecular docking method was done by redocking the crystallographic ligand of the target enzyme, gainst 5NN8, which testified the validation of the docking calculations. Alignment of the best pose of acarbose in the active site of α-glucosidase and crystallographic ligand recorded an RMSD value of 1.45 Å via ChemScore fitness functions. Results of ChemScore fitness values were reported in Table 4. It should be noted that fitness scores in gold software are dimensionless and the higher scale value showed better interactions with the active site [54, 55].

Table 4 Docking scores and interactions of compounds against α-glucosidase

Images were created by Discovery Studio Visualizer v20.1.0.19295.

As can be seen in Fig. 7, compound 1 demonstrated the best docking score (29.49) compared to other derivatives. The chromen-2-one moiety of compound 1 made two pi-alkyl interactions with Leu650 and Leu678 as well as two pi-pi T-shaped interactions with Phe649. On the other side of the molecule, 4-hydroxy chromen-2-one showed a hydrogen bond with Arg600 plus a pi–alkyl interaction with Phe525. Another H-bounding interaction was also seen between the oxygen of 4-hydroxyphenyl ring with Asp404. 4-hydroxyphenyl also recorded additional pi-pi stacked interaction with Phe649.

Fig. 7
figure 7figure 7

Docking conformations of compounds 1–14 (orange stick) in the α-glucosidase active site. Hydrogen bonds are depicted in green dashed lines, π-π stacked interactions are depicted in dark pink dashed lines, π-aryl interactions are depicted in pink dashed lines pi-sulfur interactions are depicted in pale orange dashed lines and and pi-anion interactions are depicted in dark orange dashed lines

MD simulations

Molecular docking studies alone can be misleading, as they are performed in static conditions. To validate these results and better understand the potential for interaction of ligand 1 and their ability to coordinate with α-glucosidase active site MD simulations were performed. To study the steadiness of the protein–ligand complex, the root mean square deviation (RMSD) of the complexed backbone was investigated in MD simulation. As can be seen in Fig. 8, the simulation period was adequate to reach a balanced ligand-complex structure over the simulation time and the RMSD values stabilize around a fixed value of 1.60 Å. Changes of the order of 1–3 Å are perfectly acceptable for small, globular proteins. Changes much larger indicate the protein is undergoing a large conformational change during the simulation.

Fig. 8
figure 8

RMSD plot of the enzyme in complexed compound 1 in the MD simulation time. RMSD values of the Ca atoms of the protein are depicted in blue, and ligand-complex values are exhibited in red

The root mean square fluctuation (RMSF) is useful for characterizing local changes along the protein chain and the flexibility of the protein (Fig. 9).

Fig. 9
figure 9

RMSF plot of the α-glucosidase residue in complexed with compound 1

Type and percent of protein interactions with the ligand throughout the simulation type are exhibited in Fig. 10. As can be seen interaction with Asp616 can be seen in almost 100% frame followed by Leu677, Asp518, Arg411, Trp376, and Phe649.

Fig. 10
figure 10

Protein–ligand contacts during the whole simulation time in α-glucosidase complexed with compound 1

Interactions that occur more than 30.0% of the simulation time in the selected trajectory (0.00 ns through 100.02 ns) are shown in Fig. 11. Asp616 participated in impotent hydrogen bond interaction with 4-hydroxy phenyl in more than 90% of cases. Also, Asp518 demonstrated interaction (90%) with OH of 2H-chromen-2-one through water bridge. Derivative 1 also exhibited interactions with Trp481, Phe649, Leu677, and His674. Another water bridge was also recorded between Arg411 and C=O of 2H-chromen-2-one (Fig. 11).

Fig. 11
figure 11

2D representation of ligand-residue interactions

Overall, it can be seen that synthesized biscoumarin derivatives due to containing special functionality (OH and C=O and O) can participate in several interactions with the residue of the active site. Also, the planner structure of the aromatic ring in biscoumarin provides several pi-interactions with the active site which stabilizes the binding site of the enzyme. Also, the binding energy of compound 1 calculated using the MM/GBSA method is presented in Table 5.

Table 5 Binding free energy estimates via MM-GBSA after 100 ns MD simulations. All values are in Kcal/mol


In summary, an easy and efficient protocol for the synthesis of biscoumarin derivatives in the presence of MoO3 nanoparticles was explained. This methodology showed considerable synthetic advantages in terms of product diversity, simplicity of the reaction procedure, mild reaction conditions, and good to excellent yields. This procedure could be classified as green chemistry due to the elimination of any hazardous organic solvent. Similarity search analysis proposed the biscoumarin pharmacophore as ideal α-glucosidase inhibitors and molecular docking studies exhibited that these derivatives effectively fitted within the α-glucosidase active site with Chemscore value in the range of 17.56 to 29.49. MD assessments of compound 1 recorded critical interactions with Asp616, Leu677, Asp518, Arg411, Trp376, and Phe649 which prove the potential of this chemical family as promising drug candidates for further experimental analysis.

These findings demonstrated the applicability of similarity search followed by molecular docking and MD assessments as a promising method in the early stage of drug discovery.

Availability of data and materials

The datasets analyzed during the current study are available in the Worldwide Protein Data Bank (wwPDB) repository (


  1. Mahmoodi NO, Ghanbari Pirbasti F, Jalalifard Z. Recent advances in the synthesis of biscoumarin derivatives. J Chin Chem Soc. 2018;65(4):383–94.

    CAS  Article  Google Scholar 

  2. Carotti A, Carrieri A, Chimichi S, Boccalini M, Cosimelli B, Gnerre C, et al. Natural and synthetic geiparvarins are strong and selective MAO-B inhibitors. Synthesis and SAR studies. Bioorg Med Chem Lett. 2002;12(24):3551–5.

    CAS  PubMed  Article  Google Scholar 

  3. Riaz M, Malik A. Daphsaifnin, a dimeric coumairn glucoside from Daphne oleoides. Heterocycles. 2001;55(4):769–73.

    CAS  Article  Google Scholar 

  4. Riaz M, Malik A. Structure determination of daphjamilin, a new bicoumarin glycoside, by NMR spectroscopy. Magn Reson Chem. 2001;39(10):641–2.

    CAS  Article  Google Scholar 

  5. Borges F, Roleira F, Milhazes N, Santana L, Uriarte E. Simple coumarins and analogues in medicinal chemistry: occurrence, synthesis and biological activity. Curr Med Chem. 2005;12(8):887–916.

    CAS  PubMed  Article  Google Scholar 

  6. Schalekamp T, Brassé BP, Roijers JF, Chahid Y, van Geest-Daalderop JH, de Vries-Goldschmeding H, et al. VKORC1 and CYP2C9 genotypes and acenocoumarol anticoagulation status: interaction between both genotypes affects overanticoagulation. Clin Pharmacol Ther. 2006;80(1):13–22.

    CAS  PubMed  Article  Google Scholar 

  7. Kurt BZ, Dag A, Doğan B, Durdagi S, Angeli A, Nocentini A, et al. Synthesis, biological activity and multiscale molecular modeling studies of bis-coumarins as selective carbonic anhydrase IX and XII inhibitors with effective cytotoxicity against hepatocellular carcinoma. Bioorg Chem. 2019;87:838–50.

    CAS  PubMed  Article  Google Scholar 

  8. Bavandi H, Habibi Z, Yousefi M. Porcine pancreas lipase as a green catalyst for synthesis of bis-4-hydroxy coumarins. Bioorg Chem. 2020;103:104139.

    CAS  PubMed  Article  Google Scholar 

  9. Zeynizadeh B, Gilanizadeh M. Synthesis and characterization of a magnetic graphene oxide/Zn–Ni–Fe layered double hydroxide nanocomposite: an efficient mesoporous catalyst for the green preparation of biscoumarins. New J Chem. 2019;43(47):18794–804.

    CAS  Article  Google Scholar 

  10. Hagiwara H, Miya S, Suzuki T, Ando M, Yamamoto I, Kato M. Et2AlCl-promoted coupling reactions of 4-Hydroxy-2-pyrone or 4-hydroxycoumarin with aldehydes: synthesis of methylenebis (4-hydroxy-2-pyrone) or methylenebis (4-hydroxycoumarin) derivatives. ChemInform. 1999.

    Article  Google Scholar 

  11. Kidwai M, Bansal V, Mothsra P, Saxena S, Somvanshi RK, Dey S, et al. Molecular iodine: a versatile catalyst for the synthesis of bis (4-hydroxycoumarin) methanes in water. J Mol Catal A: Chem. 2007;268(1–2):76–81.

    CAS  Article  Google Scholar 

  12. Khurana JM, Kumar S. Tetrabutylammonium bromide (TBAB): a neutral and efficient catalyst for the synthesis of biscoumarin and 3, 4-dihydropyrano [c] chromene derivatives in water and solvent-free conditions. Tetrahedron Lett. 2009;50(28):4125–7.

    CAS  Article  Google Scholar 

  13. Heravi MM, Sadjadi S, Haj NM, Oskooie HA, Bamoharram FF. Role of various heteropolyacids in the reaction of 4-hydroxycoumarin, aldehydes and ethylcyanoacetate. Catal Commun. 2009;10(13):1643–6.

    CAS  Article  Google Scholar 

  14. Singh P, Kumar P, Katyal A, Kalra R, Dass SK, Prakash S, et al. Phosphotungstic acid: an efficient catalyst for the aqueous phase synthesis of bis-(4-hydroxycoumarin-3-yl) methanes. Catal Lett. 2010;134(3–4):303–8.

    CAS  Article  Google Scholar 

  15. Mehrabi H, Abusaidi H. Synthesis of biscoumarin and 3, 4-dihydropyrano [c] chromene derivatives catalysed by sodium dodecyl sulfate (SDS) in neat water. J Iran Chem Soc. 2010;7(4):890–4.

    CAS  Article  Google Scholar 

  16. Sheldon RA. Selective catalytic synthesis of fine chemicals: opportunities and trends. J Mol Catal A Chem. 1996;107(1–3):75–83.

    CAS  Article  Google Scholar 

  17. Amrute AP, De Bellis J, Felderhoff M, Schüth F. Mechanochemical synthesis of catalytic materials. Chemistry (Weinheim an der Bergstrasse, Germany). 2021;27(23):6819.

    CAS  Google Scholar 

  18. Toda F, Takumi H, Yamaguchi H. Grignard reactions in the solid state. ChemInform. 1989.

    Article  Google Scholar 

  19. Toda F, Kiyoshige K, Yagi M. NaBH4 reduction of ketones in the solid state. Angew Chem, Int Ed Engl. 1989;28(3):320–1.

    Article  Google Scholar 

  20. Thorwirth R, Stolle A, Ondruschka B, Wild A, Schubert US. Fast, ligand-and solvent-free copper-catalyzed click reactions in a ball mill. Chem Commun. 2011;47(15):4370–2.

    CAS  Article  Google Scholar 

  21. Ren Z, Cao W, Tong W. The Knoevenagel condensation reaction of aromatic aldehydes with malononitrile by grinding in the absence of solvents and catalysts. Synth Commun. 2002;32(22):3475–9.

    CAS  Article  Google Scholar 

  22. Chavali MS, Nikolova MP. Metal oxide nanoparticles and their applications in nanotechnology. SN Appl Sci. 2019;1(6):607.

    CAS  Article  Google Scholar 

  23. Moghaddam FM, Saeidian H. Controlled microwave-assisted synthesis of ZnO nanopowder and its catalytic activity for O-acylation of alcohol and phenol. Mater Sci Eng B. 2007;139(2–3):265–9.

    CAS  Article  Google Scholar 

  24. Manivel A, Lee G-J, Chen C-Y, Chen J-H, Ma S-H, Horng T-L, et al. Synthesis of MoO3 nanoparticles for azo dye degradation by catalytic ozonation. Mater Res Bull. 2015;62:184–91.

    CAS  Article  Google Scholar 

  25. Reddy BM, Reddy VR. Influence of SO42−, Cr2O3, MoO3, and WO3 on the stability of ZrO2-tetragonal phase. J Mater Sci Lett. 2000;19(9):763–5.

    CAS  Article  Google Scholar 

  26. Zarenezhad E, Rad MNS, Mosslemin MH, Tabatabaee M, Behrouz S. Nano-MoO3 as a highly efficient heterogeneous catalyst for a one-pot synthesis of tetrahydropyrimidine derivatives in water. J Chem Res. 2014;38(10):607–10.

    CAS  Article  Google Scholar 

  27. Stahura FL, Bajorath J. New methodologies for ligand-based virtual screening. Curr Pharm Des. 2005;11(9):1189–202.

    CAS  PubMed  Article  Google Scholar 

  28. Willett PJC. Fusing similarity rankings in ligand-based virtual screening. Comput Struct Biotechnol J. 2013;5(6): e201302002.

    PubMed  PubMed Central  Article  Google Scholar 

  29. Lyne PD. Structure-based virtual screening: an overview. Drug Discov Today. 2002;7(20):1047–55.

    CAS  PubMed  Article  Google Scholar 

  30. Cavasotto CN, W Orry AJ. Ligand docking and structure-based virtual screening in drug discovery. Curr Top Med Chem. 2007;7(10):1006–14.

    CAS  PubMed  Article  Google Scholar 

  31. Varela-Rial A, Majewski M, De Fabritiis G. Structure based virtual screening: Fast and slow. Wiley Interdiscip Rev Comput Mol. 2021;12:e1544.

    Google Scholar 

  32. Sohrabi M, Binaeizadeh MR, Iraji A, Larijani B, Saeedi M, Mahdavi M. A review on α-glucosidase inhibitory activity of first row transition metal complexes: a futuristic strategy for treatment of type 2 diabetes. RSC Adv. 2022;12(19):12011–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. Ciccotti G, Ferrario M, Schuette CJE. Molecular dynamics simulation. 2014;16(233):1

  34. Nasli Esfahani A, Iraji A, Alamir A, Moradi S, Asgari MS, Hosseini S, et al. Design and synthesis of phenoxymethybenzoimidazole incorporating different aryl thiazole-triazole acetamide derivatives as α-glycosidase inhibitors. Mol Divers. 2021.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Al-Warhi T, Sabt A, Elkaeed EB, Eldehna WM. Recent advancements of coumarin-based anticancer agents: an up-to-date review. Bioorg Chem. 2020;103: 104163.

    CAS  PubMed  Article  Google Scholar 

  36. Srikrishna D, Godugu C, Dubey PK. A review on pharmacological properties of coumarins. Mini Rev Med Chem. 2018;18(2):113–41.

    CAS  PubMed  Article  Google Scholar 

  37. Upadhyay HC. Coumarin-1, 2, 3-triazole hybrid molecules: an emerging scaffold for combating drug resistance. Curr Top Med Chem. 2021;21(8):737–52.

    CAS  PubMed  Article  Google Scholar 

  38. Roig-Zamboni V, Cobucci-Ponzano B, Iacono R, Ferrara MC, Germany S, Bourne Y, et al. Structure of human lysosomal acid α-glucosidase–a guide for the treatment of Pompe disease. Nat Commun. 2017;8(1):1111.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  39. Roy K, Kar S, Das RN. Chapter 5—computational chemistry. In: Roy K, Kar S, Das RN, editors. Understanding the basics of QSAR for applications in pharmaceutical sciences and risk assessment. Boston: Academic Press; 2015. p. 151–89.

    Chapter  Google Scholar 

  40. Adcock SA, McCammon JA. Molecular dynamics: survey of methods for simulating the activity of proteins. Chem Rev. 2006;106(5):1589–615.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. Hongmao S. Chapter 2—structure-based ligand design i: with structures of protein/lead compound complex available. In: Hongmao S, editor. A practical guide to rational drug design. Woodhead Publishing; 2016. p. 15–60.

    Chapter  Google Scholar 

  42. Moghaddam FM, Daneshfar M, Daneshfar Z, Iraji A, Samandari-Najafabad A, Faramarzi MA, et al. Synthesis and characterization of 1-amidino-O-alkylureas metal complexes as α- glucosidase Inhibitors: structure-activity relationship, molecular docking, and kinetic studies. J Mol Struct. 2022;1250: 131726.

    CAS  Article  Google Scholar 

  43. Karami M, Hasaninejad A, Mahdavi H, Iraji A, Mojtabavi S, Faramarzi Mohammad A, et al. One-pot multi-component synthesis of novel chromeno[4,3-b]pyrrol-3-yl derivatives as alpha-glucosidase inhibitors. Mol Divers. 2021.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Shareghi-Boroujeni D, Iraji A, Mojtabavi S, Faramarzi MA, Akbarzadeh T, Saeedi M. Synthesis, in vitro evaluation, and molecular docking studies of novel hydrazineylideneindolinone linked to phenoxymethyl-1,2,3-triazole derivatives as potential α-glucosidase inhibitors. Bioorg Chem. 2021;111: 104869.

    CAS  PubMed  Article  Google Scholar 

  45. Maestro, Schrödinger, LLC, New York, NY; 2021.

  46. Bowers KJ, Chow DE, Xu H, Dror RO, Eastwood MP, Gregersen BA, et al., editors. Scalable algorithms for molecular dynamics simulations on commodity clusters. In: SC'06: proceedings of the 2006 ACM/IEEE conference on supercomputing. IEEE; 2006.

  47. Sastry GM, Adzhigirey M, Day T, Annabhimoju R, Sherman W. Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. J Comput Aided Mol Des. 2013;27(3):221–34.

    PubMed  Article  CAS  Google Scholar 

  48. Soltani Rad MN, Behrouz S, Mohammadtaghi-Nezhad J, Zarenezhad E, Agholi M. Silica-tethered cuprous acetophenone thiosemicarbazone (STCATSC) as a novel hybrid nano-catalyst for highly efficient synthesis of new 1, 2, 3-triazolyl-based metronidazole hybrid analogues having potent antigiardial activity. Appl Organomet Chem. 2019;33(4): e4799.

    Article  CAS  Google Scholar 

  49. Sadeghi B, Ziya T. A fast, highly efficient, and green protocol for synthesis of biscoumarins catalyzed by silica sulfuric acid nanoparticles as a reusable catalyst. J Chem. 2013;2013:1–5.

    Google Scholar 

  50. Tabatabaeian K, Heidari H, Khorshidi A, Mamaghani M, Mahmoodi NO. Synthesis of biscoumarin derivatives by the reaction of aldehydes and 4-hydroxycoumarin using ruthenium (III) chloride hydrate as a versatile homogeneous catalyst. J Serb Chem Soc. 2012;77(4):407–13.

    CAS  Article  Google Scholar 

  51. Sohrabi M, Nazari Montazer M, Farid SM, Tanideh N, Dianatpour M, Moazzam A, et al. Design and synthesis of novel nitrothiazolacetamide conjugated to different thioquinazolinone derivatives as anti-urease agents. Sci Rep. 2022;12(1):2003.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. Moghaddam FM, Daneshfar M, Daneshfar Z, Iraji A, Samandari-Najafabad A, Faramarzi MA, et al. Synthesis and characterization of 1-amidino-O-alkylureas metal complexes as α-glucosidase Inhibitors: Structure-activity relationship, molecular docking, and kinetic studies. J Mol Struct. 2022;1250: 131726.

    CAS  Article  Google Scholar 

  53. Iraji A, Shareghi-Brojeni D, Mojtabavi S, Faramarzi MA, Akbarzadeh T, Saeedi M. Cyanoacetohydrazide linked to 1,2,3-triazole derivatives: a new class of α-glucosidase inhibitors. Sci Rep. 2022;12(1):8647.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  54. Allen FH, Bellard S, Brice M, Cartwright BA, Doubleday A, Higgs H, et al. The Cambridge Crystallographic Data Centre: computer-based search, retrieval, analysis and display of information. Acta Crystallogr Sect B Struct Crystallogr Cryst Chem. 1979;35(10):2331–9.

    Article  Google Scholar 

  55. Groom CR, Allen FH. The Cambridge Structural Database in retrospect and prospect. Angew Chem Int Ed. 2014;53(3):662–71.

    CAS  Article  Google Scholar 

Download references


Not applicable.


The authors wish to thank the support of Noncommunicable Diseases Research Center, Fasa University of Medical Sciences (Grant No, 400003) and Also, this research was ethically approved, IR.FUMS.REC.1400.036.

Author information

Authors and Affiliations



EZ contributed to the designing of the work; interpretation and supervising this study, as well as the synthesis, purification, and characterization of compounds and provided the initial draft and revising of the article; MT carried out the synthesis, purification, and characterization of catalyst. MNM contributed to the acquisition, and analysis of molecular dynamics simulation data and provided the initial draft of the article. CI contributed to the acquisition, analysis, and interpretation of in silico data and provide the initial draft of the article. AI contributes to the designing of the work, the interpretation and supervising in silico part of the study as well as revising the manuscript. All authors read and approved the article.

Corresponding authors

Correspondence to Cambyz Irajie or Aida Iraji.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zarenezhad, E., Montazer, M.N., Tabatabaee, M. et al. New solid phase methodology for the synthesis of biscoumarin derivatives: experimental and in silico approaches. BMC Chemistry 16, 53 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Molecular dynamics simulations
  • Biscoumarin
  • MoO3-nanoparticle