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

Molecular mechanisms of drug resistance and compensation in SARS-CoV-2 main protease: the interplay between E166 and L50

The SARS-CoV-2 main protease (Mpro) is essential for viral replication and is a primary target for COVID-19 antivirals. Direct-acting antivirals such as nirmatrelvir, the active component of Paxlovid, target the Mpro active site to block viral polyprotein cleavage and thus replication. However, drug resistance mutations at the active site residue Glu166 (E166) have emerged during in vitro selection studies, raising concerns about the durability of current antiviral strategies. Here, we investigate the molecular basis of drug resistance conferred by E166A and E166V mutations against nirmatrelvir and the related PF-00835231, individually and in combination with the distal mutation L50F. We found that E166 mutations reduce nirmatrelvir potency by up to 3,000-fold while preserving substrate cleavage, with catalytic efficiency reduced by only up to twofold. This loss of catalytic efficiency was compensated for by the addition of L50F in the double-mutant variants. We have determined three cocrystal structures of the E166 variants (E166A, E166V, and E166V/L50F) bound to PF-00835231. Comparison of these structures with wild-type enzyme demonstrated that E166 is crucial for dimerization and for shaping the substrate-binding S1 pocket. Our findings highlight the mutability of E166, a prime site for resistance for inhibitors that leverage direct interactions with this position, and the potential emergence of highly resistant and active variants in combination with the compensatory mutation L50F. These insights support the design of inhibitors that target conserved protease features and avoid E166 side-chain interactions to minimize susceptibility to resistance.

IMPORTANCE

Drug resistance remains a great challenge to modern medicine. This study investigates SARS-CoV-2 main protease variants E166A and E166V which confer nirmatrelvir resistance. These variants can retain considerable enzymatic activity through combination with the compensatory mutation L50F. For single- and double-mutant variant enzymes, we assessed catalytic efficiency, measured loss in potency for nirmatrelvir and its analog PF-00835231, and cocrystallized with inhibitors to investigate drug resistance caused by these mutations. Our results contribute toward understanding of molecular mechanisms of resistance and combinations of mutations, which pushes toward resistance-thwarting inhibitor design. These principles also apply broadly to many quickly evolving drug targets in infectious diseases.

Structural Analysis of Inhibitor Binding to Enterovirus-D68 3C Protease

Enterovirus-D68 (EV68) continues to present as a global health issue causing respiratory illness and outbreaks associated with long-lasting neurological disease, with no antivirals or specific treatment options. The development of antiviral therapeutics, such as small-molecule inhibitors that target conserved proteins like the enteroviral 3C protease, remains to be achieved. While various 3C inhibitors have been investigated, their design does not consider the potential emergence of drug resistance mutations. For other antivirals where resistance has been a challenge, we have demonstrated that the likelihood of resistance can be minimized by designing inhibitors that leverage the evolutionary constraints of the target. Here, we characterize a series of 3C inhibitors against EV68-3C protease through enzyme inhibition, protein crystallography, and structural analysis. We have determined and analyzed three high-resolution inhibitor-bound crystal structures of EV68-3C protease, which revealed possible sites of resistance mutations, a key structural water molecule conserved during ligand binding, and the conformational flexibility of the catalytic histidine H40. This structural analysis combined with enzymatic assays provides insights for the rational design of inhibitors that are robust against resistance toward developing antiviral treatments for EV68 infections.

Viral DNA polymerase structures reveal mechanisms of antiviral drug resistance

DNA polymerases are important drug targets, and many structural studies have captured them in distinct conformations. However, a detailed understanding of the impact of polymerase conformational dynamics on drug resistance is lacking. We determined cryoelectron microscopy (cryo-EM) structures of DNA-bound herpes simplex virus polymerase holoenzyme in multiple conformations and interacting with antivirals in clinical use. These structures reveal how the catalytic subunit Pol and the processivity factor UL42 bind DNA to promote processive DNA synthesis. Unexpectedly, in the absence of an incoming nucleotide, we observed Pol in multiple conformations with the closed state sampled by the fingers domain. Drug-bound structures reveal how antivirals may selectively bind enzymes that more readily adopt the closed conformation. Molecular dynamics simulations and the cryo-EM structure of a drug-resistant mutant indicate that some resistance mutations modulate conformational dynamics rather than directly impacting drug binding, thus clarifying mechanisms that drive drug selectivity.

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.

