The hepatitis C virus (HCV) is a major global health challenge and a leading cause of morbidity and mortality. Many direct-acting antivirals (DAAs) target essential macromolecules involved in the virus' life cycle. Although such DAAs achieve great success in reducing the viral load in genotype 1 infections, other genotypes demonstrate different levels of response. This study focused on mutation sites associated with patients with genotype 4a infections that failed to respond to treatment with sofosbuvir. The genotyping of HCV samples from patients with virological failure, and responder patients, was conducted using Geno2Pheno webserver-based full NS5B sequences. We constructed 3D structural models for all the samples and used structural analysis to investigate the effect of amino acid substitution on the observed resistance to SOF-based treatment, and the docking of sofosbuvir into the active sites of the 10 models was performed. Finally, 10 molecular dynamic (MD) simulation experiments were conducted to compare the stability of the 3D models of the resistant samples against the stability of the 3D models of the responder samples. The results highlighted the presence of HCV subtype 4a in all ten samples; in addition, an amino acid (aa) substitution in the palm region may hinder HCV polymerase activity. In this study, we provide evidence that a mutation in the NS5B gene that induces resistance to sofosbuvir in patients with the S282T/C/R mutant virus is present in the Egyptian population. Overall, the docking and MD results support our findings and highlight the significant impact of the identified mutations on the resistance of HCV NS5B RNA-dependent RNA polymerase to direct-acting antivirals (DAAs).
Keywords: HCV; RNA sequencing; RNA-dependent RNA polymerase; molecular dynamics; resistance.