BACKGROUND
herpes simplex virus type  <dig> and type  <dig>  are two members of herpesviridae family, which infect almost 85% of the world population  <cit> .
they are the causative agents of a gamut of diseases ranging from mild ones like cold sores in mouth, eye cornea and genitals to more severe life threatening ones like the fatal herpes encephalitis  <cit> .
people with suppressed immune system like those suffering from aids are more prone to get infections from hsv than others  <cit> .
drugs for hsv
there is no permanent cure for these infections till date.
present day treatment involves the use of antiviral drugs to reduce the physical severity of outbreak-associated lesions and viral shedding, though this helps decreasing the chances of transmission to others only by maximum 50%  <cit> .
there are two types of drugs that are clinically useful against hsv infections.
the first category consists of nucleoside analogs like acyclovir and its prodrug valacyclovir, ganciclovir, penciclovir and its prodrug famciclovir, sorivudine and brivudine.
these require phosphorylation by viral thymidine kinase to form triphosphates that are active inhibitors of viral dna polymerase.
the second category consists of direct viral dna polymerase inhibitors like vidarbine, foscarnet and cidofovir  <cit> .
thus, both types of drugs target in dysfunctioning the replication centre i.e.
dna polymerase of the viral genome  <cit> .
development of alternative treatments for hsv
however, in past few years, a number of acyclovir drug resistant viral strains have been isolated especially from immuno-compromised patients  <cit> .
in this era, where the number of immuno-suppressed patients like those suffering from hiv is continuously increasing, there is an immediate need to find new drugs to treat hsv infections which have a higher efficacy or have an alternative mode of action  <cit> .
resistance to acyclovir is mainly due to mutations in the viral thymidine kinase  gene which impair the initial drug phosphorylation  <cit> .
these drug resistant strains have been of significant clinical attention  <cit> , indicating the need for alternative anti- hsv drugs.
previous antiviral studies on herpes viruses have focussed on developing nucleoside analogues that will inhibit viral polymerase and terminate the replicating viral dna.
a number of new anti-viral drugs against hsv dna polymerase are currently under research and development; these focus other domains of the polymerase than those targeted by the commercially available drugs  <cit> .
one such novel class of compounds is that of 4-oxo-dhqs belonging to the non-nucleoside anti-herpetic drugs family  <cit> .
traditional medicines
medicinal plants products have been used over centuries as traditional remedies for different kinds of diseases including viral diseases.
recently, there have been studies which report anti-viral activities of extracts from plants like swertia chirata, aloe forex and withania somnifera against hsv  <cit> .
these plant extracts inhibit the formation of hsv- <dig> plaque above a certain minimum concentration and their activities can be compared to the commercial drugs like acyclovir.
withania somnifera or winter cherry or indian ginseng is a proud herb of ayurveda, classified as rasayan   <cit> .
it is held in high repute in traditional indian medicine mainly because of its constituents called withanolides  <cit> .
they are built on an ergostane framework, which is oxidized at c- <dig> and c- <dig> to form a six-member lactone ring.
withaferin a , the first withanolide to be isolated and the major withanolide present in indian variety of plants has been widely researched for its pharamacological activities including anti-inflammatory, anti-cancer, anti-stress and immunomodulatory, adaptogenic, central nervous system, endocrine and cardiovascular activities  <cit> .
leaves of withania somnifera have been reported to have the highest content of wa   <cit> .
in this study, we report a possible mode of action of withaferin a against hsv by inhibition of its dna polymerase.
molecular docking studies have been used to identify the binding modes.
dynamic structural patterns were studied using molecular dynamics simulations.
methods
ligand and receptors
the crystal structure of hsv dna polymerase  was obtained from the protein data bank   <cit> .
the structure of this enzyme has been described elsewhere  <cit> .
this structure was then subjected to certain modifications which would make it suitable for docking as described elsewhere  <cit> .
the structure of the ligand molecule withaferin a  retrieved from ncbi-pubchem compound database  <cit>  is shown in figure  <dig>  amber  <dig>  <cit> , was used to minimise the energy of the ligand and the receptor using steepest descent and conjugate gradient methods.
semi-flexible docking
semi-flexible molecular docking of the hsv dna polymerase along with ligand wa was implemented using autodock  <dig>   <cit> .
the general procedure for performing docking is described elsewhere  <cit> .
the outputs from autodock were rendered with pymol  <cit>  and accelerys viewerlite  <dig> .
confirmation of the results were achieved using pardock  <cit> .
md simulations in water
the energy minimization and md simulations of hsv pol and its complex with wa were carried out using amber package as fully described elsewhere by the authors  <cit> .
RESULTS
docking of withaferin a into hsv pol
we performed docking of wa into x-ray crystal structure of hsv pol.
using binding pocket analysis, a cavity around the active site/subsites of hsv pol was predicted where wa was found binding to the critical active site residues of hsv pol with a binding energy of - <dig>  kcal/mol.
the ligand formed several h-bond interactions with the crucial residues of the polymerase .
ring amino hydrogen of his  <dig> formed bond with terminal hydroxyl group of wa.
the same hydroxyl of wa bonded up with amino group of asn  <dig> also.
the other terminal hydroxyl of wa also interacted with side chain carbonyl of gln  <dig>  gln  <dig> is one of the active site residues of hsv pol, being a part of the hydrophobic pocket.
apart from h-bonds, wa was also found forming hydrophobic interactions.
as shown in figure 2b, wa formed van der waals interactions with the highly conserved residues tyr  <dig> and tyr  <dig>  which lie in close proximity to the catalytic site.
the properties of the docked ligand are shown in table  <dig>
md simulations in water
since molecular docking provides only a static view of the protein-ligand interactions, we performed molecular dynamics simulations on hsv pol in complex with wa to study the interactions in motion.
the overall goal of these simulation steps was to account for protein flexibility and movement that however could not be achieved in the docking simulations alone.
figure 3a shows the h-bond interactions of the docked wa with hsv pol in a  <dig> ns simulated snapshot.
wa was found forming h-bond interactions with gln  <dig>  phe  <dig>  asn  <dig> and asp  <dig> of hsv pol.
strong hydrophobic interactions by the residues phe  <dig>  pro  <dig>  tyr  <dig> and tyr  <dig> of hsv pol were also observed in this simulated structure .
to probe the dynamic flexibility changes in the protein, due to the inhibitor binding in its critical regions, b-factors  for the cα atoms were calculated using the following relation:
bn = 8/ <dig> π rn2
where rn is the root-mean-square atomic fluctuation of the cα atom of residue n. as shown in figure 4a, the calculated bn values for the undocked protein show a major peak for the residues 615- <dig> close to the active site that possibly determines its function through large-amplitude motion.
however, on docking with wa, a decrease in bn values of the protein was observed that shows its conformational rigidity.
further, this points out that the large-amplitude motion of the flexible residues 615- <dig>  which encompasses two consecutive catalytically-active glutamines  <dig> and  <dig>  is retarded by the binding of wa to active site of hsv pol.
this suggests the inhibitory activity of wa against hsv pol, as in our docked structures wa established hydrogen bonds with these residues.
in another instance, for the region 800- <dig> in which most of the catalytic site residues are located, bn values for wa bound polymerase were seen attenuated in comparison to those of the enzyme alone.
the rmsds for the trajectory of wa bound hsv pol reveals the stability of the complex .
the rmsd of the complex has achieved a stationary phase during the later stage of the simulation and is always less than  <dig>  Å for the entire simulation length.
the thermodynamic stability of the complex was evident as the energy of the complex  was found always lower than that of the protein  throughout the simulation .
analysis of pre- and post-md simulated structures
from the interactions shown in figure 5a and 5b as ligplots  <cit> , a significant difference in the number and nature of interactions was observed.
the binding energy based on the physico-chemical properties of the active site as well as that of the ligand were calculated as - <dig>  and - <dig>  kcal/mol for the post- and pre-md simulated structures respectively.
the comparative parametric values of the two structures calculated using raspd server  <cit>  are shown n table  <dig>  close visual inspection of the md results permitted us to observe that post-md, wa was able to establish a better structurally stable conformation.
as evident from figure 5c, in the pre-md structure wa was not able to anchor appropriately its tail inside the cavity of hsv pol.
on the other hand, the behaviour of wa was much more stable in the post-md simulated structure where its tail is well anchored inside the deep narrow gorge of the polymerase.
moreover, the ligand substantially moved towards the residues of hsv pol to form additional h-bonds.
the residues phe  <dig> and asp  <dig> with which wa established h-bonds in the post-simulated structure participate in metal-ion coordination essential for the catalysis of polymerase.
another wa h-bonding residue asn  <dig> of hsv pol, is a highly conserved residue which besets on the prolonged second α helix of the finger domain.
the continuous nature of this structure is responsible for bringing the other highly conserved residues tyr  <dig> and gly  <dig> into close proximity of the catalytic site for facilitating proper functioning of the polymerase  <cit> .
in the post-md simulated structure, the intensity of hydrophobic interactions was also found enhanced as compared to those present before md.
the residues phe  <dig>  tyr  <dig>  pro  <dig> and tyr  <dig> provided significant hydrophobic contacts for stabilization of wa inside the cavity.
the ribose of the incoming nucleotide has been earlier report to interact with the conserved tyr  <dig> residue and induces a steric access effect against inclusion of the nucleotides  <cit> .
non-nucleosidic inhibition of hsv pol by wa
it has been earlier reported that 4-oxo-dihydroquinolines  have shown broad anti-hsv activity  <cit> .
this class of compounds inhibit most human herpes viruses, which is associated with dna polymerase inhibition.
from the binding assays of pnu- <dig> , a radiolabelled dhq on hsv pol, it was observed that this inhibitor binds only to hsv pol in complex with the dna duplex; while no binding was observed either with hsv pol or with dna duplex alone  <cit> .
however, in our case we observed that wa is able to bind to the segregated hsv pol itself.
from visual inspection of the structures of wa and the 4-oxo-dhq inhibitor, it was observed that both these ligands contain  <dig> ketonic groups, an oxygen containing heterocyclic ring and a tail containing an electronegative atom.
in the present study, the binding mode of wa was found similar to that of the radiolabelled 4-oxo-dhq in which the residues gln  <dig>  gln  <dig>  asn  <dig> and tyr  <dig> play critical role in the stabilization of the ligand.
it can be deduced from our studies that wa can be a potent non-nucleosidic inhibitor of hsv pol whose binding would result in a conformational change in the polymerase that distorts the positioning of the residues that bind dna, inhibiting polymerization.
4-oxo-dhqs have shown high specificity index in inhibiting dna polymerases belonging to the herpesviridae family because unrelated dna and rna viruses were not susceptible to their inhibitory effect, and they also proved to have broad spectrum antiviral effects  <cit> .
the same can be expected for wa owing to its analogous mode of action.
the inhibition constant of wa found in our study  was also quite comparable to the currently established nucleosidic drugs   <cit> .
the emergence of hsv resistance to antiviral drugs is also a major concern.
three basic mechanisms have been identified: altered thymidine kinase substrate specificity, absent or partial production of viral thymidine kinase and altered viral dna polymerase  <cit> .
the most common mechanism found in clinical isolates is deficient tk activity.
for the foremost approved drug- acyclovir, resistant isolates of hsv have been observed in immuno-compromised individuals, especially aids patients  <cit> .
since wa is observed exerting its inhibitory effect via interaction with a viral dna polymerase site that is less important for the binding of deoxynucleoside triphosphates, it holds potential to exert its influence even on these resistant isolates.
the non-nucleosidic mode of action of wa holds promise for prevention of infection, as it can selectively target only the viral enzymes.
moreover being a naturally occurring herbal drug candidate, wa will also be able to address the issue of safety and bioavailability.
based on the results from wa-hsv pol complex, it appears that interactions with the residues gln  <dig>  gln  <dig>  tyr  <dig>  asn  <dig> and asp  <dig> of hsv pol are important for inhibitory activity of wa.
a comparison between the conformation obtained from docking and that from molecular dynamics simulations show that substantial changes in binding conformations have occurred.
these results indicate that the initial receptor-ligand interaction observed after docking can be limited due to the receptor rigid docking algorithm and that the conformations and interactions observed after simulation runs are more energetically favoured and should be better representations of derivative poses in receptor.
CONCLUSIONS
we have performed docking and molecular dynamics simulation studies to elucidate the binding mechanism of prospective herbal drug withaferin a onto the structure of dna polymerase of herpes simplex virus.
our docking simulations results give high binding affinity of the ligand to the receptor.
long de novo md simulations for  <dig> ns performed allowed us to evaluate the dynamic behaviour of the system studied and corroborate the docking results, as well as identify key residues in the enzyme-inhibitor interactions.
the present md simulations support the hypothesis that wa is a prospective ligand that has potential to target/inhibit dna polymerase of the herpes simplex virus.
results of these studies will also guide the design of selective inhibitors of dna pol with high specificity and potent activity in order to strengthen the therapeutic arsenal available today against the dangerous biological warfare agent represented by herpes simplex virus.
authors' contributions
ag, vsb and ds designed the methods and experimental setup.
ag carried out the implementation of the various methods.
as and va assisted ag in this process.
ag and ds wrote the manuscript.
all authors have read and approved the final manuscript.
competing interests
the authors declare that they have no competing interests

 herpes simplex virus type and type are two members of herpesviridae family, which infect almost 85% of the world population. they are the causative agents of a gamut of diseases ranging from mild ones like cold sores in mouth, eye cornea and genitals to more severe life threatening ones like the fatal herpes encephalitis. people with suppressed immune system like those suffering from aids are more prone to get infections from hsv than others. drugs for hsv there is no permanent cure for these infections till date. present day treatment involves the use of antiviral drugs to reduce the physical severity of outbreak-associated lesions and viral shedding, though this helps decreasing the chances of transmission to others only by maximum 50%. there are two types of drugs that are clinically useful against hsv infections. the first category consists of nucleoside analogs like acyclovir and its prodrug valacyclovir, ganciclovir, penciclovir and its prodrug famciclovir, sorivudine and brivudine. these require phosphorylation by viral thymidine kinase to form triphosphates that are active inhibitors of viral dna polymerase. the second category consists of direct viral dna polymerase inhibitors like vidarbine, foscarnet and cidofovir. thus, both types of drugs target in dysfunctioning the replication centre i.e. dna polymerase of the viral genome. development of alternative treatments for hsv however, in past few years, a number of acyclovir drug resistant viral strains have been isolated especially from immuno-compromised patients. in this era, where the number of immuno-suppressed patients like those suffering from hiv is continuously increasing, there is an immediate need to find new drugs to treat hsv infections which have a higher efficacy or have an alternative mode of action. resistance to acyclovir is mainly due to mutations in the viral thymidine kinase gene which impair the initial drug phosphorylation. these drug resistant strains have been of significant clinical attention, indicating the need for alternative anti- hsv drugs. previous antiviral studies on herpes viruses have focussed on developing nucleoside analogues that will inhibit viral polymerase and terminate the replicating viral dna. a number of new anti-viral drugs against hsv dna polymerase are currently under research and development; these focus other domains of the polymerase than those targeted by the commercially available drugs. one such novel class of compounds is that of 4-oxo-dhqs belonging to the non-nucleoside anti-herpetic drugs family. traditional medicines medicinal plants products have been used over centuries as traditional remedies for different kinds of diseases including viral diseases. recently, there have been studies which report anti-viral activities of extracts from plants like swertia chirata, aloe forex and withania somnifera against hsv. these plant extracts inhibit the formation of hsv- plaque above a certain minimum concentration and their activities can be compared to the commercial drugs like acyclovir. withania somnifera or winter cherry or indian ginseng is a proud herb of ayurveda, classified as rasayan. it is held in high repute in traditional indian medicine mainly because of its constituents called withanolides. they are built on an ergostane framework, which is oxidized at c- and c- to form a six-member lactone ring. withaferin a, the first withanolide to be isolated and the major withanolide present in indian variety of plants has been widely researched for its pharamacological activities including anti-inflammatory, anti-cancer, anti-stress and immunomodulatory, adaptogenic, central nervous system, endocrine and cardiovascular activities. leaves of withania somnifera have been reported to have the highest content of wa. in this study, we report a possible mode of action of withaferin a against hsv by inhibition of its dna polymerase. molecular docking studies have been used to identify the binding modes. dynamic structural patterns were studied using molecular dynamics simulations. methods ligand and receptors the crystal structure of hsv dna polymerase was obtained from the protein data bank. the structure of this enzyme has been described elsewhere. this structure was then subjected to certain modifications which would make it suitable for docking as described elsewhere. the structure of the ligand molecule withaferin a retrieved from ncbi-pubchem compound database is shown in figure amber, was used to minimise the energy of the ligand and the receptor using steepest descent and conjugate gradient methods. semi-flexible docking semi-flexible molecular docking of the hsv dna polymerase along with ligand wa was implemented using autodock. the general procedure for performing docking is described elsewhere. the outputs from autodock were rendered with pymol and accelerys viewerlite. confirmation of the results were achieved using pardock. md simulations in water the energy minimization and md simulations of hsv pol and its complex with wa were carried out using amber package as fully described elsewhere by the authors.  docking of withaferin a into hsv pol we performed docking of wa into x-ray crystal structure of hsv pol. using binding pocket analysis, a cavity around the active site/subsites of hsv pol was predicted where wa was found binding to the critical active site residues of hsv pol with a binding energy of - kcal/mol. the ligand formed several h-bond interactions with the crucial residues of the polymerase. ring amino hydrogen of his formed bond with terminal hydroxyl group of wa. the same hydroxyl of wa bonded up with amino group of asn also. the other terminal hydroxyl of wa also interacted with side chain carbonyl of gln gln is one of the active site residues of hsv pol, being a part of the hydrophobic pocket. apart from h-bonds, wa was also found forming hydrophobic interactions. as shown in figure 2b, wa formed van der waals interactions with the highly conserved residues tyr and tyr which lie in close proximity to the catalytic site. the properties of the docked ligand are shown in table md simulations in water since molecular docking provides only a static view of the protein-ligand interactions, we performed molecular dynamics simulations on hsv pol in complex with wa to study the interactions in motion. the overall goal of these simulation steps was to account for protein flexibility and movement that however could not be achieved in the docking simulations alone. figure 3a shows the h-bond interactions of the docked wa with hsv pol in a ns simulated snapshot. wa was found forming h-bond interactions with gln phe asn and asp of hsv pol. strong hydrophobic interactions by the residues phe pro tyr and tyr of hsv pol were also observed in this simulated structure. to probe the dynamic flexibility changes in the protein, due to the inhibitor binding in its critical regions, b-factors for the cα atoms were calculated using the following relation: bn = 8/ π rn2 where rn is the root-mean-square atomic fluctuation of the cα atom of residue n. as shown in figure 4a, the calculated bn values for the undocked protein show a major peak for the residues 615- close to the active site that possibly determines its function through large-amplitude motion. however, on docking with wa, a decrease in bn values of the protein was observed that shows its conformational rigidity. further, this points out that the large-amplitude motion of the flexible residues 615- which encompasses two consecutive catalytically-active glutamines and is retarded by the binding of wa to active site of hsv pol. this suggests the inhibitory activity of wa against hsv pol, as in our docked structures wa established hydrogen bonds with these residues. in another instance, for the region 800- in which most of the catalytic site residues are located, bn values for wa bound polymerase were seen attenuated in comparison to those of the enzyme alone. the rmsds for the trajectory of wa bound hsv pol reveals the stability of the complex. the rmsd of the complex has achieved a stationary phase during the later stage of the simulation and is always less than Å for the entire simulation length. the thermodynamic stability of the complex was evident as the energy of the complex was found always lower than that of the protein throughout the simulation. analysis of pre- and post-md simulated structures from the interactions shown in figure 5a and 5b as ligplots, a significant difference in the number and nature of interactions was observed. the binding energy based on the physico-chemical properties of the active site as well as that of the ligand were calculated as - and - kcal/mol for the post- and pre-md simulated structures respectively. the comparative parametric values of the two structures calculated using raspd server are shown n table close visual inspection of the md results permitted us to observe that post-md, wa was able to establish a better structurally stable conformation. as evident from figure 5c, in the pre-md structure wa was not able to anchor appropriately its tail inside the cavity of hsv pol. on the other hand, the behaviour of wa was much more stable in the post-md simulated structure where its tail is well anchored inside the deep narrow gorge of the polymerase. moreover, the ligand substantially moved towards the residues of hsv pol to form additional h-bonds. the residues phe and asp with which wa established h-bonds in the post-simulated structure participate in metal-ion coordination essential for the catalysis of polymerase. another wa h-bonding residue asn of hsv pol, is a highly conserved residue which besets on the prolonged second α helix of the finger domain. the continuous nature of this structure is responsible for bringing the other highly conserved residues tyr and gly into close proximity of the catalytic site for facilitating proper functioning of the polymerase. in the post-md simulated structure, the intensity of hydrophobic interactions was also found enhanced as compared to those present before md. the residues phe tyr pro and tyr provided significant hydrophobic contacts for stabilization of wa inside the cavity. the ribose of the incoming nucleotide has been earlier report to interact with the conserved tyr residue and induces a steric access effect against inclusion of the nucleotides. non-nucleosidic inhibition of hsv pol by wa it has been earlier reported that 4-oxo-dihydroquinolines have shown broad anti-hsv activity. this class of compounds inhibit most human herpes viruses, which is associated with dna polymerase inhibition. from the binding assays of pnu-, a radiolabelled dhq on hsv pol, it was observed that this inhibitor binds only to hsv pol in complex with the dna duplex; while no binding was observed either with hsv pol or with dna duplex alone. however, in our case we observed that wa is able to bind to the segregated hsv pol itself. from visual inspection of the structures of wa and the 4-oxo-dhq inhibitor, it was observed that both these ligands contain ketonic groups, an oxygen containing heterocyclic ring and a tail containing an electronegative atom. in the present study, the binding mode of wa was found similar to that of the radiolabelled 4-oxo-dhq in which the residues gln gln asn and tyr play critical role in the stabilization of the ligand. it can be deduced from our studies that wa can be a potent non-nucleosidic inhibitor of hsv pol whose binding would result in a conformational change in the polymerase that distorts the positioning of the residues that bind dna, inhibiting polymerization. 4-oxo-dhqs have shown high specificity index in inhibiting dna polymerases belonging to the herpesviridae family because unrelated dna and rna viruses were not susceptible to their inhibitory effect, and they also proved to have broad spectrum antiviral effects. the same can be expected for wa owing to its analogous mode of action. the inhibition constant of wa found in our study was also quite comparable to the currently established nucleosidic drugs. the emergence of hsv resistance to antiviral drugs is also a major concern. three basic mechanisms have been identified: altered thymidine kinase substrate specificity, absent or partial production of viral thymidine kinase and altered viral dna polymerase. the most common mechanism found in clinical isolates is deficient tk activity. for the foremost approved drug- acyclovir, resistant isolates of hsv have been observed in immuno-compromised individuals, especially aids patients. since wa is observed exerting its inhibitory effect via interaction with a viral dna polymerase site that is less important for the binding of deoxynucleoside triphosphates, it holds potential to exert its influence even on these resistant isolates. the non-nucleosidic mode of action of wa holds promise for prevention of infection, as it can selectively target only the viral enzymes. moreover being a naturally occurring herbal drug candidate, wa will also be able to address the issue of safety and bioavailability. based on the results from wa-hsv pol complex, it appears that interactions with the residues gln gln tyr asn and asp of hsv pol are important for inhibitory activity of wa. a comparison between the conformation obtained from docking and that from molecular dynamics simulations show that substantial changes in binding conformations have occurred. these results indicate that the initial receptor-ligand interaction observed after docking can be limited due to the receptor rigid docking algorithm and that the conformations and interactions observed after simulation runs are more energetically favoured and should be better representations of derivative poses in receptor.  we have performed docking and molecular dynamics simulation studies to elucidate the binding mechanism of prospective herbal drug withaferin a onto the structure of dna polymerase of herpes simplex virus. our docking simulations results give high binding affinity of the ligand to the receptor. long de novo md simulations for ns performed allowed us to evaluate the dynamic behaviour of the system studied and corroborate the docking results, as well as identify key residues in the enzyme-inhibitor interactions. the present md simulations support the hypothesis that wa is a prospective ligand that has potential to target/inhibit dna polymerase of the herpes simplex virus. results of these studies will also guide the design of selective inhibitors of dna pol with high specificity and potent activity in order to strengthen the therapeutic arsenal available today against the dangerous biological warfare agent represented by herpes simplex virus. authors' contributions ag, vsb and ds designed the methods and experimental setup. ag carried out the implementation of the various methods. as and va assisted ag in this process. ag and ds wrote the manuscript. all authors have read and approved the final manuscript. competing interests the authors declare that they have no competing interests 
