Inhibition of autoprocessing of natural variants and multidrug resistant mutant precursors of HIV-1 protease by clinical inhibitors

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro^WT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373) and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decreases in free energy (ΔG) and thermal stability (T_m), respectively, relative to MPro^WT. The H41A mutation results in a 20-fold increase in the dissociation constant (K_d), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T_m, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro^H41A, whereas no significant change is observed in binding to MPro^WT. Structures of the four inhibitor complexes with MPro^1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro^WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro^WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro^C145A also results in a significant, but smaller, increase in K_d and decrease in ΔG and T_m, relative to NMV.

Unmasking the Conformational Stability and Inhibitor Binding to SARS-CoV-2 Main Protease Active Site Mutants and Miniprecursor

We recently demonstrated that inhibitor binding reorganizes the oxyanion loop of a monomeric catalytic domain of SARS CoV-2 main protease (MPro) from an unwound (E) to a wound (active, E*) conformation, independent of dimerization. Here we assess the effect of the flanking N-terminal residues, to imitate the MPro precursor prior to its autoprocessing, on conformational equilibria rendering stability and inhibitor binding. Thermal denaturation (T_m) of C145A mutant, unlike H41A, increases by 6.8 °C, relative to wild-type mature dimer. An inactivating H41A mutation to maintain a miniprecursor containing TSAVL[Q or E] of the flanking nsp4 sequence in an intact form [^(-6)MPro^H41A and ^(-6*)MPro^H41A, respectively], and its corresponding mature MPro^H41A were systematically examined. While the H41A mutation exerts negligible effect on T_m and dimer dissociation constant (K_dimer) of MPro^H41A, relative to the wild type MPro, both miniprecursors show a 4-5 °C decrease in T_m and > 85-fold increase in K_dimer as compared to MPro^H41A. The K_d for the binding of the covalent inhibitor GC373 to ^(-6*)MPro^H41A increases ∼12-fold, relative to MPro^H41A, concomitant with its dimerization. While the inhibitor-free dimer exhibits a state in transit from E to E* with a conformational asymmetry of the protomers' oxyanion loops and helical domains, inhibitor binding restores the asymmetry to mature-like oxyanion loop conformations (E*) but not of the helical domains. Disorder of the terminal residues 1-2 and 302-306 observed in both structures suggest that N-terminal autoprocessing is tightly coupled to the E-E* equilibrium and stable dimer formation.

Viral proteases as therapeutic targets

Some medically important viruses-including retroviruses, flaviviruses, coronaviruses, and herpesviruses-code for a protease, which is indispensable for viral maturation and pathogenesis. Viral protease inhibitors have become an important class of antiviral drugs. Development of the first-in-class viral protease inhibitor saquinavir, which targets HIV protease, started a new era in the treatment of chronic viral diseases. Combining several drugs that target different steps of the viral life cycle enables use of lower doses of individual drugs (and thereby reduction of potential side effects, which frequently occur during long term therapy) and reduces drug-resistance development. Currently, several HIV and HCV protease inhibitors are routinely used in clinical practice. In addition, a drug including an inhibitor of SARS-CoV-2 main protease, nirmatrelvir (co-administered with a pharmacokinetic booster ritonavir as Paxlovid®), was recently authorized for emergency use. This review summarizes the basic features of the proteases of human immunodeficiency virus (HIV), hepatitis C virus (HCV), and SARS-CoV-2 and discusses the properties of their inhibitors in clinical use, as well as development of compounds in the pipeline.

Cryo-EM structure of the HIV-1 Pol polyprotein provides insights into virion maturation

Key proteins of retroviruses and other RNA viruses are translated and subsequently processed from polyprotein precursors by the viral protease (PR). Processing of the HIV Gag-Pol polyprotein yields the HIV structural proteins and enzymes. Structures of the mature enzymes PR, reverse transcriptase (RT), and integrase (IN) aided understanding of catalysis and design of antiretrovirals, but knowledge of the Pol precursor architecture and function before PR cleavage is limited. We developed a system to produce stable HIV-1 Pol and determined its cryo-electron microscopy structure. RT in Pol has a similar arrangement to the mature RT heterodimer, and its dimerization may draw together two PR monomers to activate proteolytic processing. HIV-1 thus may leverage the dimerization interfaces in Pol to regulate assembly and maturation of polyprotein precursors.

A kind of HIV-1 protease inhibitors containing phenols with antiviral activity against DRV-resistant variants

Based upon the preliminary design of enhancing genetic barrier to drug-resistant viral mutants by maximizing hydrogen-bonding or other van der Waals contacts, we have designed, synthesized and biologically evaluated a new class of HIV-1 protease inhibitors with phenol derived P2 ligands and nitro or halogens in P2' ligands. Results indicate that a majority of inhibitors exhibit robust enzyme inhibitory with IC₅₀ values in picomolar or single digit nanomolar ranges. Among which, compound 17d displays potency with IC₅₀ value of 21 pM and high protease selectivity. Of note, 17d exhibits greater antiviral activity against the DRV-resistant variant than the efficacy against the wild type virus. Furthermore, the molecular modeling studies demonstrate important interactions between 17d and the active sites of both the wild-type and DRV-resistant HIV-1 protease, as well as furnish insights for further optimization of new inhibitors.

Precursors of Viral Proteases as Distinct Drug Targets

Viral proteases are indispensable for successful virion maturation, thus making them a prominent drug target. Their enzyme activity is tightly spatiotemporally regulated by expression in the precursor form with little or no activity, followed by activation via autoprocessing. These cleavage events are frequently triggered upon transportation to a specific compartment inside the host cell. Typically, precursor oligomerization or the presence of a co-factor is needed for activation. A detailed understanding of these mechanisms will allow ligands with non-canonical mechanisms of action to be designed, which would specifically modulate the initial irreversible steps of viral protease autoactivation. Binding sites exclusive to the precursor, including binding sites beyond the protease domain, can be exploited. Both inhibition and up-regulation of the proteolytic activity of viral proteases can be detrimental for the virus. All these possibilities are discussed using examples of medically relevant viruses including herpesviruses, adenoviruses, retroviruses, picornaviruses, caliciviruses, togaviruses, flaviviruses, and coronaviruses.

Novel HIV PR inhibitors with C4-substituted bis-THF and bis-fluoro-benzyl target the two active site mutations of highly drug resistant mutant PRS17

