A 3D structural model and dynamics of hepatitis C virus NS3/4A protease (genotype 4a, strain ED43) suggest conformational instability of the catalytic triad: implications in catalysis and drug resistivity

Hepatitis C Virus NS3/4A Inhibition and Host Immunomodulation by Tannins from Terminalia chebula: A Structural Perspective

Terminalia chebula Retz. forms a key component of traditional folk medicine and is also reported to possess antihepatitis C virus (HCV) and immunomodulatory activities. However, information on the intermolecular interactions of phytochemicals from this plant with HCV and human proteins are yet to be established. Thus, by this current study, we investigated the HCV NS3/4A inhibitory and host immune-modulatory activity of phytocompounds from T. chebula through in silico strategies involving network pharmacology and structural bioinformatics techniques. To start with, the phytochemical dataset of T. chebula was curated from biological databases and the published literature. Further, the target ability of the phytocompounds was predicted using BindingDB for both HCV NS3/4A and other probable host targets involved in the immune system. Further, the identified targets were docked to the phytochemical dataset using AutoDock Vina executed through the POAP pipeline. The resultant docked complexes with significant binding energy were subjected to 50 ns molecular dynamics (MD) simulation in order to infer the stability of complex formation. During network pharmacology analysis, the gene set pathway enrichment of host targets was performed using the STRING and Reactome pathway databases. Further, the biological network among compounds, proteins, and pathways was constructed using Cytoscape 3.6.1. Furthermore, the druglikeness, side effects, and toxicity of the phytocompounds were also predicted using the MolSoft, ADVERpred, and PreADMET methods, respectively. Out of 41 selected compounds, 10 were predicted to target HCV NS3/4A and also to possess druglike and nontoxic properties. Among these 10 molecules, Chebulagic acid and 1,2,3,4,6-Pentagalloyl glucose exhibited potent HCV NS3/4A inhibitory activity, as these scored a lowest binding energy (BE) of −8.6 kcal/mol and −7.7 kcal/mol with 11 and 20 intermolecular interactions with active site residues, respectively. These findings are highly comparable with Asunaprevir (known inhibitor of HCV NS3/4A), which scored a BE of −7.4 kcal/mol with 20 key intermolecular interactions. MD studies also strongly suggest that chebulagic acid and 1,2,3,4,6-Pentagalloyl glucose as promising leads, as these molecules showed stable binding during 50 ns of production run. Further, the gene set enrichment and network analysis of 18 protein targets prioritized 10 compounds and were predicted to potentially modulate the host immune system, hemostasis, cytokine levels, interleukins signaling pathways, and platelet aggregation. On overall analysis, this present study predicts that tannins from T. chebula have a potential HCV NS3/4A inhibitory and host immune-modulatory activity. However, further experimental studies are required to confirm the efficacies.

Toward the Identification of Potential α-Ketoamide Covalent Inhibitors for SARS-CoV-2 Main Protease: Fragment-Based Drug Design and MM-PBSA Calculations

Since December 2019, the world has been facing the outbreak of the SARS-CoV-2 pandemic that has infected more than 149 million and killed 3.1 million people by 27 April 2021, according to WHO statistics. Safety measures and precautions taken by many countries seem insufficient, especially with no specific approved drugs against the virus. This has created an urgent need to fast track the development of new medication against the virus in order to alleviate the problem and meet public expectations. The SARS-CoV-2 3CL main protease (Mpro) is one of the most attractive targets in the virus life cycle, which is responsible for the processing of the viral polyprotein and is a key for the ribosomal translation of the SARS-CoV-2 genome. In this work, we targeted this enzyme through a structure-based drug design (SBDD) protocol, which aimed at the design of a new potential inhibitor for Mpro. The protocol involves three major steps: fragment-based drug design (FBDD), covalent docking and molecular dynamics (MD) simulation with the calculation of the designed molecule binding free energy at a high level of theory. The FBDD step identified five molecular fragments, which were linked via a suitable carbon linker, to construct our designed compound RMH148. The mode of binding and initial interactions between RMH148 and the enzyme active site was established in the second step of our protocol via covalent docking. The final step involved the use of MD simulations to test for the stability of the docked RMH148 into the Mpro active site and included precise calculations for potential interactions with active site residues and binding free energies. The results introduced RMH148 as a potential inhibitor for the SARS-CoV-2 Mpro enzyme, which was able to achieve various interactions with the enzyme and forms a highly stable complex at the active site even better than the co-crystalized reference.

Investigating into the molecular interactions of flavonoids targeting NS2B-NS3 protease from ZIKA virus through in-silico approaches

