Improving the Resistance Profile of Hepatitis C NS3/4A Inhibitors: Dynamic Substrate Envelope Guided Design

Elucidating the Substrate Envelope of Enterovirus 68-3C Protease: Structural Basis of Specificity and Potential Resistance

Enterovirus-D68 (EV68) has emerged as a global health concern over the last decade with severe symptomatic infections resulting in long-lasting neurological deficits and death. Unfortunately, there are currently no FDA-approved antiviral drugs for EV68 or any other non-polio enterovirus. One particularly attractive class of potential drugs are small molecules inhibitors, which can target the conserved active site of EV68-3C protease. For other viral proteases, we have demonstrated that the emergence of drug resistance can be minimized by designing inhibitors that leverage the evolutionary constraints of substrate specificity. However, the structural characterization of EV68-3C protease bound to its substrates has been lacking. Here, we have determined the substrate specificity of EV68-3C protease through molecular modeling, molecular dynamics (MD) simulations, and co-crystal structures. Molecular models enabled us to successfully characterize the conserved hydrogen-bond networks between EV68-3C protease and the peptides corresponding to the viral cleavage sites. In addition, co-crystal structures we determined have revealed substrate-induced conformational changes of the protease which involved new interactions, primarily surrounding the S1 pocket. We calculated the substrate envelope, the three-dimensional consensus volume occupied by the substrates within the active site. With the elucidation of the EV68-3C protease substrate envelope, we evaluated how 3C protease inhibitors, AG7088 and SG-85, fit within the active site to predict potential resistance mutations.

Pharmacophore-Assisted Covalent Docking Identifies a Potential Covalent Inhibitor for Drug-Resistant Genotype 3 Variants of Hepatitis C Viral NS3/4A Serine Protease

The emergence of drug-resistance-inducing mutations in Hepatitis C virus (HCV) coupled with genotypic heterogeneity has made targeting NS3/4A serine protease difficult. In this work, we investigated the mutagenic variations in the binding pocket of Genotype 3 (G3) HCV NS3/4A and evaluated ligands for efficacious inhibition. We report mutations at 14 positions within the ligand-binding residues of HCV NS3/4A, including H57R and S139P within the catalytic triad. We then modelled each mutational variant for pharmacophore-based virtual screening (PBVS) followed by covalent docking towards identifying a potential covalent inhibitor, i.e., cpd-217. The binding stability of cpd-217 was then supported by molecular dynamic simulation followed by MM/GBSA binding free energy calculation. The free energy decomposition analysis indicated that the resistant mutants alter the HCV NS3/4A-ligand interaction, resulting in unbalanced energy distribution within the binding site, leading to drug resistance. Cpd-217 was identified as interacting with all NS3/4A G3 variants with significant covalent docking scores. In conclusion, cpd-217 emerges as a potential inhibitor of HCV NS3/4A G3 variants that warrants further in vitro and in vivo studies. This study provides a theoretical foundation for drug design and development targeting HCV G3 NS3/4A.

In Silico-Based Identification of Natural Inhibitors from Traditionally Used Medicinal Plants that can Inhibit Dengue Infection

Dengue fever (DF) is an endemic disease that has become a public health concern around the globe. The NS3 protease-helicase enzyme is an important target for the development of antiviral drugs against DENV (dengue virus) due to its impact on viral replication. Inhibition of the activity of the NS3 protease-helicase enzyme complex significantly inhibits the infection associated with DENV. Unfortunately, there are no scientifically approved antiviral drugs for its prevention. However, this study has been developed to find natural bioactive molecules that can block the activity of the NS3 protease-helicase enzyme complex associated with DENV infection through molecular docking, MM-GBSA (molecular mechanics-generalized born surface area), and molecular dynamics (MD) simulations. Three hundred forty-two (342) compounds selected from twenty traditional medicinal plants were retrieved and screened against the NS3 protease-helicase protein by molecular docking and MM-GBSA studies, where the top six phytochemicals have been identified based on binding affinities. The six compounds were then subjected to pharmacokinetics and toxicity analysis, and we conducted molecular dynamics simulations on three protein-ligand complexes to validate their stability. Through computational analysis, this study revealed the potential of the two selected natural bioactive inhibitors (CID-440015 and CID-7424) as novel anti-dengue agents.