The emergence of multidrug resistant (MDR) HIV strains severely reduces the effectiveness of antiretroviral therapy. Clinical inhibitor darunavir (1) has picomolar binding affinity for HIV-1 protease (PR), however, drug resistant variants like PR^S17 show poor inhibition by 1, despite the presence of only two mutated residues in the inhibitor-binding site. Antiviral inhibitors that target MDR proteases like PR^S17 would be valuable as therapeutic agents. Inhibitors 2 and 3 derived from 1 through substitutions at P1, P2 and P2' positions exhibit 3.4- to 500-fold better inhibition than clinical inhibitors for PR^S17 with the exception of amprenavir. Crystal structures of PR^S17/2 and PR^S17/3 reveal how these inhibitors target the two active site mutations of PR^S17. The substituted methoxy P2 group of 2 forms new interactions with G48V mutation, while the modified bis-fluoro-benzyl P1 group of 3 forms a halogen interaction with V82S mutation, contributing to improved inhibition of PR^S17.

HIV Protease: Historical Perspective and Current Research

The retroviral protease of human immunodeficiency virus (HIV) is an excellent target for antiviral inhibitors for treating HIV/AIDS. Despite the efficacy of therapy, current efforts to control the disease are undermined by the growing threat posed by drug resistance. This review covers the historical background of studies on the structure and function of HIV protease, the subsequent development of antiviral inhibitors, and recent studies on drug-resistant protease variants. We highlight the important contributions of Dr. Stephen Oroszlan to fundamental knowledge about the function of the HIV protease and other retroviral proteases. These studies, along with those of his colleagues, laid the foundations for the design of clinical inhibitors of HIV protease. The drug-resistant protease variants also provide an excellent model for investigating the molecular mechanisms and evolution of resistance.

Nanoscale flow cytometry reveals interpatient variability in HIV protease activity that correlates with viral infectivity and identifies drug-resistant viruses

HIV encodes an aspartyl protease that is activated during, or shortly after, budding of viral particles from the surface of infected cells. Protease-mediated cleavage of viral polyproteins is essential to generating infectious viruses, a process known as ‘maturation’ that is the target of FDA-approved antiretroviral drugs. Most assays to monitor protease activity rely on bulk analysis of millions of viruses and obscure potential heterogeneity of protease activation within individual particles. In this study we used nanoscale flow cytometry in conjunction with an engineered FRET reporter called VIral ProteasE Reporter (VIPER) to investigate heterogeneity of protease activation in individual, patient-derived viruses. We demonstrate previously unappreciated interpatient variation in HIV protease processing efficiency that impacts viral infectivity. Additionally, monitoring of protease activity in individual virions distinguishes between drug sensitivity or resistance to protease inhibitors in patient-derived samples. These findings demonstrate the feasibility of monitoring enzymatic processes using nanoscale flow cytometry and highlight the potential of this technology for translational clinical discovery, not only for viruses but also other submicron particles including exosomes, microvesicles, and bacteria.

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.

Visualizing Tetrahedral Oxyanion Bound in HIV-1 Protease Using Neutrons: Implications for the Catalytic Mechanism and Drug Design

HIV-1 protease is indispensable for virus propagation and an important therapeutic target for antiviral inhibitors to treat AIDS. As such inhibitors are transition-state mimics, a detailed understanding of the enzyme mechanism is crucial for the development of better anti-HIV drugs. Here, we used room-temperature joint X-ray/neutron crystallography to directly visualize hydrogen atoms and map hydrogen bonding interactions in a protease complex with peptidomimetic inhibitor KVS-1 containing a reactive nonhydrolyzable ketomethylene isostere, which, upon reacting with the catalytic water molecule, is converted into a tetrahedral intermediate state, KVS-1_TI. We unambiguously determined that the resulting tetrahedral intermediate is an oxyanion, rather than the gem-diol, and both catalytic aspartic acid residues are protonated. The oxyanion tetrahedral intermediate appears to be unstable, even though the negative charge on the oxyanion is delocalized through a strong n → π* hyperconjugative interaction into the nearby peptidic carbonyl group of the inhibitor. To better understand the influence of the ketomethylene isostere as a protease inhibitor, we have also examined the protease structure and binding affinity with keto-darunavir (keto-DRV), which similar to KVS-1 includes the ketomethylene isostere. We show that keto-DRV is a significantly less potent protease inhibitor than DRV. These findings shed light on the reaction mechanism of peptide hydrolysis catalyzed by HIV-1 protease and provide valuable insights into further improvements in the design of protease inhibitors.

Protease Inhibitors for the Treatment of HIV/AIDS: Recent Advances and Future Challenges

Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS.

Potent antiviral HIV-1 protease inhibitor combats highly drug resistant mutant PR20

Drug-resistance threatens effective treatment of HIV/AIDS. Clinical inhibitors, including darunavir (1), are ineffective for highly resistant protease mutant PR20, however, antiviral compound 2 derived from 1 with fused tricyclic group at P2, extended amino-benzothiazole P2' ligand and two fluorine atoms on P1 shows 16-fold better inhibition of PR20 enzyme activity. Crystal structures of PR20 and wild-type PR complexes reveal how the extra groups of 2 counteract the expanded ligand-binding pocket, dynamic flaps, and faster dimer dissociation of PR20.

Diverse Folding Pathways of HIV-1 Protease Monomer on a Rugged Energy Landscape

The modern energy landscape theory of protein folding predicts multiple folding pathways connecting a myriad of unfolded conformations and a well-defined folded state. However, direct experimental observation of heterogeneous folding pathways is difficult. Naturally evolved proteins typically exhibit a smooth folding energy landscape for fast and efficient folding by avoiding unfavorable kinetic traps. In this case, rapid fluctuations between unfolded conformations result in apparent two-state behavior and make different pathways indistinguishable. However, the landscape roughness can be different, depending on the selection pressures during evolution. Here, we characterize the unusually rugged folding energy landscape of human immunodeficiency virus-1 protease monomer using single-molecule Förster resonance energy transfer spectroscopy. Our data show that fluctuations between unfolded conformations are slow, which enables the experimental observation of heterogeneous folding pathways as predicted by the landscape theory. Although the landscape ruggedness is sensitive to the mutations and fluorophore locations, the folding rate is similar for various protease constructs. The natural evolution of the protease to have a rugged energy landscape likely results from intrinsic pressures to maintain robust folding when human immunodeficiency virus-1 mutates frequently, which is essential for its survival.

