Well-designed viral protease inhibitors (PIs) potently inhibit replication as well as create a high genetic barrier for resistance. Through in vivo selective pressure, we have generated high-level resistance against ten HIV-1 PIs and their precursor, the FDA-approved drug darunavir (DRV), achieving 1,000-fold resistance over the starting EC50. The accumulation of mutations revealed two pathways to high-level resistance, resulting in protease variants with up to 14 mutations in and outside of the active site. The two pathways demonstrate the interplay between drug resistance and viral fitness. Replicate selections showed that one inhibitor could select for resistance through either pathway, although subtle changes in chemical structure of the inhibitors led to preferential use of one pathway over the other. 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 these fourth-generation protease inhibitors, and the interplay between selection of mutations to confer resistance while maintaining viral fitness.

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 analogs to primary S1' mutations. The analogs had modifications at the hydrophobic P1' moiety 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 analog with the largest P1' moiety, indicating interdependency between the S1' subsite and the flap region. Collective analysis of protease-inhibitor van der Waals (vdW) interactions in the crystal structures using principle component analysis indicated I84V mutation underlying the largest variation in the vdW contacts. Interestingly, the principle components were able to distinguish inhibitor identity and relative potency solely based on vdW interactions of active site residues in the crystal structures. Our results reveal the interplay between inhibitor P1' moiety and primary S1' mutations, as well as suggesting a novel method for distinguishing the interdependence of resistance through principle component analyses.

Hepatitis C virus (HCV), 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 (DAAs) 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 (GZR) indicated the mechanism of GLEs 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.

Abstract The evolutionarily conserved AAA + ATPases Rvb1 and Rvb2 proteins form a heteromeric complex (Rvb1/2) required for assembly or remodeling of macromolecular complexes in essential cellular processes ranging from chromatin remodeling to ribosome biogenesis. Rvb1 and Rvb2 have a high degree of sequence and structural similarity, and both contain the classical features of ATPases of their clade, including an N-terminal AAA + subdomain with the Walker A motif, an insertion domain that typically interacts with various binding partners, and a C-terminal AAA + subdomain containing a Walker B motif, the Sensor I and II motifs, and an arginine finger. In this study, we find that despite the high degree of structural similarity, Rvb1 and Rvb2 have distinct active sites that impact their activities and regulation within the Rvb1/2 complex. Using a combination of biochemical and genetic approaches, we show that replacing the homologous arginine fingers of Rvb1 and Rvb2 with different amino acids not only has distinct effects on the catalytic activity of the complex, but also impacts cell growth, and the Rvb1/2 interactions with binding partners. Using molecular dynamics simulations, we find that changes near the active site of Rvb1 and Rvb2 cause long-range effects on the protein dynamics in the insertion domain, suggesting a molecular basis for how enzymatic activity within the catalytic site of ATP hydrolysis can be relayed to other domains of the Rvb1/2 complex to modulate its function. Further, we show the impact that the arginine finger variants have on snoRNP biogenesis and validate the findings from molecular dynamics simulations using a targeted genetic screen. Together, our results reveal new aspects of the regulation of the Rvb1/2 complex by identifying a relay of long-range molecular communication from the ATPase active site of the complex to the binding site of cofactors. Most importantly, our findings suggest that despite high similarity and cooperation within the same protein complex, the two proteins have evolved with unique properties critical for the regulation and function of the Rvb1/2 complex. Significance AAA ATPases constitute a large family of proteins involved in various essential cellular functions in living organisms in all kingdoms of life. Members of this family typically form homo or hetero multimers that convert the energy from ATP hydrolysis to mechanical work. How the conserved features of AAA ATPases relay the energy from ATP hydrolysis to other functional domains of the complex remains largely unknown. Here, using arginine finger variants of Rvb1 and Rvb2, two evolutionarily conserved closely related AAA + ATPases that form a heterohexameric complex, we reveal how individual protomers in a heteromeric complex can uniquely contribute to the overall function of the complex and how changes in the ATP binding site can be relayed to distal functional domains.

ABSTRACT 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.