Open science discovery of potent noncovalent SARS-CoV-2 main protease inhibitors

We report the results of the COVID Moonshot, a fully open-science, crowdsourced, and structure-enabled drug discovery campaign targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease. We discovered a noncovalent, nonpeptidic inhibitor scaffold with lead-like properties that is differentiated from current main protease inhibitors. Our approach leveraged crowdsourcing, machine learning, exascale molecular simulations, and high-throughput structural biology and chemistry. We generated a detailed map of the structural plasticity of the SARS-CoV-2 main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. All compound designs (>18,000 designs), crystallographic data (>490 ligand-bound x-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2400 compounds) for this campaign were shared rapidly and openly, creating a rich, open, and intellectual property-free knowledge base for future anticoronavirus drug discovery.

Viral proteases: Structure, mechanism and inhibition

Viral proteases are diverse in structure, oligomeric state, catalytic mechanism, and substrate specificity. This chapter focuses on proteases from viruses that are relevant to human health: human immunodeficiency virus subtype 1 (HIV-1), hepatitis C (HCV), human T-cell leukemia virus type 1 (HTLV-1), flaviviruses, enteroviruses, and coronaviruses. The proteases of HIV-1 and HCV have been successfully targeted for therapeutics, with picomolar FDA-approved drugs currently used in the clinic. The proteases of HTLV-1 and the other virus families remain emerging therapeutic targets at different stages of the drug development process. This chapter provides an overview of the current knowledge on viral protease structure, mechanism, substrate recognition, and inhibition. Particular focus is placed on recent advances in understanding the molecular basis of diverse substrate recognition and resistance, which is essential toward designing novel protease inhibitors as antivirals.

Drug Design Strategies to Avoid Resistance in Direct-Acting Antivirals and Beyond

Drug resistance is prevalent across many diseases, rendering therapies ineffective with severe financial and health consequences. Rather than accepting resistance after the fact, proactive strategies need to be incorporated into the drug design and development process to minimize the impact of drug resistance. These strategies can be derived from our experience with viral disease targets where multiple generations of drugs had to be developed to combat resistance and avoid antiviral failure. Significant efforts including experimental and computational structural biology, medicinal chemistry, and machine learning have focused on understanding the mechanisms and structural basis of resistance against direct-acting antiviral (DAA) drugs. Integrated methods show promise for being predictive of resistance and potency. In this review, we give an overview of this research for human immunodeficiency virus type 1, hepatitis C virus, and influenza virus and the lessons learned from resistance mechanisms of DAAs. These lessons translate into rational strategies to avoid resistance in drug design, which can be generalized and applied beyond viral targets. While resistance may not be completely avoidable, rational drug design can and should incorporate strategies at the outset of drug development to decrease the prevalence of drug resistance.

Avoiding Drug Resistance in HIV Reverse Transcriptase

HIV reverse transcriptase (RT) is an enzyme that plays a major role in the replication cycle of HIV and has been a key target of anti-HIV drug development efforts. Because of the high genetic diversity of the virus, mutations in RT can impart resistance to various RT inhibitors. As the prevalence of drug resistance mutations is on the rise, it is necessary to design strategies that will lead to drugs less susceptible to resistance. Here we provide an in-depth review of HIV reverse transcriptase, current RT inhibitors, novel RT inhibitors, and mechanisms of drug resistance. We also present novel strategies that can be useful to overcome RT's ability to escape therapies through drug resistance. While resistance may not be completely avoidable, designing drugs based on the strategies and principles discussed in this review could decrease the prevalence of drug resistance.