Zika virus (ZIKV), belongs to the flavivirus genus and Flaviviridae family that associated with serious diseased conditions like microcephaly and other neurological disorders (Guillan-Barré syndrome). As there is no vaccine or therapies available against ZIKV to date. Hence, it is an unmet need to find potential drug candidates and target sites against Zika virus infection. NS2B-NS3 protease making an attractive target for therapeutic intervention in ZIKV infections because of its critical role in hydrolysis of a single polyprotein encoded by Zika virus. Recently, there are some experimental evidence about the flavonoids as Zika virus NS2B-NS3 protease inhibitors. However, molecular interaction between protease complex and inhibitors at atomic levels has not been explored. Here, we have taken the experimentally validated thirty-eight flavonoids inhibitors against NS2B-NS3 protease to examine the molecular interaction using molecular docking and molecular dynamics simulations. We found out few flavonoids such as EGCG and its two derivatives, isoquercetin, rutin and sanggenon O showing interaction with catalytic triad (His51, Asp75, and Ser135) of the active site of NS2B-NS3 protease and found to be stable throughout the simulation. Therefore it is evident that interaction with the catalytic triad playing a vital role in the inhibition of the enzyme activity as a result inhibition of the virus propagation. However these compounds can be explored further for understanding the mechanism of action of these compounds targeting NS2B-NS3 protease for inhibition of Zika virus.

A novel molecular mechanism to explain mutations of the HCV protease associated with resistance against covalently bound inhibitors

NS3 is an important therapeutic target for direct-acting antiviral (DAA) drugs. However, many patients treated with DAAs have unsustained virologic response (UVR) due to the high mutation rate of HCV. The aim of this work was to shed some light on the puzzling molecular mechanisms of the virus's of patients who showed high viral loads even under treatment with DAA. Bioinformatics tools, molecular modelling analyses were employed to identify mutations associated with HCV resistance to boceprevir and possible structural features related to this phenomenon. We identified two mutations of NS3 that may be associated with HCV resistance: D168N and L153I. The substitution D168N was previously reported in the literature as related with drug failure. Additionally, we identified that its molecular resistance mechanism can be explained by the destabilization of receptor-ligand hydrogen bonds. For the L153I mutation, the resistance mechanism is different from previous models reported in the literature. The L153I substitution decreases the S139 deprotonation susceptibility, and consequently, this mutation impairs the covalent binding between the residue S139 from NS3 and the electrophilic trap on boceprevir, which can induce drug failure. These results were supported by the time course analysis of the mutations of the NS3 protease, which showed that boceprevir was designed for enzymes with an L residue at position 153; however, the sequences with I153 are predominant nowadays. The results presented here could be used to infer about resistance in others DAA, mainly protease inhibitors.

The NS4A Cofactor Dependent Enhancement of HCV NS3 Protease Activity Correlates with a 4D Geometrical Measure of the Catalytic Triad Region

We are developing a 4D computational methodology, based on 3D structure modeling and molecular dynamics simulation, to analyze the active site of HCV NS3 proteases, in relation to their catalytic activity. In our previous work, the 4D analyses of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) yielded divergent, gradual and genotype-dependent, 4D conformational instability measures, which strongly correlate with the known disparate catalytic activities among genotypes. Here, the correlation of our 4D geometrical measure is extended to intra-genotypic alterations in NS3 protease activity, due to sequence variations in the NS4A activating cofactor. The correlation between the 4D measure and the enzymatic activity is qualitatively evident, which further validates our methodology, leading to the development of an accurate quantitative metric to predict protease activity in silico. The results suggest plausible “communication” pathways for conformational propagation from the activation subunit (the NS4A cofactor binding site) to the catalytic subunit (the catalytic triad). The results also strongly suggest that the well-sampled (via convergence quantification) structural dynamics are more connected to the divergent catalytic activity observed in HCV NS3 proteases than to rigid structures. The method could also be applicable to predict patients’ responses to interferon therapy and better understand the innate interferon activation pathway.

Hepatitis C Virus Genotype 1 to 6 Protease Inhibitor Escape Variants: In Vitro Selection, Fitness, and Resistance Patterns in the Context of the Infectious Viral Life Cycle

Hepatitis C virus (HCV) NS3 protease inhibitors (PIs) are important components of novel HCV therapy regimens. Studies of PI resistance initially focused on genotype 1. Therefore, knowledge about the determinants of PI resistance for the highly prevalent genotypes 2 to 6 remains limited. Using Huh7.5 cell culture-infectious HCV recombinants with genotype 1 to 6 NS3 protease, we identified protease positions 54, 155, and 156 as hot spots for the selection of resistance substitutions under treatment with the first licensed PIs, telaprevir and boceprevir. Treatment of a genotype 2 isolate with the newer PIs vaniprevir, faldaprevir, simeprevir, grazoprevir, paritaprevir, and deldeprevir identified positions 156 and 168 as hot spots for resistance; the Y56H substitution emerged for three newer PIs. Substitution selection also depended on the specific recombinant. The substitutions identified conferred cross-resistance to several PIs; however, most substitutions selected under telaprevir or boceprevir treatment conferred less resistance to certain newer PIs. In a single-cycle production assay, across genotypes, PI treatment primarily decreased viral replication, which was rescued by PI resistance substitutions. The substitutions identified resulted in differential effects on viral fitness, depending on the original recombinant and the substitution. Across genotypes, fitness impairment induced by resistance substitutions was due primarily to decreased replication. Most combinations of substitutions that were identified increased resistance or fitness. Combinations of resistance substitutions with fitness-compensating substitutions either rescued replication or compensated for decreased replication by increasing assembly. This comprehensive study provides insight into the selection patterns and effects of PI resistance substitutions for HCV genotypes 1 to 6 in the context of the infectious viral life cycle, which is of interest for clinical and virological HCV research.

