In silico mechanistic analysis of IRF3 inactivation and high-risk HPV E6 species-dependent drug response

The high-risk human papillomavirus E6 (hrHPV E6) protein has been widely studied due to its implication in cervical cancer. In response to viral threat, activated kinases phosphorylate the IRF3 autoinhibitory domain, inducing type1 interferon production. HPV circumvents the antiviral response through the possible E6 interaction with IRF3 and abrogates p53’s apoptotic activity by recruiting E6-associated protein. However, the molecular mechanism of IRF3 inactivation by hrHPV E6 has not yet been delineated. Therefore, we explored this mechanism through in silico examination of protein-protein and protein-ligand docking, binding energy differences and computational alanine mutagenesis. Our results suggested that the LxxLL motifs of IRF3 binds within the hydrophobic pocket of E6, precluding Ser-patch phosphorylation, necessary for IRF3 activation and interferon induction. This model was further supported by molecular dynamics simulation. Furthermore, protein-ligand docking and drug resistance modeling revealed that the polar patches in the pocket of E6, which are crucial for complex stability and ligand binding, are inconsistent among hrHPV species. Such variabilities pose a risk of treatment failure owing to point mutations that might render drugs ineffective and allude to multi-drug therapy. Overall, this study reveals a novel perspective of innate immune suppression in HPV infections and suggests a plausible therapeutic intervention.