In silico tools to study molecular targets of neglected diseases: inhibition of TcSir2rp3, an epigenetic enzyme of Trypanosoma cruzi

There is a growing interest to study and address neglected tropical diseases (NTD). To this end, in silico methods can serve as the bridge that connects academy and industry, encouraging the development of future treatments against these diseases. This chapter discusses current challenges in the development of new therapies, available computational methods and successful cases in computer-aided design with particular focus on human trypanosomiasis. Novel targets are also discussed. As a case study, we identify amentoflavone as a potential inhibitor of TcSir2rp3 (sirtuine) from Trypanosoma cruzi (20.03 μM) with a workflow that integrates chemoinformatic approaches, molecular modeling, and theoretical affinity calculations, as well as in vitro assays.

Avoiding Drug Resistance by Substrate Envelope-Guided Design: Toward Potent and Robust HCV NS3/4A Protease Inhibitors

Hepatitis C virus (HCV) infects millions of people worldwide, causing chronic liver disease that can lead to cirrhosis, hepatocellular carcinoma, and liver transplant. In the last several years, the advent of direct-acting antivirals, including NS3/4A protease inhibitors (PIs), has remarkably improved treatment outcomes of HCV-infected patients. However, selection of resistance-associated substitutions and polymorphisms among genotypes can lead to drug resistance and in some cases treatment failure. A proactive strategy to combat resistance is to constrain PIs within evolutionarily conserved regions in the protease active site. Designing PIs using the substrate envelope is a rational strategy to decrease the susceptibility to resistance by using the constraints of substrate recognition. We successfully designed two series of HCV NS3/4A PIs to leverage unexploited areas in the substrate envelope to improve potency, specifically against resistance-associated substitutions at D168. Our design strategy achieved better resistance profiles over both the FDA-approved NS3/4A PI grazoprevir and the parent compound against the clinically relevant D168A substitution. Crystallographic structural analysis and inhibition assays confirmed that optimally filling the substrate envelope is critical to improve inhibitor potency while avoiding resistance. Specifically, inhibitors that enhanced hydrophobic packing in the S4 pocket and avoided an energetically frustrated pocket performed the best. Thus, the HCV substrate envelope proved to be a powerful tool to design robust PIs, offering a strategy that can be translated to other targets for rational design of inhibitors with improved potency and resistance profiles.IMPORTANCE Despite significant progress, hepatitis C virus (HCV) continues to be a major health problem with millions of people infected worldwide and thousands dying annually due to resulting complications. Recent antiviral combinations can achieve >95% cure, but late diagnosis, low access to treatment, and treatment failure due to drug resistance continue to be roadblocks against eradication of the virus. We report the rational design of two series of HCV NS3/4A protease inhibitors with improved resistance profiles by exploiting evolutionarily constrained regions of the active site using the substrate envelope model. Optimally filling the S4 pocket is critical to avoid resistance and improve potency. Our results provide drug design strategies to avoid resistance that are applicable to other quickly evolving viral drug targets.

Metronidazole and Secnidazole Carbamates: Synthesis, Antiprotozoal Activity, and Molecular Dynamics Studies

We prepared a series of 10 carbamates derivatives based on two common antiprotozoal drugs: metronidazole (1–5) and secnidazole (6–10). The compounds were tested in vitro against a set of two amitochondriate protozoa: Giardia duodenalis and Trichomonas vaginalis. Compounds 1–10 showed strong antiprotozoal activities, with potency values in the low micromolar-to-nanomolar range, being more active than their parent drugs. Metronidazole carbamate (1) was the most active of the series, with nanomolar activities against G. duodenalis (IC₅₀ = 460 nM) and T. vaginalis (IC₅₀ = 60 nM). The potency of compound 1 was 10 times greater than that of metronidazole against both parasites. None of compounds showed in vitro cytotoxicity against VERO cells tested at 100 µM. Molecular dynamics of compounds 1–10, secnidazole, and metronidazole onto the ligand binding site of pyruvate–ferredoxin oxidoreductase of T. vaginalis and the modeled β-tubulin of G. duodenalis revealed putative molecular interactions with key residues in the binding site of both proteins implicated in the mode of action of the parent drugs.

