HIV-1 protease inhibitors with a P1 phosphonate modification maintain potency against drug-resistant variants by increased interactions with flap residues.

Protease inhibitors are the most potent antivirals against HIV-1, but they still lose efficacy against resistant variants. Improving the resistance profile is key to developing more robust inhibitors, which may be promising candidates for simplified next-generation antiretroviral therapies. In this study, we explored analogs of darunavir with a P1 phosphonate modification in combination with increasing size of the P1' hydrophobic group and various P2' moieties to improve potency against resistant variants. The phosphonate moiety substantially improved potency against highly mutated and resistant HIV-1 protease variants, but only when combined with more hydrophobic moieties at the P1' and P2' positions. Phosphonate analogs with a larger hydrophobic P1' moiety maintained excellent antiviral potency against a panel of highly resistant HIV-1 variants, with significantly improved resistance profiles. The cocrystal structures indicate that the phosphonate moiety makes extensive hydrophobic interactions with the protease, especially with the flap residues. Many residues involved in these protease-inhibitor interactions are conserved, enabling the inhibitors to maintain potency against highly resistant variants. These results highlight the need to balance inhibitor physicochemical properties by simultaneous modification of chemical groups to further improve resistance profiles.

Selection of HIV-1 for resistance to fifth-generation protease inhibitors reveals two independent pathways to high-level resistance.

Darunavir (DRV) is exceptional among potent HIV-1 protease inhibitors (PIs) in high drug concentrations that are achieved in vivo. Little is known about the de novo resistance pathway for DRV. We selected for resistance to high drug concentrations against 10 PIs and their structural precursor DRV. Mutations accumulated through two pathways (anchored by protease mutations I50V or I84V). Small changes in the inhibitor P1'-equivalent position led to preferential use of one pathway over the other. Changes in the inhibitor P2'-equivalent position determined differences in potency that were retained in the resistant viruses and that impacted the selected mutations. Viral variants from the two pathways showed differential selection of compensatory mutations in Gag cleavage sites. These results reveal the high level of selective pressure that is attainable with fifth-generation PIs and how features of the inhibitor affect both the resistance pathway and the residual potency in the face of resistance.

Deciphering the Molecular Mechanism of HCV Protease Inhibitor Fluorination as a General Approach to Avoid Drug Resistance.

Third generation Hepatitis C virus (HCV) NS3/4A protease inhibitors (PIs), glecaprevir and voxilaprevir, are highly effective across genotypes and against many resistant variants. Unlike earlier PIs, these compounds have fluorine substitutions on the P2-P4 macrocycle and P1 moieties. Fluorination has long been used in medicinal chemistry as a strategy to improve physicochemical properties and potency. However, the molecular basis by which fluorination improves potency and resistance profile of HCV NS3/4A PIs is not well understood. To systematically analyze the contribution of fluorine substitutions to inhibitor potency and resistance profile, we used a multi-disciplinary approach involving inhibitor design and synthesis, enzyme inhibition assays, co-crystallography, and structural analysis. A panel of inhibitors in matched pairs were designed with and without P4 cap fluorination, tested against WT protease and the D168A resistant variant, and a total of 22 high-resolution co-crystal structures were determined. While fluorination did not significantly improve potency against the WT protease, PIs with fluorinated P4 caps retained much better potency against the D168A protease variant. Detailed analysis of the co-crystal structures revealed that PIs with fluorinated P4 caps can sample alternate binding conformations that enable adapting to structural changes induced by the D168A substitution. Our results elucidate molecular mechanisms of fluorine-specific inhibitor interactions that can be leveraged in avoiding drug resistance.

Inhibiting HTLV-1 Protease: A Viable Antiviral Target.

Human T-cell lymphotropic virus type 1 (HTLV-1) is a retrovirus that can cause severe paralytic neurologic disease and immune disorders as well as cancer. An estimated 20 million people worldwide are infected with HTLV-1, with prevalence reaching 30% in some parts of the world. In stark contrast to HIV-1, no direct acting antivirals (DAAs) exist against HTLV-1. The aspartyl protease of HTLV-1 is a dimer similar to that of HIV-1 and processes the viral polyprotein to permit viral maturation. We report that the FDA-approved HIV-1 protease inhibitor darunavir (DRV) inhibits the enzyme with 0.8 µM potency and provides a scaffold for drug design against HTLV-1. Analogs of DRV that we designed and synthesized achieved submicromolar inhibition against HTLV-1 protease and inhibited Gag processing in viral maturation assays and in a chronically HTLV-1 infected cell line. Cocrystal structures of these inhibitors with HTLV-1 protease highlight opportunities for future inhibitor design. Our results show promise toward developing highly potent HTLV-1 protease inhibitors as therapeutic agents against HTLV-1 infections.

Structural Analysis of Potent Hybrid HIV-1 Protease Inhibitors Containing Bis-tetrahydrofuran in a Pseudosymmetric Dipeptide Isostere.