Target-based drug design strategies to overcome resistance to antiviral agents: opportunities and challenges

Viral infections have a major impact in human health. Ongoing viral transmission and escalating selective pressure have the potential to favor the emergence of vaccine- and antiviral drug-resistant viruses. Target-based approaches for the design of antiviral drugs can play a pivotal role in combating drug-resistant challenges. Drug design computational tools facilitate the discovery of novel drugs. This review provides a comprehensive overview of current drug design strategies employed in the field of antiviral drug resistance, illustrated through the description of a series of successful applications. These strategies include technologies that enhance compound-target affinity while minimizing interactions with mutated binding pockets. Furthermore, emerging approaches such as virtual screening, targeted protein/RNA degradation, and resistance analysis during drug design have been harnessed to curtail the emergence of drug resistance. Additionally, host targeting antiviral drugs offer a promising avenue for circumventing viral mutation. The widespread adoption of these refined drug design strategies will effectively address the prevailing challenge posed by antiviral drug resistance.

Modeling and Analysis of HIV-1 Pol Polyprotein as a Case Study for Predicting Large Polyprotein Structures

Acquired immunodeficiency syndrome (AIDS) is caused by human immunodeficiency virus (HIV). HIV protease, reverse transcriptase, and integrase are targets of current drugs to treat the disease. However, anti-viral drug-resistant strains have emerged quickly due to the high mutation rate of the virus, leading to the demand for the development of new drugs. One attractive target is Gag-Pol polyprotein, which plays a key role in the life cycle of HIV. Recently, we found that a combination of M50I and V151I mutations in HIV-1 integrase can suppress virus release and inhibit the initiation of Gag-Pol autoprocessing and maturation without interfering with the dimerization of Gag-Pol. Additional mutations in integrase or RNase H domain in reverse transcriptase can compensate for the defect. However, the molecular mechanism is unknown. There is no tertiary structure of the full-length HIV-1 Pol protein available for further study. Therefore, we developed a workflow to predict the tertiary structure of HIV-1 NL4.3 Pol polyprotein. The modeled structure has comparable quality compared with the recently published partial HIV-1 Pol structure (PDB ID: 7SJX). Our HIV-1 NL4.3 Pol dimer model is the first full-length Pol tertiary structure. It can provide a structural platform for studying the autoprocessing mechanism of HIV-1 Pol and for developing new potent drugs. Moreover, the workflow can be used to predict other large protein structures that cannot be resolved via conventional experimental methods.

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.

Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.

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.

Activity-based CRISPR scanning uncovers allostery in DNA methylation maintenance machinery

Allostery enables dynamic control of protein function. A paradigmatic example is the tightly orchestrated process of DNA methylation maintenance. Despite the fundamental importance of allosteric sites, their identification remains highly challenging. Here, we perform CRISPR scanning on the essential maintenance methylation machinery-DNMT1 and its partner UHRF1-with the activity-based inhibitor decitabine to uncover allosteric mechanisms regulating DNMT1. In contrast to non-covalent DNMT1 inhibition, activity-based selection implicates numerous regions outside the catalytic domain in DNMT1 function. Through computational analyses, we identify putative mutational hotspots in DNMT1 distal from the active site that encompass mutations spanning a multi-domain autoinhibitory interface and the uncharacterized BAH2 domain. We biochemically characterize these mutations as gain-of-function, exhibiting increased DNMT1 activity. Extrapolating our analysis to UHRF1, we discern putative gain-of-function mutations in multiple domains, including key residues across the autoinhibitory TTD-PBR interface. Collectively, our study highlights the utility of activity-based CRISPR scanning for nominating candidate allosteric sites, and more broadly, introduces new analytical tools that further refine the CRISPR scanning framework.

Contingency and Entrenchment of Drug-Resistance Mutations in HIV Viral Proteins

