Available online on 15.11.2024 at http://jddtonline.info

Journal of Drug Delivery and Therapeutics

Open Access to Pharmaceutical and Medical Research

Copyright  © 2024 The   Author(s): This is an open-access article distributed under the terms of the CC BY-NC 4.0 which permits unrestricted use, distribution, and reproduction in any medium for non-commercial use provided the original author and source are credited

Open Access  Full Text Article                                                                                                                                                    Research Article

Potential Inhibitor of DENV-2 Virus Protease (NS2B-NS3): An In-Silico Studies of Anti-Viral Plants

Mhashevolu Rhakho *

Department of Botany, St. Joseph’s College (A) Jakhama, Nagaland, India.

Article Info:

_______________________________________________

Article History:

Received 12 Aug 2024  

Reviewed 27 Sep 2024  

Accepted 23 Oct 2024  

Published 15 Nov 2024  

_______________________________________________

Cite this article as: 

Rhakho M, Potential Inhibitor of DENV-2 Virus Protease (NS2B-NS3): An In-Silico Studies of Anti-Viral Plants, Journal of Drug Delivery and Therapeutics. 2024; 14(11):10-23            DOI: http://dx.doi.org/10.22270/jddt.v14i11.6870        _______________________________________________

*Address for Correspondence:  

Mhashevolu, Assistant professor, Department of Botany, St. Joseph’s College (A) Jakhama, Nagaland, India. 

Abstract

_______________________________________________________________________________________________________________

Dengue virus (DENV) is a mosquito-borne pathogen that affects millions of people worldwide. The DENV-2 protease is a vital enzyme responsible for viral replication and is a promising target for antiviral therapy.  The objective of the study is to identify potential inhibitors of DENV-2 protease using In-Silico approaches with phytocompounds from ten antiviral plants. Initially, 133 phytoconstituents were collected with anti-dengue properties from previously reported studies which were virtually screened using SWISS ADME for ADME properties. The DENV-2 protease structure (2FOM) was obtained from the Protein Data Bank and molecular docking was performed using AutoDock Vina. The best-scoring compounds were evaluated and top five potential inhibitors with high binding affinity and stability were selected. The top-scoring compounds were Ligand-91 (Terchebin, -8.1 kcal/mol), Ligand-13 (7-desacetyl-7-benzoylgedunin, -7.8 kcal/mol), Ligand-100 (Triterpenoid, -7.8 kcal/mol), Ligand-12 (7-desacetyl-7-benzoylazadiradione, -7.7 kcal/mol), Ligand-20 (Azadirolic acid, -7.7 kcal/mol), Ref.1 (Doxycycline, -6.6 kcal/mol), Ref.2(Monosdenvir, -7.5 kcal/mol), and Ref.3 (Zanamivir, -5.6 kcal/mol). The result of the study shows that 7-desacetyl-7-benzoylazadiradione and 7-desacetyl-7-benzoylgeduninas compounds with high binding affinity for the target protein. These compounds are found in Azadirachta indica making it a promising candidate for further experimental validation and development of antiviral agents against DENV-2.

Keywords: Molecular docking, Anti-dengue, Anti-viral, ADME analysis

 


 
  1. INTRODUCTION

Medicinal plants have been used in healthcare in traditional medicine since time immemorial. It has been a cornerstone of traditional medicine for centuries, providing natural remedies for various ailments. These plants are a rich source of bioactive compounds that have evolved to protect them from environmental stresses, pathogens, and insects. Humans have tapped into this chemical arsenal to develop medicines that prevent and treat diseases. Studies have been carried out globally to focus on plant-based chemicals for drug discovery which have led to the production of plant-based medicines 1. Medicinal plants have also played a crucial role in the development of modern medicine. Many conventional drugs are derived from plants, such as aspirin from willow bark, quinine from cinchona trees, and morphine from opium poppies. Additionally, plants are a rich source of leads for drug discovery, with many pharmaceutical companies exploring plant-based compounds for new medicines 2. Herbal medicines have gained its popularity in the last few decades for its therapeutic treatments like Ayurveda, Unani, Homeopathy, Sidha, etc, which use herbs as a major constituent. Plants consist of multiple alkaloids/compounds that occur naturally (as against single extracts) exhibit synergistic actions such as antiviral, antibacterial, anti-protozoa, and antioxidant 3. Medicinal plants used in traditional medicine and indigenous knowledge to treat a variety of ailments include bulbs, shrubs, ferns, and trees4.Among the different plant parts, the leaves were most frequently used for the treatment of diseases followed by root, fruit, flower, tuber, rhizome, bark, stem5. Herbal medicines can be administered as a paste, decoction, infusion, juice or poultice, or taken orally with no preparation. Medicinal plants contain a wide variety of secondary metabolites or compounds such as tannins terpenoids, alkaloids, flavonoids and possess significant antibacterial, antifungal, anticancer, antidiuretic, anti-inflammatory and anti-diabetic. Pharmacognosy provides a vast understanding of potential medicinal plants applications for the prevention and treatment of numerous conditions and diseases including Dengue6. There is an immense ongoing interest among the researchers in developing new antiviral therapeutic agents from natural products or the extract of plants. The numbers of studies that have been undertaken to document and preserve medicinal plant knowledge in Nagaland are few. The numbers of scientific research into the efficacy and safety of phytomedicines and scientific investigations are carried out to confirm the anti- viral properties which can be used to combat various viral diseases. The following antiviral plants are considered for the study, viz., Azadirachta indicaBrassica campestrisCarica papayaCatharanthus roseusEuphorbia hirtaOcimum sanctumPsidium guajavaAllium sativumHouttuynia cordata, and Mimosa catechu7.

