QSAR Study of Thieno [2, 3-d] Pyrimidine as a Promising Scaffold Using HQSAR, CoMFA and CoMSIA

1. INTRODUCTION

Podocarpusnagi (P. nagi, named Zhubai in Chinese) is widely distributed in south districts of Yangtze River, such as Fujian, Hunan, Guangxi, Guangdong, etc. This plant contains different kinds of biological compounds (such as volatile oil, flavonoids, steroids, sugar and glycosides, lactones and so on) and exhibits a wide spectrum of biological activities like hemostasis, bone setting, anti-bacterial, anti-tumor, antiviral, antioxidant and detume-scence activities^[1]. According to the folk records of the Yao Nationality, P. nagi has ever been used to treat trauma, stop-bleeding, fractures, knife wounds, gunshot wounds, body odor, eye diseases, colds, and so forth. The fresh bark or root of P. nagi was also used to treat the rheumatoid arthritis^[2-4]. Some work about the chemical components and biological activities of P.nagi has been reported: Ye Yang and XuYaming's groups isolated Podocarpus nagilactones from P. nagiplanted in Guangdong province and evaluated their biological activity. The results showed that most of them exhibited higher antitumor activity^{[5, 6]}. Chen Yegao's group isolated several bioflavo- noids and few steroids from the leaves of P. nagi grown in Yunnan^[7]. However, P. nagi was also distributed in Nan- ping of Fujian province. In recent years, a large scale of P. nagi was planted in Yangli town of Fujian province. Our research group has extracted the essential oil from its fruits, confirmed its chemical components, and evaluated its biological activities. The results showed this oil contains many active components: such as abundant unsaturated fatty acids, flavonoids, β-vanillin, vitamin E and essential microelements for human body; the biological evaluation results showed that it exhibited higher anti-oxidant^[8]. We checked the published papers and some old records and found that there are some differences about the P. nagi grown in Guangdong, Yunnan and Fujian. To the best of our knowledge, there are no reports about the chemical components of P. nagi planted in Fujian, so our research group took the lead to isolate the chemical components from the leaves of P. nagi planted in Fujian.

In this work, three sterols, (24R)-3β, 5α-dihydroxy-24- ethyl-5α-cholestan-6-one (1), 26, 27-dinorcholest-5-en-3-β- ol (2), and β-sitosterol (3), were isolated using the silica gel column chromatography. The preparative thin layer chromatography (PTLC) together with the recrystallization from the leaves of P.nagi and their structures was confirmed by NMR and XRD methods. Compound 1 was isolated from nature source for the first time and its crystal structure has not been reported. The crystal structure of 1 clearly explained its absolute configuration, and provided the reference for the assign of this kind of compounds. The crystal structures of 2 and 3 have been reported^[9-17]. The chemical structures of 1~3 are listed in Fig. 1. Compounds 1 and 2 are rare compounds with the same skeleton as the β-sitosterol, and were isolated for the first time from the leaves of P.nagi grown universally in Fujian. Compound 2 showed good to moderate in vitro anticancer activity against gastric cancer (NCI-N87), breast cancer MCF-7 (HTB-22), lung cancer A549 (CCL-185) and Hela (CCL-2)cell lines using the cell counting kit-8 (CCK-8) method^[18].

Figure 1

Figure 1.  Chemical structures of compounds 1~3

DownLoad: Full-Size Img PowerPoint

2. EXPERIMENTAL

2.1 Materials and instruments

Plant materials The leaves of Podocarpusnagi were collected in September of 2018 from the Yangli town of Fujian province, China and identified by one of the authors (J.P. Yong).

Instruments NMR spectra were recorded on a Bruker AV-400 spectrometer. Column chromatography (CC) was carried out onsilica gel (100~200 mesh, Qingdao Marine Chemical Inc., Qingdao, China). Melting points were determined on a XT-4 apparatus equipped with a microscope and uncorrected. Crystallography data were obtained from Rigaku SuperNova, with CCD detector and X-ray source of CuKα radiation (λ = 1.54184 Å). The structure was solved by direct methods with Olex2 Crystallographic Software.

2.2 Extraction and isolation

