Hepatitis C virus genotype 1a (HCV-1a) comprises clades I and II. The Q80K polymorphism is found predominantly in clade I but rarely in clade II. Here, we investigated whether natural polymorphisms in HCV-1a clade II entailed structural protein changes when occurrence of the Q80K variant was simulated. Based on HCV-1a clade I and II protein sequences, the structure of the HCV-1a Q80K mutant NS3-4A was obtained by comparative modeling. Its physicochemical properties were studied by molecular dynamics simulations and network analysis. Results demonstrate that, in the presence of the K80 variant, clade II protease polymorphisms A91 and S/G174 led to variations in hydrogen bond occupancies. Structural analyses revealed differences in (i) flexibility of the H57 catalytic residue on the NS3 protease and (ii) correlations between amino acids on the NS3 protease and the NS4A cofactor. The latter indicated possible destabilization of interactions, resulting in increased separation of these proteins. The present findings describe how the relationships between different HCV-1a NS3 protease amino acid residues could affect the appearance of viral variants and the existence of distinct genetic barriers to HCV-1a isolates.
Keywords: NS3 protease; Q80K variant; hepatitis C virus; molecular dynamics; network analysis; polymorphisms; simulation.