Identification of novel proteolytically inactive mutations in coronavirus 3C-like protease using a combined approach

Coronaviruses (CoVs) infect humans and multiple other animal species, causing highly prevalent and severe diseases. 3C-like proteases (3CLpros) from CoVs (also called main proteases) are essential for viral replication and are also involved in polyprotein cleavage and immune regulation, making them attractive and effective targets for the development of antiviral drugs. Herein, the 3CLpro from the porcine epidemic diarrhea virus, an enteropathogenic CoV, was used as a model to identify novel crucial residues for enzyme activity. First, we established a rapid, sensitive, and efficient luciferase-based biosensor to monitor the activity of PDEV 3CLproin vivo. Using this luciferase biosensor, along with confirming the well-known catalytic residues (His41 and Cys144), we identified 4 novel proteolytically inactivated mutants of PDEV 3CLpro, which was also confirmed in mammalian cells by biochemical experiments. Our molecular dynamics (MD) simulations showed that the hydrogen bonding interactions occurring within and outside of the protease's active site and the dynamic fluctuations of the substrate, especially the van der Waals contacts, were drastically altered, a situation related to the loss of 3CLpro activity. These data suggest that changing the intermolecular dynamics in protein-substrate complexes eliminates the mechanism underlying the protease activity. The discovery of novel crucial residues for enzyme activity in the binding pocket could potentially provide more druggable sites for the design of protease inhibitors. In addition, our in-depth study of the dynamic substrate's envelope model using MD simulations is an approach that could augment the discovery of new inhibitors against 3CLpro in CoVs and other viral 3C proteases.-Zhou, J., Fang, L., Yang, Z., Xu, S., Lv, M., Sun, Z., Chen, J., Wang, D., Gao, J., Xiao, S. Identification of novel proteolytically inactive mutations in coronavirus 3C-like protease using a combined approach.

Picomolar to Micromolar: Elucidating the Role of Distal Mutations in HIV-1 Protease in Conferring Drug Resistance

Drug resistance continues to be a growing global problem. The efficacy of small molecule inhibitors is threatened by pools of genetic diversity in all systems, including antibacterials, antifungals, cancer therapeutics, and antivirals. Resistant variants often include combinations of active site mutations and distal "secondary" mutations, which are thought to compensate for losses in enzymatic activity. HIV-1 protease is the ideal model system to investigate these combinations and underlying molecular mechanisms of resistance. Darunavir (DRV) binds wild-type (WT) HIV-1 protease with a potency of <5 pM, but we have identified a protease variant that loses potency to DRV 150 000-fold, with 11 mutations in and outside the active site. To elucidate the roles of these mutations in DRV resistance, we used a multidisciplinary approach, combining enzymatic assays, crystallography, and molecular dynamics simulations. Analysis of protease variants with 1, 2, 4, 8, 9, 10, and 11 mutations showed that the primary active site mutations caused ∼50-fold loss in potency (2 mutations), while distal mutations outside the active site further decreased DRV potency from 13 nM (8 mutations) to 0.76 μM (11 mutations). Crystal structures and simulations revealed that distal mutations induce subtle changes that are dynamically propagated through the protease. Our results reveal that changes remote from the active site directly and dramatically impact the potency of the inhibitor. Moreover, we find interdependent effects of mutations in conferring high levels of resistance. These mechanisms of resistance are likely applicable to many other quickly evolving drug targets, and the insights may have implications for the design of more robust inhibitors.

Identification of potential inhibitors for HCV NS5b of genotype 4a by combining dynamic simulation, protein-ligand interaction fingerprint, 3D pharmacophore, docking and 3D QSAR

