Introduction
To bring out a new drug from the laboratory to the market, it takes approximately 10-12 years with an average cost US $1.2 to $1.4 billion or more per drug.1 Thus, there is a need for some fast processes and methods for the designing, discovery and development of a new chemical entity (NCE). There are many computational and Insilico tools such as genomics, proteomics, bioinformatics, and excellent technologies like, combinatorial chemistry, high throughput screening (HTS), virtual screening in vitro, in silico ADMET screening, de novo and structure-based drug design that plays an important role to speed up the modern era drug discovery and development process. 2 Computational drug design approaches are mainly concentrating on the design of ligands/molecules for active or target sites with identified three-dimensional structure. Computational methods are also helpful to check the Drug- likeness properties of the molecule, after that molecule is docked with the target, selected according to their binding affinities. These molecules are further optimized to enhance binding characteristics and the toxicity is predicted by using different online web servers. Computational techniques integrate biological, mathematical and computer-based models to predict the most favourable binding conformation of ligands in the active site of the particular receptor. The prediction power of these techniques is increasing day by day due to advancement in the field of Molecular Biology, Biotechnology, Bioinformatics, Mathematics, and Chemistry.3 The insilico drug discovery process is described in figure below-
To select a novel drug candidate, in-silico or computational methods play key role to design and develop new drugs and target proteins in the fields of pharmaceutics and biotechnology. These methods are used to investigate molecular modeling of gene, gene sequence, gene analysis and 3D structure of proteins which play a great role in target identification and prediction of new chemical entity.
Methods used in in-silico drug designing
Some important in-silico techniques which are helpful in drug design are follows.
Homology modelling
Homology modelling allows making an unidentified model of a target protein according to its amino acid arrangement and experimental three-dimensional structure of a related homologous protein and displays similarity with the sample sequence. 4
Molecular docking
Molecular docking includes the interaction of two or more molecules to form a stable complex. Stability of the complex depends upon the binding properties of ligand and target protein. 5 The most appropriate complexes are selected on the basis of scoring function in the software. The molecules are docked against the active site and then scored to identify the one which binds more tightly to the target protein. 6
Virtual screening
Virtual screening allows testing of large libraries of molecules and compounds for their potential to interact with specific site of target protein which can screen more than thousands of possible molecules to a practicable number that can be easily synthesized, purchased and tested in the laboratory. 7
Quantitative structure-activity relationship
Quantitative structure-activity relationship is widely used to demonstrate the relationship between the structural properties of molecules with their biological activities. 8
Comparative molecular field analysis
Comparative molecular field analysis is a well-recognized 3D QSAR method which gives the values of ClogP which shows the solvent repellent constraints of the ligands and also describes the steric and electrostatic values of the ligands. 9
Comparative molecular similarity indices analysis
Comparative Molecular Similarity Indices Analysis is used in the drug discovery process to recognize the common properties, which are important for the suitable biological receptor binding with ligands. 10
3D Pharmacophore mapping
3D Pharmacophore mapping quickly predicts lead compounds along with a preferred target which has been used to make it most powerful and successful insilico and computational method.11
Conformational analysis
Conformations of molecules play a key role in the prediction of the physio-chemical properties but also helpful in prediction of biological activity of the compound. 12
Monte carlo simulation
Monte Carlo simulation helps in the generation of suitable conformations of a system by using computer simulation to allow thermodynamic, structural, and numerical properties to be calculated as a weighted average of these properties over these conformations. 13
In this review we describe the different medicinal compounds which can be screened by the use of computational methods.
Materials and Methods
Selection and refinement of various targets of JEV protein
The reported three-dimensional (3D) crystal structures of different microbial proteins were retrieved from Protein Data Bank (PDB) available at (http://www.rcsb.org/pdb/home/home.do). Retrieved structures were downloaded in PDB format and exported into the Maestro software (Schrodinger, LLC, Cambridge, USA) and refined by mediating insertion of processing, optimization, and minimization steps indicated in protein preparation wizard.
Binding site detection
Site map tool from Schrodinger was used to find the active sites in protein which have not internal ligand in their structure. The site score and draggability scores were used to select potential binding sites. The sites having the draggability score near to 1 were considered for docking study. 14
Preparation of ligands
Chemical structures of the herbal compound were obtained from the PubChem chemistry database, and the structures drawn and optimized in ChemDraw (Perkin Elmer Informatics, Waltham, MA, USA). The low energy 3D conformation of ligand was executed with the support of LigPrep OPLS3e force field. The ionization states were also adjusted prior to docking simulations.
Molecular docking
Maestro suites (Schrodinger, LLC, Cambridge, USA) was used to dock herbal compounds against the target proteins of microbes. The receptor grid was generated by specifying the internal ligand and by doing the site mapping for the proteins which have not internal ligand in their structure. 15
3.5. Molecular dynamics simulations were performed for herbal compounds towards the target proteins by using Desmond module of Schrödinger with 100ns simulation time. The root mean square deviation (RMSD) for the protein backbone and root mean square fluctuation (RMSF) of the protein residues were plotted to examine the convergence of the ligands to equilibrium.
Review of Literature
In-silico study or computational screening of compounds against SARS
Shen et al.16 have screened 290 compounds against HCoV-OC43 strain and found 27 compounds to be effective against SARS.
Table 1
Table 2
Jo et al.17 have screened various compounds through induced fit docking studies and proteolytic assay and found that compounds such as Herbacetin, rhoifolin pectolinarin were found to be effective against SARS CoV.
Nguyen et al.23 have screened various flavonoids through molecular docking studies against 3CLpro of SARS CoV and found gallocatechin gallate (with binding energy -14kCal/mol) to be the best compound effective against SARS-CoV. Ryu et al.28 reported flavones such as apigenin, luteolin and quercetin were able to inhibit 3CLpro activity with IC50 values of 280.8, 20.2, and 23.8 µM respectively.
Chen et al. 56 performed molecular docking studies and SPR/FRET-based bioassay on quercetin-3-β-galactoside which might be effective against SARS CoV through inhibition of 3CLproactivity. Ho et al. 35 have performed computational screening of 312 Chinese medicinal herbs and found 3 compounds belonging to the family Polygonaceae to be effective against SARS virus.
Wang et al.57 screened 11 molecules from TCMD primarily based on template strain MDL 28170 and found MOL736 (derived from Artemisia annua) with the lowest binding energy and highest similarity reported to be effective against SARS CoV by inhibiting endosomal protease activity of SARS virus.
Toney et al. 58 have screened 14000 molecules and reported 1853 molecules to serve as 3CLpro inhibitors. Further molecular docking studies showed that the 10 compounds have the lowest docking energy (-11.6 kCal/mol) with desirable chemical properties as therapeutics.
Herbals used against SARS (in-vitro study)
Many of the natural compounds like lycorine, emetine, herbacetin, pectolinarin, rhoifolin, glycyrrhiza, lycoris radiata, aloe-emodin, baicalin and isatis indigotica are found to be effective against in SARS infection by inhibiting viral replication and nucleocapsid protein against SARS CoV strain. The relevant studies have been extracted from Pub Med and compiled in the below mentioned Table 1.
Correlation of In-silico data to Experimental Data
In order to screen the thousands of compounds in library the in-silico predictions are the best tools to select the compound which has potential affinity to inhibit the JE infection. However, the experimental data also must determine the activity of compound against specific target and the data obtained in the in-silico may not be reliable as compared to in-vitro or in-vivo.
The drugs show different activity in-vivo as compared to in-silico and in-vitro, this may be due to poor pharmacokinetic properties or physicochemical properties or inability to interact with the target at particular environment etc. So, in current review we have correlated the in-silico and experimental (in-vitro or in-vivo) activity of the molecule against JEV infection. However, we have found valid and similar results against JE infection under in-silico and also in-vitro and in-vivo except for some compounds.
By computer simulations, Chemdiv-3 was found to be effective against JEV, and at a concentration of 20µM has potential to inhibit the viral growth in animal (Balb-2 strain mice), it has shown moderate activity against JE infection and the mice survived 3 more days compared to control mice (which JE infected mice). 46 In case of Phytolacca-mother tincture, excellent results were shown in-silico as well as in-vivo (in JE infected human), at the concentration of 26.5µM upon daily administration. 55
Kaempferol and daidzein have also shown great activity by computational simulations. However, the experimental data has shown that moderate activity (in-vitro) towards JE infected BHK-21 cells which inhibitted 70% of viral replication at the concentration of 25.7 µM kaempferol and 29.7 µM for daidzein. 48 Doxycycline and kanamycin are capable of inhibiting JE infection and inflammation under in-vivo conditions similar to computational study. But, doxycycline has shown more activity than kanamycin intraperitoneally into Swiss albino mice. 59 CW-33 has shown similar interaction under in-silico and in-vitro by inhibiting the RNA replication against JEV. 50
Niclosamide, Cilnidipine and FGIN-1-27 and Rosameric acid have shown very effective against JEV by computational simulation. However, Niclosamide, Clinidipine have shown poor activity in-vivo (mice) due to poor penetration into blood brain barrier, whereas Rosameric acid and FGIN-1-27 have shown similar activity as in-silico which prevented increasing of viral load in JE infected mice. 51, 52 Fang et al. 53 have discovered 2 azole compounds which have shown similar interaction in-silico whereas, compound 1 has shown 71% inhibition and compound 2 has shown 95% inhibition in-vitro.
Other drugs like 4-hydroxypandurantin, Aminotetrahydroquinazoline derivative were tested only by in-silico approaches, further investigation is required for the exploration of those drugs for the treatment of JEV infection. 54
Conclusion
Studies conducted on treatment of viruses like SARS CoV and Japanese Encephalitits Virus with various herbal compounds. The efficacy of the drugs or compounds was tested by different In silico and In vitro methods. The drugs have shown varied binding efficacy targeting different viral proteins and few of them interacting with viral RNA. The inhibitory action of the drug also varied with the concentration. However, the behavior of the drug is different In silico and In vitro.