1. Introduction
Dengue flavivirus (DENV) is a virus transmitted by mosquitoes and occurs in more than 100 countries. In the vast majority of cases, DENV presents in humans as a mild disease characterized by fever, myalgia, arthralgia and retroocular headache, which improve without intervention or any complication; however, severe forms can develop due to plasma leakage, which can quickly be fatal. DENV is a tropical disease; its highest amount of growth has occurred over the last three decades. DENV has an annual incidence of approximately 50 million cases. Worldwide, it has threatened over 50% of the population. Specifically, its mortality rate in the Americas is approximately 1.2%. Therefore, in the absence of an effective vaccine or therapeutic approach, there is an urgent need to develop new therapies [1,2,3].
Dengue virus has an 11 kb positive-sense single-stranded genomic RNA encoding a precursor polyprotein (5′-C-prME-NS1-NS2A-NS2B-NS3-NS4A-NS4BNS5-3′) that is cleaved to produce three structural proteins (capsid protein C, membrane protein prM and envelope protein E) and seven nonstructural proteins (NS1, NS2a, NS2b, NS3, NS4a, NS4b and NS5). Of these, the NS2B/NS3 protease complex inhibits the host type I interferon (IFN) pathway by cleaving a key adaptor protein, the human mediator of IRF3 activation (MITA). The NS2B/NS3 protease is thus required for viral replication and promotes the progression of DENV infection [4,5,6,7]. Therefore, it has been selected as an important target for the design of drugs against dengue infection [8,9].
A promising alternative is the design of small molecules directed to the allosteric site. Allosteric sites are defined as regions of the protein that, when linked to a small-molecule ligand, undergo a change in conformation or an alteration in their conformational balance that affects enzyme function. In the literature, several articles are found exploring the allosteric sites of DENV NS2B/NS3pro in search of more potent inhibitors [10,11,12,13,14,15]. In DENV NS2B/NS3pro, this binding site is located behind the active site and is formed mainly by the residues NS3-Asp71, NS3-Lys73, NS3-Lys74, NS3-Trp83, NS3-Leu85, NS3-Gly87, NS3-Glu88, NS3-Trp89, NS3- Glu91, NS3-Thr118, NS3-Thr120, NS3-Val147, NS3-Leu149, NS3-Asn152, NS3-Val155, NS3-Ala164, NS3-Ile165 and NS3-Asn167 [12].
The fact that the DENV NS3pro active site is planar and preferentially binds to peptide-based substrates, such as the competitive peptidomimetics inhibitors of NS2B/NS3pro DENV found in the literature, has conferred limited oral bioavailability and difficulties in in vivo assays in studies of potential inhibitors targeting the active site of this enzyme [13,16,17]. However, this strengthens the strategy of exploring allosteric sites that are known to be successful targets of small molecule inhibitors, as presented in the works of [13,14,15,16].
Because of the need to develop more effective therapeutic agents and insurance against DENV, the aim of this work is to carry out the virtual screening of inhibitor-like compounds, analogous to tiaguanine, proposed by Hariono et al. [18]. In addition, we plan to target these compounds to the allosteric site in order to design DENV NS2B/NS3pro non-competitive inhibitors. Thus, we performed molecular docking, molecular dynamic simulations and free-energy calculations. We hope to contribute to the efforts in discovering novel and potent anti-DENV agents.
2. Results and Discussion
2.1. Molecular Docking
In this work, 71 ligands were selected and submitted to molecular docking simulations where the ranking was performed according to the binding affinity energy values obtained by the PLPchem scoring function of the GOLD program. The four best ligands (highest energy) were selected to assess the interactions between the residues of the NS2B/NS3pro (PDB 2FOM) allosteric site and the ligand (Table 1).
From the docking results, several interactions were identified regarding the residues of the NS2B/NS3pro allosteric site with ZINC000001680989, ZINC000001679427, ZINC000006694490 and ZINC000001684485. For the ligand ZINC000001680989, hydrogen bond interactions were observed with the Leu179(149), Asn182(152) and Ala194(164) residues, and hydrophobic interactions were observed with the Leu104(74), Leu106(76) and Trp113(83) residues; the latter also resulted in π–π interactions. For ZINC000001679427, hydrogen bond interactions with the residues Lys104(74), Leu179(149) and Ala194(164) were observed; hydrophobic contacts were also reported with the residues Trp113(83), Gly178(148), Val184(154) and Ala194(164).
For ZINC000006694490, hydrogen bond interactions with the residues Asn182(152), Thr150(120) and Val185(155)—along with hydrophobic interactions with Leu106(76) and Val184(154)—were observed. ZINC000001684485 interacted with hydrogen bonds in Thr150(120), Asn182(152) and Asn197(167), and hydrophobic interactions were observed with Val184(154) and Thr150(120). These interactions are shown in Figure 1. The numbers in the parentheses indicate the numbering within the Erbel model [19].
Notably, all of the ligands interacted with residues within the DENV NS2BNS3pro allosteric cavity [12]. The ZINC000001680989, ZINC000001679427 and ZINC000001684485 binders interacted with the Asn182(152) residue, which is considered essential for the interaction with NS2B/NS3pro allosteric inhibitors and is in accordance with other molecular coupling studies that involve DENV NS2B/NS3pro allosteric inhibitors [12,13,20]. This suggests that these ligands may show inhibitory activity toward the NS2B/NS3pro allosteric site. In particular, it was observed that the ligands ZINC000001680989 and ZINC000001679427 formed hydrogen bond interactions with the residues Asn182(152), Leu179(149) and Ala194(164), which is in agreement with the results reported by Mukhametov et al. [12] and Millies et al. [13] for DENV NS2B/NS3pro allosteric inhibitors in their studies in silico and in vitro. Another fact that corroborated this result is that these binders presented the best PLPChem score energies according to the PLPchem scoring function.
2.2. Molecular Dynamics
To prove the stability and conformational changes, MD simulation was performed to understand the characteristics of the NS2B/NS3pro-ZINC000001680989, NS2B/NS3pro-ZINC000001679427, NS2B/NS3pro-ZINC000006694490 and NS2B/NS3pro-ZINC000001684485 ratios in relation to the nanosecond scale. According to Figure 2, the backbone RMSDs for the NS2B/NS3pro-ligand complexes showed that the selected ligands remained stable at the protein allosteric site, converging at an RMSD value below 2 Å, while in the native NS2B/NS3 protein, the RMSD fluctuation is significantly greater, suggesting the inhibitory effect of the ligands used in this work. Fluctuations between 1 and 3 Å within a reference protein structure are perfectly acceptable and reflect complex stability [21].
For a more refined visualization, the average RMSD values were extracted with their error bars. The results suggested that the protein/ZINC0000016799427 complex had a better stabilization in relation to the other complexes (Figure 3).
In addition, an analysis of Figure 3 showed that the RMSD values of the NS2B/NS3pro-ZINC000001679427 and NS2B/NS3pro-ZINC000001680989 complexes fluctuated within a range of ±1.2 Å over time, thus exhibiting the greater stability of the protein complex link during the simulation [22]. Notably, these binders presented the best PLPChem score and interacted with the essential residues for inhibition at the DENV NS2B/NS3 allosteric site [12,13,16].
Local fluctuations in the protein are measured by the mean quadratic fluctuation (RMSF) and are useful for characterizing local changes along the protein chain [23]. The RMSF values of NS2B/NS3pro after binding to the four ligands are shown in Figure 4. The peaks in the graph indicate the amino acid residues in the protein that fluctuate the most during the MD simulation. Typically, the tails (N- and C-terminals) fluctuate more than any other part of the protein. Secondary structure elements such as alpha helices and beta strands are usually more rigid than the unstructured part of the protein and thus fluctuate less than the loop regions [16]. In the NS2B/NS3pro-ZINC000001680989, NS2B/NS3pro-ZINC000001684485 and NS2B/NS3pro-ZINC000006694490 complexes, the rate of fluctuation is higher at residue No. 32 compared to that in the NS2B/NS3pro-ZINC000001679427 complex. In the other residues, the fluctuations of the four complexes remained almost the same. All of the protein–ligand interactions had a fluctuation rate of less than 2 Å, which is completely acceptable. The complexes were considered stable throughout the MD simulation period. These plotted data indicate that all of the complexes showed active interaction with the protein and remained stable at the allosteric site during the simulation. It is interesting to note that, in the native NS2B/NS3 protein (black line), the RMSD fluctuation is significantly greater, showing the importance of the inhibitors effect at the system’s stabilization, which suggests NS2B/NS3 enzyme inactivation.
The calculation of the radius of gyration (Rg) was performed to evaluate the protein–ligand systems’ stability by calculating the structural compactness along the MD trajectories [24]. After the MD simulation calculation, the calculation of Rg was also used to determine the stably folded or unfolded protein and complexes system. The average values of RG for all of the complexes were calculated and were around 16.032, 16.111, 16.124 and 16,093 Å for NS2B/NS3Pro-ZINC000001679427, NS2B/NS3Pro-ZINC000001680989, NS2B/NS3Pro-ZINC0001684485 and NS2B/NS3Pro-ZINC00000694490, respectively (Figure 5). All of the complexes exhibited relatively similar and consistent values of Rg, indicating systems stability.
The solvent accessible surface area (SASA) was also determined to evaluate the constancy of simulated systems and to determine how much of the receptor area was exposed to the solvents. A larger SASA value indicates the growth of protein volume during the MD simulation. The interaction of the ligands in the active site can change the SASA and protein folding [25]. In the case of the complexes NS2B/NS3pro-ZINC000001679427 and NS2B/NS3pro-ZINC000001680989, it was found to be around 8556.41 A2 and 8645.09 A2, respectively (Figure 6). Likewise, the complexes NS2B/NS3pro-ZINC000001684485 and NS2B/NS3pro-ZINC000006694490 showed average values of SASA around 8916.33 A2 and 8711.76 A2. These calculations suggest that all of the complexes were least exposed to the water solvent during the 100 ns MD simulations, which indicates the relatively stable nature of these complexes.
To better understand the mode of interaction of these new inhibitors, it is important to estimate their binding affinities (ΔGbind) and the values of the other energy contributions. Hence, the MM-GBSA and MM-PBSA methods were used to calculate the free energy. These methods have been used successfully in the reproduction of experimental results [26] and in studies of the allosteric inhibitors of DENV NS2B/NS3pro [12].
The ΔGbind and the values of the energy contributions of van der Waals (ΔEvdW), electrostatic (ΔEele), polar (ΔGGB) and nonpolar (ΔGnonpol) forces are summarized in Table 2 (all results in Kcal/mol). Using the MMGBSA approach, the ΔGbind contributions ordered the binding affinity of the ligands to the DENV NS2B/NS3pro allosteric site as ZINC000001680989 > ZINC000001679427 > ZINC000006694490 > ZINC000001684485. The MMPBSA approach ordered them as ZINC000001680989 > ZINC000006694490 > ZINC000001679427 > ZINC000001684485 (Table 2). This classification is in accordance with the values from the molecular docking, which suggests the same classification for these compounds.
The results obtained in this study demonstrate that the ligands coupled to the NS2B/NS3pro allosteric site had greater affinity compared to the results obtained using the free energy by Hariono et al. [18], where the matrix ligand presented an energy of −16.37 ± 3.22 kcal/mol by the MMPBSA method. This suggests the potential of these compounds to inhibit the allosteric site in study. Mukhametov et al. [12], in their in silico and in vitro studies of DENV NS2B/NS3pro allosteric inhibitors, predicted the ΔGbind values of the compound 2-hydroxy3-methoxybenzaldehyde N-adenosinylhydrazine (compound 10) protease (65.1%) to be −24.02 kcal/mol and −14.45 kcal/mol, respectively, by the MMGBSA and MMPBSA methods. The ΔGbind results for the ligands ZINC000001680989, ZINC000001679427 and ZINC000006694490 obtained in this study indicate the better affinity of these ligands with DENV NS2B/NS3pro (Table 2), a fact that refocuses our hypothesis regarding allosteric inhibitors. In general, the protein–ligand interaction is stabilized by several types of weak interactions, and the hydrogen bond interaction is potentially one of the most important. Table 3 shows the hydrogen bonds between the ligands ZINC000001680989, ZINC000001679427, ZINC000006694490 and ZINC000001684485, with residues from the NS2B/NS3pro allosteric site.
The Leu179(149), Ala 194(164) and Ans182(152) residues in NS2B/NS3pro played an important role in the formation of hydrogen bonds that remained occupied for more than 100%, 99% and 23% of the simulation time with ZINC000001680989, respectively. The ligand ZINC000001679427 also formed hydrogen bonds with the residues Ala194(164), Leu179(149) and Ans182(152) and was occupied for more than 73%, 50% and 14% of the simulation time, respectively. The ligand ZINC000006694490 interacted with the Asn182(152) and Lys104(74) residues and was occupied for more than 62% and 17% of the simulation time, respectively. For the ligand ZINC000001684485, the most significant interaction was with the residue Lys104(74), and it was occupied for 100% of the simulation time.
It is important to highlight that, in the study by Hariono et al. [18], the matrix compound presented the highest occupancy of hydrogen bonds during the simulation time, 60%, a lower number than the numbers the ligands ZINC000001680989, ZINC000001679427, ZINC000006694490 and ZINC000001684485 presented with the residues of the DENV NS2B/NS3pro allosteric site.
Therefore, the greater occupation of the hydrogen bonds of ZINC000001680989 with the NS2B/NS3pro allosteric residues may explain the better bond-free energy observed for the DENV NS2B/NS3pro-ZINC000001680989 complex (−40.65 ± 0.09 kcal/mol) compared to those observed for the DENV NS2B/NS3pro-ZINC000001679427 (−35.53 ± 0.11 kcal/mol), DENV NS2B/NS3pro-ZINC000006694490 (−27.76 ± 0.09 kcal/mol, and DENV NS2B/NS3pro- ZINC000001684485 (−11.7 ± 0.14 kcal/mol) complexes.
It is important to highlight that the interactions of ZINC000001680989 and ZINC000001679427 with the residues Leu179(149), Asn182(152) and Ala194(164) of the NS2B/NS3pro allosteric site had already been observed in the results of the molecular docking. Mukhametov et al. [12], in their in silico and in vitro studies on DENV NS2B/NS3pro allosteric inhibitors, reported that the most potent inhibitor (compound 10) formed hydrogen bonds with the residues Leu179(149), Asn182(152) and Ala194(164), which proved the high binding potential of the two ligands to the NS2B/NS3pro allosteric site. As was already mentioned, the hydrogen bond interaction with the Asn182(152) residue is considered essential in the inhibition process and corroborates the in silico and in vivo studies of [17,27]. The residues Leu149, Asn152 and Ala164 in this protease play key roles in the allosteric inhibition of the active site when binding with the ligand by bringing out significant conformational changes and disabling the protein to make its function [12].
In order to assess which residues contributed to the stabilization and maintenance of the simulated compounds in the NS2B/NS3pro allosteric site, an analysis of the free energy decomposition by the residue was conducted (Figure 7). In the graph, positive values correspond to unfavorable interactions, while negative values correspond to favorable interactions, which are those that contributed to the stabilization of the ligand in the complex.
An analysis of the interaction energy graph for every residue revealed that the residues that contributed most significantly to the total interaction energy and, therefore, to the stabilization of the complex were Lys73, Lys74, Leu76, Trp83, Leu85, Ile123, Val147, Gly148, Leu149, Asn152, Val154, Ala164, Ile165, Ala166 and Asn167 for ZINC000001680989. Likewise, Lys73, Lys74, Leu76, Trp83, Leu85, Ile123, Asn152, Val154, Ala164, Ile165 and Ala166 contributed most significantly for ZINC000001679427. For ZINC000006694490, the residues Lys74, Asp75, Leu76, Trp83, Ile123 and Asn152 were observed, as also evinced by the hydrogen bond (Table 3), Val154, Ala164 and Asn167. The ligand ZINC000001684485 showed an unfavorable interaction with the Glu88 residue and a favorable interaction with the residues Lys74, Leu76, Gly148, Leu149, Asn152, Ile165 and Ala166.
The per-residue interaction profile indicates that the four ligands present favorable interactions with the residues from the NS2B/NS3pro allosteric site, which is crucial for the protein–ligand complex’s stability. The results of the MD simulation and the MM/GBSA analysis validate the molecular docking results. The negative binding free energy of all four complexes demonstrates their stable configuration, which indicates that these compounds have enough affinity with NS2B/NS3pro DENV to be considered as allosteric inhibitors.
The greater number of residues with significant contributions from the ZINC000001680989 binder also corroborated the hypothesis of it having greater binding affinity. The hydrogen bond interactions, together with the hydrophobic clusters, contributed to the greater affinity of ZINC000001680989 bonding, with special emphasis on the hydrogen bond formed with the Leu179 residue, which also contributed to ZINC000001679427.
In particular, by analyzing the ligand ZINC000001680989, which showed the best binding affinities (ΔGbind) based on the MM/GBSA method, it was possible to observe interactions with the residues Lys73, Lys74, Leu149 and Asn152, which, according to experimental and theoretical studies, are essential for protease inhibition [12,13,17,27]. Othman and collaborators [28] observed that interactions with the residues Lys74 and Leu149 explained the differences in the inhibition activities of non-competitive inhibitors in their studies. This interaction with Lys74 is directly linked to Asp75, inducing a change in the region of the catalytic triad. This, presumably, could interrupt the electron transfer process necessary for the binding of the substrate at the active site, thus affecting the activity of the protease. Thus, the results obtained reinforce our strategy of using this class of compounds as allosteric inhibitors of DENV NS2B/NS3pro.
Another highlight was that the residues of Lys74, Leu76 and Asn152 appeared in all of the simulated systems. This indicated that the interaction of these residues with the ligands in the allosteric site was essential for their maintenance and stability. It is possible to explore these interactions in future studies.
It is interesting to note that, when we have conducted visualization of the frames corresponding to the times of 0, 25, 50, 75 and 100 ns, we have noted that the ligands ZINC000001680989, ZINC000006694490 and ZINC000001679427, according to the MMGBSA energies, remain in the enzyme active site (Figure 8).
2.3. SwissADME Prediction of Selected Compounds
The data predicted for the lipophilicity, druglikeness, solubility, physicochemical characteristics, pharmacokinetics and medicinal chemistry of the selected compounds evaluated by SwissADME are given in Table 4. According to Lipinski′s rule of five, the molecular weights of the selected compounds reached between 313.42 Da (ZINC000001679427) and 369.44 Da (ZINC000001684485), respectively, and within the limits up to 500 Da. The log Po/w values were less than 5 and were in the range of 2.33 (ZINC000001684485) up to 3.50 (ZINC000001680989), suggesting a good potential for lipophilicity. The HBA number of all the selected products was between 3 and 5 and less than 10; the number of HBD atoms was 3 (ZINC000006694490) and 1 for all other hits, respectively—both less than 5 [29,30].
For physicochemical properties such as molecular refractivity, the hit ZINC000001680989 showed a value of 101.23, indicating a greater size and polarizability. The surface sum over all of the polar atoms or molecules, primarily oxygen and nitrogen, including their attached hydrogen atoms, was evaluated from the topological polar surface area (TPSA) [31], in which the hit ZINC000006694490 exhibited an index of 120.89, the highest among the compounds selected.
The potential of the selected hits for water solubility was evaluated from Log S (SILICOS-IT), a hybrid method relying on 27 fragments and 7 topological descriptors. The solubility class: Log S Scale: insoluble < −10 poorly < −6, moderately <−4 soluble < −2. Only the hit ZINC000001680989 showed a low potential for solubility in water, with a Log S = −6.02. All of the other compounds exhibited potential values for moderate solubility [32].
For the druglikeness properties, the evaluated hits did not show any variation from Linpinski′s rules, and all of the hits had a bioavailability score of 0.55. The bioavailability score is similar but seeks to predict the probability of a hit having at least 10% oral bioavailability in rats or measurable Caco-2 permeability. This semi-quantitative rule-based score relying on total charge, TPSA and the violation of the Lipinski filter defines four classes of compounds with probabilities of 11%, 17%, 56% and 85% [30,33].
The skin permeation coefficient (kp) in pharmacokinetics calculates the conductance of skin to a particular chemical from a specific aqueous vehicle using three different models. The ZINC000001684485 presented the lowest value (−6.91) in the series of compounds, while the index of −5.17 (ZINC000001679427) can be considered the highest in the series. It is worth noting the little statistical significance between the values of all the selected hits.
All of the selected compounds showed negative indices for the potential penetration of the blood brain barrier. On the other hand, all of the hits studied expressed high human intestinal absorption (GI). All of the compounds can be promising agents that can very easily be absorbed by the gastrointestinal tract without potential BBB permeability. Since these compounds cannot cross the BBB, they do not cause central nervous system depression or drowsiness as side effects [30].
Regarding the parameters of SwissADME medicinal chemistry, the synthetic accessibility (SA) score is basically based on the assumption that the frequency of molecular fragments in ‘really’ obtainable compounds correlates with the ease of synthesis. The fragmental contribution to SA should be favorable for frequent chemical moieties and unfavorable for infrequent moieties. The scale score is normalized to range from 1 (very easy) to 10 (very difficult to synthesize) [34]. All of the compounds selected in this study showed indexes lower than 4, especially ZINC000001679427, which had an index lower than 3, corroborating the presence of chemical groups that are not rare in the structure of the selected molecules and suggesting the ease of new syntheses.
2.4. Toxicological Analysis by Derek NexusTM
The compound ZINC000006694490 exhibited a warning with a plausible level of toxicity for the purine base chemical moiety. In this analysis, reports were identified in the literature for hepatotoxicity and chromosome damage in vivo in mouse, rat and human organisms.
This alert describes the hepatotoxicity of purine derivatives. Purines have been shown to cause severe liver toxicity in humans and animals. Steatosis, cholestatic jaundice, necrosis and veno-occlusive disease were reported in patients treated with purine derivatives during immunosuppression or antiretroviral therapies [35]. In animals, elevations of the liver enzymes and damage of the centrilobular hepatocytes were observed [36]. The mechanism of toxicity of purines has dose-dependent and accumulative characteristics, but hypersensitivity reactions have also been reported after a short time of treatment [37].
The prolonged exposure of HIV-positive patients to didanosine (daily dose: 400 mg) during antiretroviral therapy led to the elevation of liver enzymes, various degrees of micro or macrovesicular steatosis and mild cholestasis [38,39]. Centrilobular Mallory bodies have also been observed [40]. In children, therapy with didanosine (120–270 mg/m2) caused severe hepatic injury such as necrotising hepatitis, which in some cases led to fatal hepatic failure [41]. Patients undergoing immunosuppression therapy with azathioprine (daily dose: 50–2500 mg) for the management of autoimmune diseases such as inflammatory bowel disease, rheumatoid arthritis, lupus erythematous and chronic hepatitis developed liver injury manifested with the elevation of liver enzymes, icterus and, in some cases, fever, diarrhea and vomiting [42,43,44]. The liver biopsies of these patients showed various degrees of liver damage from canalicular cholestasis, with ballooning centrilobular hepatocytes, necrosis with portal infiltrates and cirrhosis [43,44].
2.5. Chromosome Damage: In Vivo Chromosome Aberration Test; In Vivo Micronucleus Test
This alert describes the activity of exogenous purine base, nucleoside and nucleotide analogues in the in vivo chromosome aberration and micronucleus tests. Such compounds are widely used in medicine as antiviral or antineoplastic therapies and have been shown to induce chromosome damage in vivo. For example, positive results in the mouse bone marrow micronucleus test have been reported for acyclovir, ganciclovir and penciclovir when administered intravenously [45] and for 6-mercaptopurine [46] and 6-thioguanine [NTP] following the administration by oral gavage.
In addition, 6-mercaptopurine demonstrated activity when tested in the mouse bone marrow chromosome aberration test [47]. Didanosine (2′,3′-dideoxyinosine) failed to induce micronuclei in mouse bone marrow when administered orally but demonstrated activity following the administration by intraperitoneal injection [48]. This may be explained by the decomposition of didanosine under the acidic conditions of the stomach [49].
There are several mechanisms that may contribute to the clastogenicity of these compounds. Exogenous base and nucleoside analogues may become incorporated into mammalian DNA [50]. Once incorporated, the type of mechanism that leads to cytogenetic damage may be dependent upon the structural modification in the analogue. For example, ring-modified purines may pair with an incorrect nucleotide on the complementary DNA strand and induce mutations. By contrast, compounds with modified sugar moieties may distort DNA, causing repair and replication errors, or induce the premature termination of growing DNA chains [50].
3. Materials and Methods
3.1. Virtual Screening for Ligand Selection
Initially, we selected the S-[2-(pentanoylamino)-9H-purin-6-yl]pentanethioate(1), a thioguanine analogue (Figure 9) that was designed as a DENV NS2B/NS3pro inhibitor (IC50 of 0.38 μM) in a study conducted by Hariono and collaborators [18]. This compound served as a matrix in the present study.
After that, compounds structurally similar to 1 were searched for using the Swiss Similarity web tool (
Swiss Similarity is based on the LBVS hypothesis, and similar molecules are prone to exhibit similar biological activity. This tool allows for the rapid screening of large-scale drug libraries in more than 30 chemical databases including more than 205 million virtual compounds synthesizable from commercially obtained reagents [51]. This tool has been used for in silico studies in drug planning [52].
3.2. Molecular Docking
The 71 ligands selected in the previous step were subjected to molecular coupling. The DENV NS2B/NS3pro crystallographic model was retrieved from the Protein Data Bank with the code PDB 2FOM and a resolution of 1.50 Å [19]. To perform the molecular docking simulations, the CSDGOLD program [53] was used, and the docking grids with a 10 Å radius were positioned on the specific allosteric binding site of the non-competitive inhibitor DENV NS2BNS3pro according to [12,13]. For the CSD-GOLD protocol, we set the following parameters for the ChemPLP algorithm: all water molecules and ions have been removed. The CSDGOLD program uses the empirical aptitude function called ChemPLP, which consists of the application of terms of hydrogen and metal bonding and linear potential by parts (PLP) to model the steric complementarity between the protein and the ligand. The dimensionless scoring scale measures the success of the pose; higher scores indicate better nesting positions [54].
3.3. Molecular Dynamics Simulation
The best docking results will be subjected to molecular dynamics (MD) simulation in order to monitor the behavior of the systems in the function of time [55]. The impact of MD simulations on molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules completely at the atomic level and in very fine temporal resolutions. Major improvements in terms of simulation speed, accuracy and accessibility, along with the proliferation of structural and experimental data, have increased the search for biomolecular simulations for experimentalists [56].
To analyze the conformational changes of the protein and ligand structures, as well as the stability of the ligand-receptor complexes, we performed MD simulations, accelerating the graphics processing unit (GPU) in Amber18 [57]. The Amber ff14SB [58] and the general amber force field (GAFF) [59] were selected as the parameters to describe the protein and ligands, respectively. The restrained electrostatic potential (RESP) charges for the ligands were calculated at the HF/6-31G(d) level of theory using the Gaussian 09 program [60]. The systems were solvated in a cubic water-box with the explicit solvation model TIP3P using a distance of 12.0 Å with periodic boundary conditions. The Na+ ions were added to maintain the electroneutrality of the systems. To eliminate the steric clashes, each system was subjected to three minimization steps. First, all hydrogen atoms were minimized by 2000 steps of the steepest descent, followed by 3000 steps of the conjugate gradient algorithm. Next, using the same protocol, the positions of the water molecules were relaxed. Lastly, the whole system was energy-minimized for 5000 steps of the steepest descent plus 5000 steps of the conjugate gradients. Afterward, we started the heat of the system from 0 to 300 K, running 200 ps of MD and then 300 ps of density balance with position restraints on the protein–ligand atoms at a constant volume. Before performing the production step, all of the protein–ligand systems were balanced with 500 ps of MD without positional restraints at a constant pressure. The temperature was maintained at 300 K by coupling to a Langevin thermostat using a collision frequency of 2ps-1. A cutoff of 10 Å was employed for the nonbonded interactions, and the Particle Mesh Ewald (PME) [61] method and the SHAKE [62] algorithm were used to restrict the bond lengths involving the hydrogen atoms. The MD simulations (production) were performed using 100 ns at a temperature of 300 K. The structural analysis of each system was performed by using the root means square deviation (RMSD) values of the heavy atoms of the protein′s backbones, the root mean square fluctuation (RMSF), the radius of gyration (RG) and the solvent accessible surface area (SASA) obtained by the CPPTRAJ module of Amber18.
3.4. Binding Free Energy Calculations MM/PBSA and MM/GBSA
The values of free binding energy (ΔGbind) for each system were calculated using the MM/PBSA and MM/GBSA approaches implemented in the MMPBSA.py module. For this, the last 10 ns of the DM simulation were used.
The equations that describe the calculations of the energy are:
∆Gbind = ∆H − T∆ ≈ ∆EMM + ∆Gsol + T∆S,(1)
∆EMM = ∆Einternal + E∆electrostatic + ∆Evdw (2)
∆Gsol = ∆GPB/GB + ∆GSASA (3)
As shown in Equation (1), the enthalpy part is expressed as the summation of the molecular mechanical energy (∆EMM) and the solvation free energy (∆Gsol), where ∆EMM is composed of the intra-molecular energy (∆Einternal, including the bond, angle and dihedral energies of the system), the electrostatic energy (∆Eele) and the van der Waals interactions (∆EvdW) (Equation (2)). The solvation free energy (∆Gsol) is also composed of two parts, namely, the polar part (∆Gpol) and the non-polar (∆Gnp) part (Equation (3)), where ∆Gpol is usually computed by the Generalized Born (GB) model or by solving the Poisson–Boltzmann (PB) equation, while ∆Gnp is estimated by the solvent accessible surface area (SASA)-based approach [26,63].
3.5. Decomposition Energy by Residue
The method of the decomposition of energy by residue was used to determine the total energy contribution of each residue to the drug–receptor interaction and also to investigate the chemical nature of its interactions [64].
The interaction energy of the residues can be described from four terms: the electrostatic contribution of van der Waals (ΔEvdw), the electrostatic contribution (ΔEele) in the gas phase, the contribution of polar solvation (ΔGpol) and the contribution of non-polar solvation (ΔGnonpol), according to the equation:
∆Ginhibitor = ∆Evdw +∆Eele + ∆Gpol + ∆Gnonpol(4)
3.6. Pharmacokinetics Analysis by SwissADME
The computational studies of the compounds from the ZINC database 1–3 were selected to predict the molecular properties using the SwissADME online server (
3.7. Toxicity Analysis by DerekTM
The predictions described below were made on a computer with six logical processors of an Intel® Core TM i5-9400, 4.10 GHz processor using the Windows 10 Professional operating system.
Derek Nexus is an ‘expert rule-based’ system for the prediction of toxicity [68,69,70] that uses a virtual base of molecules with reported toxicity alerts—for example, mutagenicity, genotoxicity, teratogenicity, carcinogenicity and hepatotoxicity—for comparison with the potential toxic group present in each query molecule [6]. The program generates an alert based on the literature evidence and describes the potential toxicity for the complete structure [33].
4. Conclusions
In the present study, in order to investigate the binding affinity, selectivity and stability of candidates for allosteric inhibitors of the DENV NS2B-NS3pro enzyme, the Swiss Similarity web tool was utilized in order to search for compounds structurally similar to ligand 1 described by Hariono et al. Within this web server, 71 compounds were selected for molecular docking studies using the GOLD program. The four best-scored ligands according to the scoring function energies were then subjected to an MD protocol involving 100 ns of simulation in the AMBER package. The free bond energies, hydrogen bonds and energy contributions per residue were also calculated for the simulated systems. The computational approach used here proved to be useful for planning new inhibitors to combat dengue. In the simulations, the residues of Lys74, Leu176 and Asn152 contributed to the stability of the ligand interactions at the allosteric site in all of the simulated systems. Meanwhile, the Leu149 and Ala166 residues, in addition to contributing significantly to the energy decomposition by the residue, showed high rates of hydrogen bond permanence. This indicated that these interactions are essential for maintenance and stability.
The selected compounds showed negative indices for the potential penetration of the blood brain barrier and expressed high human intestinal absorption. They can be promising agents that can very easily be absorbed by the gastrointestinal tract without potential BBB permeability, indicating a low risk of central nervous system depression or drowsiness as side effects. All of the compounds selected in this study showed synthetic accessibility (SA) scores lower than 4, especially ZINC000001679427, with an index lower than 3, corroborating the presence of chemical groups that are not rare in the structure of the selected molecules and suggesting the ease of new syntheses.
In the toxicological analysis, the ZINC000006694490 compound exhibited an alert with a plausible level of toxicity for the purine base chemical moiety, indicating hepatotoxicity and chromosome damage in vivo in mouse, rat and human organisms. These results agree with those of previous studies. Hence, there is scope for further work related to the selective planning of candidates for DENV NS2B-NS3pro inhibitors.
Conceptualization, R.A.d.C., J.A.P.d.R. and D.d.S.B.B.; methodology, R.A.d.C., J.A.P.d.R., A.S.P., R.C.S. and D.d.S.B.B.; software, R.A.d.C., J.A.P.d.R. and A.d.S.G.; validation, R.A.d.C., A.d.S.S.d.C., E.C.M.d.R. and D.d.S.B.B.; formal analysis, R.A.d.C., A.S.P., A.d.S.S.d.C., E.C.M.d.R. and D.d.S.B.B.; investigation, R.A.d.C., J.A.P.d.R., A.d.S.G. and D.d.S.B.B.; resources, R.A.d.C. and D.d.S.B.B.; data curation, R.A.d.C., J.A.P.d.R., R.C.S., A.d.S.G. and D.d.S.B.B.; writing—original draft preparation, R.A.d.C., C.B.R.S. and D.d.S.B.B.; writing—review and editing, J.A.P.d.R., A.S.P. and A.d.S.G.; visualization, R.A.d.C. and C.B.R.S.; supervision, D.d.S.B.B.; project administration, R.A.d.C. and D.d.S.B.B.; funding acquisition, R.A.d.C. and D.d.S.B.B. All authors have read and agreed to the published version of the manuscript.
Not applicable.
The authors would like to acknowledge the PROPESP/UFPA (PAPQ) for the financial assistance.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 1. Predicted interactions in the docking with the NS2B/NS3pro enzyme of the ligands ZINC000001680989 (A), ZINC000001679427 (B), ZINC000006694490 (C) and ZINC000001684485 (D) in the GOLD program. The dashed black lines represent the hydrogen bonds and the green lines represent the hydrophobic interactions. The images were obtained through the PoseView Web 1.97 server.
Figure 2. . Graphical representation of the RMSD values versus the simulation time to the backbone Apo NS2B/NS3 (black line) ZINC000001679427 (red line), ZINC000001680989 (green line), ZINC000001684485 (blue line) and ZINC000006694490 (magenta line) studied.
Figure 3. Average RMSD values between protein/protein (orange) and ligand/ligand (cyan).
Figure 4. Graphical representation of the RMSF values versus the simulation time of the NS2B/NS3pro Apo backbone (black line) and its complex with the ZINC000001679427 (red line), ZINC000001680989 (green line), ZINC000001684485 (blue line) and ZINC000006694490 (magenta line) complexes studied.
Figure 5. Radius of gyration (Rg) of the NS2B/NS3pro ZINC000001679427 (black line), ZINC000001680989 (red line), ZINC000001684485 (green line) and ZINC000006694490 (blue line) complexes studied.
Figure 6. Solvent surface accessible area (SASA) of the NS2B/NS3pro ZINC000001679427 (black line), ZINC000001680989 (red line), ZINC000001684485 (green line) and ZINC000006694490 (blue line) complexes studied.
PLPChem score values obtained in the GOLD program and the structures of the four selected compounds.
Ligands | PLPchem | Structures |
---|---|---|
ZINC000001680989 | 62.05 | [Image omitted. Please see PDF.] |
ZINC000001679427 | 59.69 | [Image omitted. Please see PDF.] |
ZINC000006694490 | 58.41 | [Image omitted. Please see PDF.] |
ZINC000001684485 | 58.36 | [Image omitted. Please see PDF.] |
Estimated binding affinities and energetic contributions for the systems under study.
Ligand ID | Terms | |||||
---|---|---|---|---|---|---|
ΔEvdW | ΔEele | ΔGGB | ΔGnonpol | ΔGMMGBSA | ΔGMMPBSA | |
ZINC000001680989 | −44.83 ± 0.10 | −21.84 ± 0.11 | 31.28 ± 0.09 | −5.15 ± 0.01 | −40.65 ± 0.09 | −39.43 ± 0.12 |
ZINC000001679427 | −42.58 ± 0.10 | −21.53 ± 0.19 | 33.98 ± 0.16 | −5.40 ± 0.01 | −35.53 ± 0.11 | −25.72 ± 0.13 |
ZINC000006694490 | −38.22 ± 0.09 | −33.08 ± 0.25 | 48.50 ± 0.20 | −4.96 ± 0.01 | −27.76 ± 0.09 | −30.30 ± 0.12 |
ZINC000001684485 | −22.14 ± 0.07 | −31 11± 0.15 | 45.84 ± 0.35 | −4.31 ± 0.01 | −11.71 ± 0.14 | −19.67 ± 0.22 |
Hydrogen bonds between the compounds ZINC000001680989, ZINC000001679427, ZINC000006694490 and ZINC000001684485 with NS2B/NS3pro, which found at least 2.0% occupancy throughout 100 ns of simulation time.
Complex | Hydrogen Bond Formation | Distance (Å) | Occupancy (%) |
---|---|---|---|
NS2B/NS3pro-ZINC000001680989 | LEU179@O-LIG@H16-LIG@N3 | 2.82 | 86.39 |
ALA194@O-LIG@H15-LIG@N3 | 2.83 | 82.74 | |
LIG@N2-ASN182@HD21-ASN_182@ND2 | 2.91 | 20.09 | |
ALA194@O-LIG@H16-LIG@N3 | 2.82 | 16.96 | |
LEU179@O-LIG@H15-LIG@N3 | 2.82 | 16.85 | |
LIG@N4-LEU179@H-LEU179@N | 2.94 | 8.91 | |
LIG@N3-ASN182@HD21-ASN182@ND2 | 2.93 | 2.49 | |
NS2B/NS3pro-ZINC000001679427 | ALA194@O-LIG@H15-LIG@N3 | 2.83 | 46.96 |
LEU179@O-LIG@H14-LIG@N3 | 2.83 | 32.68 | |
ALA194@O-LIG@H14-LIG@N3 | 2.83 | 26.39 | |
LEU179@O-LIG@H15-LIG@N3 | 2.82 | 15.65 | |
LIG@N3-ASN182@HD21-ASN182@ND2 | 2.92 | 11.13 | |
LIG@N4-ASN182@HD21- ASN182@ND2 | 2.93 | 3.37 | |
LIG@N2-LEU179@H-LEU179@N | 2.94 | 2.58 | |
LYS104@O-LIG@H14-LIG@N3 | 2.82 | 2.37 | |
NS2B/NS3pro-ZINC000006694490 | ASN182@OD1- LIG@H23-LIG@N4 | 2.82 | 62.53 |
LYS104@O-LIG@H23-LIG@N4 | 2.83 | 8.58 | |
ILE195@O-LIG@H11-LIG@N5 | 2.88 | 6.83 | |
ILE_195@O-LIG@H-LIG@N | 2.88 | 5.58 | |
LIG@N1-ASN197@HD21-ASN197@ND2 | 2.91 | 4.47 | |
LIG@O-LYS104@HZ1-LYS104@NZ | 2.82 | 3.27 | |
LIG@OLYS104@HZ2-LYS104@NZ | 2.81 | 3.19 | |
LIG@O-LYS104@HZ3-LYS104@NZ | 2.81 | 2.86 | |
NS2B/NS3pro-ZINC000001684485 | LYS104@NZ-LIG@H-LIG@N | 2.80 | 95.48 |
LIG@N-LYS104@HZ1-LYS104@NZ | 2.80 | 55.70 | |
LIG@N-LYS104@HZ3-LYS104@NZ | 2.80 | 20.84 | |
LIG@N-LYS104@HZ2-LYS104@NZ | 2.79 | 18.93 | |
LIG@N4-ASN197@HD21-ASN197@ND2 | 2.91 | 2.04 |
Lipophilicity, druglikeness, solubility, physicochemical characteristics, pharmacokinetics and medicinal chemistry of the selected compounds evaluated by SwissADME.
Ligands | MR | TPSA (Å2) | Log Po/w | Log S | Log Kp (cm/s) | BS | SA |
---|---|---|---|---|---|---|---|
ZINC000001684485 | 100.5 | 107.23 | 2.33 | Moderately soluble | −6.91 | 0.55 | 3.69 |
ZINC000006694490 | 95.55 | 120.89 | 2.40 | Moderately soluble | −6.27 | 0.55 | 3.19 |
ZINC000001680989 | 101.23 | 94.92 | 3.50 | Poorly soluble | −5.18 | 0.55 | 3.06 |
ZINC000001679427 | 91.62 | 94.92 | 2.87 | Moderately soluble | −5.17 | 0.55 | 2.86 |
References
1. Clapham, H. Zika Virus Increases Risk of Dengue Disease. Science; 2020; 369, pp. 1055-1056. [DOI: https://dx.doi.org/10.1126/science.abd5922]
2. Tovar, J.H.; Pinzón, M.A.; Rincón, D.; Jiménez-Canizalez, C.E.; Arrieta-Mendoza, M.E.; Mondragón-Cardona, Á. Consideraciones Anestésicas En El Paciente Con Enfermedad Por Virus Dengue. Rev. Chil. Anest.; 2018; 47, pp. 20-26. [DOI: https://dx.doi.org/10.25237/revchilanestv47n01.04]
3. Vial, T.; Tan, W.-L.; Deharo, E.; Missé, D.; Marti, G.; Pompon, J. Mosquito Metabolomics Reveal That Dengue Virus Replication Requires Phospholipid Reconfiguration via the Remodeling Cycle. Proc. Natl. Acad. Sci. USA; 2020; 117, pp. 27627-27636. [DOI: https://dx.doi.org/10.1073/pnas.2015095117]
4. Kuhn, R.J.; Zhang, W.; Rossmann, M.G.; Pletnev, S.V.; Corver, J.; Lenches, E.; Jones, C.T.; Mukhopadhyay, S.; Chipman, P.R.; Strauss, E.G. et al. Structure of Dengue Virus. Cell; 2002; 108, pp. 717-725. [DOI: https://dx.doi.org/10.1016/S0092-8674(02)00660-8]
5. Yu, C.-Y.; Chang, T.-H.; Liang, J.-J.; Chiang, R.-L.; Lee, Y.-L.; Liao, C.-L.; Lin, Y.-L. Dengue Virus Targets the Adaptor Protein MITA to Subvert Host Innate Immunity. PLoS Pathog.; 2012; 8, e1002780. [DOI: https://dx.doi.org/10.1371/journal.ppat.1002780]
6. Green, A.M.; Beatty, P.R.; Hadjilaou, A.; Harris, E. Innate Immunity to Dengue Virus Infection and Subversion of Antiviral Responses. J. Mol. Biol.; 2014; 426, pp. 1148-1160. [DOI: https://dx.doi.org/10.1016/j.jmb.2013.11.023]
7. Su, Y.; Huang, Y.; Wu, Y.; Chen, H.; Wu, Y.; Hsu, C.; Hsu, Y.; Lee, J. MicroRNA-155 Inhibits Dengue Virus Replication by Inducing Heme Oxygenase-1-mediated Antiviral Interferon Responses. FASEB J.; 2020; 34, pp. 7283-7294. [DOI: https://dx.doi.org/10.1096/fj.201902878R] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32277848]
8. Brecher, M.; Zhang, J.; Li, H. The Flavivirus Protease as a Target for Drug Discovery. Virol. Sin.; 2013; 28, pp. 326-336. [DOI: https://dx.doi.org/10.1007/s12250-013-3390-x]
9. Luo, D.; Vasudevan, S.G.; Lescar, J. The Flavivirus NS2B–NS3 Protease–Helicase as a Target for Antiviral Drug Development. Antivir. Res.; 2015; 118, pp. 148-158. [DOI: https://dx.doi.org/10.1016/j.antiviral.2015.03.014]
10. Merdanovic, M.; Mönig, T.; Ehrmann, M.; Kaiser, M. Diversity of Allosteric Regulation in Proteases. ACS Chem. Biol.; 2013; 8, pp. 19-26. [DOI: https://dx.doi.org/10.1021/cb3005935]
11. Hauske, P.; Ottmann, C.; Meltzer, M.; Ehrmann, M.; Kaiser, M. Allosteric Regulation of Proteases. ChemBioChem; 2008; 9, pp. 2920-2928. [DOI: https://dx.doi.org/10.1002/cbic.200800528]
12. Mukhametov, A.; Newhouse, E.I.; Aziz, N.A.; Saito, J.A.; Alam, M. Allosteric Pocket of the Dengue Virus (Serotype 2) NS2B/NS3 Protease: In Silico Ligand Screening and Molecular Dynamics Studies of Inhibition. J. Mol. Graph. Model.; 2014; 52, pp. 103-113. [DOI: https://dx.doi.org/10.1016/j.jmgm.2014.06.008]
13. Millies, B.; von Hammerstein, F.; Gellert, A.; Hammerschmidt, S.; Barthels, F.; Göppel, U.; Immerheiser, M.; Elgner, F.; Jung, N.; Basic, M. et al. Proline-Based Allosteric Inhibitors of Zika and Dengue Virus NS2B/NS3 Proteases. J. Med. Chem.; 2019; 62, pp. 11359-11382. [DOI: https://dx.doi.org/10.1021/acs.jmedchem.9b01697]
14. Lim, L.; Dang, M.; Roy, A.; Kang, J.; Song, J. Curcumin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting Its Active Conformation. ACS Omega; 2020; 5, pp. 25677-25686. [DOI: https://dx.doi.org/10.1021/acsomega.0c00039] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33073093]
15. Yildiz, M.; Ghosh, S.; Bell, J.A.; Sherman, W.; Hardy, J.A. Allosteric Inhibition of the NS2B-NS3 Protease from Dengue Virus. ACS Chem. Biol.; 2013; 8, pp. 2744-2752. [DOI: https://dx.doi.org/10.1021/cb400612h]
16. Uday, R.V.S.; Misra, R.; Harika, A.; Dolui, S.; Saha, A.; Pal, U.; Ravichandiran, V.; Maiti, N.C. Dabrafenib, Idelalisib and Nintedanib Act as Significant Allosteric Modulator for Dengue NS3 Protease. PLoS ONE; 2021; 16, e0257206. [DOI: https://dx.doi.org/10.1371/journal.pone.0257206]
17. Brecher, M.; Li, Z.; Liu, B.; Zhang, J.; Koetzner, C.A.; Alifarag, A.; Jones, S.A.; Lin, Q.; Kramer, L.D.; Li, H. A Conformational Switch High-Throughput Screening Assay and Allosteric Inhibition of the Flavivirus NS2B-NS3 Protease. PLoS Pathog.; 2017; 13, e1006411. [DOI: https://dx.doi.org/10.1371/journal.ppat.1006411]
18. Hariono, M.; Choi, S.B.; Roslim, R.F.; Nawi, M.S.; Tan, M.L.; Kamarulzaman, E.E.; Mohamed, N.; Yusof, R.; Othman, S.; Abd Rahman, N. et al. Thioguanine-Based DENV-2 NS2B/NS3 Protease Inhibitors: Virtual Screening, Synthesis, Biological Evaluation and Molecular Modelling. PLoS ONE; 2019; 14, e0210869. [DOI: https://dx.doi.org/10.1371/journal.pone.0210869]
19. Erbel, P.; Schiering, N.; D’Arcy, A.; Renatus, M.; Kroemer, M.; Lim, S.P.; Yin, Z.; Keller, T.H.; Vasudevan, S.G.; Hommel, U. Structural Basis for the Activation of Flaviviral NS3 Proteases from Dengue and West Nile Virus. Nat. Struct. Mol. Biol.; 2006; 13, pp. 372-373. [DOI: https://dx.doi.org/10.1038/nsmb1073]
20. da Costa, R.A.; da Rocha, J.A.P.; Pinheiro, A.S.; da Costa, A.S.S.; da Rocha, E.C.M.; Josino, L.P.C.; da Gonçalves, A.S.; Lima, A.H.L.; Brasil, D.S.B. In Silico Identification of Novel Allosteric Inhibitors of Dengue Virus NS2B/NS3 Serine Protease. J. Serb. Chem. Soc.; 2022; 11. [DOI: https://dx.doi.org/10.2298/JSC210929011D]
21. Opo, F.A.D.M.; Rahman, M.M.; Ahammad, F.; Ahmed, I.; Bhuiyan, M.A.; Asiri, A.M. Structure Based Pharmacophore Modeling, Virtual Screening, Molecular Docking and ADMET Approaches for Identification of Natural Anti-Cancer Agents Targeting XIAP Protein. Sci. Rep.; 2021; 11, 4049. [DOI: https://dx.doi.org/10.1038/s41598-021-83626-x]
22. Kato, K.; Nakayoshi, T.; Fukuyoshi, S.; Kurimoto, E.; Oda, A. Validation of Molecular Dynamics Simulations for Prediction of Three-Dimensional Structures of Small Proteins. Molecules; 2017; 22, 1716. [DOI: https://dx.doi.org/10.3390/molecules22101716]
23. Pal, U.; Sen, S.; Maiti, N.C. Cα–H Carries Information of a Hydrogen Bond Involving the Geminal Hydroxyl Group: A Case Study with a Hydrogen-Bonded Complex of 1,1,1,3,3,3-Hexafluoro-2-Propanol and Tertiary Amines. J. Phys. Chem. A; 2014; 118, pp. 1024-1030. [DOI: https://dx.doi.org/10.1021/jp411488a]
24. Lobanov, M.Y.; Bogatyreva, N.S.; Galzitskaya, O.V. Radius of Gyration as an Indicator of Protein Structure Compactness. Mol. Biol.; 2008; 42, pp. 623-628. [DOI: https://dx.doi.org/10.1134/S0026893308040195]
25. Shahroz, M.M.; Sharma, H.K.; Altamimi, A.S.A.; Alamri, M.A.; Ali, A.; Ali, A.; Alqahtani, S.; Altharawi, A.; Alabbas, A.B.; Alossaimi, M.A. et al. Novel and Potential Small Molecule Scaffolds as DYRK1A Inhibitors by Integrated Molecular Docking-Based Virtual Screening and Dynamics Simulation Study. Molecules; 2022; 27, 1159. [DOI: https://dx.doi.org/10.3390/molecules27041159]
26. Sun, H.; Duan, L.; Chen, F.; Liu, H.; Wang, Z.; Pan, P.; Zhu, F.; Zhang, J.Z.H.; Hou, T. Assessing the Performance of MM/PBSA and MM/GBSA Methods. 7. Entropy Effects on the Performance of End-Point Binding Free Energy Calculation Approaches. Phys. Chem. Chem. Phys.; 2018; 20, pp. 14450-14460. [DOI: https://dx.doi.org/10.1039/C7CP07623A]
27. Wu, H.; Bock, S.; Snitko, M.; Berger, T.; Weidner, T.; Holloway, S.; Kanitz, M.; Diederich, W.E.; Steuber, H.; Walter, C. et al. Novel Dengue Virus NS2B/NS3 Protease Inhibitors. Antimicrob. Agents Chemother.; 2015; 59, pp. 1100-1109. [DOI: https://dx.doi.org/10.1128/AAC.03543-14]
28. Othman, R.; Kiat, T.S.; Khalid, N.; Yusof, R.; Irene Newhouse, E.; Newhouse, J.S.; Alam, M.; Rahman, N.A. Docking of Noncompetitive Inhibitors into Dengue Virus Type 2 Protease: Understanding the Interactions with Allosteric Binding Sites. J. Chem. Inf. Model.; 2008; 48, pp. 1582-1591. [DOI: https://dx.doi.org/10.1021/ci700388k]
29. Mannhold, R.; Poda, G.I.; Ostermann, C.; Tetko, I.V. Calculation of Molecular Lipophilicity: State-of-the-Art and Comparison of LogP Methods on More than 96,000 Compounds. J. Pharm. Sci.; 2009; 98, pp. 861-893. [DOI: https://dx.doi.org/10.1002/jps.21494]
30. Arnott, J.A.; Planey, S.L. The Influence of Lipophilicity in Drug Discovery and Design. Expert Opin. Drug Dis.; 2012; 7, pp. 863-875. [DOI: https://dx.doi.org/10.1517/17460441.2012.714363]
31. Ji, D.; Xu, M.; Udenigwe, C.C.; Agyei, D. Physicochemical Characterisation, Molecular Docking, and Drug-Likeness Evaluation of Hypotensive Peptides Encrypted in Flaxseed Proteome. Curr. Res. Food Sci.; 2020; 3, pp. 41-50. [DOI: https://dx.doi.org/10.1016/j.crfs.2020.03.001]
32. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep.; 2017; 7, 42717. [DOI: https://dx.doi.org/10.1038/srep42717]
33. Silva, R.C.; Freitas, H.F.; Campos, J.M.; Kimani, N.M.; Silva, C.H.T.P.; Borges, R.S.; Pita, S.S.R.; Santos, C.B.R. Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro. Int. J. Mol. Sci.; 2021; 22, 11739. [DOI: https://dx.doi.org/10.3390/ijms222111739]
34. Silva, R.; Poiani, J.; Ramos, R.; Costa, J.; Silva, C.; Brasil, D.; Santos, C. Ligand and Structure-Based Virtual Screening from 16-((Diisobutylamino)Methyl)-6α-Hydroxyivouacapane-7β,17β-Lactone a Compound with Potential Anti-Prostate Cancer Activity. J. Serb. Chem. Soc.; 2019; 84, pp. 153-174. [DOI: https://dx.doi.org/10.2298/JSC180129047S]
35. Kontorinis, N.; Agarwal, K.; Gondolesi, G.; Fiel, M.I.; O’Rourke, M.; Schiano, T.D. Diagnosis of 6 Mercaptopurine Hepatotoxicity Post Liver Transplantation Utilizing Metabolite Assays. Am. J. Transp.; 2004; 4, pp. 1539-1542. [DOI: https://dx.doi.org/10.1111/j.1600-6143.2004.00543.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15307844]
36. Clarke, D.A.; Philips, F.S.; Sternberg, S.S.; Stock, C.C.; Elion, G.B.; Hitchings, G.H. 6-Mercaptopurine: Effects in Mouse Sarcoma 180 and in Normal Animals. Cancer Res.; 1953; 13, pp. 593-604. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/13082494]
37. Gisbert, J.P.; González-Lama, Y.; Maté, J. Thiopurine-Induced Liver Injury in Patients with Inflammatory Bowel Disease: A Systematic Review. Am. J. Gastroenterol.; 2007; 102, pp. 1518-1527. [DOI: https://dx.doi.org/10.1111/j.1572-0241.2007.01187.x]
38. Bissuel, F.; Bruneel, F.; Habersetzer, F.; Chassard, D.; Cotte, L.; Chevallier, M.; Bernuau, J.; Lucet, J.-C.; Trepo, C. Fulminant Hepatitis with Severe Lactate Acidosis in HIV-Infected Patients on Didanosine Therapy. J. Int. Med.; 1994; 235, pp. 367-372. [DOI: https://dx.doi.org/10.1111/j.1365-2796.1994.tb01088.x]
39. Maida, I.; Núñez, M.; Ríos, M.J.; Martín-Carbonero, L.; Sotgiu, G.; Toro, C.; Rivas, P.; Barreiro, P.; Mura, M.S.; Babudieri, S. et al. Severe Liver Disease Associated with Prolonged Exposure to Antiretroviral Drugs. JAIDS J. Acquir. Immune Defic. Syndr.; 2006; 42, pp. 177-182. [DOI: https://dx.doi.org/10.1097/01.qai.0000221683.44940.62]
40. Hu, B.; French, S.W. 2′,3′-Dideoxyinosine–Induced Mallory Bodies in Patients With HIV. Am. J. Clin. Pathol.; 1997; 108, pp. 280-283. [DOI: https://dx.doi.org/10.1093/ajcp/108.3.280]
41. Lacaille, F.; Ortigao, M.B.; Debré, M.; Rouzioux, C.; Brousse, N.; Blanche, S. Hepatic Toxicity Associated with 2′-3ʼDideoxyinosine in Children with AIDS. J. Pediatr. Gastroenterol. Nutr.; 1995; 20, pp. 287-290. [DOI: https://dx.doi.org/10.1097/00005176-199504000-00004]
42. Bastida, G.; Nos, P.; Aguas, M.; Beltran, B.; Rubin, A.; Dasi, F.; Ponce, J. Incidence, Risk Factors and Clinical Course of Thiopurine-Induced Liver Injury in Patients with Inflammatory Bowel Disease. Aliment. Pharmacol. Ther.; 2005; 22, pp. 775-782. [DOI: https://dx.doi.org/10.1111/j.1365-2036.2005.02636.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16225485]
43. Davis, M.; Eddleston, A.L.; Williams, R. Hypersensitivity and Jaundice Due to Azathioprine. Postgrad. Med. J.; 1980; 56, pp. 274-275. [DOI: https://dx.doi.org/10.1136/pgmj.56.654.274]
44. Jeurissen, M.E.; Boerbooms, A.M.; van de Putte, L.B.; Kruijsen, M.W. Azathioprine Induced Fever, Chills, Rash, and Hepatotoxicity in Rheumatoid Arthritis. Ann. Rheum. Dis.; 1990; 49, pp. 25-27. [DOI: https://dx.doi.org/10.1136/ard.49.1.25] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/2138007]
45. Haynes, P.; Lambert, T.R.; de G. Mitchell, I. Comparative In-Vivo Genotoxicity of Antiviral Nucleoside Analogues; Penciclovir, Acyclovir, Ganciclovir and the Xanthine Analogue, Caffeine, in the Mouse Bone Marrow Micronucleus Assay. Mutat. Res. Genet. Toxicol.; 1996; 369, pp. 65-74. [DOI: https://dx.doi.org/10.1016/S0165-1218(96)90049-X]
46. Hara, T.; Makita, T.; Horiya, N.; Ozawa, S.; Ohba, M.; Naito, J.; Shibuya, T. Micronucleus Test with 6-Mercaptopurine Monohydrate Administered Intraperitoneally and Orally. Mutat. Res. Genet. Toxicol.; 1989; 223, pp. 349-352. [DOI: https://dx.doi.org/10.1016/0165-1218(89)90083-9]
47. Holden, H.E.; Ray, V.A.; Wahrenburg, M.G.; Zelenski, J.D. Mutagenicity Studies with 6-Mercaptopurine: I. Cytogenetic Activity in Vivo. Mutat. Res. Fundam. Mol. Mech. Mutagen.; 1973; 20, pp. 257-263. [DOI: https://dx.doi.org/10.1016/0027-5107(73)90195-4]
48. Phillips, M.D.; Nascimbeni, B.; Tice, R.R.; Shelby, M.D.; van Zeeland, A.A. Induction of Micronuclei in Mouse Bone Marrow Cells: An Evaluation of Nucleoside Analogues Used in the Treatment of AIDS. Environ. Mol. Mutagen.; 1991; 18, pp. 168-183. [DOI: https://dx.doi.org/10.1002/em.2850180305] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/1915312]
49. Yarchoan, R.; Mitsuya, H.; Thomas, R.V.; Pluda, J.M.; Hartman, N.R.; Perno, C.-F.; Marczyk, K.S.; Allain, J.-P.; Johns, D.G.; Broder, S. In Vivo Activity Against HIV and Favorable Toxicity Profile of 2′,3′-Dideoxyinosine. Science; 1989; 245, pp. 412-415. [DOI: https://dx.doi.org/10.1126/science.2502840]
50. Wutzler, P.; Thust, R. Genetic Risks of Antiviral Nucleoside Analogues–a Survey. Antivir. Res.; 2001; 49, pp. 55-74. [DOI: https://dx.doi.org/10.1016/S0166-3542(00)00139-X]
51. Zoete, V.; Daina, A.; Bovigny, C.; Michielin, O. SwissSimilarity: A Web Tool for Low to Ultra High Throughput Ligand-Based Virtual Screening. J. Chem. Inf. Model.; 2016; 56, pp. 1399-1404. [DOI: https://dx.doi.org/10.1021/acs.jcim.6b00174]
52. Majumder, R.; Mandal, M. Screening of Plant-Based Natural Compounds as a Potential COVID-19 Main Protease Inhibitor: An in Silico Docking and Molecular Dynamics Simulation Approach. J. Biomol. Struct. Dyn.; 2020; 40, pp. 696-711. [DOI: https://dx.doi.org/10.1080/07391102.2020.1817787]
53. Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and Validation of a Genetic Algorithm for Flexible Docking 1 1Edited by F. E. Cohen. J. Mol. Biol.; 1997; 267, pp. 727-748. [DOI: https://dx.doi.org/10.1006/jmbi.1996.0897] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9126849]
54. Cole, J.; Willem, M.; Nissink, J.; Taylor, R. Protein-Ligand Docking Virtual Screening with GOLD. Virtual Screening in Drug Discovery; Taylor & Francis CRC Press: Boca Raton, FL, USA, 2005; pp. 379-415.
55. Werner, T.; Morris, M.B.; Dastmalchi, S.; Church, W.B. Structural Modelling and Dynamics of Proteins for Insights into Drug Interactions. Adv. Drug Deliv. Rev.; 2012; 64, pp. 323-343. [DOI: https://dx.doi.org/10.1016/j.addr.2011.11.011]
56. Hollingsworth, S.A.; Dror, R.O. Molecular Dynamics Simulation for All. Neuron; 2018; 99, pp. 1129-1143. [DOI: https://dx.doi.org/10.1016/j.neuron.2018.08.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30236283]
57. Cerutti, D.S.; Mermelstein, D.; Lin, C.; LeGrand, S.; Giese, T.J.; Roitberg, A.; Case, D.A.; Walker, R.C.; York, D.M. GPU-Accelerated Molecular Dynamics and Free Energy Methods in Amber18: Performance Enhancements and New Features. J. Chem. Inf. Model.; 2018; 58, pp. 2043-2050. [DOI: https://dx.doi.org/10.1021/acs.jcim.8b00462]
58. Song, L.F.; Lee, T.-S.; Zhu, C.; York, D.M.; Merz, K.M. Using AMBER18 for Relative Free Energy Calculations. J. Chem. Inf. Model.; 2019; 59, pp. 3128-3135. [DOI: https://dx.doi.org/10.1021/acs.jcim.9b00105] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31244091]
59. Wang, J.M.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General Amber Force Field. J. Comput. Chem.; 2004; 25, pp. 1157-1174. [DOI: https://dx.doi.org/10.1002/jcc.20035]
60. Frisch, M.; Trucks, G.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. Gaussian 09; Revision d. 01 Gaussian, Inc.: Wallingford, CT, USA, 2009; 201.
61. Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N·log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys.; 1993; 98, pp. 10089-10092. [DOI: https://dx.doi.org/10.1063/1.464397]
62. Ryckaert, J.P.; Ciccotti, G.; Berendsen, H.J.C. Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes. J. Comput. Phys.; 1977; 23, pp. 327-341. [DOI: https://dx.doi.org/10.1016/0021-9991(77)90098-5]
63. Costa, R.A.; Cruz, J.N.; Nascimento, F.C.A.; Silva, S.G.; Silva, S.O.; Martelli, M.C.; Carvalho, S.M.L.; Santos, C.B.R.; Neto, A.M.J.C.; Brasil, D.S.B. Studies of NMR, Molecular Docking, and Molecular Dynamics Simulation of New Promising Inhibitors of Cruzaine from the Parasite Trypanosoma Cruzi. Med. Chem. Res.; 2019; 28, pp. 246-259. [DOI: https://dx.doi.org/10.1007/s00044-018-2280-z]
64. Gohlke, H.; Kiel, C.; Case, D.A. Insights into Protein–Protein Binding by Binding Free Energy Calculation and Free Energy Decomposition for the Ras–Raf and Ras–RalGDS Complexes. J. Mol. Biol.; 2003; 330, pp. 891-913. [DOI: https://dx.doi.org/10.1016/S0022-2836(03)00610-7]
65. Szakács, G.; Váradi, A.; Özvegy-Laczka, C.; Sarkadi, B. The Role of ABC Transporters in Drug Absorption, Distribution, Metabolism, Excretion and Toxicity (ADME–Tox). Drug Dis. Today; 2008; 13, pp. 379-393. [DOI: https://dx.doi.org/10.1016/j.drudis.2007.12.010]
66. Costa, E.B.; Silva, R.C.; Espejo-Román, J.M.; Neto, M.F.d.A.; Cruz, J.N.; Leite, F.H.A.; Silva, C.H.T.P.; Pinheiro, J.C.; Macêdo, W.J.C.; Santos, C.B.R. Chemometric Methods in Antimalarial Drug Design from 1,2,4,5-Tetraoxanes Analogues. SAR QSAR Environ. Res.; 2020; 31, pp. 677-695. [DOI: https://dx.doi.org/10.1080/1062936X.2020.1803961] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32854545]
67. Kirchmair, J.; Göller, A.H.; Lang, D.; Kunze, J.; Testa, B.; Wilson, I.D.; Glen, R.C.; Schneider, G. Predicting Drug Metabolism: Experiment and/or Computation?. Nat. Rev. Drug Dis.; 2015; 14, pp. 387-404. [DOI: https://dx.doi.org/10.1038/nrd4581] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25907346]
68. Greene, N.; Judson, P.N.; Langowski, J.J.; Marchant, C.A. Knowledge-Based Expert Systems for Toxicity and Metabolism Prediction: DEREK, StAR and METEOR. SAR QSAR Environ. Res.; 1999; 10, pp. 299-314. [DOI: https://dx.doi.org/10.1080/10629369908039182]
69. Braggio, S.; Corsi, M.; Feriani, A.; Fontana, S.; Marocchio, L.; Virginio, C. CHAPTER 15. Discovery Toxicology in Lead Optimisation. The Handbook of Medicinal Chemistry; Royal Society of Chemistry: Cambridge, UK, 2015; pp. 364-412.
70. Sutter, A.; Amberg, A.; Boyer, S.; Brigo, A.; Contrera, J.F.; Custer, L.L.; Dobo, K.L.; Gervais, V.; Glowienke, S.; van Gompel, J. et al. Use of in Silico Systems and Expert Knowledge for Structure-Based Assessment of Potentially Mutagenic Impurities. Regul. Toxicol. Pharmacol.; 2013; 67, pp. 39-52. [DOI: https://dx.doi.org/10.1016/j.yrtph.2013.05.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23669331]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Dengue virus (DENV) is a danger to more than 400 million people in the world, and there is no specific treatment. Thus, there is an urgent need to develop an effective method to combat this pathology. NS2B/NS3 protease is an important biological target due it being necessary for viral replication and the fact that it promotes the spread of the infection. Thus, this study aimed to design DENV NS2B/NS3pro allosteric inhibitors from a matrix compound. The search was conducted using the Swiss Similarity tool. The compounds were subjected to molecular docking calculations, molecular dynamics simulations (MD) and free energy calculations. The molecular docking results showed that two compounds, ZINC000001680989 and ZINC000001679427, were promising and performed important hydrogen interactions with the Asn152, Leu149 and Ala164 residues, showing the same interactions obtained in the literature. In the MD, the results indicated that five residues, Lys74, Leu76, Asn152, Leu149 and Ala166, contribute to the stability of the ligand at the allosteric site for all of the simulated systems. Hydrophobic, electrostatic and van der Waals interactions had significant effects on binding affinity. Physicochemical properties, lipophilicity, water solubility, pharmacokinetics, druglikeness and medicinal chemistry were evaluated for four compounds that were more promising, showed negative indices for the potential penetration of the Blood Brain Barrier and expressed high human intestinal absorption, indicating a low risk of central nervous system depression or drowsiness as the the side effects. The compound ZINC000006694490 exhibited an alert with a plausible level of toxicity for the purine base chemical moiety, indicating hepatotoxicity and chromosome damage in vivo in mouse, rat and human organisms. All of the compounds selected in this study showed a synthetic accessibility (SA) score lower than 4, suggesting the ease of new syntheses. The results corroborate with other studies in the literature, and the computational approach used here can contribute to the discovery of new and potent anti-dengue agents.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details



1 Graduate Program in Science and Environment, Institute of Exact and Natural Sciences, Federal University of Pará (UFPA), Belém 66075-110, PA, Brazil;
2 Federal Institute of Education, Science and Technology of Pará—Campus Bragança, Bragança 68600-000, PA, Brazil;
3 Graduate Program in Science and Environment, Institute of Exact and Natural Sciences, Federal University of Pará (UFPA), Belém 66075-110, PA, Brazil;
4 Departamento de Química, Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto 14040-901, SP, Brazil;
5 Federal Institute of Education, Science and Technology of Espírito Santo, Vila Velha 29106-010, ES, Brazil;
6 Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapá 68902-280, AP, Brazil;