Abstract
Free full text
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.
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.
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.
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).
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
Articles from Bioinformation are provided here courtesy of Biomedical Informatics Publishing Group
Full text links
Read article at publisher's site: https://doi.org/10.6026/97320630015233
Read article for free, from open access legal sources, via Unpaywall: http://www.bioinformation.net/015/97320630015233.pdf
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
In silico testing of flavonoids as potential inhibitors of protease and helicase domains of dengue and Zika viruses.
PeerJ, 10:e13650, 04 Aug 2022
Cited by: 2 articles | PMID: 35945938 | PMCID: PMC9357371
High Throughput Screening Targeting the Dengue NS3-NS5 Interface Identifies Antivirals against Dengue, Zika and West Nile Viruses.
Cells, 11(4):730, 18 Feb 2022
Cited by: 5 articles | PMID: 35203378 | PMCID: PMC8870125
Dronabinol as an answer to flavivirus infections: an in-silico investigation.
J Biomol Struct Dyn, 41(20):11219-11230, 28 Dec 2022
Cited by: 1 article | PMID: 36576139
Peptide derivatives as inhibitors of NS2B-NS3 protease from Dengue, West Nile, and Zika flaviviruses.
Bioorg Med Chem, 27(18):3963-3978, 23 Jul 2019
Cited by: 16 articles | PMID: 31351847
Review