Structural and binding insights into HIV-1 protease and P2-ligand interactions through molecular dynamics simulations, binding free energy and principal component analysis

HIV-1 protease (HIV-1-pr) plays an important role in viral replication and maturation, making it one of the most attractive targets for anti-retroviral therapy. To design new effective inhibitors able to combat drug resistance in mutant HIV-1-pr variants, it is essential to gain further understanding about the mechanisms by which the recently proposed inhibitors deactivate the mutant HIV-1-pr variants. In the present work, we explored the interactions between two P2-ligands (DRV, and one new derivative, 4UY) with wild type (WT) and two multiple mutant HIV-1-pr variants (p20 and p51) with all atom molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA). The trajectories of MD simulations show that both 4UY and DRV primarily bind with the active sites, flap and 80s loop regions of HIV-1-pr variants through either hydrogen bonds or hydrophobic interactions. More hydrogen bonds and hydrophobic interactions were located for 4UY/HIV-1-pr complexes than for DRV/HIV-1-pr counterparts. More importantly, 4UY was found to have an extra hydrogen bond with the backbone of Gly48' in the flap region of the HIV-1-prs. The flap tip-tip distance (I50-I50') and flap tip-active site distance (I50-D25 and I50'-D25') indicate that the flaps turn more closed in 4UY bound HIV-1-prs than DRV bound ones, and the former also have more compact hydrophobic cavities than the latter. Further, the vector projections from PCA indicate that 4UY/DRV inhibitor binding projects the closing of flap in HIV-1-pr variants. In line with the above trajectory analysis, the thermodynamics calculation with MM-PBSA method suggests much stronger binding affinity for 4UY/HIV-1-pr than DRV/HIV-1-pr by 4.3-6.4 kcal/mol. Although p20 and p51 also induce weaker binding due to multiple mutants for 4UY inhibitor by 1.9-1.8 kcal/mol, their bindings to the new P2 ligand (4UY) are indeed significantly enhanced as compared to DRV. The thermodynamic components responsible for the binding differences and the contribution from key residues to the binding were also discussed in detail.

Inhibition of HIV Maturation via Selective Unfolding and Cross-Linking of Gag Polyprotein by a Mercaptobenzamide Acetylator

For HIV to become infectious, any new virion produced from an infected cell must undergo a maturation process that involves the assembly of viral polyproteins Gag and Gag-Pol at the membrane surface. The self-assembly of these viral proteins drives formation of a new viral particle as well as the activation of HIV protease, which is needed to cleave the polyproteins so that the final core structure of the virus will properly form. Molecules that interfere with HIV maturation will prevent any new virions from infecting additional cells. In this manuscript, we characterize the unique mechanism by which a mercaptobenzamide thioester small molecule (SAMT-247) interferes with HIV maturation via a series of selective acetylations at highly conserved cysteine and lysine residues in Gag and Gag-Pol polyproteins. The results provide the first insights into how acetylation can be utilized to perturb the process of HIV maturation and reveal a new strategy to limit the infectivity of HIV.

Exploiting the unique features of Zika and Dengue proteases for inhibitor design

Zika and Dengue viruses have attracted substantial attention from researchers in light of recent outbreaks of Dengue fever and increases in cases of congenital microcephaly in areas with Zika incidence. This review summarizes the current state of knowledge about Zika and Dengue proteases. These enzymes have several interesting features: 1) NS3 serine protease requires the activating co-factor NS2B, which is anchored in the membrane of the endoplasmic reticulum; 2) NS2B displays extensive conformational dynamics; 3) NS3 is a multidomain protein with proteolytic, NTPase, RNA 5' triphosphatase and helicase activity and has many protein-protein interaction partners; 4) NS3 is autoproteolytically released from its precursor. Attempts to design tight-binding and specific active-site inhibitors are complicated by the facts that the substrate pocket of the NS2B-NS3 protease is flat and the active-site ligands are charged. The ionic character of potential active-site inhibitors negatively influences their cell permeability. Possibilities to block cis-autoprocessing of the protease precursor have recently been considered. Additionally, potential allosteric sites on NS2B-NS3 proteases have been identified and allosteric compounds have been designed to impair substrate binding and/or block the NS2B-NS3 interaction. Such compounds could be specific to viral proteases, without off-target effects on host serine proteases, and could have favorable pharmacokinetic profiles. This review discusses various groups of inhibitors of these proteases according to their mechanisms of action and chemical structures.

Targeting HIV-1 Protease Autoprocessing for High-throughput Drug Discovery and Drug Resistance Assessment