Dengue is a virus-based infection caused by the dengue virus (DENV), a member of the genus Flavivirus in the family Flaviviridae which is transmitted to humans through the bites of infected female mosquitoes, primarily the Aedes aegypti mosquito. The virus in other species within the Aedes genus such as Aedes albopictus, Aedes polynesiensis, Aedes scutellaris, Aedes (Stegomyia) aegypti, Aedes albopictus, etc can also act as vectors, but their contribution is secondary to Aedes aegypti8-9. Along with the dengue virus, Flavivirus includes Japanese Encephalitis Virus (JEV), Yellow fever virus (YFV), West Nile virus (WNV) and tick-borne Encephalitis virus (TBEV). Zika virus outbreak has been reported in various countries which is a prominent human-pathogenic flaviviruses 10-11. Dengue is endemic in more than 125 countries, a significant cause of mortality in tropical and subtropical regions. Infections range from asymptomatic self-limited, acute, febrile disease called dengue fever (DF) to life-threatening dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS). There has been a drastic increase case of Dengue and over 2.5 billion people, in more than 100 countries, are at risk of dengue infection, with millions of cases occurring around the world every year. 12-13. According to WHO, an estimation of 390 million infections is reported annually worldwide, of which 96 million are clinically apparent. The largest number of dengue cases reported was in 2023. The WHO Region of the Americas reported 4.5 million cases, with 2300 deaths. A high number of cases were reported in Asia: Bangladesh (321 000), Malaysia (111 400), Thailand (150 000), and Vietnam (369 000). Dengue is considered to be ubiquitous throughout the tropics, with local spatial variations in risk influenced strongly by rainfall, temperature and the degree of urbanization. Dengue fever can be better framed for health interventions in terms of specific environmental features and assemblages of high-risk spaces. 14-15. Half of the world’s population is at risk of dengue infection as there is no specific antiviral treatment that act against DENV is currently available despite considerable efforts to find potent inhibitors for the Dengue protease. The only dengue vaccine available on the market is CYD-TDV, which was developed by Sanofi Pasteur (marketed as Dengvaxia and is currently approved in 20 countries in Latin America, Asia, and Australia. While CYD-TDV is effective for the prevention of severe infection in previously infected people, it also increases the risk of severe dengue in individuals who have not been previously infected by dengue. The results highlighted were important in directing future vaccine development to identify immune correlates of protection 16-17

The causative agent for dengue fever (DF) is dengue virus (DENV), an RNA virus from the Flavivirus genus belonging to the Flaviviridae family. It is roughly spherical in shape, an enveloped single-stranded positive-sense RNA virus with an icosahedral nucleocapsid covered by the lipid bilayer. It has a positive-strand RNA genome inside a protein capsid also known as nucleocapsid and 11 kb long DENV genome can function as mRNA, and, similar to that in eukaryotes, there are untranslated regions (UTRs) at both the 5′ and 3′ end flanking the open reading frame (ORF). DENV has the ability to infect a wide range of cell types including cells of the human immune system ranging from dendritic cells, monocytes, B- and T-cells, hepatocytes, endothelial cells. Dengue infections are caused by 4 different serotypes, which are: DENV 1, DENV2, DENV3, and DENV4. These serotypes share approximately 65% of their genome similarity, but there is some genetic variation in each serotype. Despite these variations, infection with each of the dengue serotypes results in the same disease and range of clinical symptoms13, 18.

The pathogenic female Aedes aegypti mosquito releases the Dengue virus via saliva into the skin of the mammalian host. Salivary components of Aedes aegypti mosquitoes have been identified as significant contributors to increased viral replication. A salivary protein named “34 kDa protein” increases DENV viral titer in human keratinocytes, reducing the expression of the antimicrobial peptides LL-37 and S100A7 and type I interferons. Viral replication starts in the salivary glands of the vector which leads to virion release in the saliva that initiates binding of the virus to host cell receptors for entry via receptor-mediated endocytosis 8,19,20. DENV after receptor-mediated endocytosis, the positive-strand viral RNAis released into the cytoplasm and translated into a polyprotein. The translated polyprotein are further cleaved by viral and host proteases into three structural (capsid: C, membrane: M, and envelope: E) and seven nonstructural (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) proteins which helps in replication. The RNA genome and the capsid protein interact to form a complex, while other structural proteins form part of the virion envelope. Although the NS proteins are absent within the virion, they assist in virus replication and evasion of the immune system within an infected cell. 21,22,23,26. C protein localized in the cytoplasm and nuclei is a foremost structural component of DENV which is thought to be crucial for its well-organized replication. The lipid bilayer of virions is formed by lipid between the nucleocapsid core and E/M outer shell. 

During the replication of DENV, a membrane-bound replication complex formation helps to incorporate host factors, viral proteins, and genomic RNA. The positive-strand (+) DENV genomic RNA acts as a template to synthesize complementary negative-strand (-) RNA, which acts as a template for the synthesis of new positive-strand viral RNAs21,24,25. Viral replication occurs on the endoplasmic reticulum by NS proteins followed by viral assembly and trafficking of immature viral particles to the trans-Golgi network (TGN) where it undergoes pH-dependent maturation. The acidic environment of the TGN facilitates viral maturation where the pr domain of precursor membrane protein (prM) is cleaved by the host protease furin. prM shields the envelope proteins from premature fusion and pH-induced reorganization during viral secretion. Subsequently, the fully matured infectious virus particle are then released from the hostcell.8,25,26.

Among the four serotypes of Dengue -DENV 1, DENV2, DENV3, and DENV4, DENV2 is most prevalent globally which contains a single-stranded RNA of positive polarity. The Dengue virus (DENV) is roughly spherical in shape with a diameter of approximately 50 nm, which enveloped single-stranded positive-sense RNA virus of 11 kb long which is translated into a large polyprotein during the infectious life cycle. This polyprotein is processed by cellular and viral proteases into three mature structural proteins: the capsid (C), envelope (E), and membrane (M) proteins along with other seven nonstructural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5). These nonstructural proteins play roles in viral replication, virion assembly and attenuation of the host antiviral response 8,27,28.


 

 

 image

Figure 1: Dengue protease (NS2B-NS3) with catalytic site


 

