Molecular docking study
In order to identify the potential drug candidate for managing dengue fever, Molecular Docking analysis was performed for Fucoidan against active site of DENVNS3 proteases of various DEN V Strains. This analysis provides to investigate the nature of binding mode between the protein-ligand complexes.Table 1 reveals the Docking Score ofFucoidan-DENV2 NS2B NS3 protease, Fucoidan-DENV3NS3 protease and Fucoidan-DENV4NS3 protease. From the table it is suggested thatthe fucoidan phytocompounds could be a potent lead compound with upon binding with target site of DENV NS3 proteases.
Table 1. Docking result of Fucoidan phytocompound against DENV NS3 proteases.
S.No
|
Target protein
|
Compound Name
|
Docking Score kcal/mol
|
Glide Energy
kcal/mol
|
Glide emodel
kcal/mol
|
ΔGBind
kcal/mol
|
1
|
DENV2/NS2B NS3 protease
|
Fucoidan
|
-6.458
|
-32.287
|
-48.241
|
-41.25
|
2
|
DENV3NS3 protease
|
-6.987
|
-39.891
|
-41.311
|
-44.38
|
3
|
DENV4 NS3 protease
|
-7.483
|
-40.202
|
-44.291
|
-47.22
|
The fig.1 (A) exhibits the 2D and 3D interactions diagram of fucoidan-DENV2/NS2B NS3 protease complex. The binding conformation forfucoidan-DENV2/NS2B NS3 protease complex showed docking score is -6.458 kcal/mol and forms two hydrogen bonds with HIS51 and GLY151 amino acids with hydroxyl groups of fucoidan. The bond distances for His51 with OH and Gly151 with OH are found to be for 2.58 Å and 1.76 Å respectively. In addition, numbers of hydrophobic interactions are formed With Met49, Val52, Asp75, Leu128, Phe130, Ser131, Pro132, Gly133, Thr134, Ser135, Tyr150, Asn152, Gly153, Val154, Tyr161, Val163 and Ser163 amino acid residues. The docking score for fucoidan-DENV3/NS3 protease complex is -6.987 kcal/mol which has formedfour hydrogen bonds with Phe130 (bond distance 1.72 Å), Gly133 (bond distance 1.89 Å), Thr134 (bond distance 2.73 Å) and Ser135 (bond distance 1.83 Å) with OH group. The complexis also attaching with other amino acids namely His51, Asp129, Lys131, Pro132, Thr150, Gly151 and Tyr161 are from hydrophobic interactions (Fig 1B). The binding mode of fucoidan-DENV4 NS3 protease complex exhibited a docking score of - 7.483 kcal/mol and hydroxyl group of fucoidan build a three hydrogen bonds with amino acids of Val36, His51 and Gly133. The bond distance between these complexes is found to be 2.01 Å, 2.00 Å and 2.26 Å respectively. In addition, this fucoidan phytocompound are connected with six other amino acids such as Gln35, Gly37, Val52, Pro132, Ser135 and Asn152 by developing a hydrophobic interaction (Fig 1C).
The literature reports suggested that the NS3 protease is a trypsin-like serine protease which made up of 6 β-strands that can forms two β-barrels with the catalytic triad of His51, Asp75 and Ser135 amino acids [22, 23]. There are two mechanisms involved for inhibition of protease; firstly ligand can either disturb electronic density or produce steric effect relating to the catalytic triad (His51, Asp75 and Ser135) amino acid residues. Second, it may disturb the movement of the C-terminal of the NS2B necessary for the transition between the "open" and "closed" formation of the protease [24]. In addition, interaction with one of the catalytic triad amino acid residues of protease can help to form the hydrogen bond which can disturb the electron transfer between the carboxyl groups containing Asp75 amino acid and nitrogen atom on the imidazole group containing His51. This happening leads to disruption of residue His51 ability to trigger the nucleophilic attack of the OH group of Ser135, basically required in the initiation of NS3 protease proteolytic activity [25]. In the present study, similar binding interaction with important catalytic triad of His51, Asp75 and Ser135 amino acids was observed and strongly formed hydrogen bond with OH group of fucoidan and hydrophobic interactions. The observed catalytic triad amino acids are actively participated with fucoidan phytocompound which may inhibit the protease activity. Therefore, fucoidan bioactive compound could be explored as a potential inhibitor for the further development of effective DENV drugs.
Binding Free Energy Calculation
Further, the MM/GBSA approach was applied to estimate the Binding Free Energy for fucoidan against NS3 proteases. A negative value of MM/GBSA suggested superior binding of the protein and ligand molecule. The Binding Free Energiesfor fucoidan with DENV2/NS2B NS3, DENV3 NS3 and DENV4 NS3 proteases complexes are found to be 41.25 kcal/mol, -44.38 kcal/mol and -47.22 kcal/mol respectively (Table 1). All the complexes exhibited better binding energy [26]. Therefore, the fucoidan phytocompound is deduced with substantiate binding affinity that can be strongly bound with catalytic triad of DENV NS3 proteases which may inhibit the viral proteolytic activity and prevent DENV fever.
Molecular Dynamic Simulation
Further, to substantiate the structure stability and conformational changes, MD simulation was performedfor docked complexes. Fig 2 (A-C) indicates the RMSD plot of Fucoidan-DENV2/NS2B NS3, Fucoidan-DENV3 NS3 and Fucoidan-DENV4 NS3 proteasecomplexes for 50 ns. Figure A shows the RMSD plot of backbone atoms and Fucoidan-DENV2/NS2B NS3 of protease complex. It seen that the average RMSD value for backbone (black colour) and complex (red colour) are 4.5 ± 0.7 and 5.3 ± 0.85nm respectively. In figure B, during the simulation period, there is no significant change observed between backbone for DENV3 NS3 protease and Fucoidan-DENV3 NS3 protease complex except equilibration stage from 0-5 ns. The average RMSD value for backbone and fucoidan-DENV3 NS3 protease complex are found to be 0.35 ± 0.5 and 0.3 ± 0.8 nm respectively. In the fucoidan-DENV4 NS3 protease complex, no significant structural changes are obtained and similar RMSD value is noticed (0.45±1.0nm).The set of RMSD results indicated that the fucoidan phytocompound interacted with the all proteases are constant stable and no significant deviation observed whole simulation period. Thus, all the complexes are considered for further analysis.
Similar dynamic characteristics for the backbone atoms are also reflected by the RMSF in the MDS. To further assess the protein stability and to realize the effects of fucoidan phytocompound on the NS3 protease structure’s flexibility, the mean per residue fluctuation of the protein backbone was determined using RMSF analysis for the period of 50ns. The RMSF of backbone atom of all the NS3 proteases indicated that the protein exhibited low fluctuation when binding of with fucoidan phytocompound. The average RMSF values for DENV2/NS2B NS3, DENV3 NS3 and DENV4 NS3 proteases are maintained at 0.24 nm, 0.20nm and 0.32nm respectively (Fig 3. A-C). In figure 3 A, the following amino acid residues such as Leu18 (0.5 nm), Glu19 (0.48 nm), Trp50 (0.4 nm), Pro138 (0.36 nm), Gln110 (0.35 nm), Lys128 (0.35 nm), Thr111 (0.34 nm), Lys112 (0.34 nm), and Ser163 (0.34 nm) are in high fluctuations. The fucoidan-DENV3 NS3 protease complex has shown higher fluctuation with the amino acid residues of Val4 (0.76nm), Pro8 (0.74nm), Ser9 (0.59nm), Pro10 (0.46nm), Thr120 (0.45nm), LYS157 (0.44nm) and Asn158 (0.43nm) (Figure 3 B). The fucoidan-DENV4 NS3 protease complex exhibited higher fluctuations with amino acid residues such as Ser19 (0.70 nm), Asp520 (0.61 nm) Lys117 (0.56nm), Leu119 (0. 56nm), Glu20 (0.54nm), Arg28 (0.54nm), Lys157 (0.54nm), Ser159 (0.52nm) and Lys104 (0.51nm) (Figure 3C). From the RMSF plot it reveals that most of the amino acid residues are stable and high fluctuation of amino acid residues are lying in the loop region; those fluctuations did not affect the fucoidan-NS3 protease complexes.
Hydrogen Bond analysis
The Hydrogen Bond interaction analysis is the key tool to determining the molecular recognition and binding strength of the complexes. Fig4 (A-C) depicts the Hydrogen Bond interaction between the fucoidan-DENV2/NS2B NS3, fucoidan-DENV3 NS3 and fucoidan-DENV4 NS3 protease complexes. During the simulation time, all the complexes have indicated the strong and constant binding stability due to the non-conformational changes observed. Therefore, the fucoidan phytocompound could be a potential candidate for inhibiting the NS3 proteases activity.
DFT study
The DFT analysis gives important information about the global and local indices of the fucoidan phytocompound which can helps to examine the chemical reactivity of the fucoidan phytocompound. According to the frontier molecular orbital theory, the electron donor and acceptor properties are determined. Fig5 depicts the HOMO and LUMO distribution plot of fucoidan phytocompound. The figure shows that the energy values of HOMO and LUMO are found to be -0.249 eV and 0.060 eV respectively. In the fucoidan phytocompound, HOMO region is located in the (2S, 3S, 4S, 5R)-4, 5-dihydroxy-2-methyltetrahydro-2H-pyran-3-yl hydrogen sulfate group while the LUMO region is occupied in the (2S, 3S)-1, 3-dihydroxybutan-2-yl hydrogen sulfate group. Form the analysis in infers that the fucoidan phytocompound has high HOMO energy value when compared with LUMO value, suggesting that the phytocompound mainly involved in the nucleophilic reaction [27, 28]. The HOMO-LUMO energy gap (HLG) is obtained at -0.189 eV, confirming that the phytocompound may demonstrate a higher inhibitory activity against the DENV Infection.
ADME properties
This study evaluates the fucoidan phytocompound for pharmacological relevance to ensure the drug likeness and drug’s pharmacokinetics properties [28]. Table 2 gives the ADME properties of the fucoidan phytocompound. The table divulges that all the pharmacokinetic properties for fucoidan has acceptable range due to obey the Lipinski’s rule of five and no violation obtained. Further physicochemical descriptors including QPlogPo/w, QPlogS, QPlogHERG, QPlogBB, and % of Human Oral Absorption are fall in the acceptable ADME range. The overall results reveal that the fucoidan phytocompound is obeying the drug-likeness potteries for human use. Therefore this phytocompound could be considered as a potential antiviral agent against DENV.
Table 2. ADME properties of the fucoidan compound.
S.
No
|
Compound
Name
|
MW1
|
DonorHB2
|
AccepterHB3
|
QPlogPw4
|
QPlogS5
|
%HOA6
|
QPlog
HERG7
|
1
|
Fucoidan
|
242.243
|
3.00
|
9.00
|
14.343
|
-1.011
|
65.387
|
-1.368
|
BOILED-Egg plot study
For the drug discovery process, gastrointestinal absorption and BBB penetration play an imperative role. This plot helps in the computation of polarity and lipophilicity of fucoidan phytocompound. In the analysis, the white region is indicates the gastrointestinal tract and yellow region notices the high polarity of the BBB. In addition, the blue color indicator reveals that the compound is actively effluxes by P-glycoprotein (PGP+), while the red color indicator exhibited the non substrate of Pgp (PGP−) [30]. In the present analysis, the fucoidan phytocompound is located in the white region, suggesting that the fucoidan phytocompound has highly gastrointestinal absorption which is easily absorbed by the gastrointestinal tract and do not exhibit BBB permeation (Fig6).