The ability of HIV-1 to rapidly mutate leads to antiretroviral therapy (ART) failure among infected patients. Drug-resistance mutations (DRMs), which cause a fitness penalty to intrinsic viral fitness, are compensated by accessory mutations with favorable epistatic interactions which cause an evolutionary trapping effect, but the kinetics of this overall process has not been well characterized. Here, using a Potts Hamiltonian model describing epistasis combined with kinetic Monte Carlo simulations of evolutionary trajectories, we explore how epistasis modulates the evolutionary dynamics of HIV DRMs. We show how the occurrence of a drug-resistance mutation is contingent on favorable epistatic interactions with many other residues of the sequence background and that subsequent mutations entrench DRMs. We measure the time-autocorrelation of fluctuations in the likelihood of DRMs due to epistatic coupling with the sequence background, which reveals the presence of two evolutionary processes controlling DRM kinetics with two distinct time scales. Further analysis of waiting times for the evolutionary trapping effect to reverse reveals that the sequences which entrench (trap) a DRM are responsible for the slower time scale. We also quantify the overall strength of epistatic effects on the evolutionary kinetics for different mutations and show these are much larger for DRM positions than polymorphic positions, and we also show that trapping of a DRM is often caused by the collective effect of many accessory mutations, rather than a few strongly coupled ones, suggesting the importance of multiresidue sequence variations in HIV evolution. The analysis presented here provides a framework to explore the kinetic pathways through which viral proteins like HIV evolve under drug-selection pressure.

The evolution of the HIV-1 protease folding stability

The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4⁺ T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.

HIV-1 protease with 10 lopinavir and darunavir resistance mutations exhibits altered inhibition, structural rearrangements and extreme dynamics

Antiretroviral drug resistance is a therapeutic obstacle for people with HIV. HIV protease inhibitors darunavir and lopinavir are recommended for resistant infections. We characterized a protease mutant (PR10x) derived from a highly resistant clinical isolate including 10 mutations associated with resistance to lopinavir and darunavir. Compared to the wild-type protease, PR10x exhibits ∼3-fold decrease in catalytic efficiency and Ki values of 2-3 orders of magnitude worse for darunavir, lopinavir, and potent investigational inhibitor GRL-519. Crystal structures of the mutant were solved in a ligand-free form and in complex with GRL-519. The structures show altered interactions in the active site, flap-core interface, hydrophobic core, hinge region, and 80s loop compared to the corresponding wild-type protease structures. The ligand-free crystal structure exhibits a highly curled flap conformation which may amplify drug resistance. Molecular dynamics simulations performed for 1 μs on ligand-free dimers showed extremely large fluctuations in the flaps for PR10x compared to equivalent simulations on PR with a single L76V mutation or wild-type protease. This analysis offers insight about the synergistic effects of mutations in highly resistant variants.

Prediction and molecular field view of drug resistance in HIV-1 protease mutants

Conquering the mutational drug resistance is a great challenge in anti-HIV drug development and therapy. Quantitatively predicting the mutational drug resistance in molecular level and elucidating the three dimensional structure-resistance relationships for anti-HIV drug targets will help to improve the understanding of the drug resistance mechanism and aid the design of resistance evading inhibitors. Here the MB-QSAR (Mutation-dependent Biomacromolecular Quantitative Structure Activity Relationship) method was employed to predict the molecular drug resistance of HIV-1 protease mutants towards six drugs, and to depict the structure resistance relationships in HIV-1 protease mutants. MB-QSAR models were constructed based on a published data set of K_i values for HIV-1 protease mutants against drugs. Reliable MB-QSAR models were achieved and these models display both well internal and external prediction abilities. Interpreting the MB-QSAR models supplied structural information related to the drug resistance as well as the guidance for the design of resistance evading drugs. This work showed that MB-QSAR method can be employed to predict the resistance of HIV-1 protease caused by polymorphic mutations, which offer a fast and accurate method for the prediction of other drug target within the context of 3D structures.

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.

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Reviewing Antiviral Research Against Viruses Causing Human Diseases - A Structure Guided Approach

The littlest of all the pathogens, viruses have continuously been the foremost strange microorganisms to consider. Viral Infections can cause extreme sicknesses as archived by the HIV/AIDS widespread or the later Ebola or Zika episodes. Apprehensive framework distortions are too regularly watched results of numerous viral contaminations. Besides, numerous infections are oncoviruses, which can trigger different sorts of cancer. Nearly every year a modern infection species rises debilitating the world populace with an annihilating episode. Subsequently, the need of creating antivirals to combat such rising infections. In any case, from the innovation of to begin with antiviral medicate Idoxuridine in 1962 to the revelation of Baloxavir marboxil (Xofluza) that was FDA-approved in 2018, the hone of creating antivirals has changed significantly. In this article, different auxiliary science strategies have been described that can be referral for therapeutics innovation.

Structural basis of substrate specificity in human cytidine deaminase family APOBEC3s

