key: cord-0027167-krjink6a
authors: Shimu, Mst. Sharmin Sultana; Mahmud, Shafi; Tallei, Trina Ekwati; Sami, Saad Ahmed; Adam, Ahmad Akroman; Acharjee, Uzzal Kumar; Paul, Gobindo Kumar; Emran, Talha Bin; Zaman, Shahriar; Uddin, Md. Salah; Saleh, Md. Abu; Alshehri, Sultan; Ghoneim, Mohammed M; Alruwali, Maha; Obaidullah, Ahmad J.; Jui, Nabilah Rahman; Kim, Junghwan; Kim, Bonglee
title: Phytochemical Compound Screening to Identify Novel Small Molecules against Dengue Virus: A Docking and Dynamics Study
date: 2022-01-20
journal: Molecules
DOI: 10.3390/molecules27030653
sha: 378c3e5783ccffc6b9541ce23b7c9fad6135988c
doc_id: 27167
cord_uid: krjink6a

The spread of the Dengue virus over the world, as well as multiple outbreaks of different serotypes, has resulted in a large number of deaths and a medical emergency, as no viable medications to treat Dengue virus patients have yet been found. In this paper, we provide an in silico virtual screening and molecular dynamics-based analysis to uncover efficient Dengue infection inhibitors. Based on a Google search and literature mining, a large phytochemical library was generated and employed as ligand molecules. In this investigation, the protein target NS2B/NS3 from Dengue was employed, and around 27 compounds were evaluated in a docking study. Phellodendroside (−63 kcal/mole), quercimeritrin (−59.5 kcal/mole), and quercetin-7-O-rutinoside (−54.1 kcal/mole) were chosen based on their binding free energy in MM-GBSA. The tested compounds generated numerous interactions at Lys74, Asn152, and Gln167 residues in the active regions of NS2B/NS3, which is needed for the protein’s inhibition. As a result, the stable mode of docked complexes is defined by various descriptors from molecular dynamics simulations, such as RMSD, SASA, Rg, RMSF, and hydrogen bond. The pharmacological properties of the compounds were also investigated, and no toxicity was found in computational ADMET properties calculations. As a result, this computational analysis may aid fellow researchers in developing innovative Dengue virus inhibitors.

The Dengue virus is spread by the mosquito Aedes aegypti, which is mostly found in tropical and subtropical areas [1] . Currently, 2.5 billion people around the world are at danger of contracting Dengue fever, with hundreds of millions of people suffering from the disease [2] . Climate change, population increase, and deforestation all contribute to the spread of Dengue fever, resulting in recurring outbreaks and more intense transmission of diverse serotypes [3] . The Dengue virus has four different serotypes, and infection with one of them does not confer protection against the other serotypes [4] . As a result, Dengue infections caused by two different serotypes carry the risk of Dengue hemorrhagic fever [5] . Despite the fact that many studies are still ongoing and that a large number of people are affected by Dengue fever on a daily basis, no effective inhibitors have been identified. As a result, a novel technique for the identification of a powerful Dengue inhibitor is required [6, 7] .

As a result, the Dengue virus is classified as an RNA virus with a single-stranded positive-sense RNA genome of 10.7 kb [8] [9] [10] . The Dengue virus RNA genome is transcribed as a single polyprotein, and subsequently co-and post-translationally cleaved into three structural proteins (Capsid (C), pre-membrane (prM), envelope (E) and seven non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) [11] [12] [13] . These structural proteins, on the other hand, mediate the proteolytic processing that is mediated by the endoplasmic reticulum-based host cell signal peptidase [14] . Cleavage occurs in the NS2A/NS2B, NS2B/NS3, NS3/NS4A, and NS4B/NS5 sites by the NS2B/NS3 protease. As a result, this viral protease from Dengue also cleaves the internal side of C, NS2A, NS3, NS4A, and C/prM. As a result, the NS2B/NS3 function as a key viral protease, which is significant for the development of antiviral drugs against the Dengue virus that can interfere with the virus's functions [15, 16] .