The detailed isolation processes are listed below: 10 kilograms of the air-dried and powdered leaves were added into a 25 L container and the material was dipped in 20 L70 % ethanol-water solution for one month and then filtered. The solution was concentrated under the reduced pressure, and the residue was dispersed in 5 L water and extracted with 1 L ethyl acetate for three times. The ethyl acetate layers were combined and concentrated under the reduced pressure to obtain another residue, which was rechromatographed over a column of silica gel with petroleum ether, petroleum ether-ethyl acetate (V_{petroleumether}: V_{ethyl acetate}, 10:1 to 0:1) as eluents to obtain some fractions: 10 fractions using petroleum ether as eluent; 10 fractions using V_{petroleumether}: V_{ethyl acetate}, 10:1 as eluent; 18 fractions using V_{petroleumether}: V_{ethyl acetate}, 5:1 as eluent; 20 fractions using V_{petroleumether}: V_{ethyl acetate}, 2:1 as eluent; 48 fractions using V_{petroleumether}: V_{ethyl acetate}, 1:1 as eluent; and 21 fractions using ethyl acetate as eluent. After the simple TLC analysis, we selected some fractions and combined to obtain another 6 fractions for further isolation: fraction 1 (petroleum ether as eluent); fraction 2 (V_{petroleum ether}: V_{ethyl acetate}, 10:1 as eluent); fraction 3 (V_{petroleumether}: V_{ethyl acetate}, 5:1 as eluent); fraction 4 (V_{petroleumether}: V_{ethyl acetate}, 2:1 as eluent); fraction 5 (V_{petroleumether}: V_{ethyl acetate}, 1:1 as eluent) and fraction 6 (ethyl acetate as eluent).

Compounds 1, 2 and 3 were isolated from the fraction 3 using silica gel column separation, preparative thin layer chromatrgraphy (PTLC) together with recrystallization. We checked the crystals under microscope and found that the appearance of the crystals was very different. We selected different crystals and analyzed their structures by XRD method, obtaining three different structural compounds.

2.3 Structure characterization

The isolated compounds were characterized using NMR and XRD methods. NMR was recorded on a 400 MHz Bruker AVANCE III spectrometer in CDCl₃. The chemical shifts were expressed in ppm relative to tetramethylsilane (TMS) as the internal standard. XRD were recorded on a SuperNova, Dual, Cu at zero, Atlas diffractometer equipped with graphite-monochromated CuKα radiation (λ = 1.54184 Å).

Compound 1: white lamellar single crystal, m.p.: 253~256 ℃; HR-MS for C₂₉H₅₀O₃Na, [M+Na]⁺: Calcd. 469.3652, found: 469.3652. This compound was confirmed by XRD analysis. A white lamellar single crystal of compound 1 with dimensions of 0.18mm × 0.18mm × 0.06mm was used for X-ray diffraction analysis. A total of 10599 reflections were collected at 100.01(16) K in the range of 5.24≤2θ≤149.60º by using an ω-scan mode, of which 6989 were unique with R_int = 0.0419 and R_sigma = 0.0595 and 6989 were observed with I > 2σ(I). The final R = 0.0561 and wR = 0.1209. The structure was solved by direct methods with SHELXS-2014 and refined by full-matrix least-squares methods with SHELXL-2014 program package^[19]. All of the non-hydrogen atoms were located with successive difference Fourier synthesis. Hydrogen atoms were added in idealized positions. The non-hydrogen atoms were refined anisotropically. Selected bond lengths and bond angles from XRD data are listed in Table 1. The XRD data are ideal and physical data agree well with (24R)-3β, 5α-dihydroxy-24- ethyl-5α-cholestan-6-one^[9]. The HR-MS result was also consistent well with its molecular weight.

Table 1

Table 1.  Selected Bond Lengths (Å) and Bond Angles (°) of Compound 1

DownLoad: CSV

Bond	Dist.		Angle	(°)
C(3)–O(1)	1.433(5)		O(1)–C(3)–C(2)	110.1(3)
C(5)–O(2)	1.442(4)		O(1)–C(3)–C(4)	108.8(3)
C(6)=O(3)	1.217(5)		O(2)–C(5)–C(4)	107.5(3)
C(5)–C(6)	1.540(4)		O(2)–C(5)–C(6)	105.2(3)
			O(2)–C(5)–C(10)	109.6(3)
			O(3)=C(6)–C(5)	123.3(3)
			O(3)=C(6)–C(7)	122.1(4)