The human cytidine deaminase family of APOBEC3s (A3s) plays critical roles in both innate immunity and the development of cancers. A3s comprise seven functionally overlapping but distinct members that can be exploited as nucleotide base editors for treating genetic diseases. Although overall structurally similar, A3s have vastly varying deamination activity and substrate preferences. Recent crystal structures of ssDNA-bound A3s together with experimental studies have provided some insights into distinct substrate specificities among the family members. However, the molecular interactions responsible for their distinct biological functions and how structure regulates substrate specificity are not clear. In this study, we identified the structural basis of substrate specificities in three catalytically active A3 domains whose crystal structures have been previously characterized: A3A, A3B- CTD, and A3G-CTD. Through molecular modeling and dynamic simulations, we found an interdependency between ssDNA substrate binding conformation and nucleotide sequence specificity. In addition to the U-shaped conformation seen in the crystal structure with the CTC0 motif, A3A can accommodate the CCC0 motif when ssDNA is in a more linear (L) conformation. A3B can also bind both U- and L-shaped ssDNA, unlike A3G, which can stably recognize only linear ssDNA. These varied conformations are stabilized by sequence-specific interactions with active site loops 1 and 7, which are highly variable among A3s. Our results explain the molecular basis of previously observed substrate specificities in A3s and have implications for designing A3-specific inhibitors for cancer therapy as well as engineering base-editing systems for gene therapy.

Report of the National Institutes of Health SARS-CoV-2 Antiviral Therapeutics Summit

The NIH Virtual SARS-CoV-2 Antiviral Summit, held on 6 November 2020, was organized to provide an overview on the status and challenges in developing antiviral therapeutics for coronavirus disease 2019 (COVID-19), including combinations of antivirals. Scientific experts from the public and private sectors convened virtually during a live videocast to discuss severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) targets for drug discovery as well as the preclinical tools needed to develop and evaluate effective small-molecule antivirals. The goals of the Summit were to review the current state of the science, identify unmet research needs, share insights and lessons learned from treating other infectious diseases, identify opportunities for public-private partnerships, and assist the research community in designing and developing antiviral therapeutics. This report includes an overview of therapeutic approaches, individual panel summaries, and a summary of the discussions and perspectives on the challenges ahead for antiviral development.

Deciphering Complex Mechanisms of Resistance and Loss of Potency through Coupled Molecular Dynamics and Machine Learning

Drug resistance threatens many critical therapeutics through mutations in the drug target. The molecular mechanisms by which combinations of mutations, especially those remote from the active site, alter drug binding to confer resistance are poorly understood and thus difficult to counteract. A machine learning strategy was developed that coupled parallel molecular dynamics simulations with experimental potency to identify specific conserved mechanisms underlying resistance. Physical features were extracted from the simulations, analyzed, and integrated into one consistent and interpretable elastic network model. To rigorously test this strategy, HIV-1 protease variants with diverse mutations were used, with potencies ranging from picomolar to micromolar to the drug darunavir. Feature reduction resulted in a model with four specific features that predicts for both the training and test sets inhibitor binding free energy within 1 kcal/mol of the experimental value over this entire range of potency. These predictive features are physically interpretable, as they vary specifically with affinity and diagonally transverse across the protease homodimer. This physics-based strategy of parallel molecular dynamics and machine learning captures mechanisms by which complex combinations of mutations confer resistance and identify critical features that serve as bellwethers of affinity, which will be critical in future drug design.

Acquired HIV-1 Protease Conformational Flexibility Associated with Lopinavir Failure May Shape the Outcome of Darunavir Therapy after Antiretroviral Therapy Switch

Understanding the underlying molecular interaction during a therapy switch from lopinavir (LPV) to darunavir (DRV) is essential to achieve long-term virological suppression. We investigated the kinetic and structural characteristics of multidrug-resistant South African HIV-1 subtype C protease (HIV-1 PR) during therapy switch from LPV to DRV using enzyme activity and inhibition assay, fluorescence spectroscopy, and molecular dynamic simulation. The HIV-1 protease variants were from clinical isolates with a combination of drug resistance mutations; MUT-1 (M46I, I54V, V82A, and L10F), MUT-2 (M46I, I54V, L76V, V82A, L10F, and L33F), and MUT-3 (M46I, I54V, L76V, V82A, L90M, and F53L). Enzyme kinetics analysis shows an association between increased relative resistance to LPV and DRV with the progressive decrease in the mutant HIV-1 PR variants' catalytic efficiency. A direct relationship between high-level resistance to LPV and intermediate resistance to DRV with intrinsic changes in the three-dimensional structure of the mutant HIV-1 PR as a function of the multidrug-resistance mutation was observed. In silico analysis attributed these structural adjustments to the multidrug-resistance mutations affecting the LPV and DRV binding landscape. Though DRV showed superiority to LPV, as a lower concentration was needed to inhibit the HIV-1 PR variants, the inherent structural changes resulting from mutations selected during LPV therapy may dynamically shape the DRV treatment outcome after the therapy switch.

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.

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.