HCV NS5B polymerase has been one of the most attractive targets for developing new drugs for HCV infection and many drugs were successfully developed, but all of them were designed for targeting Hepatitis C Virus genotype 1 (HCV GT1). Hepatitis C virus genotype 4a (HCV GT4a) dominant in Egypt has paid less attention. Here, we describe our protocol of virtual screening in identification of novel potential potent inhibitors for HCV NS5B polymerase of GT4a using homology modeling, protein-ligand interaction fingerprint (PLIF), docking, pharmacophore, and 3D CoMFA quantitative structure activity relationship (QSAR). Firstly, a high-quality 3D model of HCV NS5B polymerase of GT4a was constructed using crystal structure of HCV NS5B polymerase of GT1 (PDB ID: 3hkw) as a template. Then, both the model and the template were simulated to compare conformational stability. PLIF was generated using five crystal structures of HCV NS5B (PDB ID: 4mia, 4mib, 4mk9, 4mka, and 4mkb), which revealed the most important residues and their interactions with the co-crystalized ligands. After that, a 3D pharmacophore model was developed from the generated PLIF data and then used as a screening filter for 17000328 drug-like zinc database compounds. 900 compounds passed the pharmacophore filter and entered the docking-based virtual screening stage. Finally, a 3D CoMFA QSAR was developed using 42 compounds as a training and 19 compounds as a test set. The 3D CoMFA QSAR was used to design and screen some potential inhibitors, these compounds were further evaluated by the docking stage. The highest ranked five hits from docking result (compounds (p1-p4) and compound q1) were selected for further analysis.Communicated by Ramaswamy H. Sarma.

Resistance outside the substrate envelope: hepatitis C NS3/4A protease inhibitors

Direct acting antivirals have dramatically increased the efficacy and tolerability of hepatitis C treatment, but drug resistance has emerged with some of these inhibitors, including nonstructural protein 3/4 A protease inhibitors (PIs). Although many co-crystal structures of PIs with the NS3/4A protease have been reported, a systematic review of these crystal structures in the context of the rapidly emerging drug resistance especially for early PIs has not been performed. To provide a framework for designing better inhibitors with higher barriers to resistance, we performed a quantitative structural analysis using co-crystal structures and models of HCV NS3/4A protease in complex with natural substrates and inhibitors. By comparing substrate structural motifs and active site interactions with inhibitor recognition, we observed that the selection of drug resistance mutations correlates with how inhibitors deviate from viral substrates in molecular recognition. Based on this observation, we conclude that guiding the design process with native substrate recognition features is likely to lead to more robust small molecule inhibitors with decreased susceptibility to resistance.

Elucidating the Interdependence of Drug Resistance from Combinations of Mutations

HIV-1 protease is responsible for the cleavage of 12 nonhomologous sites within the Gag and Gag-Pro-Pol polyproteins in the viral genome. Under the selective pressure of protease inhibition, the virus evolves mutations within (primary) and outside of (secondary) the active site, allowing the protease to process substrates while simultaneously countering inhibition. The primary protease mutations impede inhibitor binding directly, while the secondary mutations are considered accessory mutations that compensate for a loss in fitness. However, the role of secondary mutations in conferring drug resistance remains a largely unresolved topic. We have shown previously that mutations distal to the active site are able to perturb binding of darunavir (DRV) via the protein's internal hydrogen-bonding network. In this study, we show that mutations distal to the active site, regardless of context, can play an interdependent role in drug resistance. Applying eigenvalue decomposition to collections of hydrogen bonding and van der Waals interactions from a series of molecular dynamics simulations of 15 diverse HIV-1 protease variants, we identify sites in the protease where amino acid substitutions lead to perturbations in nonbonded interactions with DRV and/or the hydrogen-bonding network of the protease itself. While primary mutations are known to drive resistance in HIV-1 protease, these findings delineate the significant contributions of accessory mutations to resistance. Identifying the variable positions in the protease that have the greatest impact on drug resistance may aid in future structure-based design of inhibitors.