Multiple efforts are currently underway to find efficient Dengue virus inhibitors and vaccines by targeting various structural and non-structural proteins. As a result, Dengvaxia, a Dengue vaccine, has been studied in two independent phase 3 clinical trials in Latin America and Asia, with efficacy ranging from 31.3 percent to 79 percent, demonstrating diversity in vaccine efficacy across different geographies and serotypes [17, 18] . Despite various efforts by experts, there are currently no antiviral medications available on the market for the Dengue virus. Furthermore, Zingiberaceae chalcone compounds have been shown to inhibit the DENV2 protease. Various isolated substances, such as panduratin A and 4-hydroxypanduratin, also have inhibitory activity [19] . Quercetin and agathisflavone [20] , dehydronaphthalene [21] , benzimidazole [22] , thiadiazoloacrylamide [23] , and drug designing process [24, 25] are some of the powerful compounds that target the NSB2/NS3 from phytochemicals.

Curdlan sulfate [26] , rolitetracycline [26] , fucoidan [27] , balapiravir [28] , suramin [29] , carnosine [30] , policresulen [31] , dasatinib [32] , duramycin [33] , luteolin [34, 35] , and imatinib [36] are examples of DAA or direct acting antivirals that efficiently interact by targeting NS and structural proteins. DAA, on the other hand, has the ability to target viral proteins and usually has a lower level of toxicity, but it comes with the acknowledged drawback of a higher likelihood of resistance development [37] . Furthermore, numerous DAA candidates for Dengue failed in mice experiments despite showing promise in in vitro assays, and only balapiravir has been studied in human trials among many DAA candidates [38] . Furthermore, due to the lower chance of acquiring resistance, host directed antivirals, or HDAs, have a lot of potential in a broader range, but only a few candidates have made it to the clinical trial stage [39] . Cleavages of florescence energy transfer (FRET), HPLC-based detection, isothermal titration calorimetry, and surface plasmon resonance have also been used to identify various protease inhibitors [37, [40] [41] [42] .

We conducted a virtual screening working flow of phytochemicals against the viral protein target of NS2B/NS3 in this investigation. The binding affinity and MM-GBSA techniques aid in the screening of phytochemical libraries and the selection of molecules with the highest binding energy. A molecular dynamics simulation was also carried out to further evaluate the docked complexes' binding conformations and stability.

For improved energy against viral NS2B/NS3 protease (PDB ID: 2FOM), a total of 27 compounds were chosen, with docking scores ranging from −6.5 kcal/mol to −3.4 kcal/mol. The binding affinity of choerospondin diaplayed (−6.576 kcal/mol) to the target receptor was the highest whereas the control exhibit −5. Therefore, the XP score was found higher for cherospondin (−5.1 kcal/mol), angustidine (−6 kcal/mol), forsythoside A (−5.8 kcal/mol), luteon (−5.9 kcal/mol), quercimeritrin (−6 kcal/mol), hydroxytanshinone (−5.4 kcal/mol), phellodendroside (−6.1 kcal/mol), quercetin-7-O-rutinoside (−5.9 kcal/mol), phellamuretin (−5.8 kcal/mol), and control (−5.4 kcal/mol), respectively. The majority of the compounds interacted with our target receptor's Lys74, Leu76, Trp83, Leu85, Asn152, and Ile165 active site residues via hydrogen and hydrophobic contacts. Table 1 and Figure S1 show the docking score and amino acid interactions between the selected drugs and the Dengue virus NS2B/NS3 protease. 

Lys74, Leu76, Ala166 Ile165 -

Leu128 --

We used Prime MM-GBSA attributes to compute the binding energy of the ligands with the Dengue virus NS2B/NS3 protease to assess the clarity of the molecular docking simulation. The binding free energy of a ligand molecule with a lower negative value suggests that it has better binding abilities. The results showed that most of the ligands' binding energies imply a strong binding engagement with the target receptor. Phellodendroside (−63 kcal/mole), quercimeritrin (−59.5 kcal/mole), and quercetin-7-O-rutinoside (−54.1 kcal/mole) are the top ligand-receptor complexes, and the control exhibits a higher binding energy of −51.3 Kcal/mole. Table 1 and Figures 1-4 show the results of the binding free energies, three compounds (Figures 1-3) , and a control ( Figure 4 ) with improved amino acid interactions at the active site of the target protein, respectively. 

Pharmacological and toxicity evaluation are critical in assuring the efficacy and safety of lead compounds. The screened compounds were evaluated using the parameters molecular weight, h bond acceptor, h bond donor, CYP2D6 substrate, CYP3A4 substrate, CYP1A2 inhibitor, CYP2C19 inhibitor, CYP2C9 inhibitor, CYP2D6 inhibitor, CYP3A4 inhibitor, AMES toxicity, and oral acute toxicity (Table 2) . Phellodendroside, quercimeritrin, and quercetin-7-O-rutinoside had molecular weights of 518.51, 464.38, and 610.52g/mol, respectively. There is no CYP inhibition in any of the drugs. In AMES toxicity and oral rat acute toxicity profile, the chemical displays no toxicity. 

To further understand the stability and stiffness of the docked complexes, molecular dynamics simulations were used. The simulated trajectories were used to calculate the root mean square deviations from the C-alpha atoms of the complexes. Figure 5A demonstrates that at the start of the simulations, phellodendroside, quercetin-7-O-rutinoside, and quercimeritin have an upper trend in RMSD, indicating that the complexes were still flexible. The complexes, on the other hand, had established a stable state after 20ns and had maintained their integrity until the simulations' final images. The quercimeritin displayed a greater RMSD trend and more fluctuations at 40-55 ns periods, indicating the complex's more flexible character; however, the complexes eventually settled into a stable state. RMSD was less than 2.5 Å for all complexes until the end of the simulations, indicating that the complexes were overall stable.

Furthermore, the solvent accessible surface area (SASA) of the simulated complexes was examined, with a larger trend in SASA indicating the expansion of the protein's surface area and a lower SASA indicating the complex's truncated nature. Figure 5B shows that phellodendroside had a larger SASA trend and a higher degree of departures from 60-70 and 80-100 ns compared to other complexes, indicating that when ligand binds to the protein, the surface area expands and the stability decreases. The radius of gyrations from the simulated trajectories was also determined, with the larger Rg indicating a more mobile character of the complexes and the lower Rg indicating a stiffer structure. From 30 ns to the rest of the simulation durations, the complexes of phellodendroside had a greater Rg than the other complexes, indicating that the complexes are more labile ( Figure 5C ). In simulation trajectories, the other complexes show smaller degrees of divergence. In determining the stable state of the complexes, the hydrogen bond plays a crucial role. Figure 5D shows that all complexes showed lower degrees of hydrogen bond patterning aberrations. The root mean square fluctuations (RMSF) of the complexes, which determine the flexibility across the amino acid residues, were also calculated. Figure 6 shows that the maximum residues had an RMSF of less than 2.5, indicating that the complexes were overall stable. 

Molecular docking is one of the most widely used methods in computational drug design for predicting the orientation of tiny molecules attached to the binding pocket of an enzyme or receptor [43] . The binding affinities of chosen drugs against Dengue virus NS2B/NS3 protease were investigated using molecular docking simulation in this work. In conjugation with the NS2B cofactor, NS3, a trypsin-like viral protease, is a crucial regulator of viral replication [44] . Potential protease inhibitors derived from Carica papaya leaves interacted with the NS3 protease active site's Leu149 and conserved substrate-binding residue Asn152. Bioflavonoids extracted from Azadirachta indica, such as rutin, hyperoside, kaempferol-3-O-rutinoside, and epicatechin, interacted with the Asn152 residue of the NS3 protease active site [45] . Lys74, Leu76, Trp83, Leu85, Glu88, Ile165, Ala166, and Asn167 active site residues of the dengue NS2B/NS3 protease were found to be involved in hydrophobic and hydrogen bond formation in a recent study based on antiviral phytochemicals [46] . Similar interactions were also found for a few amino acid residues in the dengue NS2B/NS3 protease active site (Lys74, Asn152, and Gln167) in several experimentally tested synthetic compounds [47] . Three compounds were chosen for this study: quercimeritrin, phellodendroside, and quercetin-7-O-rutinoside, based on their amino acid interactions, docking scores, and binding free energies. Quercimeritrin interacted with the NS3 protease active site residues Lys74, Leu76, Trp83, Leu85, Glu88, Asn152, Asn167, whereas phellodendroside interacted with Trp69, Lys74, Trp83, Trp89, Val146, Val147, Gly148, Ile165. Quercetin-7-O-rutinoside established hydrogen bonds with the NS2B/NS3 protease residues Asp71, Lys74, Trp83, Gly87, Val147, and Ile165. Multiple hydrogen bonds were formed between the control ligand molecules and the target protein at Val72, His51, Arg54, and Phe130, as well as one hydrophobic interaction at Leu128.

Multiple interactions were also observed at Asn152 and Lys74 by benzothiazole derivatives in molecular modeling studies, and their inhibition of replications was validated in vitro and in a cell-based PR assay for distinct Dengue serotypes [48] . Benzimide also generated hydrogen-bond interactions at Asn152, and the hydrophobic side of the benziminazole ring immediately attached to engage with Lys74 and Leu76, which are identical to the interactions of the top three compounds in this study [49] . In recent fragment-based drug design efforts against Dengue virus NS2B/NS3 revealed that Asn167, Leu85, and Glu88 were important in generating non-covalent interactions with the NS2B/NS3 protein [50] . Our screened compounds' interaction patterns were also shown to be similar for papraine A, robustune, eryvarin, sigmodin, osaji, and laburnetin [51] .

Based on the findings of the previous investigations, these amino acid interactions showed that these compounds could be effective as dengue NS2B/NS3 protease inhibitors. Although molecular docking is a reliable method for predicting a small molecule's binding poses, it still requires additional validation to appropriately rank the ligands' affinities to the target receptor. We used Prime MM/GBSA analysis to verify the docking experiments, which is a revolutionary QM/MM property for calculating relative binding affinities using the best poses from ligand-receptor interactions [52] . Within the binding site of the Dengue virus NS2B/NS3 protease, phellodendroside has the highest MM-GBSA binding energy of approximately −63 kcal/mole. The binding free energies of quercimeritrin and quercetin-7-O-rutinoside were also acceptable, confirming the results of our docking experiment.

In addition, a molecular dynamics simulation investigation was carried out for further docking conformation validations [53] [54] [55] . For a better understanding of the docked complexes' stable nature, numerous characteristics from simulated trajectories were evaluated. The RMSD, RMSF, Rg, hydrogen bond, and SASA define the complexes and show that they are more stable. As a result, the complexes' pharmacological properties were investigated to assure their safety and toxicity level [56] [57] [58] . To be used as drug candidates, substances must have particular drug-like qualities and have a low level of toxicity [59] [60] [61] . In ADMET calculations, the top three compounds show no toxicity profiling and have favorable drug similarity features.

In addition, the MM-GBSA binding free energy following molecular dynamics simulations shows that the complex; quercimeritrin, phellodendroside, and quercetin-7-O-rutinoside exhibited binding free energies of −51.3Kcal/mol, −54.3Kcal/mol, and −52.7Kcal/mol, respectively. The superimpositions of the pre-and post-MD docked structures revealed (Figure 7 ) that the phellodendroside, quercetin-7-O-rutinoside, and quercimeritrin complexes had RMSD of 1.186, 1.01, and 1.04, respectively, indicating that the docked complexes had less structural deviations during molecular dynamics simulations. The top three phytochemicals may have the ability to interfere with the operation of Dengue virus NS2B/NS3 protease, according to combinatorial bioinformatics methodologies such as molecular dynamics and molecular dynamics simulations. This work also tabulated fresh phytochemical information, which will aid researchers in using additional targeted viral proteins or signaling molecules. Although the data in this work were generated exclusively through computational methods, they still have to be confirmed in a wet lab setting using several enzymatic assays.

Based on the literature research, about 2000 phytochemical substances were collected from the Pubchem database [62] (Table S1 , Supplementary Materials). The structures were neutralized using Epik 2.2 at pH 7.0 ± 2.0, and the force field OPLS 2005 contained in Maestro, version 10.1 (Schrödinger suite) [63] was used to reduce them. Per ligand, up to 32 stereoisomers were preserved.

The RCSB PDB [64, 65] provided three-dimensional crystallographic structures of the proteins utilized in this study: Dengue virus NS2B/NS3 protease (PDB ID: 2FOM) [66] . To prepare and modify crystallographic structures for molecular docking experiments, the Protein Preparation Wizard of Maestro version 10.1 (Schrödinger suite) was used. Hydrogens were added to the heavy atoms, selenomethionines were transformed to methionines, and all waters were removed. Minimization was carried out using the force field OPLS 2005, with the maximum heavy atom RMSD (root mean square deviation) set to 0.30 [67, 68] .

For virtual screening, the Glide software from the Schrodinger working flow was used, which offers different docking protocols, HTVS (high-throughput virtual screening), standard precision, and SP. The best compounds from the pool were identified using the prior techniques. Every ligand was docked against the receptor using HTVS, which results in a single posture. About half of all plant-derived chemicals were shifted from HTVS to SP, reducing the number of false positives [69, 70] . Here, 1-(4-{5-[(piperidin-4-yl)methoxy]-3-[4-(1H-pyrazol-4-yl)phenyl]pyrazin-2-yl}phenyl)methanamine was used a control ligand molecules [71] . [72] . Here, the Glide pose viewer file of the best conformation was given as the source in Prime MM-GBSA simulation [73] . The total free energy of binding:

∆G bind = G complex − (G protein + G ligand ), where G = EMM + GSGB + GNP (1)

SwissADME [74] , admetSAR [75] , and pkCSM [76] were used to calculate the pharmacological and pharmacokinetic parameters of the three compounds. As an entry system, the ligand molecules' canonical simplified molecular input line entry system (SMILES) was used.

Molecular dynamics simulations of docked complexes were carried out using YASARA dynamics [77] and the AMBER14 force field [78] . Initial cleaning, optimization, and orientation of hydrogen bond networks were performed on docked complexes. In a cubic simulation cell with periodic boundary conditions, the TIP3P water solvation model was used [79] . The simulation cell's physiological conditions were set at 298K, pH 7.4, and 0.9% NaCl. The particle mesh Ewalds (PME) approach was used to calculate long-range electrostatic interactions with a cutoff radius of 8.0 [80] [81] [82] . The steepest gradient techniques were used to implement the initial energy minimization procedure using the simulated annealing approach (5000 cycles). The simulation's time step was set to 2.0 fs [83] . After every 100 ps period, the simulation trajectories were saved. The simulation was extended for 100 ns by using a constant pressure and Berendsen thermostat, and trajectories were used to compute the root mean square deviation (RMSD), root mean square fluctuations (RMSF), solvent accessible surface area (SASA), radius of gyration (Rg), and hydrogen bond [84] [85] [86] [87] .

To select the effective inhibitors from the compounds, this study used structurebased screening of phytochemicals. As a result, hundreds of phytochemicals have been identified in plants, particularly in Asian species. Based on greater binding affinity and MM-GBSA binding free energy, the top three compounds were chosen: phellodendroside (−63 kcal/mole), quercimeritrin (−59.5 kcal/mole), and quercetin-7-O-rutinoside (−54.1 kcal/mole). As a result, a molecular dynamics analysis was carried out to confirm the structural stability and binding pose. The toxicity and carcinogenicity of the tested molecules shows no possible adverse effect of the compounds. This study relies primarily on computational screening and pipelines; nonetheless, these studies must be confirmed in the wet lab using several enzymatic assays. 

Dengue and dengue hemorrhagic fever

The pathology of dengue hemorrhagic fever

The global distribution and burden of dengue

Dengue virus life cycle: Viral and host factors modulating infectivity

Secondary infection as a risk factor for dengue hemorrhagic fever/dengue shock syndrome: An historical perspective and role of antibody-dependent enhancement of infection

Global traffic and disease vector dispersal

The first licensed dengue vaccine: An important tool for integrated preventive strategies against dengue virus infection

Design of New Competitive Dengue Ns2b/Ns3 Protease Inhibitors-A Computational Approach

Higher catalytic efficiency of N-7-methylation is responsible for processive N-7 and 2 -O methyltransferase activity in dengue virus

Anthracene-based inhibitors of dengue virus NS2B-NS3 protease

Flavivirus Genome Organization, Expression, And Replication

In vitro processing of dengue virus type 2 nonstructural proteins NS2A, NS2B, and NS3

Both nonstructural proteins NS2B and NS3 are required for the proteolytic processing of dengue virus nonstructural proteins

In vitro processing of dengue virus structural proteins: Cleavage of the pre-membrane protein

Discovery of novel cyclic peptide inhibitors of dengue virus NS2B-NS3 protease with antiviral activity

Deletion analysis of dengue virus type 4 nonstructural protein NS2B: Identification of a domain required for NS2B-NS3 protease activity

Clinical efficacy and safety of a novel tetravalent dengue vaccine in healthy children in Asia: A phase 3, randomised, observer-masked, placebo-controlled trial

Inhibitory activity of cyclohexenyl chalcone derivatives and flavonoids of fingerroot, Boesenbergia rotunda (L.), towards dengue-2 virus NS3 protease

Flavonoids as noncompetitive inhibitors of Dengue virus NS2B-NS3 protease: Inhibition kinetics and docking studies

Identification of a Novel Inhibitor of Dengue Virus Protease through Use of a Virtual Screening Drug Discovery Web Portal

A small molecule inhibitor of dengue virus type 2 protease inhibits the replication of all four dengue virus serotypes in cell culture

Identification of novel thiadiazoloacrylamide analogues as inhibitors of dengue-2 virus NS2B/NS3 protease

Discovery of antiviral molecules for dengue: In silico search and biological evaluation

Nonsubstrate Based Inhibitors of Dengue Virus Serine Prote-ase: A Molecular Docking Approach to Study Binding Interactions between Protease and Inhibitors

Sulfated Polysaccharide, Curdlan Sulfate, Efficiently Prevents Entry/Fusion and Restricts Antibody-Dependent Enhancement of Dengue Virus Infection In Vitro: A Possible Candidate for Clinical Application

Structure and anti-dengue virus activity of sulfated polysaccharide from a marine alga

Activation of Peripheral Blood Mononuclear Cells by Dengue Virus Infection Depotentiates Balapiravir

Suramin inhibits helicase activity of NS3 protein of dengue virus in a fluorescence-based high throughput assay format

Carnosine exhibits significant antiviral activity against Dengue and Zika virus

Policresulen, a novel NS2B/NS3 protease inhibitor, effectively inhibits the replication of DENV2 virus in BHK-21 cells

The Small Molecules AZD0530 and Dasatinib Inhibit Dengue Virus RNA Replication via Fyn Kinase

Virion-associated phosphatidylethanolamine promotes TIM1-mediated infection by Ebola, dengue, and West Nile viruses

Luteolin restricts dengue virus replication through inhibition of the proprotein convertase furin

Investigative studies on the inhibition of amyloid-like fibrils formation by the extracts of Vernonia amygdalina Del. leaf

GNF-2 Inhibits Dengue Virus by Targeting Abl Kinases and the Viral E

Broad-spectrum agents for flaviviral infections: Dengue, Zika and beyond

Drug repurposing approaches to fight Dengue virus infection and related diseases

Dengue Antiviral Development: A Continuing Journey

Activity of Recombinant Dengue 2 Virus NS3 Protease in the Presence of a Truncated NS2B Co-factor, Small Peptide Substrates, and Inhibitors

Purified NS2B/NS3 Serine Protease of Dengue Virus Type 2 Exhibits Cofactor NS2B Dependence for Cleavage of Substrates with Dibasic Amino Acids in Vitro

Biochemistry and Medicinal Chemistry of the Dengue Virus Protease

Knowledge-based scoring function to predict protein-ligand interactions

Progress for dengue virus diseases: Towards the NS2B-NS3pro inhibition for a therapeutic-based ap-proach

Anti-dengue infectivity evaluation of bioflavonoid from Azadirachta indica by dengue virus serine protease inhibition

Antiviral phytochemicals as potent inhibitors against NS3 protease of dengue virus

Structure-guided screening of chemical database to identify NS3-NS2B inhibitors for effective therapeutic application in dengue infection

Novel Dengue Virus NS2B/NS3 Protease Inhibitors

Computational insight into dengue virus NS2B-NS3 protease inhibition: A combined ligand-and structure-based approach

Structureguided fragment-based in silico drug design of dengue protease inhibitors

Computational exploration of antiviral activity of phytochemicals against NS2B/NS3 proteases from dengue virus

Ur-Rahman, A. In silico identification of potential inhibitors of key SARS-CoV-2 3CL hydrolase (Mpro) via molecular docking, MMGBSA predictive binding energy calculations, and molecular dynamics simulation

Bridging Molecular Docking to Molecular Dynamics in Exploring Ligand-Protein Recognition Process: An Overview

Endocrine Disruption at the Androgen Receptor: Employing Molecular Dynamics and Docking for Improved Virtual Screening and Toxicity Prediction

Structure-function relationships in ABCG2: Insights from molecular dynamics simulations and molecular docking studies

Molecular Docking of Ferulic Acid Derivatives on P2Y12 Receptor and their ADMET Prediction

Open access in silico tools to predict the ADMET profiling of drug candidates

Homology Modelling and Molecular Docking Studies of Selected Substituted Benzo[d]imidazol-1-yl)methyl)benzimidamide Scaffolds on Plasmodium falciparum Adenylosuccinate Lyase Receptor