Management of Antiretroviral Therapy with Boosted Protease Inhibitors-Darunavir/Ritonavir or Darunavir/Cobicistat

A major challenge in the management of antiretroviral therapy (ART) is to improve the patient's adherence, reducing the burden caused by the high number of drugs that compose the treatment regimens for human immunodeficiency virus positive (HIV+) patients. Selection of the most appropriate treatment regimen is responsible for therapeutic success and aims to reduce viremia, increase the immune system response capacity, and reduce the incidence rate and intensity of adverse reactions. In general, protease inhibitor (PI) is one of the pillars of regimens, and darunavir (DRV), in particular, is frequently recommended, along with low doses of enzyme inhibitors as cobicistat (COBI) or ritonavir (RTV), by the international guidelines. The potential of clinically significant drug interactions in patients taking COBI or RTV is high due to the potent inhibitory effect on cytochrome CYP 450, which attracts significant changes in the pharmacokinetics of PIs. Regardless of the patient or type of virus, the combined regimens of DRV/COBI or DRV/RTV are available to clinicians, proving their effectiveness, with a major impact on HIV mortality/morbidity. This study presents current information on the pharmacokinetics, pharmacology, drug interactions, and adverse reactions of DRV; it not only compares the bioavailability, pharmacokinetic parameters, immunological and virological responses, but also the efficacy, advantages, and therapeutic disadvantages of DRV/COBI or DRV/RTV combinations.

Crystal Structure of SARS-CoV-2 Main Protease in Complex with the Non-Covalent Inhibitor ML188

Viral proteases are critical enzymes for the maturation of many human pathogenic viruses and thus are key targets for direct acting antivirals (DAAs). The current viral pandemic caused by SARS-CoV-2 is in dire need of DAAs. The Main protease (Mpro) is the focus of extensive structure-based drug design efforts which are mostly covalent inhibitors targeting the catalytic cysteine. ML188 is a non-covalent inhibitor designed to target SARS-CoV-1 Mpro, and provides an initial scaffold for the creation of effective pan-coronavirus inhibitors. In the current study, we found that ML188 inhibits SARS-CoV-2 Mpro at 2.5 µM, which is more potent than against SAR-CoV-1 Mpro. We determined the crystal structure of ML188 in complex with SARS-CoV-2 Mpro to 2.39 Å resolution. Sharing 96% sequence identity, structural comparison of the two complexes only shows subtle differences. Non-covalent protease inhibitors complement the design of covalent inhibitors against SARS-CoV-2 main protease and are critical initial steps in the design of DAAs to treat CoVID 19.

Darunavir-Resistant HIV-1 Protease Constructs Uphold a Conformational Selection Hypothesis for Drug Resistance

Multidrug resistance continues to be a barrier to the effectiveness of highly active antiretroviral therapy in the treatment of human immunodeficiency virus 1 (HIV-1) infection. Darunavir (DRV) is a highly potent protease inhibitor (PI) that is oftentimes effective when drug resistance has emerged against first-generation inhibitors. Resistance to darunavir does evolve and requires 10–20 amino acid substitutions. The conformational landscapes of six highly characterized HIV-1 protease (PR) constructs that harbor up to 19 DRV-associated mutations were characterized by distance measurements with pulsed electron double resonance (PELDOR) paramagnetic resonance spectroscopy, namely double electron–electron resonance (DEER). The results show that the accumulated substitutions alter the conformational landscape compared to PI-naïve protease where the semi-open conformation is destabilized as the dominant population with open-like states becoming prevalent in many cases. A linear correlation is found between values of the DRV inhibition parameter K_i and the open-like to closed-state population ratio determined from DEER. The nearly 50% decrease in occupancy of the semi-open conformation is associated with reduced enzymatic activity, characterized previously in the literature.

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-C₂-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.