long de novo md simulations for ns performed allowed us to evaluate the dynamic behaviour of the system studied and corroborate the docking results, as well as identify key residues in the enzyme-inhibitor interactions. our docking simulations results give high binding affinity of the ligand to the receptor. we have performed docking and molecular dynamics simulation studies to elucidate the binding mechanism of prospective herbal drug withaferin a onto the structure of dna polymerase of herpes simplex virus. these results indicate that the initial receptor-ligand interaction observed after docking can be limited due to the receptor rigid docking algorithm and that the conformations and interactions observed after simulation runs are more energetically favoured. results of these studies will also guide the design of selective inhibitors of dna pol with high specificity and potent activity in order to strengthen the therapeutic arsenal available today against the dangerous biological warfare agent represented by herpes simplex virus. withaferin a was found to bind very similarly to that of the previously reported 4-oxo-dhq inhibitor. the present md simulations support the hypothesis that withaferin a is a potential ligand to target/inhibit dna polymerase of the herpes simplex virus. previous antiviral studies on herpes viruses have focussed on developing nucleoside analogues that can inhibit viral polymerase and terminate the replicating viral dna. the present study is an attempt to elucidate the action mechanism of naturally occurring withaferin a in inhibiting viral dna polymerase, thus providing an evidence for its development as a novel anti-herpetic drug. herpes simplex virus and causes several infections in humans including cold sores and encephalitis. a comparison of the conformation obtained from docking and the molecular dynamics simulations shows that substantial changes in the binding conformations have occurred.
 herpes simplex virus and causes several infections in humans including cold sores and encephalitis. previous antiviral studies on herpes viruses have focussed on developing nucleoside analogues that can inhibit viral polymerase and terminate the replicating viral dna. however, these drugs bear an intrinsic non-specificity as they can also inhibit cellular polymerase apart from the viral one. the present study is an attempt to elucidate the action mechanism of naturally occurring withaferin a in inhibiting viral dna polymerase, thus providing an evidence for its development as a novel anti-herpetic drug.  withaferin a was found to bind very similarly to that of the previously reported 4-oxo-dhq inhibitor. withaferin a was observed binding to the residues gln gln asn and tyr all of which are crucial to the proper functioning of the polymerase. a comparison of the conformation obtained from docking and the molecular dynamics simulations shows that substantial changes in the binding conformations have occurred. these results indicate that the initial receptor-ligand interaction observed after docking can be limited due to the receptor rigid docking algorithm and that the conformations and interactions observed after simulation runs are more energetically favoured.  we have performed docking and molecular dynamics simulation studies to elucidate the binding mechanism of prospective herbal drug withaferin a onto the structure of dna polymerase of herpes simplex virus. our docking simulations results give high binding affinity of the ligand to the receptor. long de novo md simulations for ns performed allowed us to evaluate the dynamic behaviour of the system studied and corroborate the docking results, as well as identify key residues in the enzyme-inhibitor interactions. the present md simulations support the hypothesis that withaferin a is a potential ligand to target/inhibit dna polymerase of the herpes simplex virus. results of these studies will also guide the design of selective inhibitors of dna pol with high specificity and potent activity in order to strengthen the therapeutic arsenal available today against the dangerous biological warfare agent represented by herpes simplex virus. november- december asia pacific bioinformatics network tenth international conference on bioinformatics – first iscb asia joint conference kuala lumpur, malaysia 