Computational analysis of calculated physicochemical and ADMET properties of protein-protein interaction inhibitors

Molecular docking and ADMET studies of Allium cepa, Azadirachta indica and Xylopia aethiopica isolates as potential anti-viral drugs for Covid-19

PubChem substance and compound databases

Schrödinger Suite; Schrödinger Release

The Protein Data Bank

RCSB Protein Data Bank: Enabling biomedical research and drug discovery

Structural basis for the activation of flaviviral NS3 proteases from dengue and West Nile virus

Schrödinger Release 2015-1; LLC Schrödinger Suite

Protein Preparation Wizard

Virtual Screening for Ligand Discovery at the σ1 Receptor

Targeting a hidden site on class A beta-lactamases

X-ray Crystallography and Antiviral Activity of Allosteric Inhibitors of Flavivirus NS2B-NS3 Protease

Performance of the MM/GBSA scoring using a binding site hydrogen bond network-based frame selection: The protein kinase case

Binding-affinity predictions of HSP90 in the D3R Grand Challenge 2015 with docking, MM/GBSA, QM/MM, and free-energy simulations

SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

admetSAR: A Comprehensive Source and Free Tool for Assessment of Chemical ADMET Properties

Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures

A tool to obtain structural guidance in biocatalytic investigations

Development and testing of a general amber force field

Structure and dynamics of TIP3P, TIP4P, and TIP5P water near smooth and atomistic walls of different hydroaffinity

A smooth particle mesh Ewald method

An Implementation of the Smooth Particle Mesh Ewald Method on GPU Hardware

Fast empirical pKa prediction by Ewald summation

New ways to boost molecular dynamics simulations

In silico design of an epitope-based vaccine against choline binding protein A of Streptococcus pneumoniae

Identification and in silico molecular modelling study of newly isolated Bacillus subtilis SI-18 strain against S9 protein of Rhizoctonia solani

Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: A molecular docking study

Extensive immunoinformatics study for the prediction of novel peptide-based epitope vaccine with docking confirmation against envelope protein of Chikungunya virus: A computational biology approach

The authors declare no conflict of interest.Sample Availability: Samples of the compounds are available from the authors.