Herpes simplex virus type 1 (HSV-1) infections affect about two-thirds of the world population, and the standard treatment consists of acyclovir (ACV) and its analogs, which interact with thymidine kinase (TK) blocking viral replication. Lately, the emergence of ACV-resistant strains has been reported, especially associated with TK mutations. In this context, ACV therapy fails against isolates encoding Y172C and Y53H/R163H TK mutants, but the molecular mechanism of drug resistance remains unclear. Thus, we examined the effects of these mutations on ACV and the cofactor ATP binding through molecular modeling approaches. We showed that Y172C prevents the anchoring of the aromatic ring of ACV through π-π stacking interactions, leading to an inversed binding mode and different interactions. On the other hand, Y53H/R163H remarkably affected the cofactor binding mode which shifted away from its binding site and also influenced the interaction network of ACV. This is likely due to the loss of polar interactions with R163 residue. Unlike what was observed in the wild-type complex, both drug and cofactor binding poses were not well positioned to allow the phosphorylation reaction which explains the resistance observed. Moreover, energy analysis corroborated the experimental data and showed lower theoretical affinity of ACV with mutant enzymes resulted from energetic loss in polar solvation in Y172C and electrostatic terms in Y53H/R163H mutant enzyme. Therefore, our study shed light on the resistance mechanism toward ACV of two mutant TKs identified in clinical HSV-1 strains and may further support the development of new anti-herpetic drugs to treat resistant infections. [Formula: see text] Communicated by Ramaswamy H. Sarma.
Keywords: Herpes; acyclovir; erpes simplex virus; molecular modeling; resistance.