HIV-1 protease autoprocessing liberates the free mature protease from its Gag-Pol polyprotein precursor through a series of highly regulated autoproteolysis reactions. Herein, we report the development and validation (Z' ≥ 0.50) of a cell-based functional assay for high-throughput screening (HTS) of autoprocessing inhibitors using fusion precursors in combination with AlphaLISA (amplified luminescent proximity homogeneous assay ELISA). Through pilot screening of a collection of 130 known protease inhibitors, the AlphaLISA assay confirmed all 11 HIV protease inhibitors in the library capable of suppressing precursor autoprocessing at low micromolar concentrations. Meanwhile, other protease inhibitors had no impact on precursor autoprocessing. We next conducted HTS of ~23,000 compounds but found no positive hits. Such high selectivity is advantageous for large-scale HTS campaigns and as anticipated based on assay design because a positive hit needs simultaneously to be nontoxic, cell permeable, and inhibiting precursor autoprocessing. Furthermore, AlphaLISA quantification of fusion precursors carrying mutations known to cause resistance to HIV protease inhibitors faithfully recapitulated the reported resistance, suggesting that precursor autoprocessing is a critical step contributing to drug resistance. Taken together, this reported AlphaLISA platform will provide a useful tool for drug discovery targeting HIV-1 protease autoprocessing and for quantification of PI resistance.

Inhibition of the precursor and mature forms of HIV-1 protease as a tool for drug evaluation

HIV-1 protease (PR) is a homodimeric enzyme that is autocatalytically cleaved from the Gag-Pol precursor. Known PR inhibitors bind the mature enzyme several orders of magnitude more strongly than the PR precursor. Inhibition of PR at the precursor level, however, may stop the process at its rate-limiting step before the proteolytic cascade is initiated. Due to its structural heterogeneity, limited solubility and autoprocessing, the PR precursor is difficult to access by classical methods, and limited knowledge regarding precursor inhibition is available. Here, we describe a cell-based assay addressing precursor inhibition. We used a reporter molecule containing the transframe (TFP) and p6* peptides, PR, and N-terminal fragment of reverse transcriptase flanked by the fluorescent proteins mCherry and EGFP on its N- and C- termini, respectively. The level of FRET between EGFP and mCherry indicates the amount of unprocessed reporter, allowing specific monitoring of precursor inhibition. The inhibition can be quantified by flow cytometry. Additionally, two microscopy techniques confirmed that the reporter remains unprocessed within individual cells upon inhibition. We tested darunavir, atazanavir and nelfinavir and their combinations against wild-type PR. Shedding light on an inhibitor’s ability to act on non-mature forms of PR may aid novel strategies for next-generation drug design.

Context-dependent autoprocessing of human immunodeficiency virus type 1 protease precursors

HIV-1 protease autoprocessing is responsible for liberation of free mature protease (PR) from the Gag-Pol polyprotein precursor. A cell-based model system was previously developed to examine the autoprocessing mechanism of fusion precursors carrying the p6*-PR miniprecursor sandwiched between various proteins or epitopes. We here report that precursor autoprocessing is context-dependent as its activity and outcomes can be modulated by sequences upstream of p6*-PR. This was exemplified by the 26aa maltose binding protein (MBP) signal peptide (SigP) when placed at the N-terminus of a fusion precursor. The mature PRs released from SigP-carrying precursors are resistant to self-degradation whereas those released from SigP-lacking fusion precursors are prone to self-degradation. A H69D mutation in PR abolished autoprocessing of SigP-containing fusion precursors whereas it only partially suppressed autoprocessing of fusion precursors lacking SigP. An autoprocessing deficient GFP fusion precursor with SigP exhibited a subcellular distribution pattern distinct from the one without it in transfected HeLa cells. Furthermore, a SigP fusion precursor carrying a substitution at the P1 position released the mature PR and PR-containing fragments that were different from those released from the precursor carrying the same mutation but lacking SigP. We also examined autoprocessing outcomes in viral particles produced by a NL4-3 derived proviral construct and demonstrated the existence of several PR-containing fragments along with the mature PR. Some of these resembled the SigP precursor autoprocessing outcomes. This finding of context-dependent modulation reveals the complexity of precursor autoprocessing regulation that most likely accompanies sequence variation imposed by the evolution of the upstream Gag moiety.

Binding kinetics and substrate selectivity in HIV-1 protease-Gag interactions probed at atomic resolution by chemical exchange NMR

The conversion of immature noninfectious HIV-1 particles to infectious virions is dependent upon the sequential cleavage of the precursor group-specific antigen (Gag) polyprotein by HIV-1 protease. The precise mechanism whereby protease recognizes distinct Gag cleavage sites, located in the intrinsically disordered linkers connecting the globular domains of Gag, remains unclear. Here, we probe the dynamics of the interaction of large fragments of Gag and various variants of protease (including a drug resistant construct) using Carr-Purcell-Meiboom-Gill relaxation dispersion and chemical exchange saturation transfer NMR experiments. We show that the conformational dynamics within the flaps of HIV-1 protease that form the lid over the catalytic cleft play a significant role in substrate specificity and ordered Gag processing. Rapid interconversion between closed and open protease flap conformations facilitates the formation of a transient, sparsely populated productive complex between protease and Gag substrates. Flap closure traps the Gag cleavage sites within the catalytic cleft of protease. Modulation of flap opening through protease-Gag interactions fine-tunes the lifetime of the productive complex and hence the likelihood of Gag proteolysis. A productive complex can also be formed in the presence of a noncognate substrate but is short-lived owing to lack of optimal complementarity between the active site cleft of protease and the substrate, resulting in rapid flap opening and substrate release, thereby allowing protease to differentiate between cognate and noncognate substrates.

Decoding HIV resistance: from genotype to therapy

Genetic variation in HIV poses a major challenge for prevention and treatment of the AIDS pandemic. Resistance occurs by mutations in the target proteins that lower affinity for the drug or alter the protein dynamics, thereby enabling viral replication in the presence of the drug. Due to the prevalence of drug-resistant strains, monitoring the genotype of the infecting virus is recommended. Computational approaches for predicting resistance from genotype data and guiding therapy are discussed. Many prediction methods rely on rules derived from known resistance-associated mutations, however, statistical or machine learning can improve the classification accuracy and assess unknown mutations. Adding classifiers such as information on the atomic structure of the protein can further enhance the predictions.

Inference of Epistatic Effects Leading to Entrenchment and Drug Resistance in HIV-1 Protease

Understanding the complex mutation patterns that give rise to drug resistant viral strains provides a foundation for developing more effective treatment strategies for HIV/AIDS. Multiple sequence alignments of drug-experienced HIV-1 protease sequences contain networks of many pair correlations which can be used to build a (Potts) Hamiltonian model of these mutation patterns. Using this Hamiltonian model, we translate HIV-1 protease sequence covariation data into quantitative predictions for the probability of observing specific mutation patterns which are in agreement with the observed sequence statistics. We find that the statistical energies of the Potts model are correlated with the fitness of individual proteins containing therapy-associated mutations as estimated by in vitro measurements of protein stability and viral infectivity. We show that the penalty for acquiring primary resistance mutations depends on the epistatic interactions with the sequence background. Primary mutations which lead to drug resistance can become highly advantageous (or entrenched) by the complex mutation patterns which arise in response to drug therapy despite being destabilizing in the wildtype background. Anticipating epistatic effects is important for the design of future protease inhibitor therapies.

Structural Studies of a Rationally Selected Multi-Drug Resistant HIV-1 Protease Reveal Synergistic Effect of Distal Mutations on Flap Dynamics

We report structural analysis of HIV protease variant PR^S17 which was rationally selected by machine learning to represent wide classes of highly drug-resistant variants. Crystal structures were solved of PR^S17 in the inhibitor-free form and in complex with antiviral inhibitor, darunavir. Despite its 17 mutations, PR^S17 has only one mutation (V82S) in the inhibitor/substrate binding cavity, yet exhibits high resistance to all clinical inhibitors. PR^S17 has none of the major mutations (I47V, I50V, I54ML, L76V and I84V) associated with darunavir resistance, but has 10,000-fold weaker binding affinity relative to the wild type PR. Comparable binding affinity of 8000-fold weaker than PR is seen for drug resistant mutant PR20, which bears 3 mutations associated with major resistance to darunavir (I47V, I54L and I84V). Inhibitor-free PR^S17 shows an open flap conformation with a curled tip correlating with G48V flap mutation. NMR studies on inactive PR^S17_D25N unambiguously confirm that the flaps adopt mainly an open conformation in solution very similar to that in the inhibitor-free crystal structure. In PR^S17, the hinge loop cluster of mutations, E35D, M36I and S37D, contributes to the altered flap dynamics by a mechanism similar to that of PR20. An additional K20R mutation anchors an altered conformation of the hinge loop. Flap mutations M46L and G48V in PR^S17/DRV complex alter the Phe53 conformation by steric hindrance between the side chains. Unlike the L10F mutation in PR20, L10I in PR^S17 does not break the inter-subunit ion pair or diminish the dimer stability, consistent with a very low dimer dissociation constant comparable to that of wild type PR. Distal mutations A71V, L90M and I93L propagate alterations to the catalytic site of PR^S17. PR^S17 exhibits a molecular mechanism whereby mutations act synergistically to alter the flap dynamics resulting in significantly weaker binding yet maintaining active site contacts with darunavir.

Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution

Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.

Discovery of novel anti-HIV agents via Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry-based approach

INTRODUCTION

In recent years, a variety of new synthetic methodologies and concepts have been proposed in the search for new pharmaceutical lead structures and optimization. Notably, the Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry approach has drawn great attention and has become a powerful tool for the generation of privileged medicinal skeletons in the discovery of anti-HIV agents. This is due to the high degree of reliability, complete specificity (chemoselectivity and regioselectivity), mild conditions, and the biocompatibility of the reactants.

AREAS COVERED

Herein, the authors describe the progress thus far on the discovery of novel anti-HIV agents via the CuAAC click chemistry-based approach.

EXPERT OPINION

CuAAC click chemistry is a proven protocol for synthesizing triazole products which could serve as basic pharmacophores, act as replacements of traditional scaffold or substituent modification, be a linker of dual-target or dual-site inhibitors and more for the discovery of novel anti-HIV agents. What's more, it also provides convenience and feasibility for dynamic combinatorial chemistry and in situ screening. It is envisioned that click chemistry will draw more attention and make more contributions in anti-HIV drug discovery in the future.

Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease

Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.

Binding of Clinical Inhibitors to a Model Precursor of a Rationally Selected Multidrug Resistant HIV-1 Protease Is Significantly Weaker Than That to the Released Mature Enzyme

We have systematically validated the activity and inhibition of a HIV-1 protease (PR) variant bearing 17 mutations (PR(S17)), selected to represent high resistance by machine learning on genotype-phenotype data. Three of five mutations in PR(S17) correlating with major drug resistance, M46L, G48V, and V82S, and five of 11 natural variations differ from the mutations in two clinically derived extreme mutants, PR20 and PR22 bearing 19 and 22 mutations, respectively. PR(S17), which forms a stable dimer (<10 nM), is ∼10- and 2-fold less efficient in processing the Gag polyprotein than the wild type and PR20, respectively, but maintains the same cleavage order. Isolation of a model precursor of PR(S17) flanked by the 56-amino acid transframe region (TFP-p6pol) at its N-terminus, which is impossible upon expression of an analogous PR20 precursor, allowed systematic comparison of inhibition of TFP-p6pol-PR(S17) and mature PR(S17). Resistance of PR(S17) to eight protease inhibitors (PIs) relative to PR (Ki) increases by 1.5-5 orders of magnitude from 0.01 to 8.4 μM. Amprenavir, darunavir, atazanavir, and lopinavir, the most effective of the eight PIs, inhibit precursor autoprocessing at the p6pol/PR site with IC50 values ranging from ∼7.5 to 60 μM. Thus, this process, crucial for stable dimer formation, shows inhibition ∼200-800-fold weaker than that of the mature PR(S17). TFP/p6pol cleavage, which occurs faster, is inhibited even more weakly by all PIs except darunavir (IC50 = 15 μM); amprenavir shows a 2-fold increase in IC50 (∼15 μM), and atazanavir and lopinavir show increased IC50 values of >42 and >70 μM, respectively.

Highly resistant HIV-1 proteases and strategies for their inhibition

The virally encoded protease is an important drug target for AIDS therapy. Despite the potency of the current drugs, infections with resistant viral strains limit the long-term effectiveness of therapy. Highly resistant variants of HIV protease from clinical isolates have different combinations of about 20 mutations and several orders of magnitude worse binding affinity for clinical inhibitors. Strategies are being explored to inhibit these highly resistant mutants. The existing inhibitors can be modified by introducing groups with the potential to form new interactions with conserved protease residues, and the flexible flaps. Alternative strategies are discussed, including designing inhibitors to bind to the open conformation of the protease dimer, and inhibition of the protease-catalyzed processing of the Gag-Pol precursor.

Flap Dynamics in Aspartic Proteases: A Computational Perspective

Recent advances in biochemistry and drug design have placed proteases as one of the critical target groups for developing novel small-molecule inhibitors. Among all proteases, aspartic proteases have gained significant attention due to their role in HIV/AIDS, malaria, Alzheimer's disease, etc. The binding cleft is covered by one or two β-hairpins (flaps) which need to be opened before a ligand can bind. After binding, the flaps close to retain the ligand in the active site. Development of computational tools has improved our understanding of flap dynamics and its role in ligand recognition. In the past decade, several computational approaches, for example molecular dynamics (MD) simulations, coarse-grained simulations, replica-exchange molecular dynamics (REMD) and metadynamics, have been used to understand flap dynamics and conformational motions associated with flap movements. This review is intended to summarize the computational progress towards understanding the flap dynamics of proteases and to be a reference for future studies in this field.

Dynamical Network of HIV-1 Protease Mutants Reveals the Mechanism of Drug Resistance and Unhindered Activity

HIV-1 protease variants resist drugs by active and non-active-site mutations. The active-site mutations, which are the primary or first set of mutations, hamper the stability of the enzyme and resist the drugs minimally. As a result, secondary mutations that not only increase protein stability for unhindered catalytic activity but also resist drugs very effectively arise. While the mechanism of drug resistance of the active-site mutations is through modulating the active-site pocket volume, the mechanism of drug resistance of the non-active-site mutations is unclear. Moreover, how these allosteric mutations, which are 8-21 Å distant, communicate to the active site for drug efflux is completely unexplored. Results from molecular dynamics simulations suggest that the primary mechanism of drug resistance of the secondary mutations involves opening of the flexible protease flaps. Results from both residue- and community-based network analyses reveal that this precise action of protease is accomplished by the presence of robust communication paths between the mutational sites and the functionally relevant regions: active site and flaps. While the communication is more direct in the wild type, it traverses across multiple intermediate residues in mutants, leading to weak signaling and unregulated motions of flaps. The global integrity of the protease network is, however, maintained through the neighboring residues, which exhibit high degrees of conservation, consistent with clinical data and mutagenesis studies.

Identifying representative drug resistant mutants of HIV

BACKGROUND

Drug resistance is one of the most important causes for failure of anti-AIDS treatment. During therapy, multiple mutations accumulate in the HIV genome, eventually rendering the drugs ineffective in blocking replication of the mutant virus. The huge number of possible mutants precludes experimental analysis to explore the molecular mechanisms of resistance and develop improved antiviral drugs.

RESULTS

In order to solve this problem, we have developed a new algorithm to reveal the most representative mutants from the whole drug resistant mutant database based on our newly proposed unified protein sequence and 3D structure encoding method. Mean shift clustering and multiple regression analysis were applied on genotype-resistance data for mutants of HIV protease and reverse transcriptase. This approach successfully chooses less than 100 mutants with the highest resistance to each drug out of about 10K in the whole database. When considering high level resistance to multiple drugs, the numbers reduce to one or two representative mutants.

CONCLUSION

This approach for predicting the most representative mutants for each drug has major importance for experimental verification since the results provide a small number of representative sequences, which will be amenable for in vitro testing and characterization of the expressed mutant proteins.

Pressure-induced structural transition of mature HIV-1 protease from a combined NMR/MD simulation approach

We investigate the pressure-induced structural changes in the mature human immunodeficiency virus type 1 protease dimer, using residual dipolar coupling (RDC) measurements in a weakly oriented solution. (1)DNH RDCs were measured under high-pressure conditions for an inhibitor-free PR and an inhibitor-bound complex, as well as for an inhibitor-free multidrug resistant protease bearing 20 mutations (PR20). While PR20 and the inhibitor-bound PR were little affected by pressure, inhibitor-free PR showed significant differences in the RDCs measured at 600 bar compared with 1 bar. The structural basis of such changes was investigated by MD simulations using the experimental RDC restraints, revealing substantial conformational perturbations, specifically a partial opening of the flaps and the penetration of water molecules into the hydrophobic core of the subunits at high pressure. This study highlights the exquisite sensitivity of RDCs to pressure-induced conformational changes and illustrates how RDCs combined with MD simulations can be used to determine the structural properties of metastable intermediate states on the folding energy landscape.

Conformational variation of an extreme drug resistant mutant of HIV protease

Molecular mechanisms leading to high level drug resistance have been analyzed for the clinical variant of HIV-1 protease bearing 20 mutations (PR20); which has several orders of magnitude worse affinity for tested drugs. Two crystal structures of ligand-free PR20 with the D25N mutation of the catalytic aspartate (PR20D25N) revealed three dimers with different flap conformations. The diverse conformations of PR20D25N included a dimer with one flap in a unique "tucked" conformation; directed into the active site. Analysis of molecular dynamics (MD) simulations of the ligand-free PR20 and wild-type enzymes showed that the mutations in PR20 alter the correlated interactions between two monomers in the dimer. The two flaps tend to fluctuate more independently in PR20 than in the wild type enzyme. Combining the results of structural analysis by X-ray crystallography and MD simulations; unusual flap conformations and weakly correlated inter-subunit motions may contribute to the high level resistance of PR20.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.

Substituted Bis-THF Protease Inhibitors with Improved Potency against Highly Resistant Mature HIV-1 Protease PR20

An extremely drug resistant mutant of HIV-1 protease (PR) bearing 20 mutations (PR20) has been studied with two potent antiviral investigational inhibitors. GRL-5010A and GRL-4410A were designed to introduce hydrogen bond interactions with the flexible flaps of the PR by incorporating gem-difluorines and alkoxy, respectively, at the C4 position of the bis-THF of darunavir. PR20 provides an excellent model for high level resistance, since clinical inhibitors are >1000-fold less active on PR20 than on wild-type enzyme. GRL-5010A and GRL-4410A show inhibition constants of 4.3 ± 7.0 and 1.7 ± 1.8 nM, respectively, for PR20, compared to the binding affinity of 41 ± 1 nM measured for darunavir. Crystal structures of PR20 in complexes with the two inhibitors confirmed the new hydrogen bond interactions with Gly 48 in the flap of the enzyme. The two new compounds are more effective than darunavir in inhibiting mature PR20 and show promise for further development of antiviral agents targeting highly resistant PR mutants.

A Functional Interplay between Human Immunodeficiency Virus Type 1 Protease Residues 77 and 93 Involved in Differential Regulation of Precursor Autoprocessing and Mature Protease Activity

HIV-1 protease (PR) is a viral enzyme vital to the production of infectious virions. It is initially synthesized as part of the Gag-Pol polyprotein precursor in the infected cell. The free mature PR is liberated as a result of precursor autoprocessing upon virion release. We previously described a model system to examine autoprocessing in transfected mammalian cells. Here, we report that a covariance analysis of miniprecursor (p6*-PR) sequences derived from drug naïve patients identified a series of amino acid pairs that vary together across independent viral isolates. These covariance pairs were used to build the first topology map of the miniprecursor that suggests high levels of interaction between the p6* peptide and the mature PR. Additionally, several PR-PR covariance pairs are located far from each other (>12 Å Cα to Cα) relative to their positions in the mature PR structure. Biochemical characterization of one such covariance pair (77-93) revealed that each residue shows distinct preference for one of three alkyl amino acids (V, I, and L) and that a polar or charged amino acid at either of these two positions abolishes precursor autoprocessing. The most commonly observed 77V is preferred by the most commonly observed 93I, but the 77I variant is preferred by other 93 variances (L, V, or M) in supporting precursor autoprocessing. Furthermore, the 77I93V covariant enhanced precursor autoprocessing and Gag polyprotein processing but decreased the mature PR activity. Therefore, both covariance and biochemical analyses support a functional association between residues 77 and 93, which are spatially distant from each other in the mature PR structure. Our data also suggests that these covariance pairs differentially regulate precursor autoprocessing and the mature protease activity.

The triple threat of HIV-1 protease inhibitors

Newly released human immunodeficiency virus type 1 (HIV-1) particles obligatorily undergo a maturation process to become infectious. The HIV-1 protease (PR) initiates this step, catalyzing the cleavage of the Gag and Gag-Pro-Pol structural polyproteins. Proper organization of the mature virus core requires that cleavage of these polyprotein substrates proceeds in a highly regulated, specific series of events. The vital role the HIV-1 PR plays in the viral life cycle has made it an extremely attractive target for inhibition and has accordingly fostered the development of a number of highly potent substrate-analog inhibitors. Though the PR inhibitors (PIs) inhibit only the HIV-1 PR, their effects manifest at multiple different stages in the life cycle due to the critical importance of the PR in preparing the virus for these subsequent events. Effectively, PIs masquerade as entry inhibitors, reverse transcription inhibitors, and potentially even inhibitors of post-reverse transcription steps. In this chapter, we review the triple threat of PIs: the intermolecular cooperativity in the form of a cooperative dose-response for inhibition in which the apparent potency increases with increasing inhibition; the pleiotropic effects of HIV-1 PR inhibition on entry, reverse transcription, and post-reverse transcription steps; and their potency as transition state analogs that have the potential for further improvement that could lead to an inability of the virus to evolve resistance in the context of single drug therapy.

Conformation of inhibitor-free HIV-1 protease derived from NMR spectroscopy in a weakly oriented solution

Flexibility of the glycine-rich flaps is known to be essential for catalytic activity of the HIV-1 protease, but their exact conformations at the different stages of the enzymatic pathway remain subject to much debate. Although hundreds of crystal structures of protease-inhibitor complexes have been solved, only about a dozen inhibitor-free protease structures have been reported. These latter structures reveal a large diversity of flap conformations, ranging from closed to semi-open to wide open. To evaluate the average structure in solution, we measured residual dipolar couplings (RDCs) and compared these to values calculated for crystal structures representative of the closed, semi-open, and wide-open states. The RDC data clearly indicate that the inhibitor-free protease, on average, adopts a closed conformation in solution that is very similar to the inhibitor-bound state. By contrast, a highly drug-resistant protease mutant, PR20, adopts the wide-open flap conformation.

RC1339/APRc from Rickettsia conorii is a novel aspartic protease with properties of retropepsin-like enzymes

Members of the species Rickettsia are obligate intracellular, gram-negative, arthropod-borne pathogens of humans and other mammals. The life-threatening character of diseases caused by many Rickettsia species and the lack of reliable protective vaccine against rickettsioses strengthens the importance of identifying new protein factors for the potential development of innovative therapeutic tools. Herein, we report the identification and characterization of a novel membrane-embedded retropepsin-like homologue, highly conserved in 55 Rickettsia genomes. Using R. conorii gene homologue RC1339 as our working model, we demonstrate that, despite the low overall sequence similarity to retropepsins, the gene product of rc1339 APRc (for Aspartic Protease from Rickettsia conorii) is an active enzyme with features highly reminiscent of this family of aspartic proteases, such as autolytic activity impaired by mutation of the catalytic aspartate, accumulation in the dimeric form, optimal activity at pH 6, and inhibition by specific HIV-1 protease inhibitors. Moreover, specificity preferences determined by a high-throughput profiling approach confirmed common preferences between this novel rickettsial enzyme and other aspartic proteases, both retropepsins and pepsin-like. This is the first report on a retropepsin-like protease in gram-negative intracellular bacteria such as Rickettsia, contributing to the analysis of the evolutionary relationships between the two types of aspartic proteases. Additionally, we have also shown that APRc is transcribed and translated in R. conorii and R. rickettsii and is integrated into the outer membrane of both species. Finally, we demonstrated that APRc is sufficient to catalyze the in vitro processing of two conserved high molecular weight autotransporter adhesin/invasion proteins, Sca5/OmpB and Sca0/OmpA, thereby suggesting the participation of this enzyme in a relevant proteolytic pathway in rickettsial life-cycle. As a novel bona fide member of the retropepsin family of aspartic proteases, APRc emerges as an intriguing target for therapeutic intervention against fatal rickettsioses.

Several rickettsiae are pathogenic to humans by causing severe infections, including epidemic typhus (Rickettsia prowazekii), Rocky Mountain spotted fever (Rickettsia rickettsii), and Mediterranean spotted fever (Rickettsia conorii). Progress in correlating rickettsial genes and gene functions has been greatly hampered by the intrinsic difficulty in working with these obligate intracellular bacteria, despite the increasing insights into the mechanisms of pathogenesis of and the immune response to rickettsioses. Therefore, comparison of the multiple available genomes of Rickettsia is proving to be the most practical method to identify new factors that may play a role in pathogenicity. Here, we identified and characterized a novel retropepsin-like enzyme, APRc, that is expressed by at least two pathogenic rickettsial species, R. conorii and R. rickettsii. We have also established that APRc acts to process two major surface antigen/virulence determinants (OmpB/Sca5, OmpA/Sca0) in vitro and we suggest that this processing event is important for protein function. We demonstrate that APRc is specifically inhibited by drugs clinically used to treat HIV infections, providing the exciting possibility of targeting this enzyme for therapeutic intervention. With this work, we demonstrate that retropepsin-type aspartic proteases are indeed present in prokaryotes, suggesting that these enzymes may represent an ancestral form of these proteases.

Dimerization of HIV-1 protease occurs through two steps relating to the mechanism of protease dimerization inhibition by darunavir

Dimerization of HIV-1 protease (PR) subunits is an essential process for PR's acquisition of proteolytic activity, which plays a critical role in the maturation of HIV-1. Recombinant wild-type PR (PR(WT)) proved to dimerize, as examined with electrospray ionization mass spectrometry; however, two active site interface PR mutants (PR(T26A) and PR(R87K)) remained monomeric. On the other hand, two termini interface PR mutants (PR(1-C95A) and PR(97/99)) took both monomeric and dimeric forms. Differential scanning fluorimetry indicated that PR(1-C95A) and PR(97/99) dimers were substantially less stable than PR(WT) dimers. These data indicate that intermolecular interactions of two monomers occur first at the active site interface, generating unstable or transient dimers, and interactions at the termini interface subsequently occur, generating stable dimers. Darunavir (DRV), an HIV-1 protease inhibitor, inhibits not only proteolytic activity but also PR dimerization. DRV bound to protease monomers in a one-to-one molar ratio, inhibiting the first step of PR dimerization, whereas conventional protease inhibitors (such as saquinavir) that inhibit enzymatic activity but not dimerization failed to bind to monomers. DRV also bound to mutant PRs containing the transframe region-added PR (TFR-PR(D25N) and TFR-PR(D25N-7AA)), whereas saquinavir did not bind to TFR-PR(D25N) or TFR-PR(D25N-7AA). Notably, DRV failed to bind to mutant PR containing four amino acid substitutions (V32I, L33F, I54M, and I84V) that confer resistance to DRV on HIV-1. To our knowledge, the present report represents the first demonstration of the two-step PR dimerization dynamics and the mechanism of dimerization inhibition by DRV, which should help design further, more potent novel PIs.

Structures of darunavir-resistant HIV-1 protease mutant reveal atypical binding of darunavir to wide open flaps

The molecular basis for high resistance to clinical inhibitors of HIV-1 protease (PR) was examined for the variant designated PR_P51 that was selected for resistance to darunavir (DRV). High resolution crystal structures of PR_P51 with the active site D25N mutation revealed a ligand-free form and an inhibitor-bound form showing a unique binding site and orientation for DRV. This inactivating mutation is known to increase the dimer dissociation constant and decrease DRV affinity of PR. The PR_P51-D25N dimers were in the open conformation with widely separated flaps, as reported for other highly resistant variants. PR_P51-D25N dimer bound two DRV molecules and showed larger separation of 8.7 Å between the closest atoms of the two flaps compared with 4.4 Å for the ligand-free structure of this mutant. The ligand-free structure, however, lacked van der Waals contacts between Ile50 and Pro81′ from the other subunit in the dimer, unlike the majority of PR structures. DRV is bound inside the active site cavity; however, the inhibitor is oriented almost perpendicular to its typical position and exhibits only 2 direct hydrogen bond and two water-mediated interactions with atoms of PR_P51-D25N compared with 11 hydrogen bond interactions seen for DRV bound in the typical position in wild-type enzyme. The atypical location of DRV may provide opportunities for design of novel inhibitors targeting the open conformation of PR drug-resistant mutants.

Understanding HIV-1 protease autoprocessing for novel therapeutic development

In the infected cell, HIV-1 protease (PR) is initially synthesized as part of the GagPol polyprotein. PR autoprocessing is a virus-specific process by which the PR domain embedded in the precursor catalyzes proteolytic reactions responsible for liberation of free mature PRs, which then recognize and cleave at least ten different peptide sequences in the Gag and GagPol polyproteins. Despite extensive structure and function studies of the mature PRs as well as the successful development of ten US FDA-approved catalytic-site inhibitors, the precursor autoprocessing mechanism remains an intriguing yet-to-be-solved puzzle. This article discusses current understanding of the autoprocessing mechanism, in an effort to prompt the development of novel anti-HIV drugs that selectively target precursor autoprocessing.

The maturation of HIV-1 protease precursor studied by discrete molecular dynamics

The equilibrium properties of a HIV-1-protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native-like, but their binding can be non-native.

Enhanced stability of monomer fold correlates with extreme drug resistance of HIV-1 protease

During treatment, mutations in HIV-1 protease (PR) are selected rapidly that confer resistance by decreasing affinity to clinical protease inhibitors (PIs). As these unique drug resistance mutations can compromise the fitness of the virus to replicate, mutations that restore conformational stability and activity while retaining drug resistance are selected on further evolution. Here we identify several compensating mechanisms by which an extreme drug-resistant mutant bearing 20 mutations (PR20) with >5-fold increased Kd and >4000-fold decreased affinity to the PI darunavir functions. (1) PR20 cleaves, albeit poorly, Gag polyprotein substrates essential for viral maturation. (2) PR20 dimer, which exhibits distinctly enhanced thermal stability, has highly attenuated autoproteolysis, thus likely prolonging its lifetime in vivo. (3) The enhanced stability of PR20 results from stabilization of the monomer fold. Both monomeric PR20(T26A) and dimeric PR20 exhibit Tm values 6-7.5 °C higher than those for their PR counterparts. Two specific mutations in PR20, L33F and L63P at sites of autoproteolysis, increase the Tm of monomeric PR(T26A) by ~8 °C, similar to PR20(T26A). However, without other compensatory mutations as seen in PR20, L33F and L63P substitutions, together, neither restrict autoproteolysis nor significantly reduce binding affinity to darunavir. To determine whether dimer stability contributes to binding affinity for inhibitors, we examined single-chain dimers of PR and PR(D25N) in which the corresponding identical monomer units were covalently linked by GGSSG sequence. Linking of the subunits did not appreciably change the ΔTm on inhibitor binding; thus stabilization by tethering appears to have little direct effect on enhancing inhibitor affinity.

Extreme multidrug resistant HIV-1 protease with 20 mutations is resistant to novel protease inhibitors with P1'-pyrrolidinone or P2-tris-tetrahydrofuran

Extreme drug resistant mutant of HIV-1 protease (PR) bearing 20 mutations (PR20) has been studied with the clinical inhibitor amprenavir (1) and two potent antiviral investigational inhibitors GRL-02031 (2) and GRL-0519 (3). Clinical inhibitors are >1000-fold less active on PR20 than on wild-type enzyme, which is consistent with dissociation constants (KL) from isothermal titration calorimetry of 40 nM for 3, 178 nM for amprenavir, and 960 nM for 2. High resolution crystal structures of PR20-inhibitor complexes revealed altered interactions compared with the corresponding wild-type PR complexes in agreement with relative inhibition. Amprenavir lacks interactions due to PR20 mutations in the S2/S2' subsites relative to PR. Inhibitors 2 and 3 lose interactions with Arg8' in PR20 relative to the wild-type enzyme because Arg8' shifts to interact with mutated L10F side chain. Overall, inhibitor 3 compares favorably with darunavir in affinity for PR20 and shows promise for further development.