Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


Dengue, West Nile and Zika virus belongs to the family flaviviridae and genus flavivirus. It is of interest to design and develop inhibitors with improved activity against these diseases. We used the helicases target to screen for potential inhibitors against these viruses using molecular docking analysis. NS3 helicases of flavivirus family of viruses such as Dengue, West Nile and Zika are prime targets for drug development. The computer aided molecular docking analysis of netropsin and novobiocin with the viral protein targets HABD, MTD and RCD is reported for further consideration.

Free full text 


Logo of bioinformLink to Publisher's site
Bioinformation. 2019; 15(4): 233–239.
Published online 2019 Mar 30. https://doi.org/10.6026/97320630015233
PMCID: PMC6599439
PMID: 31285639

Molecular docking analysis of netropsin and novobiocin with the viral protein targets HABD, MTD and RCD

Abstract

Dengue, West Nile and Zika virus belongs to the family flaviviridae and genus flavivirus. It is of interest to design and develop inhibitors with improved activity against these diseases. We used the helicases target to screen for potential inhibitors against these viruses using molecular docking analysis. NS3 helicases of flavivirus family of viruses such as Dengue, West Nile and Zika are prime targets for drug development. The computer aided molecular docking analysis of netropsin and novobiocin with the viral protein targets HABD, MTD and RCD is reported for further consideration.

Keywords: Dengue, therapeutics, helicase ATP binding domain, methyl transferase, RNA catalytic domain

Background

Dengue infection is a mosquito-borne human pathogen affecting by and large in tropical regions where 3.9 billion individuals live with around 100 million sicknesses showing clinical symptoms every year 1. Transmission of the disease occurs through mosquito bite of the Aedes class (for the most part Aedesaegypti and A. albopictus) 2. These mosquito vectors are commonly found in the middle of tropical regions 3. Most contamination is either asymptomatic or it causes smooth reactions (fever, joint agony or rashes) for a couple of days. However, available data on dengue hemorrhagic fever and dengue stun disorder raise to 20,000 cases in recent years 4. The dengue contamination has a spot with the flavi viruses family in addition to West Nile infection (WNV), dengue infection and the rising Zika infection (ZIKV) 5. Moreover, arbo viruses and Zika disease are also found to erupt lately in these regions 6. Zika virus contaminates are found in South America, Central America, and the Caribbean in addition to arbo viruses in the Western Hemisphere during the last two decades 7, 8. Without doubt, Zika diseases seeks attention after West Nile infection, dengue, which appeared in 1999 and chikungunya which appeared in 2013 8. The Zika contamination has a spot with the assortment flavivirus and is generally vectored by Aedes mosquitoes 9 found over the continents 10-12.

Helicases are omnipresent engine proteins that catalyze the unwinding of double-stranded DNA or RNA by ATP hydrolysis. They move along a nucleic acid phospho-di-ester spine and separate the complementary strands by utilizing energy got from ATP hydrolysis. 13 RNA helicases are important particles for RNA metabolic procedures, for example, ribosome biogenesis, joining and interpretation. It has been demonstrated that the helicases are connected to different segments of the macromolecular machines and they have critical role in numerous processes 14. Hence, the consideration of helicases as potential viral drug targets is realized.

The viral proteome is an extended polypeptide chain in the ER. It is further known that NS2A, NS2B, NS4A, and NS4B are structural layer proteins essential for viral replication through possible protein-protein and protein-lipid associations. The functions played by NS1, NS2A and NS4A in viral replication are complex to understand 15. NS1 has different oligomerization states dependent on its glycosylation status 16-18. The methyltransferase is in charge of topping the incipient genomic RNA by successively utilizing S-adenosyl methionine as the methyl donor. This is done by consecutive methylation on the N7 molecule with the top guanine and the 2'O particle of the ribose in the main adenine 19-20. Defect in capping decrease viral multiplication to cause infection that are gradually subtle to the natural resistant reaction as they stimulate higher interferon (IFN) flagging and immunizer reactions 21. Moreover, the positive sense of viral RNA is discharged in the cytoplasm during infection. NS5 protein initially deciphers it as a negative sense strand before utilizing the negative strand (with regards to a dsRNA middle of the road) to organize excess of positive sense RNA. The viral mRNA is then used to express the poly protein by host cell ribosomes. It is also known that a hetero-duplex framed by a DNA format and an RNA primer is formed 22. This component is available in the DENV NS5 full-length protein whose preference for dsDNA is like ssRNA 23. In the context of these available data it is of interest to screen HABD (helicase ATP binding domain), MTD (methyl transferase domain) and RCD (RNA catalytic domain) for potential inhibitors as drug candidates.