Contrasting effect of new HCMV pUL54 mutations on antiviral drug susceptibility: Benefits and limits of 3D analysis

Human cytomegalovirus (HCMV) resistance to antiviral drugs is a major drawback of repeated or long-duration treatment in immunocompromised patients. Resistance testing is usually performed by genotypic assays. For accurate interpretation of these assays, the role of new mutations in HCMV resistance has to be assessed. Two previously unknown UL54 single point mutations (D515Y and V787A) were characterized for phenotypic drug-resistance by marker transfer analysis using bacterial artificial chromosome (BAC) mutagenesis. Increases in 50% inhibitory concentrations of ganciclovir and foscarnet were found for both mutated recombinant strains showing that mutations D515Y and V787A induce resistance to both antivirals. Importantly, none of those impacted the viral growth kinetics. For a better understanding of their molecular resistance mechanisms, a 3D homology model was used to localize the mutated amino-acids in functional domains of UL54 and predict their impact on UL54 function and resistance. However, 3D homology model analysis has limits and phenotypic characterization using BAC-HCMV is still essential to measure the role of unknown mutations.

Comparative molecular dynamics simulation of Hepatitis C Virus NS3/4A protease (Genotypes 1b, 3a and 4b) predicts conformational instability of the catalytic triad in drug resistant strains

The protease domain of the Hepatitis C Virus (HCV) nonstructural protein 3 (NS3) has been targeted for inhibition by several direct-acting antiviral drugs. This approach has had marked success to treat infections caused by HCV genotype 1 predominant in the USA, Europe, and Japan. However, genotypes 3 and 4, dominant in developing countries, are resistant to a number of these drugs and little progress has been made towards understanding the structural basis of their drug resistivity. We have previously developed a 4D computational methodology, based on 3D structure modeling and molecular dynamics simulation, to analyze the active sites of the NS3 proteases of HCV-1b and 4a in relation to their catalytic activity and drug susceptibility. Here, we improved the methodology, extended the analysis to include genotype 3a (predominant in South Asia including Pakistan), and compared the results of the three genotypes (1b, 3a and 4a). The 4D analyses of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) indicate conformational instability of the catalytic site in HCV-3a and 4a compared to that of HCV-1b NS3 protease. The divergence is gradual and genotype-dependent, with HCV-1b being the most stable, HCV-4a being the most unstable and HCV-3a representing an intermediate state. These results suggest that the structural dynamics behavior, more than the rigid structure, could be related to the altered catalytic activity and drug susceptibility seen in NS3 proteases of HCV-3a and 4a.

Over-expression and characterization of NS3 and NS5A of Hepatitis C virus genotype 3a

Background

Hepatitis C virus (HCV) is a common and leading cause for liver cirrhosis and hepatocellular carcinoma. Current therapies to treat HCV infection are shown to be partially effective and poorly tolerated. Therefore, ample efforts are underway to rationally design therapies targeting the HCV non-structural proteins. Most of the work carried out in this direction has been focusing mainly on HCV genotype 1. Two direct-acting antiviral agents (DAAs) Telaprevir and Boceprevir are being used against genotype 1a infection in combination therapy with interferon and ribavirin. Unfortunately these DAAs are not effective against genotype 3a. Considering the wide spread infection by HCV genotype 3a in developing countries especially South Asia, we have focused on the recombinant production of antiviral drug targets NS3 and NS5A from HCV genotype 3a. These protein targets are to be used for screening of inhibitors.

Results

High-level expression of NS3 and NS5A was achieved at 25°C, using ~1 and 0.5 mM Isopropyl β-D-1-thiogalactopyranoside (IPTG), respectively. Yields of the purified NS3 and NS5A were 4 and 1 mg per liter culture volume, respectively. Although similar amounts of purified NS3 were obtained at 25 and 14°C, specificity constant (K_cat/K_m) was somewhat higher at expression temperature of 25°C. Circular dichroism (CD) and Fourier-transform infrared (FT-IR) spectroscopy revealed that both NS3 and NS5A contain a mixture of alpha-helix and beta-sheet secondary structures. For NS3 protein, percentages of secondary structures were similar to the values predicted from homology modeling.

Conclusions

NS3 and NS5A were over-expressed and using Nickel-affinity method both proteins were purified to ~ 95% purity. Yield of the purified NS3 obtained is four fold higher than previous reports. CD spectroscopy revealed that difference in activity of NS3 expressed at various temperatures is not related to changes in global structural features of the protein. Moreover, CD and FT-IR analysis showed that NS3 and NS5A contain both alpha-helical and beta-sheet structures and for NS5A, the proportion is almost equal. The production of NS3 and NS5A in milligram quantities will allow their characterization by biophysical and biochemical means that will help in designing new strategies to fight against HCV infection.