Exploring the Drug Resistance of HCV Protease

Hepatitis C virus (HCV) currently affects several million people across the globe. One of the major classes of drugs against HCV inhibits the NS3/4A protease of the polyprotein chain. Efficacy of these drugs is severely limited due to the high mutation rate that results in several genetically related quasispecies. The molecular mechanism of drug resistance is frequently deduced from structural studies and binding free energies. However, prediction of new mutations requires the evaluation of both binding free energy of the drug as well as the parameters (k_cat and K_M) for the natural substrate. The vitality values offer a good approach to investigate and predict mutations that render resistance to the inhibitor. A successful mutation should only affect the binding of the drug and not the catalytic activity and binding of the natural substrate. In this article, we have calculated the vitality values for four known drug inhibitors that are either currently in use or in clinical trials, evaluating binding free energies by the relevant PDLD/S-LRA method and activation barriers by the EVB method. The molecular details pertaining to resistance are also discussed. We show that our calculations are able to reproduce the catalytic effects and binding free energies in a good agreement with the corresponding observed values. Importantly, previous computational approaches have not been able to achieve this task. The trend for the vitality values is in accordance with experimental findings. Finally, we calculate the vitality values for mutations that have either not been studied experimentally or reported for some inhibitors.

Hepatitis C Virus NS3/4A Protease Inhibitors Incorporating Flexible P2 Quinoxalines Target Drug Resistant Viral Variants

A substrate envelope-guided design strategy is reported for improving the resistance profile of HCV NS3/4A protease inhibitors. Analogues of 5172-mcP1P3 were designed by incorporating diverse quinoxalines at the P2 position that predominantly interact with the invariant catalytic triad of the protease. Exploration of structure-activity relationships showed that inhibitors with small hydrophobic substituents at the 3-position of P2 quinoxaline maintain better potency against drug resistant variants, likely due to reduced interactions with residues in the S2 subsite. In contrast, inhibitors with larger groups at this position were highly susceptible to mutations at Arg155, Ala156, and Asp168. Excitingly, several inhibitors exhibited exceptional potency profiles with EC₅₀ values ≤5 nM against major drug resistant HCV variants. These findings support that inhibitors designed to interact with evolutionarily constrained regions of the protease, while avoiding interactions with residues not essential for substrate recognition, are less likely to be susceptible to drug resistance.

Database-Centric Method for Automated High-Throughput Deconvolution and Analysis of Kinetic Antibody Screening Data

The state-of-the-art industrial drug discovery approach is the empirical interrogation of a library of drug candidates against a target molecule. The advantage of high-throughput kinetic measurements over equilibrium assessments is the ability to measure each of the kinetic components of binding affinity. Although high-throughput capabilities have improved with advances in instrument hardware, three bottlenecks in data processing remain: (1) intrinsic molecular properties that lead to poor biophysical quality in vitro are not accounted for in commercially available analysis models, (2) processing data through a user interface is time-consuming and not amenable to parallelized data collection, and (3) a commercial solution that includes historical kinetic data in the analysis of kinetic competition data does not exist. Herein, we describe a generally applicable method for the automated analysis, storage, and retrieval of kinetic binding data. This analysis can deconvolve poor quality data on-the-fly and store and organize historical data in a queryable format for use in future analyses. Such database-centric strategies afford greater insight into the molecular mechanisms of kinetic competition, allowing for the rapid identification of allosteric effectors and the presentation of kinetic competition data in absolute terms of percent bound to antigen on the biosensor.

Interdependence of Inhibitor Recognition in HIV-1 Protease