Highly drug-resistant HIV-1 protease reveals decreased intra-subunit interactions due to clusters of mutations

Drug-resistance is a serious problem for treatment of the HIV/AIDS pandemic. Potent clinical inhibitors of HIV-1 protease show several orders of magnitude worse inhibition of highly drug-resistant variants. Hence, the structure and enzyme activities were analyzed for HIV protease mutant HIV-1 protease (EC 3.4.23.16) (PR) with 22 mutations (PRS5B) from a clinical isolate that was selected by machine learning to represent high-level drug-resistance. PRS5B has 22 mutations including only one (I84V) in the inhibitor binding site; however, clinical inhibitors had poor inhibition of PRS5B activity with kinetic inhibition value (K_i ) values of 4-1000 nm or 18- to 8000-fold worse than for wild-type PR. High-resolution crystal structures of PRS5B complexes with the best inhibitors, amprenavir (APV) and darunavir (DRV) (K_i ~ 4 nm), revealed only minor changes in protease-inhibitor interactions. Instead, two distinct clusters of mutations in distal regions induce coordinated conformational changes that decrease favorable internal interactions across the entire protein subunit. The largest structural rearrangements are described and compared to other characterized resistant mutants. In the protease hinge region, the N83D mutation eliminates a hydrogen bond connecting the hinge and core of the protease and increases disorder compared to highly resistant mutants PR with 17 mutations and PR with 20 mutations with similar hinge mutations. In a distal β-sheet, mutations G73T and A71V coordinate with accessory mutations to bring about shifts that propagate throughout the subunit. Molecular dynamics simulations of ligand-free dimers show differences consistent with loss of interactions in mutant compared to wild-type PR. Clusters of mutations exhibit both coordinated and antagonistic effects, suggesting PRS5B may represent an intermediate stage in the evolution of more highly resistant variants. DATABASES: Structural data are available in Protein Data Bank under the accession codes 6P9A and 6P9B for PRS5B/DRV and PRS5B/APV, respectively.

Pharmacokinetic Parameters of HIV-1 Protease Inhibitors

Since the beginning of the HIV epidemic, research has been carried out to control the virus. Understanding the mechanisms of replication has given access to the various classes of drugs that over time have transformed AIDS into a manageable chronic disease. The class of protease inhibitors (PIs) gained notice in anti-retroviral therapy, once it was found that peptidomimetic molecules act by blocking the active catalytic center of the aspartic protease, which is directly related to HIV maturation. However, mutations in enzymatic internal residues are the biggest issue for these drugs, because a small change in biochemical interaction can generate resistance. Low plasma concentrations of PIs favor viral natural selection; high concentrations can inhibit even partially resistant enzymes. Food-drug/drug-drug interactions can decrease the bioavailability of PIs and are related to many side effects. Therefore, this review summarizes the pharmacokinetic properties of current PIs, the changes when pharmacological boosters are used and also lists the major mutations to help understanding of how long the continuous treatment can ensure a low viral load in patients.

Optimizing the refinement of merohedrally twinned P61 HIV-1 protease-inhibitor cocrystal structures

Twinning is a crystal-growth anomaly in which protein monomers exist in different orientations but are related in a specific way, causing diffraction reflections to overlap. Twinning imposes additional symmetry on the data, often leading to the assignment of a higher symmetry space group. Specifically, in merohedral twinning, reflections from each monomer overlap and require a twin law to model unique structural data from overlapping reflections. Neglecting twinning in the crystallographic analysis of quasi-rotationally symmetric homo-oligomeric protein structures can mask the degree of structural non-identity between monomers. In particular, any deviations from perfect symmetry will be lost if higher than appropriate symmetry is applied during crystallographic analysis. Such cases warrant choosing between the highest symmetry space group possible or determining whether the monomers have distinguishable structural asymmetries and thus require a lower symmetry space group and a twin law. Using hexagonal cocrystals of HIV-1 protease, a C2-symmetric homodimer whose symmetry is broken by bound ligand, it is shown that both assigning a lower symmetry space group and applying a twin law during refinement are critical to achieving a structural model that more accurately fits the electron density. By re-analyzing three recently published HIV-1 protease structures, improvements in nearly every crystallographic metric are demonstrated. Most importantly, a procedure is demonstrated where the inhibitor can be reliably modeled in a single orientation. This protocol may be applicable to many other homo-oligomers in the PDB.