Compound 2 is a white thin lamellar single crystal, m.p.: 127~128 ℃, and compound 3 is a white needle crystal, m.p.: 140~141 ℃. The NMR data, XRD analysis and relevant biological evaluation have been reported earlier^[9-18].

2.4 Preliminary in vitro anticancer evaluation

Compound 2 was selected to evaluate for their preliminary in vitro anticancer activity against gastric cancer (NCI-N87), breast cancer MCF-7 (HTB-22), lung cancer A549 (CCL-185) and Hela (CCL-2) cell lines using the CCK-8 method. Briefly, the cancer cell lines were seeded in 96-well plates (5000 cells/well) with 100 μL DMEM supplemented with 10% fetal bovine serum, and cultured at 37 ℃ in a humidified CO₂ incubator (95% air, 5% CO₂) for 24 h. While the cell lines grew to 90% in logarithmic growth, the culture medium was removed from each well, and 100 μL fresh DEME was added to each well. Then, 10 μL solution of compound 2 was added into each well (The experiment was repeated for 5 times) and the plates were incubated for another 48 h at 37 ℃. Subsequently, 10 μL CCK8 was added to each well, and the plates were cultured at 37 ℃ for another 4 hours. The optical density was measured at a wave-length of 450 nm on an ELISA microplate reader. DMEM and DMSO solution (V/V: 10/1) was used as a negative control. The results were expressed as the inhibition calculated at the ratio [(1-(OD₄₅₀ treated/OD₄₅₀ negative control)) × 100].

3. RESULTS AND DISCUSSION

During the isolation of fraction 3, after the silica gel column isolation, the preparative thin layer chromatography (PTLC) was used to isolate one compound. However, it was confirmed to be the mixture of compounds 1, 2 and 3 based on ¹H-NMR analysis. Because the R_f values of these compounds (1, 2 and 3) were almost the same, it is so difficult to differentiate them using TLC analysis and also very difficult to isolate them through PTLC. However, this mixture is very easy to form crystals in solution (V_{petroleum ether}: V_{ethyl acetate}, 1:1). We checked the crystals under microscope and found that their appearances are very different. We selected different crystals and analyzed their structures by XRD method, obtaining three different structural compounds: compound 1 (9 mg) as a white lamellar single crystal, compound 2 as a white thin lamellar single crystal (11 mg), and compound 3 as a colorless acicular single crystal (186 mg). It exhibited that compound 1 is in the very lower content in fraction 3 and the amount is not enough for NMR analysis, so we only finished the XRD and HR-MS analyses.

The molecular structure and ORTEP diagram of compound 1 are shown in Fig. 2. The skeleton of steroid is the same as that of compounds 2 and 3, while the hydroxyl at C(5) and carbonyl at C(6) of compound 1 are unique in comparison with compounds 2 (two hydroxyl groups at C(5) and C(6)) and 3 (double bonds between C(5) and C(6)). Besides, there are many differences of the dihedral angles of C(4)–C(10)–C(5)–C(6) between rings A and B of compounds 1 and 3. C(4)–C(5)–C(6) and C(10)–C(5)–C(6) are different planes of rings A and B, respectively. The value for 1 is 126.694º, while that for 3 is 177.734º, with the deviation to be 51.04º. The difference of dihedral angles might be caused by the type of bonds between C(5) and C(6): single bond (1.540 Å) in 1 but double bond (1.337 Å) in 3. The big groups at C(5) and C(6) of compound 1 increase the spatial effect and make rings A and B turn round accordingly. The schematic diagrams of the crystal cells and intermolecular hydrogen bonds of compound 1 are shown in Fig. 3 and Table 2. It is a supramolecular laminated structure, in which the branched alkanes of ten carbons can easily rotate and interact with the adjacent crystal cell layer. Intermolecular hydrogen bonds are formed by hydroxyl groups at C(5) between compound molecules (d2). Hydroxyl groups at C(3), C(5) and carbonyl groups at C(6) can also form hydrogen bonds with water (d1, d3, d5, d6). Effective hydrogen bonds can be formed in three-dimensional space, and strong intermolecular forces are found by calculating effective distances, the range of hydrogen bonds: 1.929~2.077 Å. The melting point (253~256 ℃) of compound 1 is much higher than that of compounds 3 (140~141 ℃) aroused by the intermolecular hydrogen bond.

Figure 2

Figure 2.  Molecular structure and ORTEP diagram of compound 1