Molecular recognition is a highly interdependent process. Subsite couplings within the active site of proteases are most often revealed through conditional amino acid preferences in substrate recognition. However, the potential effect of these couplings on inhibition and thus inhibitor design is largely unexplored. The present study examines the interdependency of subsites in HIV-1 protease using a focused library of protease inhibitors, to aid in future inhibitor design. Previously a series of darunavir (DRV) analogs was designed to systematically probe the S1′ and S2′ subsites. Co-crystal structures of these analogs with HIV-1 protease provide the ideal opportunity to probe subsite interdependency. All-atom molecular dynamics simulations starting from these structures were performed and systematically analyzed in terms of atomic fluctuations, intermolecular interactions, and water structure. These analyses reveal that the S1′ subsite highly influences other subsites: the extension of the hydrophobic P1′ moiety results in 1) reduced van der Waals contacts in the P2′ subsite, 2) more variability in the hydrogen bond frequencies with catalytic residues and the flap water, and 3) changes in the occupancy of conserved water sites both proximal and distal to the active site. In addition, one of the monomers in this homodimeric enzyme has atomic fluctuations more highly correlated with DRV than the other monomer. These relationships intricately link the HIV-1 protease subsites and are critical to understanding molecular recognition and inhibitor binding. More broadly, the interdependency of subsite recognition within an active site requires consideration in the selection of chemical moieties in drug design; this strategy is in contrast to what is traditionally done with independent optimization of chemical moieties of an inhibitor.

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.

Molecular Basis for Differential Patterns of Drug Resistance in Influenza N1 and N2 Neuraminidase

Neuraminidase (NA) inhibitors are used for the prevention and treatment of influenza A virus infections. Two subtypes of NA, N1 and N2, predominate in viruses that infect humans, but differential patterns of drug resistance have emerged in each subtype despite highly homologous active sites. To understand the molecular basis for the selection of these drug resistance mutations, structural and dynamic analyses on complexes of N1 and N2 NA with substrates and inhibitors were performed. Comparison of dynamic substrate and inhibitor envelopes and interactions at the active site revealed how differential patterns of drug resistance have emerged for specific drug resistance mutations, at residues I222, S246, and H274 in N1 and E119 in N2. Our results show that the differences in intermolecular interactions, especially van der Waals contacts, of the inhibitors versus substrates at the NA active site effectively explain the selection of resistance mutations in the two subtypes. Avoiding such contacts that render inhibitors vulnerable to resistance by better mimicking the dynamics and intermolecular interactions of substrates can lead to the development of novel inhibitors that avoid drug resistance in both subtypes.

Molecular and Dynamic Mechanism Underlying Drug Resistance in Genotype 3 Hepatitis C NS3/4A Protease

Hepatitis C virus (HCV), affecting an estimated 150 million people worldwide, is the leading cause of viral hepatitis, cirrhosis and hepatocellular carcinoma. HCV is genetically diverse with six genotypes (GTs) and multiple subtypes of different global distribution and prevalence. Recent development of direct-acting antivirals against HCV including NS3/4A protease inhibitors (PIs) has greatly improved treatment outcomes for GT-1. However, all current PIs exhibit significantly lower potency against GT-3. Lack of structural data on GT-3 protease has limited our ability to understand PI failure in GT-3. In this study the molecular basis for reduced potency of current inhibitors against GT-3 NS3/4A protease is elucidated with structure determination, molecular dynamics simulations and inhibition assays. A chimeric GT-1a3a NS3/4A protease amenable to crystallization was engineered to recapitulate decreased sensitivity of GT-3 protease to PIs. High-resolution crystal structures of this GT-1a3a bound to 3 PIs, asunaprevir, danoprevir and vaniprevir, had only subtle differences relative to GT-1 despite orders of magnitude loss in affinity. In contrast, hydrogen-bonding interactions within and with the protease active site and dynamic fluctuations of the PIs were drastically altered. The correlation between loss of intermolecular dynamics and inhibitor potency suggests a mechanism where polymorphisms between genotypes (or selected mutations) in the drug target confer resistance through altering the intermolecular dynamics of the protein-inhibitor complex.

Improving Viral Protease Inhibitors to Counter Drug Resistance