Methodology

Multiple sequence alignment (MSA) of the conserved domain:

The conserved domain sequences of dengue virus, zika virus and Japanese encephalitis virus were retrieved from the genome database in NCBI and BLASTp analysis was completed. The FASTA formats of the retrieved sequences were used for further analysis. Multiple sequence alignment (MSA) of Helicase ATP binding domain (201 amino acid residues), Methyl transferase domain (262 amino acid residues) and RNA Catalytic domain (149 amino acid residues) was completed using ClustalW3 and Clustal Omega. Various domains were manually assigned and confirmed by using Pfam, Prosite, SMART, PANTHER and InterProScan.

Three-dimensional structure prediction by I-TASSER:

The protein sequence (201) of Helicase ATP binding domain (HABD), the protein sequence (262) of Methyl transferase domain (MTD) and the protein sequence (149) of RNA Catalytic domain (RCD) were downloaded from the Swiss Pro database. The threedimensional model was produced utilizing the I-TASSER server which creates a 3D model of inquiry arrangement by different threading arrangements and iterative necessary gathering reenactment 24. We used this server because its accessibility, composite methodology of displaying an execution in CASP rivalry. I-TASSER technique incorporates general strides of threading, and auxiliary get together, display determination, refinement, and structure-based comments 25. An optional structure was developed by PSIPRED 26 utilizing the structure library LOMETS 27. Z-score checked the nature of the formal arrangement followed by threading arrangements 24 using Monte Carlo simulation 28. The reenactment incorporates Ca/side chain connection, H-securities, hydrophobicity, spatial controls from threading template 27 and arrangement based on contact expectations from SVMSEQ 29. The adaptations produced amid the refinement reenactment process were bunched by SPICKER 30. Furthermore, the normal of three-dimensional directions of all the grouped structure was determined to acquire bunch centroids. In the refinement process, the chosen bunch centroids were again used to perform further reenactment, which evacuates steric conflicts to refine the topology of the group centroids. The known structures in PDB were recognized by TM-adjust 31. The final structural models were generated by REMO 32 in which group centroids of second-round reenactment were utilized. The useful analogs were ranked based on TM-score, RMSD, arrangement character, and the inclusion of the structural arrangement. The nature of the model was dictated by C-score (certainty score), which is -5 to 2. It depends on threading arrangement and the combination of auxiliary for refinement reenactments.

Preparation of ligands and protein targets:

The 3D structures of HABD, MTD and RCD were built using ITASSER. The hydrogen atoms having polar nature were then included. The buildup structures less in numbers were removed and the fragmented side chains were later replaced by Auto Dock Tools (ADT) version 1.5.6 downloaded from the Scripps Research Institute. Further, particles having Gasteiger charges were included and the non-polar hydrogen iotas were added to the protein structure. The built structures were then stored in PDBQT format in ADT 33. 3D structures of netropsin and novobiocin were drawn using ChemBioDraw Office 12.0. The 3D co-ordinates of the ligands was then created and stored in PDBQT format using ADT 33.

Receptor grid formation:

Networks predetermine matrix maps of restricting energies in various particle types (for example, hydrogen holding oxygen, carbons, and aliphatic carbons in a macromolecule (for example a RNA/DNA, protein)) before docking 34. These network maps are then utilized in AutoDock 4.2 docking computations to characterize the absolute restricting vitality for a ligand with a macromolecule 33. Network mapping computes the vital parameters over the protein nuclear information and determines the directions of the HABD, MTD, and RCD for docking. Likewise, lattice mapping manages an appropriate surface topology for the iotas of mixes for association with the HABD, MTD and RCD dynamic sites. Network mapping is a necessity to guide netropsin, and novobiocin mixes to search for their locale for binding with the HABD, MTD and RCD dynamic sites.

The network measurements for the HABD protein was 58 x 48 x 52 matrix focused with separating 1.00 Å between the framework focused on the ligand for protein (59.686, 70.660 and 46.318 directions). The framework measurements for MTD protein was 48 x 52 x 58 matrix focused with dispersing 1.00 Å between the network focuses; however fixated on the ligand for protein (61.363, 69.710 and 54.870 co-ordinates). The lattice measurements for malate synthase protein was 52 x 52 x 55 framework focused with dividing 1.00 Å between the framework focuses yet fixated on the ligand for protein (60.820, 65.487 and 70.722 directions). The network was made for selecting promising cooperation for docking between ligands and targets 35.