DownLoad: Full-Size Img PowerPoint

Figure 3

Figure 3.  Crystal cell and hydrogen bonds of compound 1

DownLoad: Full-Size Img PowerPoint

Table 2

Table 2.  Hydrogen Bonds for Compound 1

DownLoad: CSV

D–H⋅⋅⋅A	d(D–H)/Å	d(H⋅⋅⋅A)/Å	d(D⋅⋅⋅A)/Å	∠D–H⋅⋅⋅A/°
O(2)–H(2)⋅⋅⋅O(2)1	0.82	2.08	2.899(3)	177
O(2)–H(2)⋅⋅⋅O(H2O)	0.82	2.02	2.840(3)	173
O(1)–H(1)⋅⋅⋅O(H2O)	0.82	2.00	2.757(3)	154
O(2)–H(2)⋅⋅⋅O(2)	0.82	2.08	2.872(3)	163
O(H2O)–H(H2O)⋅⋅⋅O(H2O)2	0.95	1.93	2.869(3)	170
12 – x, –1/2 + y, 1 – z; 22 – x, 1/2 + y, 1 – z

In this work, three compounds were isolated and confirmed from the leaves of P. nagi. Compound 1 is very rare compounds with the same skeleton as β-sitosterol and was isolated for the first time from the leaves of P. nagi. Its absolute configuration was confirmed using XRD. During isolation of compounds 1, 2 and 3, we used PTLC together with recrystallization methods. In addition, we checked the appearances of the crystals under microscope to distinguish one from another, and selected out one by one from the mixture. According to the reported results^{[5, 7]}, only β-sitosterol was isolated from the leaves of P. nagi grown in Yunan and Guangdong provinces. Maybe there are compounds 1 and 2 in the leaves of P. nagi grown in Yunan and Guangdong provinces (This assumption needs to be determined by comparing the chemical components planted in different provinces through studying the finger print of this plant later), but the authors did not isolate them, because these two compounds are so difficult to discover and isolate. But we obtained compounds 1 and 2 through selecting the crystals under microscope. This work provides an effective and worthy separation method for some compounds with the smaller differences of polarity.

It was reported that compound 3 could inhibit the proliferation of cancer cells and induce apoptosis^{[16, 17]}, indicating that this series of compounds are promising for in vitro anticancer. In current work, we tested compound 2 for its preliminary in vitro anticancer against gastric cancer, breast cancer (MCF-7), lung cancer (A549) and Hela cell lines. The results showed that compound 2 exhibited good to moderate inhibition against the four cancer cell lines with the inhibition of 89.16% ± 1.17, 97.02% ± 0.53, 47.20% ± 2.58 and 36.89% ± 1.22, respectively at the concentration of 1.4 × 10^-2 M. We have improved the cell viability experiment in vitro and enriched the tested cancer cell lines. These results indicated that these series of sterols maybe have good anticancer activity. It means that we found the new anticancer agent in this plant medicine. Inspired by this wok, more compounds will be isolated and their anticancer activity will be evaluated for the development of anticancer drugs.

References

1. [1]

   Pédeboscq, S.; Gravier, D.; Casadebaig, F.; Hou, G.; Gissot, A.; Rey, C.; Ichas, F.; De, G. F.; Lartigue, L.; Pometan, J. P. Synthesis and evaluation of apoptosis induction of thienopyrimidine compounds on KRAS and BRAF mutated colorectal cancer cell lines. Bioorgan. Med. Chem. 2012, 22, 6724−6731.
   

2. [2]

   El-Sayed, W. A.; Ali, O. M.; Zyada, R. A. F.; Mohamed, A. A.; Abdel-Rahman, A. A. H. Synthesis and antimicrobial activity of new substituted thienopyrimidines, their tetrazolyl and sugar derivatives. Acta Pol. Pharm. 2012, 3, 439−447.
   

3. [3]

   El-Sherbeny, M. A.; El-Ashmawy, M. B.; El-Subbagh, H. I.; El-Emam, A. A.; Badria, F. A. Synthesis, antimicrobial and antiviral evaluation of certain thienopyrimidine derivatives. Eur. J. Med. Chem. 1995, 5, 445−449.
   