NS3, a large multifunctional protein of 618 amino acids belongs to the superfamily 2 of RNA helicases/NTPases having Walker A, GK(S/T) and Walker B, DEx(D/H) motifs, along with other superfamily-characteristic conserved sequence motifs present in subdomains I and II (which possess the fold of the RecA protein) endowed with protease, helicase, nucleoside 5′-triphosphatase (NTPase), as well as 5′-terminal RNA triphosphatase activities, plays an important role in viral polyprotein processing and genome replication. The N-terminal of NS3 comprises a serine protease domain, with the protein NS2B acting as a membrane-anchoring cofactor, necessary for proteolytic activity (figure 1). Helicase activity of NS3 catalyzes the unwinding of the duplex RNA in the 3′–5′ direction releasing the single stranded RNA available to NS5 as a template for replication in the presence of a divalent cation and the energy derived from the hydrolysis of ATP 29,30, 31. NS2B is a viral serine protease and, along with other NS proteins NS1, NS2A, NS4, and NS4B, plays diverse roles in viral replication, assembly, and release. The C-terminal end of the NS3 protein has three enzymatic properties: a 5′ RNA-triphosphatase (RTP), a nucleoside triphosphatase (NTPase), and a helicase. NS3 forms a complex with NS5 and assists in replication through the unwinding of viral RNA and dephosphorylation prior to 5′-end capping 8, 29, 32. The dengue virus NS3 protease, a member of the flavivirin enzyme family (EC 3.4.21.91), is located in the N-terminal 184 residues of the multifunctional 69 kDa NS3 protein and contains a functional catalytic triad consisting of His51, Asp75 and Ser135 (figure 1). In addition to the serine protease, the NS3 protein contains enzymatic activities of a nucleoside triphosphatase, a 5' - RNA triphosphatase (RTPase) and a RNA - stimulated RNA helicase. The NS3 protease catalyses the post-translational cleavage of the viral polyprotein precursor in the non-structural region at the NS2A/NS2B, NS2B/NS3, NS3/NS4A and NS4B/NS5 sites and at additional sites within the viral capsid protein, NS2A, NS4A and within a C-terminal region of NS3 itself. The overall conformation of the dengue virus NS3 protease displays the β-barrel conformation typical for serine proteases, although the viral enzyme appears to possess higher compactness with short or absent loop structures and a relatively shallow substrate binding site 33,34. Based on a number of studies, the methyltransferase (MTase) domain of the DV non-structural protein NS5 (NS5 MTase) is thought to be a promising antiviral target35, 36] However, the closed conformation is the major form even in the unbound state, which thus represents the best model for structure-guided drug designs27, 37, 38.

  1. METHODS

2.1  Selection of Phytocompounds

After comprehensive search on databases, 10 plants having antiviral properties were selected to confirm their efficiency against DENV-2. A total of 133 phytoconstituents from the selected plants were identified.


 

Table 1: List of selected phytocompounds.

Sl. No

Phytocompounds 

Name of Plants

1

Nimbin

 

Azadirachta indica39

 

2

desacetylnimbin

3

desacetylsalannin

4

kaempferol 3-O-rutinoside

5

Epicatechin

6

nimbinene,

7

6-deacetylnimbinene

8

Nimbandiol

9

Nimbolide

10

ascorbic acid

11

n-hexacosanol

12

7-desacetyl-7-benzoylazadiradione

13

7-desacetyl-7-benzoylgedunin

14

7-hydroxyazadiradione

15

Nimbiol

16

Quercetin

17

ß-sitosterol

18

Rutin

19

Limbonin

20

Azadirolic acid

21

Nimbochalcin

22

Azadirachtin

23

Limbocidin

24

3-p-coumaroylquinic

Brassica campestris40

25

Caffeic

26

Ferulic

27

sinapic acids

28

kaempferol 3-sophoroside-7-glucosides

29

Citric

30

Aconitic

31

2-ketoglutaric

32

Maleic

33

Shikimic

34

fumaric acids

35

Zeaxanthin

36

Lutein

37