Molecular docking:

AutoDock 4.2 with standard parameters was used to dock the netropsin, and novobiocin mixes into the dynamic sites of HABD, MTD, and RCD 33, 36. The energy calculation was done using the Lamarckian hereditary estimation (LGA). The adaptations with the most significant free restricting vitality were chosen for dissecting the connections between the receptor and ligands using PyMOL and Discovery studio 33.

Results and Discussion

Sequence analysis and domain organization:

An analysis of amino acid alignment sequence of three unique domain that is Helicase ATP binding domain (201 amino acid), Methyl transferase domain (262 amino acid) and RNA Catalytic domain (149 amino acid) shown in Figure 1 was done using the Clustal Omega program. The InterProScan examination demonstrates that Dengue, Zika infection and Japanese Encephalitis viruses have conserved motifs (Figure 2 A-F). The comparative study provided useful insights on conserved domains, which are present in all three viruses. This is useful to study targets for antiviral drug design.

An external file that holds a picture, illustration, etc.
Object name is 97320630015233F1.jpg

Multiple sequence alignment (MSA) of conserved domain amino acid sequence of dengue virus, zika virus and japanese encephalitis virus. MSA was done using the CLUSTALW2 program (http://www.ebi.ac.uk).

An external file that holds a picture, illustration, etc.
Object name is 97320630015233F2.jpg

Domain organization of Dengue virus, Zika virus and Japanese encephalitis is shown. The conserved sequences of each domain are written inside the boxes. The text written in a box of various conserved domains and the numbers refer to the amino acids spanning the various domains.

Three-dimensional structure models:

The molecular model of HABD, MTD, and RCD proteins was developed using I-TASSER. The models obtained from server incorporate supporting structure with certainty score (0 to 9), anticipated dissolvable openness. We obtained five models with C-score, top ten formats from PDB used in the arrangement; top ten PDB basic analogs, useful analogs of protein, and restricting site deposits. The HABD model (Figure 3A) was chosen with Cscore 0.92, TM-score 0.60 ± 0.14 and RMSD 7.3 Å 4.2 Å. MTD (Figure 3B) was selected as the best-expected model with C-score 1.17, TM-score 0.87±0.07, and RMSD 3.6 ± 2.5 Å. MTD (Figure 3 C) was chosen as the best model with C-score 0.99, TM-score 0.85 ± 0.08, and RMSD 2.8 ± 2.1 Å. C-score with higher esteem mirror a model for better quality 18. Standardized Zscore commonly evaluates threading. A normalized Z-score > 1 esteem mirror a specific arrangement. TM-adjust distinguished 5xdrA, 2px5A and 4v0qA2 in PDB library as the best scoring models from I-TASSER with a TM-score of 0.792, 0.979 and 0.979, respectively.

An external file that holds a picture, illustration, etc.
Object name is 97320630015233F3.jpg

(A) (a) Three-dimensional structure of HABD protein anticipated by I-TASSER. (b) Alignment of question protein (Green) with auxiliary simple (Pink) 5xdrA in PDB library. (B) (c) Threedimensional structure of MTD protein developed by I-TASSER is shown. (d) Alignment of question protein (green) with basic simple (pink) 2px5A in PDB library. (C) (e) Three-dimensional structure of RCD protein anticipated by I-TASSER.(f)Alignment of inquiry protein (green) with auxiliary simple (pink) 4v0qA2 in PDB library.

Molecular docking studies:

Auto docking 4.2 was used to set up the coupling capacity of the netropsin, and novobiocin mixes with HABD, MTD, and RCD for docking interactions. We examined the interaction of HABD, MTD and RCD with the inhibitors using PyMOL. Results show that the compound netropsin is compatible with the dynamic site of HABD, MTD and RCD with least restricting vitality (ΔG) -6.7 kcal/mol, (ΔG) -7.3 kcal/mol and -6.3 kcal/mol. The compound novobiocin is pleasantly limited into the dynamic site of HABD, MTD and RCD with least restricting vitality (ΔG) -7.5 kcal/mol, (ΔG) - 8.9 kcal/mol and -7.8 kcal/mol. The compound netropsin forms four hydrogen bonds - one each with HABD buildup (Ala7; 2.2 Å, Ile87; 2.6 Å, Tyr89; 2.6 Å and Thr91; Ile87; 2.4 and 2.3 Å) (Figure 4A (b).

An external file that holds a picture, illustration, etc.
Object name is 97320630015233F4.jpg

Molecular docking of compounds with HABD: (A) (a) 2D schematic diagram showing interactions of compound netropsin. (b) Cartoon view of HABD protein with compound netropsin. (B) (c) 2D schematic diagram showing interactions of compound novobiocin. (d) Cartoon view of HABD protein with compound novobiocin. Residues involved in hydrogen bonding, van der Waals interactions, carbon hydrogen and Pi-alkyl are represented in different color indicated in the inset. (C) (e) 2D schematic outline appearing of compound netropsin. (f) Cartoon perspective on MTD protein with compound netropsin. (D) (g) 2D schematic outline appearing of compound novobiocin. (h) Cartoon perspective on MTD protein with compound novobiocin. Deposits engaged with hydrogen holding, van der Waals (vdW) cooperation, carbon hydrogen and Pi-alkyl are spoken to in various shading shown in the inset. (E) (i) 2D schematic graph appearing of compound netropsin. (j) Cartoon perspective on RCD protein with compound netropsin. (F) (k) 2D schematic graph appearing of compound novobiocin. (l) Cartoon perspective on RCD protein with compound novobiocin. Deposits engaged with hydrogen holding, van der Waals (vdW) cooperation, carbon hydrogen and pi-alkyl are spoken to in various shading demonstrated in the inset.

The compound netropsin is compatible inside the HABD binding site by connecting with different deposits appeared in Figure 4C (e). The compound novobiocin forms three hydrogen bonds with the HABD buildup (Gly43; 2.3 Å and Lys48; 2.3 and 2.6 Å) as shown in Figure 4D (h). The compound novobiocin likewise communicates with the HABD restricting site by associating with different deposits appeared in Figure 4D (g). The compound netropsin through six hydrogen bonds with MTD (Cys81; 2.3 Å, Gly82; 2.2 Å, Lys104; 2.9 Å, His109; 3.3 Å, Glu148; 3.4 Å and Asp145; 2.4 Å) as shown (Figure 4E (j). It communicates with the MTD binding site by interfacing with different deposits appeared in Figure 4E (i). Novobiocin associated well through three hydrogen bonds MTD buildup (Ser55; 2.6 Å, Gly80; 2.8 Åand Arg159; 2.2 and 2.5 Å,) appeared in Figure 4F (l). The compound novobiocin likewise associates with the MTD restricting site by collaborating with different build ups appeared in Figure 4F (k). Compound netropsin showed three hydrogen bonds with RCD build up (Asp6; 2.5 Å and Gly80; 2.1 and 2.5 Å) (Figure 4F (l). The compound netropsin likewise accepts good interaction inside the RCD restricting site by interfacing with different buildups as shown in Figure 4E (i). Novobiocin collaborated through six hydrogen securities RCD build up (Asn22; 2.0 Å, Glu23; 2.0 Å, Ile26; 2.6 Å, Tyr80; 2.2 and 2.6 Å, Gly81; 2.5 Å and Ser134; 2.5 Å) as shown in Figure 4F (l). The compound novobiocin likewise collaborates with the RCD restricting site by associating with different deposits as shown in Figure 4F (k).

Conclusion

We report the molecular docking analysis of netropsin and novobiocin against the helicases targets of Dengue, West Nile and Zika viruses. The analysis shows that netropsin and novobiocin bind to viral targets HABD, MTD and RCD with high binding ability for further in vitro and in vivo studies.

Conflict of Interest

Authors declare no conflict of interest.

Acknowledgments

Authors are grateful to King Abdulaziz City for Science and Technology (KACST, Riyadh, Saudi Arabia), General Directorate of the research grants program for funding this study with grant No. (1-18-01-009-0035). This work was also funded by the Deanship of Scientific Research (DSR) grant, King Abdulaziz University, Jeddah, Saudi Arabia, The authors, therefore, acknowledge the DSR technical and financial support.

Edited by P Kangueane

Citation: Alwabli et al. Bioinformation 15(4): 233-239 (2019)

References

1 . Wilder-Smith A, et al. The Journal of infectious diseases. 2016;214:1796. [Abstract] [Google Scholar]
2 . Gratz N. Medical and veterinary entomology . 2004;18:215. [Abstract] [Google Scholar]
3 . Machado CM, et al. Revista do Instituto de Medicina Tropical deSao Paulo . 2009;51:309. [Abstract] [Google Scholar]
4 . Kuno G. Dengue and dengue hemorrhagic fever. New York: CABInternational . 1997;23 [Google Scholar]
5 . Samarasekera U, Triunfol M. The Lancet. 2016;387:521. [Abstract] [Google Scholar]
6 . Benelli G, Mehlhorn H. Parasitol Res. 2016;115:1747. [Abstract] [Google Scholar]
7 . Attar N. Nature Reviews Microbiology. 2016;14:62. [Abstract] [Google Scholar]
8 . Fauci AS, Morens DM. New England Journal of Medicine. 2016;374:601. [Abstract] [Google Scholar]
9 . Marcondes CB, Ximenes MdFFd. Revista da SociedadeBrasileira de Medicina Tropical . 2016;49:4. [Abstract] [Google Scholar]
10 . Mehlhorn H. Arthropods as vectors of emerging diseasesSpringer Science and Business Media. 2012;3:38. [Google Scholar]
11 . Becker N, et al. Exotic mosquitoes conquer the world, 2012;31 [Google Scholar]
12 . Melaun C, et al. Parasitology Research . 2015;114:1051. [Abstract] [Google Scholar]
13 . Cordin O, et al. Gene . 2006;367:17. [Abstract] [Google Scholar]
14 . Bae S, et al. Nature Nanotechnology. 2010;5:574. [Abstract] [Google Scholar]
15 . Yu LK, et al. Virology. 2013;446:365. [Abstract] [Google Scholar]
16 . Rastogi MN, et al. Virology Journal. 2016;13:131. [Europe PMC free article] [Abstract] [Google Scholar]
17 . Somnuke P, et al. Virology . 2011;413:253. [Europe PMC free article] [Abstract] [Google Scholar]
18 . Avirutnan P, et al. PLoS Pathogens. 2007;3:183. [Europe PMC free article] [Abstract] [Google Scholar]
19 . Zhao Y, et al. Proceedings of the National Academy of Sciences. 2015;112:14834. [Google Scholar]
20 . Issur M, et al. RNA . 2009;15:2340. [Europe PMC free article] [Abstract] [Google Scholar]
21 . Decroly E, et al. .Nature Reviews Microbiology. 2012;10:51. [Europe PMC free article] [Abstract] [Google Scholar]
22 . Lu G, et al. Antimicrobial Agents and Chemotherapy. 2017;61:1967. [Google Scholar]
23 . Szymanski MR, et al. J Biol Chem . 2011;286:33095. [Europe PMC free article] [Abstract] [Google Scholar]
24 . Wu S, et al. BMC Biology . 2007;5:17. [Europe PMC free article] [Abstract] [Google Scholar]
25 . Roy AA, et al. Nature Protocols . 2010;5:725. [Europe PMC free article] [Abstract] [Google Scholar]
26 . Jones DT. Journal of Molecular Biology . 1999;292:195. [Abstract] [Google Scholar]
27 . Wu S, Zhang Y. Nucleic Acids Research . 2007;35:3375. [Europe PMC free article] [Abstract] [Google Scholar]
28 . Zhang Y, et al. Proteins: Structure, Function and Bioinformatics. 2002;48:192. [Abstract] [Google Scholar]
29 . Wu S, Zhang Y. Bioinformatics . 2008;24:924. [Europe PMC free article] [Abstract] [Google Scholar]
30 . Zhang Y, Skolnick J. Journal of computational chemistry. 2004;25:865. [Abstract] [Google Scholar]
31 . Zhang Y, Skolnick J. Nucleic Acids Research . 2005;33:2302. [Europe PMC free article] [Abstract] [Google Scholar]
32 . Li Y, Zhang Y. Proteins: Structure, Function, and Bioinformatics. 2009;76:665. [Europe PMC free article] [Abstract] [Google Scholar]
33 . Ahmad K, et al. Letters in Drug Design and Discovery . 2017;14:488. [Google Scholar]
34 . Gill BS, Kumar S. Medicinal Chemistry Research. 2015;24:3483. [Google Scholar]
35 . Repasky MP, et al. Current Protocols in Bioinformatics . 2007;18:8. [Abstract] [Google Scholar]
36 . Morris GM, et al. J. Comput. Chem. . 1998;19:1639. [Google Scholar]

Articles from Bioinformation are provided here courtesy of Biomedical Informatics Publishing Group

Similar Articles