4. [4]

   Kotaiah, Y.; Harikrishna, N.; Nagaraju, K.; Rao, C. V. Synthesis and antioxidant activity of 1, 3, 4-oxadiazole tagged thieno[2, 3-d]pyrimidine derivatives. Eur. J. Med. Chem. 2012, 12, 340−345.
   

5. [5]

   Shirole, N. L.; Shirole, K. D.; Deore, R. D.; Fursule, R. A.; Talele, G. S. Synthesis, characterization and pharmacological evaluation of 2-substituted thieno[2, 3-d]pyrimidine-4(3H)-ones. Asian J. Chem. 2007, 7, 4985−4992.
   

6. [6]

   Perrissin, M.; Favre, M.; Cuong, L. D.; Huguet, F.; Gaultier, C.; Narcisse, G. ChemInform abstract: synthesis and pharmacological activities of some substituted thienopyrimidin-4-ones. J. Cheminformatics 1989, 13, 104−105.
   

7. [7]

   Bousquet, E.; Guerrera, F.; Siracusa, M. A.; Caruso, A.; Amico-Roxas, M. Synthesis and pharmacological activity of 3-substituted pyrido[3΄, 2΄: 4, 5] thieno[3, 2-d] pyrimidin-4(3H)-ones. Farmaco 1984, 2, 110−119.
   

8. [8]

   Ashour, H. M.; Shaaban, O. G.; Rizk, O. H.; El-Ashmawy, I. M. Synthesis and biological evaluation of thieno [2΄, 3΄: 4, 5]pyrimido[1, 2-b][1, 2, 4]triazines and thieno[2, 3-d][1, 2, 4]triazolo[1, 5-a]pyrimidines as anti-inflammatory and analgesic agents. Eur. J. Med. Chem. 2013, 62, 341−351.  doi: 10.1016/j.ejmech.2012.12.003
   

9. [9]

   Tong, J.; Shan, L.; Qin, S. S.; Wang, Y. QSAR studies of TIBO derivatives as HIV-1 reverse transcriptase inhibitors using HQSAR, CoMFA and CoMSIA. J. Mol. Struct. 2018, 1168, 56−64.  doi: 10.1016/j.molstruc.2018.05.005
   

10. [10]

    Kubinyi, H. QSAR and 3D-QSAR in drug design part 1: methodology. Drug Discov. Today 1997, 11, 457−467.
    

11. [11]

    Abdizadeh, T.; Ghodsi, R.; Hadizadeh, F. 3D-QSAR (CoMFA, CoMSIA) and molecular docking studies on histone deacetylase 1 selective inhibitors. Recent Pat. Anti-Canc. 2017, 4, 365−383.
    

12. [12]

    Verma, J.; Khedkar, V. M.; Coutinho, E. C. 3D-QSAR in drug design. Curr. Top. Med. Chem. 2010, 1, 95−115.
    

13. [13]

    Kellogg, G. E.; Semus, S. F.; Abraham, D. J. HINT: a new method of empirical hydrophobic field calculation for CoMFA. J. Comput. Aid. Mol. Des. 1991, 6, 545−552.
    

14. [14]

    Bunce, J. D.; Patterson, D. E.; Frank, I. E. Crossvalidation, bootstrapping, and partial least squares compared with multiple regression in conventional QSAR studies. Quant. Struct-Act Rel. 1988, 1, 18−25.
    

15. [15]

    Klebe, G.; Abraham, U.; Mietzner, T. Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. J. Med. Chem. 2002, 24, 4130−4146.
    

16. [16]

    Salum, L. B.; Andricopulo, A. D. Fragment-based QSAR strategies in drug design. Expert. Opin. Drug. Dis. 2010, 5, 405−412.  doi: 10.1517/17460441003782277
    

17. [17]

    Andricopulo, A. D.; Salum, L. B. Fragment-based QSAR:  perspectives in drug design. Mol. Divers. 2009, 3, 277−285.
    

18. [18]

    Dunn III, W. J.; Wold, S.; Edlund, U.; Hellberg, S.; Gasteiger, J. Multivariate structure-activity relationships between data from a battery of biological tests and an ensemble of structure descriptors: the PLS method. Quant. Struct-Act Rel. 1984, 4, 131−137.
    

19. [19]

    Ali, E. M. H.; Abdel-Maksoud, M. S.; Oh, C. H. Thieno 2, 3-d pyrimidine as a promising scaffold in medicinal chemistry: recent advances. Bioorgan. Med. Chem. 2019, 7, 1159−1194.
    