The design, synthesis, and X-ray structural analysis of hybrid HIV-1 protease inhibitors (PIs) containing bis-tetrahydrofuran (bis-THF) in a pseudo-C2-symmetric dipeptide isostere are described. A series of PIs were synthesized by incorporating bis-THF of darunavir on either side of the Phe-Phe isostere of lopinavir in combination with hydrophobic amino acids on the opposite P2/P2' position. Structure-activity relationship studies indicated that the bis-THF moiety can be attached at either the P2 or P2' position without significantly affecting potency. However, the group on the opposite P2/P2' position had a dramatic effect on potency depending on the size and shape of the side chain. Cocrystal structures of inhibitors with wild-type HIV-1 protease revealed that the bis-THF moiety retained similar interactions as observed in the darunavir-protease complex regardless of the position on the Phe-Phe isostere. Analyses of cocrystal structures and molecular dynamics simulations provide insights into optimizing HIV-1 PIs containing bis-THF in non-sulfonamide dipeptide isosteres.

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.

Molecular and Structural Mechanism of Pan-Genotypic HCV NS3/4A Protease Inhibition by Glecaprevir.

Hepatitis C virus, causative agent of chronic viral hepatitis, infects 71 million people worldwide and is divided into seven genotypes and multiple subtypes with sequence identities between 68 to 82%. While older generation direct-acting antivirals had varying effectiveness against different genotypes, the newest NS3/4A protease inhibitors including glecaprevir (GLE) have pan-genotypic activity. The structural basis for pan-genotypic inhibition and effects of polymorphisms on inhibitor potency were not well-known due to lack of crystal structures of GLE-bound NS3/4A or genotypes other than 1. In this study, we determined the crystal structures of NS3/4A from genotypes 1a, 3a, 4a, and 5a in complex with GLE. Comparison with the highly similar grazoprevir indicated the mechanism of GLE's drastic improvement in potency. We found that, while GLE is highly potent against wild-type NS3/4A of all genotypes, specific resistance-associated substitutions (RASs) confer orders of magnitude loss in inhibition. Our crystal structures reveal molecular mechanisms behind pan-genotypic activity of GLE, including potency loss due to RASs at D168. Our structures permit for the first time analysis of changes due to polymorphisms among genotypes, providing insights into design principles that can aid future drug development and potentially can be extended to other proteins.

Molecular Determinants of Epistasis in HIV-1 Protease: Elucidating the Interdependence of L89V and L90M Mutations in Resistance.

Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.

HIV-1 Protease Inhibitors Incorporating Stereochemically Defined P2' Ligands To Optimize Hydrogen Bonding in the Substrate Envelope.

A structure-guided design strategy was used to improve the resistance profile of HIV-1 protease inhibitors by optimizing hydrogen bonding and van der Waals interactions with the protease while staying within the substrate envelope. Stereoisomers of 4-(1-hydroxyethyl)benzene and 4-(1,2-dihydroxyethyl)benzene moieties were explored as P2' ligands providing pairs of diastereoisomers epimeric at P2', which exhibited distinct potency profiles depending on the configuration of the hydroxyl group and size of the P1' group. While compounds with the 4-(1-hydroxyethyl)benzene P2' moiety maintained excellent antiviral potency against a panel of multidrug-resistant HIV-1 strains, analogues with the polar 4-(1,2-dihydroxyethyl)benzene moiety were less potent, and only the (R)-epimer incorporating a larger 2-ethylbutyl P1' group showed improved potency. Crystal structures of protease-inhibitor complexes revealed strong hydrogen bonding interactions of both (R)- and (S)-stereoisomers of the hydroxyethyl group with Asp30'. Notably, the (R)-dihydroxyethyl group was involved in a unique pattern of direct hydrogen bonding interactions with the backbone amides of Asp29' and Asp30'. The SAR data and analysis of crystal structures provide insights for optimizing these promising HIV-1 protease inhibitors.

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.

Structural Adaptation of Darunavir Analogues against Primary Mutations in HIV-1 Protease.

HIV-1 protease is one of the prime targets of agents used in antiretroviral therapy against HIV. However, under selective pressure of protease inhibitors, primary mutations at the active site weaken inhibitor binding to confer resistance. Darunavir (DRV) is the most potent HIV-1 protease inhibitor in clinic; resistance is limited, as DRV fits well within the substrate envelope. Nevertheless, resistance is observed due to hydrophobic changes at residues including I50, V82, and I84 that line the S1/S1' pocket within the active site. Through enzyme inhibition assays and a series of 12 crystal structures, we interrogated susceptibility of DRV and two potent analogues to primary S1' mutations. The analogues had modifications at the hydrophobic P1' moiety compared to DRV to better occupy the unexploited space in the S1' pocket where the primary mutations were located. Considerable losses of potency were observed against protease variants with I84V and I50V mutations for all three inhibitors. The crystal structures revealed an unexpected conformational change in the flap region of I50V protease bound to the analogue with the largest P1' moiety, indicating interdependency between the S1' subsite and the flap region. Collective analysis of protease-inhibitor interactions in the crystal structures using principle component analysis was able to distinguish inhibitor identity and relative potency solely based on van der Waals contacts. Our results reveal the complexity of the interplay between inhibitor P1' moiety and S1' mutations and validate principle component analyses as a useful tool for distinguishing resistance and inhibitor potency.