Drug resistance is a major problem in health care, undermining therapy outcomes and necessitating novel approaches to drug design. Extensive studies on resistance to viral protease inhibitors, particularly those of HIV-1 and hepatitis C virus (HCV) protease, revealed a plethora of information on the structural and molecular mechanisms underlying resistance. These insights led to several strategies to improve viral protease inhibitors to counter resistance, such as exploiting the essential biological function and leveraging evolutionary constraints. Incorporation of these strategies into structure-based drug design can minimize vulnerability to resistance, not only for viral proteases but for other quickly evolving drug targets as well, toward designing inhibitors one step ahead of evolution to counter resistance with more intelligent and rational design.

Structural and Thermodynamic Effects of Macrocyclization in HCV NS3/4A Inhibitor MK-5172

Recent advances in direct-acting antivirals against Hepatitis C Virus (HCV) have led to the development of potent inhibitors, including MK-5172, that target the viral NS3/4A protease with relatively low susceptibility to resistance. MK-5172 has a P2-P4 macrocycle and a unique binding mode among current protease inhibitors where the P2 quinoxaline packs against the catalytic residues H57 and D81. However, the effect of macrocyclization on this binding mode is not clear, as is the relation between macrocyclization, thermodynamic stabilization, and susceptibility to the resistance mutation A156T. We have determined high-resolution crystal structures of linear and P1-P3 macrocyclic analogs of MK-5172 bound to WT and A156T protease and compared these structures, their molecular dynamics, and experimental binding thermodynamics to the parent compound. We find that the "unique" binding mode of MK-5172 is conserved even when the P2-P4 macrocycle is removed or replaced with a P1-P3 macrocycle. While beneficial to decreasing the entropic penalty associated with binding, the constraint exerted by the P2-P4 macrocycle prevents efficient rearrangement to accommodate the A156T mutation, a deficit alleviated in the linear and P1-P3 analogs. Design of macrocyclic inhibitors against NS3/4A needs to achieve the best balance between exerting optimal conformational constraint for enhancing potency, fitting within the substrate envelope and allowing adaptability to be robust against resistance mutations.

Anti-HIV Drug Discovery and Development: Current Innovations and Future Trends

The early effectiveness of combinatorial antiretroviral therapy (cART) in the treatment of HIV infection has been compromised to some extent by rapid development of multidrug-resistant HIV strains, poor bioavailability, and cumulative toxicities, and so there is a need for alternative strategies of antiretroviral drug discovery and additional therapeutic agents with novel action modes or targets. From this perspective, we first review current strategies of antiretroviral drug discovery and optimization, with the aid of selected examples from the recent literature. We highlight the development of phosphate ester-based prodrugs as a means to improve the aqueous solubility of HIV inhibitors, and the introduction of the substrate envelope hypothesis as a new approach for overcoming HIV drug resistance. Finally, we discuss future directions for research, including opportunities for exploitation of novel antiretroviral targets, and the strategy of activation of latent HIV reservoirs as a means to eradicate the virus.

Structural analysis of asunaprevir resistance in HCV NS3/4A protease

Asunaprevir (ASV), an isoquinoline-based competitive inhibitor targeting the hepatitis C virus (HCV) NS3/4A protease, is very potent in vivo. However, the potency is significantly compromised by the drug resistance mutations R155K and D168A. In this study three crystal structures of ASV and an analogue were determined to analyze the structural basis of drug resistance susceptibility. These structures revealed that ASV makes extensive contacts with Arg155 outside the substrate envelope. Arg155 in turn is stabilized by Asp168, and thus when either residue is mutated, the enzyme’s interaction with ASV’s P_2* isoquinoline is disrupted. Adding a P₁–P₃ macrocycle to ASV enhances the inhibitor’s resistance barrier, likely due to poising the inhibitor to its bound conformation. Macrocyclic inhibitors with P_2* extension moieties avoiding interaction with the protease S₂ residues including Arg155 must be chosen for future design of more robust protease inhibitors.

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.