p-coumaric (4-hydroxycinnamic

38

Chymopapain,

Carica papaya41.42

39

Carpaine

40

Dehydrocarpaine

41

Oleic acid

42

Tocopherol

43

Squalene

44

Neophytadiene

45

butyl 9,12,15-octadecatrienoate

46

n-hexadecanoic acid

47

phytol ( tetramethyl-2-hexadecen)

48

dasycarpidan-1-methanol acetate (ester)

49

octadecenoic acid

50

D-limonene

51

bis (2-(2-chloroethoxy)ethyl) ether

52

Dimethoxydimethylsilane

53

dibenzyl ether

54

benzhydrazide,

55

o-butylisourea

56

2-chloro-5,5-dimethyl-1-phenyl-3-hexen-1-ol

57

2-methoxybenzeneacetaldehyde

58

Myricetin

59

1-methyl-2-pyrrolidinone

60

Benzonitrile

61

Nonanal

62

octanoic acid

63

1-decene

64

nonanoic acid

65

Benzene

66

Benzene, 1,3-bis(1,1-dimethylethyl)-5-methoxy

67

1-iodooctadecane

68

2-methylnaphthalene

69

2-tetradecene

70

10-undecenoic acid

71

Dodecanal

72

1,4-dimethylnaphthalene

73

9-oxononanoic acid

74

1-hentriacontane

75

2,4-di-tert-butylphenol

76

nonanedioic acid

77

dimethyl ester

78

1-octadecene

79

Ajmalicine

Catharanthus roseus43,44

80

Catharanthine

81

Tabersonine

82

Serpentine

83

Vindoline

84

Quercitol

Euphorbia hirta45

85

Alpha (α)-amyrin

86

beta (β)-amyrin

87

Friedelin

88

Taraxerol

89

Euphorbins B

90

Euphorbin E

91

Gallic

92

Geraniin

93

ellagic acid 

94

Terchebin

95

Tannic

96

tartaric acids

97

lupeol fatty acid ester

98

fatty acid  phytyl     esters

99

linoleic acid

100

Triterpenoid

101

Eugenol

Ocimum sanctum46,47

102

β-sitosterol

103

Stigmasterol

104

Campesterol

105

Reynoutrin

Psidium guajava48,49

106

Guajaverin

107

Avicularin

108

Isoquercitrin

109

Hyperoside

110

2,6-dihydroxy-3,5-dimethyl- 4-O- (6 ″ - O-galloyl-β-D-glucopyranosyl) –benzophenone

111

flavonoids (naringin)

112

flavonoids ( hesperetin)

113

flavonoids (daidzein)

114

γ-glutamyl-l-cysteine peptides

Allium sativum50,51

115

Alliin (S-allyl-l-cysteine sulfoxide)

116

diallyl thiosulfinate

117

diallyl disulfide

118

Garlicin

119

diallyl trisulfide (allitridin or DATS), 

120

Ajoene

121

vinyl-dithiins (3-Vinyl-4H-1,2-dithiin)

122

Chlorogenic acid,

Houttuynia cordata52,53

123

Glabranin

124

7-O-methyl-glabranine

125

4-methyl-heptane

Mimosa catechu54

126

carboxylic acid

127

methyl laurate

128

2-ethyl-3-methylbut-1-ene

129

tetra decanoic 

130

4-hydroxybenzoic acid

131

Afzelechin

132

Aromadendrin

133

Baicalein

 


 
  1.  ADME Screening

The three-dimensional structure of the phytoconstituents identified from the selected plants was obtained from the PUBCHEM55 and the structures that are not available were drawn using ChemDRAW. The Mol2 structure of the ligands were used as the input file for screening of the drug-like properties was performed using SWISS ADME and the parameters of Lipinski Rule of 5 i.e, Molecular Weight (MW), Partition Coefficient (LogP), Hydrogen bond donors, Hydrogen Bond acceptors, Molar Refractivity (MR) were recorded.

2.3 Molecular Docking 

The structures of the compounds were further evaluated via Molecular docking using Pyrx56,57. Structures of Ligands were minimized using Pyrx and the parameters for Energy minimization are uff for force-field and Conjugate Gradient for Optimization Algorithm. The file is then converted to pdbqt file. The structure of protease, 2FOM58 was obtained from the Protein Data Bank. The structure of Protein was prepared using Discovery Studio where the water molecules and hetero atoms were removed. Pyrx software was employed for energy minimization and conversion to pdbqt file. Molecular Docking Analysis was performed using Autodock Vina59,60 which is incorporated in Pyrx Graphical User Interface.  The Active site selected on the protease for docking are HIS51, VAL75, LYS73, LYS74, ASP75, THR120, LEU128, PRO132, SER135, GLY151, ASN152, GLY153, VAL154, and ALA164. The active sites are centered inside the grid box where the docking is performed.

  1. RESULTS AND DISCUSSION

A total of 133 phytoconstituents from 10 different plants and three reference drugs were screened using SWISS ADME. The software used for ADME Screening is SWISSADME61 where a total number of 133 ligands from the selected plants in mol2 files were fed into SWISSADME. Based on Lipinski rule of five, all the ligands were selected for docking as there were no serious violation to the rule. 

3.1 Protein Preparation

The viral Protease used for this study is Dengue Virus NS2B/NS3 Protease 2FOM. Viral Protein 2FOM structure was downloaded from PUBCHEM 62 in PDB format and prepared using MOE.  The catalytic triad His51, Asp75 and Ser 135 are present

3.2 Molecular Docking

Among the selected phytoconstituents, the top 5 phytoconstituents docking scores of the interactions with 2FOM and interacting active site amino acids are ranked based on their docking scores and shown in table 2, where Ligand91 (Terchebin) being the highest with -8.1 kcal/mol and Ligand12 (7-desacetyl-7-benzoylazadiradione) and Ligand20 (Azadirolic acid) being the lowest with a score of -7.7 kcal/mol each. Docking results shows that all the selected Phytoconstituents have higher docking scores as compared to the reference drugs used in this study.


 

 

 

 

 

Table 2: Top docking scores with important interactions

Sl. No

Ligand no

Ligands

Dock Score (kcal/mol)

Important interactions

1

91

Terchebin

-8.1

Trp50, Arg54, Asp75, Gly153, Tyr161, Leu128, Tyr150, Phe130, Pro132, Ile36, Val52, His51

2

13

7-desacetyl-7-benzoylgedunin

-7.8

Leu128, Phe130, Gly151, Gly153, Asp75, His51

3

100

Triterpenoid

- 7.8

Tyr150, Gly151

4

12

7-desacetyl-7-benzoylazadiradione

-7.7

Pro132, His51, Tyr161, Tyr150, Leu128

5

20

Azadirolic acid

-7.7

Gly151, His51, Lys73, Val72, Lys74

6

Ref.1

Doxycycline

-6.6

Val72, Trp50, Asp75, Gly151, Gly153, Val154

7

Ref.2

Monosdenvir

-7.5

Leu128, Pro132, Ser131, His51, Asp75, Asn152, Val154

8

Ref.3

Zanamivir

-5.6

Gly151, Pro132, Ser135, His51, Try150, Phe130, Leu128, Gly153

 

 

 

(A)

(B)

 

 

(C)

(D)

 

 

(E)

(F)

 

 

(G)

(H)

 

 

(I)

(J)

 

 

(K)

(L)

 

 

(M)

(N)

 

 

(O)

(P)

Figure 2: 2D and 3D images of Protein-Ligand interactions: (A) 2D image of protein-ligand91 interactions (B) 3D image of protein-ligand91 interactions (C) 2D image of protein-ligand13 interactions (D) 3D image of protein-ligand13 interactions (E) 2D image of protein-ligand100 interactions (F) 3D image of protein-ligand100 interactions (G) 2D image of protein-ligand12 interactions (H) 3D image of protein-ligand12 interactions (I) 2D image of protein-ligand20 interactions (J) 3D image of protein-ligand20 interactions (K) 2D image of protein-ligandRef.1 interactions (L) 3D image of protein-ligandRef.1 interactions (M) 2D image of protein-ligandRef.2 interactions (N) 3D image of protein-ligandRef.2 interactions (O) 2D image of protein-ligandRef.3 interactions (N) 3D image of protein-ligandRef.3 interactions.


 

The results of molecular docking analysis of the top binding ligands with the reference are shown in figure2. Result of the docking analysis of Protein-ligand91 Interactions shows 7 Conventional Hydrogen Bond with Amino acids – HIS 51, ILE36, PRO 132, PHE 130, TYR 150, TYR 161, GLY 15, Hydrogen Bond interaction with catalytic Traid –HIS 51 and One unfavourable donor interaction with - TRP 50 is observed. In Protein-ligand13 interaction, docking analysis shows 3 Cabon Hydrogen bond with PHE 130, GLY 151, GLY 153, Pi- Cation bond with the catalytic triad – HIS 51 and Pi-Anion bond with the catalytic triad –ASP 75. Additionally, one unfavourable interaction is observed. Docking analysis of Protein-ligand100 interactions shows 1 Conventional Hydrogen Bond with TYR 50, 1 Cabon Hydrogen bond with GLY 151. There is no interaction with the catalytic triad and no unfavourable interaction is observed. In case of Protein-ligand12 interactions, it shows 1Conventional Hydrogen Bond with TYR 150, 1 Cabon Hydrogen bond with LEU 128, Pi-Sigma interaction with catalytic triad- HIS 51 and No unfavourable interaction observed. A protein-ligand20 interaction shows 2 Conventional Hydrogen Bond with GLY 151 and VAL 72, 1 Cabon Hydrogen bond with LYS 74. There is no interaction with the catalytic triad, however, one unfavourable Acceptor-Acceptor interaction with LYS 73 is observed. Docking analysis of Protein-ligand Ref.1 interactions shows 5 Conventional Hydrogen Bond with VAL 72, TRP 50, ASP 75, GLY 151, GLY 153, Interaction with catalytic triad – ASP 75 and one unfavourable interaction is observed. In Protein-ligandRef.2 interactions, there is1 Cabon Hydrogen bond interaction with ASN 152, Pi- Cation bond with the catalytic triad – HIS 51 and Pi-Anion bond with the catalytic triad –ASP 75. The result shows no unfavourable interaction between the ligands and the drugs. In Protein-ligandRef.3 interactions, 8 Cabon Hydrogen bond interactions with GLY 151, PRO 132, SER 135, HIS 51, TYR 150, PHE 130, LEU 128, GLY 153, Interaction with Catalytic traid – HIS 51and SER 135 is seen. Also, one unfavourable interaction is observed.

  1. CONCLUSION

 The incidence of Dengue- mosquito-borne viral fever is increasing dramatically in recent years with   repeated outbreaks from many States and newer areas around the world in recent decades. While many dengue infections are asymptomatic or produce only mild illness, the virus can occasionally cause more severe cases, and even death. To date, there is no universal vaccine available to treat or prevent dengue. However, current research efforts in the development of effective vaccines and therapeutics against dengue are emerging to develop antivirals against dengue.

The present work is based on the use of the phytocompounds extracted from 10 medicinal plants found in Nagaland to inhibit Dengue virus. A total of 133 Phytocompounds were extracted from the selected plants as an anti-dengue agent targeting 2FOM (NS2B/NS3 Protease) of DENV. Virtual screening of the phytocompound was carried out using SWISS ADME to estimate the drug likeness and medicinal potential to inhibit the target protein of DENV. Reference drugs (Doxycycline, Monosdenvir, Zanamivirwhich were previously studied were considered to compare its binding affinity to a target. Molecular docking of the selected phytocompounds and reference drugs was performed using Pyrx. Based on scores, 5 top ligands were selected to investigate the most effective phytocompounds owing to their binding affinity and analyze their efficacy to inhibit the target. Compounds with stronger binding affinities and their mode of interactions between the compounds with highest affinities and their binding site indicate the potential inhibitors. 

The study reveals that Ligand 12 (7-desacetyl-7-benzoylazadiradione) and Ligand 13 (7-desacetyl-7-benzoylgedunin) can be potential inhibitors for DENV (NS2B/NS3 Protease). These phytocompounds are found in Azadirachta indica63,64,65 commonly known as Neem which can be considered as the most potential antiviral plant selected for these studies based on its physiochemical properties. 

The future of in silico studies of phytoconstituents holds great potential for advancing our understanding of these complex molecules, bioactivity, toxicity, prediction of phytoconstituent interactions and mechanisms to increase efficiency in research and development thus, enhancing drug discovery. More anti-viral plants can be considered for comprehensive study of phytocompunds applying increased computational power to simulate complex molecular interactions and analyze large datasets, expansion of databases containing phytoconstituent information to facilitate data sharing and research and visualization of molecular interactions and simulations to enhance understanding. Combination of in silico studies with experimental techniques can also help to validate results and improve accuracy to design new drugs using phytoconstituents with optimized bioactivity and reduced toxicity.

The combined methods of In-Silico studies with other computational methods such as Molecular Dynamic (MD) simulation, can be considered for further research to optimized and give extensive and profound results to provide deeper insights of the anti-viral activity of these natural molecules which can serves as a lead to aid in designing the therapeutic drugs against Virus infections. Besides Molecular Dynamic (MD) Simulation, other computational methods such as Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics General Born Surface Area (MM/GBSA) can be employed to understand the thermodynamics of protein-ligand interactions to predict the binding affinity of small molecules, estimate the binding free energy and calculate the free energy of binding between ligands and proteins which are essential in drug design and discovery. 

Conflict of Interest

The author declares no conflict of interest.

Acknowledgements

I would like to extend my sincere gratitude to St Joseph’s College (Autonomous), Jakhama for providing the facilities to carry out this study. Heartfelt appreciation and gratitude to Dr. Peter Solo (Dean of Science), St Joseph’s College (Autonomous), Jakhama for his constant assistance and valuable insights during the whole process.

Author Contributions: All authors have equal contribution in the preparation of manuscript and compilation.

Source of SupportNil

Funding: The authors declared that this study has received no financial support.

Data Availability Statement: The data presented in this study are available on request from the corresponding author.  

Ethics approval and consent to participate: Not applicable

REFERENCES

  1. Sofowora A, Ogunbodede E, Onayade A. The role and place of medicinal plants in the strategies for disease prevention. Afr J Tradit Complement Altern Med. 2013;10(5):210-29. https://doi.org/10.4314/ajtcam.v10i5.2. PMID: 24311829.
  2. Chaachouay Noureddine, Zidane Lahcen. Plant-Derived Natural Products: A Source for Drug Discovery and Development. Drugs Drug Candidates 2024; 3(1):184-207.  https://doi.org/10.3390/ddc3010011.
  3. Chaughule, R.S., Barve, R.S. Role of herbal medicines in the treatment of infectious diseases. Vegetos. 2024; 37, 41–51. https://doi.org/10.1007/s42535-022-00549-2.
  4. Aruku D. A brief review of medicinal plants used by the Naga tribes of Nagaland, India. International Journal of Life Science Research Archive, 2023; 04(01):143–152.https://doi.org/10.53771/ijlsra.2023.4.1.0025.
  5. Bhuyan SI, Meyiwapangla , Laskar I. Indigenous Knowledge and Traditional Use of Medicinal Plants by Four Major Tribes of Nagaland, North East India.International Journal of Innovative Science, Engineering & Technology. 2014;1(6), ISSN 2348 – 7968.
  6. Garg Amit Kumar, Faheem Mohammed, SinghSumer. Role of Medicinal Plant in Human Health Disease. Asian J. Plant Sci. Res., 2021;11(1):19-21.
  7. Meyanungsang Kichu, Teresa Malewska, Kaisarun Akter, ImchawatiImchen, David Harrington, James Kohen, Subramanyam R. Vemulpad, Joanne F. Jamie. An ethnobotanical study of medicinal plants of Chungtia village, Nagaland, India, Journal of Ethnopharmacology. 2015; 166:5-17.https://doi.org/10.1016/j.jep.2015.02.053
  8. Nanaware N, Banerjee A, Mullick Bagchi S, Bagchi P, Mukherjee A. Dengue Virus Infection: A Tale of Viral Exploitations and Host Responses. Viruses. 2021; 13(10):1967.https://doi.org/10.3390/v13101967
  9. Islam MT, Quispe C, Herrera-Bravo J, Sarkar C, Sharma R, Garg N, Fredes LI, Martorell M, Alshehri MM, Sharifi-Rad J, Daştan SD, Calina D, Alsafi R, Alghamdi S, Batiha GE, Cruz-Martins N. Production, Transmission, Pathogenesis, and Control of Dengue Virus: A Literature-Based Undivided Perspective. Biomed Res Int. 2021;4224816.https://doi.org/10.1155/2021/4224816. PMID: 34957305.
  10. Daep, C.A., Muñoz-Jordán, J.L. &Eugenin, E.A. Flaviviruses, an expanding threat in public health: Focus on dengue, West Nile, and Japanese encephalitis virus. J. Neurovirol.2014; 20, 539–560.https://doi.org/10.1007/s13365-014-0285-z.
  11. Rastogi, M., Sharma, N. & Singh, S.K. Flavivirus NS1: a multifaceted enigmatic viral protein. Virol J. 2016;13, 131.https://doi.org/10.1186/s12985-016-0590-7
  12. Yauch LE, Shresta S. Dengue virus vaccine development. Adv Virus Res. 2014;88;315-72. https://doi.org/10.1016/B978-0-12-800098-4.00007-6. PMID: 24373316.
  13.  Zhi-Shan Zhang, Yan-Sheng Yan, Yu-Wei Weng, Hai-Long Huang, Shi-Qing Li, Shi He, Jian-Ming Zhang. High-level expression of recombinant dengue virus type 2 envelope domain III protein and induction of neutralizing antibodies in BALB/C mice. Journal of Virological Methods. 2007; 143(2):125-131. https://doi.org/10.1016/j.jviromet.2007.02.012.
  14. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, Brownstein JS, Hoen AG, Sankoh O, Myers MF, George DB, Jaenisch T, Wint GR, Simmons CP, Scott TW, Farrar JJ, Hay SI. The global distribution and burden of dengue. Nature. 2013; 496(7446):504-507. https://doi.org/10.1038/nature12060. PMID: 23563266.
  15. James A. Trostle, Charlotte Robbins, Betty Corozo Angulo, Andrés Acevedo, Josefina Coloma, Joseph N.S. Eisenberg, Dengue fever is not just urban or rural: reframing its spatial categorization. Social Science & Medicine, 2024; 117384. https://doi.org/10.1016/j.socscimed.2024.117384
  16. Obi, J.O., Gutiérrez-Barbosa, H., Chua, J.V.Deredge, D.J. Current Trends and Limitations in Dengue Antiviral Research. Trop. Med. Infect. Dis. 2021; 6(180). https://doi.org/10.3390/tropicalmed6040180.
  17. Lauren E. Yauch, Sujan Shresta, Chapter Seven - Dengue Virus Vaccine Development, Editor(s): Karl Maramorosch, Frederick A. Murphy, Advances in Virus Research, Academic Press, 2014; 88:315-372. https://doi.org/10.1016/B978-0-12-800098-4.00007-6.
  18. Sirisena PDN, Mahilkar S, Sharma C, Jain J, Sunil S. Concurrent dengue infections: Epidemiology & clinical implications. Indian J Med Res. 2021;154(5):669-679. https://doi.org/10.4103/ijmr.IJMR_1219_18. PMID: 35532585.
  19. Sun, P., Nie, K., Zhu, Y. et al. A mosquito salivary protein promotes flavivirus transmission by activation of autophagy. Nat Commun. 2020; 11( 260). https://doi.org/10.1038/s41467-019-14115-z
  20. Valenzuela-Leon PC, Shrivastava G, Martin-Martin I, Cardenas JC, Londono-Renteria B, Calvo E. Multiple Salivary Proteins from Aedes aegypti Mosquito Bind to the Zika Virus Envelope Protein. Viruses. 2022;4(2), 221; https://doi.org/10.3390/v14020221. PMID: 35215815.
  21. Murugesan A, Manoharan M. Dengue Virus. Emerging and Reemerging Viral Pathogens. 2020; 281–359. https://doi.org/10.1016/B978-0-12-819400-3.00016-8. PMC7149978.
  22. Lee MF, Wu YS, Poh CL. Molecular Mechanisms of Antiviral Agents against Dengue Virus. Viruses. 2023;15(3):705. https://doi.org/10.3390/v15030705. PMID: 36992414
  23.  Sinha, S., Singh, K., Ravi Kumar, Y.S. et al. Dengue virus pathogenesis and host molecular machineries. J Biomed Sci.2024; 31, 43 .https://doi.org/10.1186/s12929-024-01030-9
  24. Balinsky CA, Schmeisser H, Ganesan S, Singh K, Pierson TC, Zoon KC. Nucleolin interacts with the dengue virus capsid protein and plays a role in formation of infectious virus particles. J Virol. 2013;87(24):13094-106.  https://doi.org/10.1128/jvi.00704-13. PMID: 24027323.
  25. Obi JO, Gutiérrez-Barbosa H, Chua JV, Deredge DJ. Current Trends and Limitations in Dengue Antiviral Research. Trop Med Infect Dis. 2021 Sep 30;6(4):180.https://doi.org/10.3390/tropicalmed6040180. PMID: 34698303.
  26. Nicholls CMR, Sevvana M, Kuhn RJ. Structure-guided paradigm shifts in flavivirus assembly and maturation mechanisms. Adv Virus Res. 2020;108:33-83. https://doi.org/10.1016/bs.aivir.2020.08.003. PMID: 3383772.
  27. Liangzhong Lim, Mei Dang, Amrita Roy, Jian Kang, and Jianxing Song. Curcumin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting Its Active Conformation. ACS Omega 2020; 5(40), 25677-25686.https://doi.org/10.1021/acsomega.0c00039
  28. Lim, S.V., Rahman, M.B.A. & Tejo, B.A. Structure-based and ligand-based virtual screening of novel methyltransferase inhibitors of the dengue virus. BMC Bioinformatics .2011;12 (24). https://doi.org/10.1186/1471-2105-12-S13-S24
  29. Xu T, Sampath A, Chao A, Wen D, Nanao M, Chene P, Vasudevan SG, Lescar J. Structure of the Dengue virus helicase/nucleoside triphosphatase catalytic domain at a resolution of 2.4 A. J Virol. 2005;79(16):10278-88. https://doi.org/10.1128/jvi.79.16.10278-10288.2005. PMID: 16051821.
  30. Chandrakala Basavannacharya, Subhash G. Vasudevan. Suramin inhibits helicase activity of NS3 protein of dengue virus in a fluorescence-based high throughput assay format. Biochemical and Biophysical Research Communications. 2014; 4539(3);539-544.https://doi.org/10.1016/j.bbrc.2014.09.113.
  31. Swarbrick CMD, Basavannacharya C, Chan KWK, Chan SA, Singh D, Wei N, Phoo WW, Luo D, Lescar J, Vasudevan SG. NS3 helicase from dengue virus specifically recognizes viral RNA sequence to ensure optimal replication. Nucleic Acids Res. 2017;45(22):12904-12920.https://doi.org/10.1093/nar/gkx1127. PMID: 29165589.
  32. Wang Y, Xie X, Shi PY. Flavivirus NS4B protein: Structure, function, and antiviral discovery. Antiviral Res. 2022;207:105423.https://doi.org/10.1016/j.antiviral.2022.105423. PMID: 36179934.
  33. Salaemae, W., Junaid, M., Angsuthanasombat, C. et al. Structure-guided mutagenesis of active site residues in the dengue virus two-component protease NS2B-NS3. J Biomed Sci. 2010;17,68. https://doi.org/10.1186/1423-0127-17-68
  34. Teo KF, Wright PJ: Internal proteolysis of the NS3 protein specified by dengue virus 2. J Gen Virol. 1997, 78 (2): 337-341.https://doi.org/10.1099/0022-1317-78-2-337
  35. El Sahili A, Lescar J. Dengue Virus Non-Structural Protein 5. Viruses. 2017;9(4):91.  https://doi.org/10.3390/v9040091. PMID: 28441781.
  36. Ray D, Shah A, Tilgner M, Guo Y, Zhao Y, Dong H, Deas TS, Zhou Y, Li H, Shi PY. West Nile virus 5'-cap structure is formed by sequential guanine N-7 and ribose 2'-O methylations by nonstructural protein 5. J Virol. 2006;80(17):8362-70.  https://doi.org/10.1128/jvi.00814-06. PMID: 16912287.
  37. Wu H, Bock S, Snitko M, Berger T, Weidner T, Holloway S, Kanitz M, Diederich WE, Steuber H, Walter C, Hofmann D, Weißbrich B, Spannaus R, Acosta EG, Bartenschlager R, Engels B, Schirmeister T, Bodem J. Novel dengue virus NS2B/NS3 protease inhibitors. Antimicrob Agents Chemother. 2015;59(2):1100-9. PMID: 25487800.
  38.  Dang M, Lim L, Roy A, Song J. Myricetin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting the Active and Locking the Inactive Conformations. ACS Omega. 2022 11;7(3):2798-2808.https://doi.org/10.1021/acsomega.1c05569. PMID: 35097276.
  39. Alzohairy MA. Therapeutics Role of Azadirachta indica (Neem) and Their Active Constituents in Diseases Prevention and Treatment. Evid Based Complement Alternat Med. 2016;2016:7382506. https://doi.org/10.1155/2016/7382506. PMID: 27034694.
  40. Nawaz, H., Shad, M. A., & Muzaffar, S. Phytochemical Composition and Antioxidant Potential of Brassica. InTech. 2018. https://doi.org/10.5772/intechopen.76120.
  41. Sharma A, Sharma R, Sharma M, Kumar M, Barbhai Mrunal D, Lorenzo M J, Sharma S, Samota Mahesh K, Atanassova M, Caruso G, Naushad Mo., Radha, Chandran D, Prakash P, Hasan M, Rais N, Dey A, Mahato D Kumar, Dhumal S, Singh S, Senapathy M, Rajalingam S, Visvanathan M, Saleena LAbdul Kalam, Mekhemar M. Carica papaya L. Leaves: Deciphering Its Antioxidant Bioactives, Biological Activities, Innovative Products, and Safety Aspects, Oxidative Medicine and Cellular Longevity. 2022; 2451733:20. https://doi.org/10.1155/2022/2451733
  42. 42. Soib HH, Ismail HF, Husin F, Abu Bakar MH, Yaakob H, Sarmidi MR. Bioassay-Guided Different Extraction Techniques of Carica papaya (Linn.) Leaves on In Vitro Wound-Healing Activities. Molecules. 2020; 25(3):517. https://doi.org/10.3390/molecules25030517
  43. 43. Noor Zarina, A. W., &Nazlina, I. (2020). Efficacy of Catharanthus roseus extract against Dengue Virus Type 2 Infection in vitro. Indian Journal of Public Health Research & Development, 11(1), 1314-1319.
  44. Wahab NZA, Azizul A, Ibrahim N. Phytochemistry, cytotoxicity and antiviral activity of Catharanthus roseus. Iran J Microbiol. 2020;12(5):460-465.https://doi.org/10.18502/ijm.v12i5.4608.PMID: 33604002.
  45. K. B. Charuvil, S. Sivanandan and Radha Kesavan LekshmiScreening for Anti-Dengue Leads from Euphorbia hirta L. through In Silico Methods.Indian J Pharm Sci 2022;84(4):950-958. doi:10.36468/pharmaceutical-sciences.989.
  46. Joshi RK, Agarwal S, Patil P, Alagarasu K, Panda K, Cherian S, Parashar D, Roy S. Anti-Dengue Activity of Lipophilic Fraction of Ocimumbasilicum L. Stem. Molecules. 2023 Feb 2;28(3):1446.https://doi.org/10.3390/molecules28031446. PMID: 36771120.
  47. Kaushik, S., Kaushik, S., Dar, L. et al. Eugenol isolated from supercritical fluid extract of Ocimum sanctum: a potent inhibitor of DENV-2. AMB Expr. 2023; 13, 105. https://doi.org/10.1186/s13568-023-01607-x
  48. Dewi BE, Taufiqqurrachman I, Desti H, Sudiro M, Fithriya, Angelina M, Inhibition mechanism of Psidium guajava leaf to dengue virus replication in vitro. IOP Conf. Ser.: Earth Environ. Sci. 2020;462 012034.https://doi.org/10.1088/1755-1315/462/1/012034
  49. Trujillo-Correa AI, Quintero-Gil DC, Diaz-Castillo F, Quiñones W, Robledo SM, Martinez-Gutierrez M. In vitro and in silico anti-dengue activity of compounds obtained from Psidium guajava through bioprospecting. BMC Complement Altern Med. 2019; 6;19(1):298. https://doi.org/10.1186/s12906-019-2695-1.   PMID: 31694638.
  50. Rouf R, Uddin SJ, Sarker DK, Islam MT, Ali ES, Shilpi JA, Nahar L, Tiralongo E, Sarker SD. Antiviral potential of garlic (Allium sativum) and its organosulfur compounds: A systematic update of pre-clinical and clinical data. Trends Food Sci Technol. 2020;104:219-234. https://doi.org/10.1016/j.tifs.2020.08.006. PMID: 32836826.
  51. Shabab Nasir, Keith F.A. Walters, Roberto M. Pereira, Muhammad Waris, Awais Ali Chatha, Munawar Hayat, Marriam Batool. Larvicidal activity of acetone extract and green synthesized silver nanoparticles from Allium sativum L. (Amaryllidaceae) against the dengue vector Aedes aegypti L. (Diptera: Culicidae). Journal of Asia-Pacific Entomology. 2022; 25(3)101937. https://doi.org/10.1016/j.aspen.2022.101937.
  52. Saleh MSM, Kamisah Y. Potential Medicinal Plants for the Treatment of Dengue Fever and Severe Acute Respiratory Syndrome-Coronavirus. Biomolecules. 2021; 11(1):42. https://doi.org/10.3390/biom11010042
  53. K.H. Chiow, M.C. Phoon, Thomas Putti, Benny K.H. Tan, Vincent T. Chow, Evaluation of antiviral activities of Houttuynia cordataThunb. extract, quercetin, quercetrin and cinanserin on murine coronavirus and dengue virus infection. Asian Pacific Journal of Tropical Medicine.2016; 9(1):1-7, https://doi.org/10.1016/j.apjtm.2015.12.002.
  54. Babbar R,  Kaur R,  Rana P,  Arora S,  Behl T,  Albratty M,  Najmi A, Meraya AM, Alhazmi HA, Singla RK. The Current Landscape of Bioactive Molecules against DENV: A Systematic Review. Evidence-based complementary and alternative medicine. 2023https://doi.org/10.1155/2023/2236210
  55. Sunghwan Kim, Jie Chen, Tiejun Cheng, Asta Gindulyte, Jia He, Siqian He, Qingliang Li, Benjamin A Shoemaker, Paul A Thiessen, Bo Yu, Leonid Zaslavsky, Jian Zhang, Evan E Bolton, PubChem 2023 update, Nucleic Acids Research.2023;51(D1);D1373–D1380.  https://doi.org/10.1093/nar/gkac956
  56. Dallakyan S, Olson AJ. Small-molecule library screening by docking with PyRx. Methods Mol Biol. 2015;1263:243-50. https://doi.org/10.1007/978-1-4939-2269-7_19. PMID: 25618350.
  57. Dallakyan, S., Olson, A.J. (2015). Small-Molecule Library Screening by Docking with PyRx. In: Hempel, J., Williams, C., Hong, C. (eds) Chemical Biology. Methods in Molecular Biology. Vol 1263. Humana Press, New York, NY. https://doi.org/10.1007/978-1-4939-2269-7_19
  58. Erbel, P., Schiering, N., D'Arcy, A. et al. Structural basis for the activation of flaviviral NS3 proteases from dengue and West Nile virus. Nat Struct Mol Biol. 2006;13:372–373. https://doi.org/10.1038/nsmb1073
  59. Eberhardt J, Santos-Martins D, Tillack AF, Forli S. AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. Journal of Chemical Information and Modeling.2021;61(8). https://pubs.acs.org/doi/10.1021/acs.jcim.1c00203.
  60. Trott O, Olson AJ, AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading, Journal of Computational Chemistry.2010;21(2): 455-461. https://doi.org/10.1002/jcc.21334
  61.  Daina, A., Michielin, O. and Zoete, V. (2017) SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Scientific Report. 7, 42717.https://doi.org/10.1038/srep42717
  62. Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., Li, Q., Shoemaker, B. A., Thiessen, P. A., Yu, B., Zaslavsky, L., Zhang, J., & Bolton, E. E.Nucleic Acids Res. 2023; 51(D1): D1373–D1380. https://doi.org/10.1093/nar/gkac956
  63. Saleh MSM, Kamisah Y. Potential Medicinal Plants for the Treatment of Dengue Fever and Severe Acute Respiratory Syndrome-Coronavirus. Biomolecules. 2020;11(1):42.  https://doi.org/10.3390/biom11010042. PMID: 33396926.
  64. Parida MM, Upadhyay C, Pandya G, Jana AM. Inhibitory potential of neem (Azadirachta indica) leaves on dengue virus type-2 replication. J Ethnopharmacol. 2002;79(2):273-8. https://doi.org/10.1016/s0378-8741(01)00395-6. PMID: 11801392
  65. Saleem S, Muhammad G, Hussain MA, Bukhari SNA. A comprehensive review of phytochemical profile, bioactives for pharmaceuticals, and pharmacological attributes of Azadirachta indica. Phytother Res. 2018;32(7):1241-1272. https://doi.org/10.1002/ptr.6076. PMID: 29671907.