20. [20]

    Sheila, A.; Malcolm, A. C.; Homer, R. W.; Tad, H.; Gregory, B. S. SYBYL Line notation (SLN): a versatile language for chemical structure representation. J. Chem. Inf. Model. 1997, 1, 71−79.
    

21. [21]

    Tong, J. B.; Zhan, P.; Wang, X. S.; Wu, Y. J. Quionolone carboxylic acid derivatives as HIV-1 integrase inhibitors: docking-based HQSAR and topomer CoMFA analyses. J. Chemometr. 2017, 12, 2934−2947.
    

22. [22]

    Damre, M. V.; Gangwal, R. P.; Dhoke, G. V.; Lalit, M.; Sharma, D.; Khandelwal, K.; Sangamwar, A. T. 3D-QSAR and molecular docking studies of amino-pyrimidine derivatives as PknB inhibitors. J. Taiwan. Inst. Chem. E 2014, 2, 354−364.
    

23. [23]

    Wang, Z. Y.; Chang, Y. Q.; Han, Y. S.; Liu, K. J.; Hou, J. S.; Dai, C. L.; Zhai, Y. H.; Guo, J. L.; Sun, P. H.; Lin, J.; Chen, W. M. 3D-QSAR and docking studies on 1-hydroxypyridin-2-one compounds as mutant isocitrate dehydrogenase 1 inhibitors. J. Mol. Struct. 2016, 6, 335−343.
    

24. [24]

    Hong, H.; Fang, H.; Xie, Q.; Perkins, R.; Sheehan, D. M.; Tong, W. Comparative molecular field analysis (CoMFA) model using a large diverse set of natural, synthetic and environmental chemicals for binding to the androgen receptor. Sar Qsar Environ. Res. 2003, 6, 373−388.
    

25. [25]

    Ghasemi, J. B.; Shiri, F. Molecular docking and 3D-QSAR studies of falcipain inhibitors using CoMFA, CoMSIA, and open 3D-QSAR. Med. Chem. Res. 2012, 10, 2788−2806.
    

26. [26]

    Bhonsle, J. B.; Venugopal, D.; Huddler, D. P.; Magill, A. J.; Hicks, R. P. Application of 3D-QSAR for identification of descriptors defining bioactivity of antimicrobial peptides. J. Med. Chem. 2008, 26, 6545−6553.
    

27. [27]

    Doddareddy, M. R.; Lee, Y. J.; Cho, Y. S.; Choi, K. I.; Koh, H. Y.; Pae, A. N. Hologram quantitative structure activity relationship studies on 5-HT6 antagonists. Bioorgan. Med. Chem. 2004, 14, 3815−3824.  doi: 10.1016/j.bmcl.2004.04.103
    

28. [28]

    Waller, C. L. A comparative QSAR study using CoMFA, HQSAR, and FRED/SKEYS paradigms for estrogen receptor binding affinities of structurally diverse compounds. J. Chem. Inf. Model. 2004, 2, 758−765.
    

29. [29]

    Patel, S.; Patel, B.; Bhatt, H. 3D-QSAR studies on 5-hydroxy-6-oxo-1, 6-dihydropyrimidine-4-carboxamide derivatives as HIV-1 integrase inhibitors. J. Taiwan. Inst. Chem. E 2016, 59, 61−68.
    

30. [30]

    ong, J. B.; Qin, S. S.; Jiang, G. Y. 3D-QSAR study of melittin and amoebapore analogues by CoMFA and CoMSIA methods. Chin. J. Struct. Chem. 2019, 2, 201−210.
    

31. [31]

    Wang, Y. J.; Zhang, G. W.; Yan, J. K.; Gong, D. M. Inhibitory effect of morin on tyrosinase: insights from spectroscopic and molecular docking studies. Food Chem. 2014, 163, 226−233.
    

32. [32]

    Jain, A. N. Effects of protein conformation in docking: improved pose prediction through protein pocket adaptation. J. Comput. Aid. Mol. Des. 2009, 6, 355−374.
    

33. [33]

    Tong, J. B.; Wang, Y.; Lei, S.; Qin, S. S. Comprehensive 3D-QSAR and binding mode of DAPY inhibitors using R-group search and molecular docking. Chin. J. Struct. Chem. 2019, 1, 25−36.