Comparison of the HIV-1 and HIV-2 proteinases using oligopeptide substrates representing cleavage sites in Gag and Gag-Pol polyproteins

Sequence dependencies and biophysical features both govern cleavage of diverse cut-sites by HIV protease

The infectivity of HIV-1 requires its protease (PR) cleave multiple cut-sites with low sequence similarity. The diversity of cleavage sites has made it challenging to investigate the underlying sequence properties that determine binding and turnover of substrates by PR. We engineered a mutational scanning approach utilizing yeast display, flow cytometry, and deep sequencing to systematically measure the impacts of all individual amino acid changes at 12 positions in three different cut-sites (MA/CA, NC/p1, and p1/p6). The resulting fitness landscapes revealed common physical features that underlie cutting of all three cut-sites at the amino acid positions closest to the scissile bond. In contrast, positions more than two amino acids away from the scissile bond exhibited a strong dependence on the sequence background of the rest of the cut-site. We observed multiple amino acid changes in cut-sites that led to faster cleavage rates, including a preference for negative charge five and six amino acids away from the scissile bond at locations where the surface of protease is positively charged. Analysis of individual cut sites using full-length matrix-capsid proteins indicate that long-distance sequence context can contribute to cutting efficiency such that analyses of peptides or shorter engineered constructs including those in this work should be considered carefully. This work provides a framework for understanding how diverse substrates interact with HIV-1 PR and can be extended to investigate other viral PRs with similar properties.

Dynamics of HIV-1 Gag Processing as Revealed by Fluorescence Lifetime Imaging Microscopy and Single Virus Tracking

The viral polyprotein Gag plays a central role for HIV-1 assembly, release and maturation. Proteolytic processing of Gag by the viral protease is essential for the structural rearrangements that mark the transition from immature to mature, infectious viruses. The timing and kinetics of Gag processing are not fully understood. Here, fluorescence lifetime imaging microscopy and single virus tracking are used to follow Gag processing in nascent HIV-1 particles in situ. Using a Gag polyprotein labelled internally with eCFP, we show that proteolytic release of the fluorophore from Gag is accompanied by an increase in its fluorescence lifetime. By tracking nascent virus particles in situ and analyzing the intensity and fluorescence lifetime of individual traces, we detect proteolytic cleavage of eCFP from Gag in a subset (6.5%) of viral particles. This suggests that for the majority of VLPs, Gag processing occurs with a delay after particle assembly.

Stephan Oroszlan and the Proteolytic Processing of Retroviral Proteins: Following A Pro

Steve Oroszlan determined the sequences at the ends of virion proteins for a number of different retroviruses. This work led to the insight that the amino-terminal amino acid of the mature viral CA protein is always proline. In this remembrance, we review Steve’s work that led to this insight and show how that insight was a necessary precursor to the work we have done in the subsequent years exploring the cleavage rate determinants of viral protease processing sites and the multiple roles the amino-terminal proline of CA plays after protease cleavage liberates it from its position in a protease processing site.

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.

Subtype-specific differences in Gag-protease replication capacity of HIV-1 isolates from East and West Africa

BACKGROUND

The HIV-1 epidemic in sub-Saharan Africa is heterogeneous with diverse unevenly distributed subtypes and regional differences in prevalence. Subtype-specific differences in disease progression rate and transmission efficiency have been reported, but the underlying biological mechanisms have not been fully characterized. Here, we tested the hypothesis that the subtypes prevalent in the East Africa, where adult prevalence rate is higher, have lower viral replication capacity (VRC) than their West African counterparts where adult prevalence rates are lower.

RESULTS

Gag-protease sequencing was performed on 213 and 160 antiretroviral-naïve chronically infected participants from West and East Africa respectively and bioinformatic tools were used to infer subtypes and recombination patterns. VRC of patient-derived gag-protease chimeric viruses from West (n = 178) and East (n = 114) Africa were determined using a green fluorescent protein reporter-based cell assay. Subtype and regional differences in VRC and amino acid variants impacting VRC were identified by statistical methods. CRF02_AG (65%, n = 139), other recombinants (14%, n = 30) and pure subtypes (21%, n = 44) were identified in West Africa. Subtypes A1 (64%, n = 103), D (22%, n = 35), or recombinants (14%, n = 22) were identified in East Africa. Viruses from West Africa had significantly higher VRC compared to those from East Africa (p < 0.0001), with subtype-specific differences found among strains within West and East Africa (p < 0.0001). Recombination patterns showed a preference for subtypes D, G or J rather than subtype A in the p6 region of gag, with evidence that subtype-specific differences in this region impact VRC. Furthermore, the Gag A83V polymorphism was associated with reduced VRC in CRF02_AG. HLA-A*23:01 (p = 0.0014) and HLA-C*07:01 (p = 0.002) were associated with lower VRC in subtype A infected individuals from East Africa.

CONCLUSIONS

Although prevalent viruses from West Africa displayed higher VRC than those from East Africa consistent with the hypothesis that lower VRC is associated with higher population prevalence, the predominant CRF02_AG strain in West Africa displayed higher VRC than other prevalent strains suggesting that VRC alone does not explain population prevalence. The study identified viral and host genetic determinants of virus replication capacity for HIV-1 CRF02_AG and subtype A respectively, which may have relevance for vaccine strategies.

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.

Evaluation of Structurally Distorted Split GFP Fluorescent Sensors for Cell-Based Detection of Viral Proteolytic Activity

Cell-based assays are essential for virus functional characterization in fundamental and applied research. Overcoming the limitations of virus-labelling strategies while allowing functional assessment of critical viral enzymes, virus-induced cell-based biosensors constitute a powerful approach. Herein, we designed and characterized different cell-based switch-on split GFP sensors reporting viral proteolytic activity and virus infection. Crucial to these sensors is the effective-yet reversible-fluorescence off-state, through protein distortion. For that, single (protein embedment or intein-mediated cyclization) or dual (coiled-coils) distortion schemes prevent split GFP self-assembly, until virus-promoted proteolysis of a cleavable sequence. All strategies showed their applicability in detecting viral proteolysis, although with different efficiencies depending on the protease. While for tobacco etch virus protease the best performing sensor was based on coiled-coils (signal-to-noise ratio, SNR, 97), for adenovirus and lentivirus proteases it was based on GFP11 cyclization (SNR 3.5) or GFP11 embedment distortion (SNR 6.0), respectively. When stably expressed, the sensors allowed live cell biosensing of adenovirus infection, with sensor fluorescence activation 24 h post-infection. The structural distortions herein studied are highly valuable in the development of cellular biosensing platforms. Additionally highlighted, selection of the best performing strategy is highly dependent on the unique properties of each viral protease.

Fine-Tuning of Sequence Specificity by Near Attack Conformations in Enzyme-Catalyzed Peptide Hydrolysis

The catalytic role of near attack conformations (NACs), molecular states that lie on the pathway between the ground state (GS) and transition state (TS) of a chemical reaction, is not understood completely. Using a computational approach that combines Bürgi–Dunitz theory with all-atom molecular dynamics simulations, the role of NACs in catalyzing the first stages of HIV-1 protease peptide hydrolysis was previously investigated using a substrate that represents the recognized SP1-NC cleavage site of the HIV-1 Gag polyprotein. NACs were found to confer no catalytic effect over the uncatalyzed reaction there ( Δ Δ G N ‡ ∼ 0 kcal/mol). Here, using the same approach, the role of NACs across multiple substrates that each represent a further recognized cleavage site is investigated. Overall rate enhancement varies by | Δ Δ G ‡ | ∼ 12–15 kcal/mol across this set, and although NACs contribute a small and approximately constant barrier to the uncatalyzed reaction (< Δ G N ‡ u > = 4.3 ± 0.3 kcal/mol), they are found to contribute little significant catalytic effect ( | Δ Δ G N ‡ | ∼ 0–2 kcal/mol). Furthermore, no correlation is exhibited between NAC contributions and the overall energy barrier ( R 2 = 0.01). However, these small differences in catalyzed NAC contributions enable rates to match those required for the kinetic order of processing. Therefore, NACs may offer an alternative and subtle mode compared to non-NAC contributions for fine-tuning reaction rates during complex evolutionary sequence selection processes—in this case across cleavable polyproteins whose constituents exhibit multiple functions during the virus life-cycle.

Circular zymogens of human ribonuclease 1

The endogenous production of enzymes as zymogens provides a means to control catalytic activities. Here, we describe the heterologous production of ribonuclease 1 (RNase 1), which is the most prevalent secretory ribonuclease in humans, as a zymogen. In folded RNase 1, the N and C termini flank the enzymic active site. By using intein-mediated cis-splicing, we created circular proteins in which access to the active site of RNase 1 is obstructed by an amino-acid sequence that is recognized by the HIV-1 protease. Installing a sequence that does not perturb the RNase 1 fold led to only modest inactivation. In contrast, the ancillary truncation of residues from each terminus led to a substantial decrease in the catalytic activity of the zymogen with the maintenance of thermostability. For optimized zymogens, activation by HIV-1 protease led to a > 104 -fold increase in ribonucleolytic activity at a rate comparable to that for the cleavage of endogenous viral substrates. Molecular modeling indicated that these zymogens are inactivated by conformational distortion in addition to substrate occlusion. Because protease levels are elevated in many disease states and ribonucleolytic activity can be cytotoxic, RNase 1 zymogens have potential as generalizable prodrugs.

Highly Drug-Resistant HIV-1 Protease Mutant PRS17 Shows Enhanced Binding to Substrate Analogues

We report the structural analysis of highly drug-resistant human immunodeficiency virus protease (PR) variant PR^S17, rationally selected by machine learning, in complex with substrate analogues. Crystal structures were solved of inhibitor-free inactive PR^S17-D25N, wild-type PR/CA-p2 complex, and PR^S17 in complex with substrate analogues, CA-p2 and p2-NC. Peptide analogues p2-NC and CA-p2 exhibit inhibition constants of 514 and 22 nM, respectively, for PR^S17 or approximately 3-fold better than for PR. CA-p2 is a better inhibitor of PR^S17 than are clinical inhibitors (K _i = 50-8390 nM) except for amprenavir (K _i = 11 nM). G48V resistance mutation induces curled flap tips in PR^S17-D25N structure. The inner P2-P2' residues of substrate analogues in PR^S17 complexes maintain similar conformations to those of wild-type complex, while significant conformational changes are observed in the peripheral residues P3, P4' of CA-p2 and P3, P4, and P3' of p2-NC. The loss of β-branched side chain by V82S mutation initiates a shift in 80's loop and reshapes the S3/S3' subsite, which enhances substrate binding with new hydrogen bonds and van der Waals interactions that are absent in the wild-type structures. The steric hindrance caused by G48V mutation in the flap of PR^S17 contributes to altered binding interactions of P3 Arg, P4' norleucine of CA-p2, and P4 and P3' of p2-NC with the addition of new hydrogen bonds and van der Waals contacts. The enhanced interaction of PR^S17 with substrate analogues agrees with their relative inhibition, suggesting that this mutant improves substrate binding while decreasing affinity for clinical inhibitors.

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

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

CCR5/CXCR4 Dual Antagonism for the Improvement of HIV Infection Therapy

HIV entry in the host cell requires the interaction with the CD4 membrane receptor, and depends on the activation of one or both co-receptors CCR5 and CXCR4. Former selective co-receptor antagonists, acting at early stages of infection, are able to impair the receptor functions, preventing the viral spread toward AIDS. Due to the capability of HIV to develop resistance by switching from CCR5 to CXCR4, dual co-receptor antagonists could represent the next generation of AIDS prophylaxis drugs. We herein present a survey on relevant results published in the last few years on compounds acting simultaneously on both co-receptors, potentially useful as preventing agents or in combination with classical anti-retroviral drugs based therapy.

An n→π* Interaction in the Bound Substrate of Aspartic Proteases Replicates the Oxyanion Hole

Aspartic proteases regulate many biological processes and are prominent targets for therapeutic intervention. Structural studies have captured intermediates along the reaction pathway, including the Michaelis complex and tetrahedral intermediate. Using a Ramachandran analysis of these structures, we discovered that residues occupying the P1 and P1' positions (which flank the scissile peptide bond) adopt the dihedral angle of an inverse γ-turn and polyproline type-II helix, respectively. Computational analyses reveal that the polyproline type-II helix engenders an n→π* interaction in which the oxygen of the scissile peptide bond is the donor. This interaction stabilizes the negative charge that develops in the tetrahedral intermediate, much like the oxyanion hole of serine proteases. The inverse γ-turn serves to twist the scissile peptide bond, vacating the carbonyl π* orbital and facilitating its hydration. These previously unappreciated interactions entail a form of substrate-assisted catalysis and offer opportunities for drug design.

HIV-1 Protease Uses Bi-Specific S2/S2' Subsites to Optimize Cleavage of Two Classes of Target Sites

Retroviral proteases (PRs) have a unique specificity that allows cleavage of sites with or without a P1' proline. A P1' proline is required at the MA/CA cleavage site due to its role in a post-cleavage conformational change in the capsid protein. However, the HIV-1 PR prefers to have large hydrophobic amino acids flanking the scissile bond, suggesting that PR recognizes two different classes of substrate sequences. We analyzed the cleavage rate of over 150 combinations of six different HIV-1 cleavage sites to explore rate determinants of cleavage. We found that cleavage rates are strongly influenced by the two amino acids flanking the amino acids at the scissile bond (P2-P1/P1'-P2'), with two complementary sets of rules. When P1' is proline, the P2 side chain interacts with a polar region in the S2 subsite of the PR, while the P2' amino acid interacts with a hydrophobic region of the S2' subsite. When P1' is not proline, the orientations of the P2 and P2' side chains with respect to the scissile bond are reversed; P2 residues interact with a hydrophobic face of the S2 subsite, while the P2' amino acid usually engages hydrophilic amino acids in the S2' subsite. These results reveal that the HIV-1 PR has evolved bi-functional S2 and S2' subsites to accommodate the steric effects imposed by a P1' proline on the orientation of P2 and P2' substrate side chains. These results also suggest a new strategy for inhibitor design to engage the multiple specificities in these subsites.

Exploring the Reasons for Decrease in Binding Affinity of HIV-2 Against HIV-1 Protease Complex Using Interaction Entropy Under Polarized Force Field

In this study, the differences of binding patterns between two type HIV (HIV-1 and HIV-2) protease and two inhibitors (darunavir and amprenavir) are analyzed and compared using the newly developed interaction entropy (IE) method for the entropy change calculation combined with the polarized force field. The functional role of protonation states in the two HIV-2 complexes is investigated and our study finds that the protonated OD1 atom of Asp25' in B chain is the optimal choice. Those calculated binding free energies obtained from the polarized force field combined with IE method are significantly consistent with the experimental observed. The bridging water W301 is favorable to the binding of HIV-1 complexes; however, it is unfavorable to the HIV-2 complexes in current study. The volume of pocket, B-factor of Cα atoms and the distance of flap tip in HIV-2 complexes are smaller than that of HIV-1 consistently. These changes may cause localized rearrangement of residues lining their surface and finally result in the different binding mode for the two types HIV. Predicated hot-spot residues (Ala28/Ala28', Ile50/Ile50', and Ile84/Ile84') are nearly same in the four systems. However, the contribution to the free energy of Asp30 residue is more favorable in HIV-1 system than in HIV-2 system. Current study, to some extent, reveals the origin for the decrease in binding affinity of inhibitors against HIV-2 compared with HIV-1 and will provides theoretical guidance for future design of potent dual inhibitors targeting two type HIV protease.

A recombinant fusion protein-based, fluorescent protease assay for high throughput-compatible substrate screening

In connection with the intensive investigation of proteases, several methods have been developed for analysis of the substrate specificity. Due to the great number of proteases and the expected target molecules to be analyzed, time- and cost-efficient high-throughput screening (HTS) methods are preferred. Here we describe the development and application of a separation-based HTS-compatible fluorescent protease assay, which is based on the use of recombinant fusion proteins as substrates of proteases. The protein substrates used in this assay consists of N-terminal (hexahistidine and maltose binding protein) fusion tags, cleavage sequences of the tobacco etch virus (TEV) and HIV-1 proteases, and a C-terminal fluorescent protein (mApple or mTurquoise2). The assay is based on the fluorimetric detection of the fluorescent proteins, which are released from the magnetic bead-attached substrates by the proteolytic cleavage. The protease assay has been applied for activity measurements of TEV and HIV-1 proteases to test the suitability of the system for enzyme kinetic measurements, inhibition studies, and determination of pH optimum. We also found that denatured fluorescent proteins can be renatured after SDS-PAGE of denaturing conditions, but showed differences in their renaturation abilities. After in-gel renaturation both substrates and cleavage products can be identified by in-gel UV detection.

Stimulated Emission Depletion Nanoscopy Reveals Time-Course of Human Immunodeficiency Virus Proteolytic Maturation

Concomitant with human immunodeficiency virus type 1 (HIV-1) budding from a host cell, cleavage of the structural Gag polyproteins by the viral protease (PR) triggers complete remodeling of virion architecture. This maturation process is essential for virus infectivity. Electron tomography provided structures of immature and mature HIV-1 with a diameter of 120-140 nm, but information about the sequence and dynamics of structural rearrangements is lacking. Here, we employed super-resolution STED (stimulated emission depletion) fluorescence nanoscopy of HIV-1 carrying labeled Gag to visualize the virion architecture. The incomplete Gag lattice of immature virions was clearly distinguishable from the condensed distribution of mature protein subunits. Synchronized activation of PR within purified particles by photocleavage of a caged PR inhibitor enabled time-resolved in situ observation of the induction of proteolysis and maturation by super-resolution microscopy. This study shows the rearrangement of subviral structures in a super-resolution light microscope over time, outwitting phototoxicity and fluorophore bleaching through synchronization of a biological process by an optical switch.

Four Amino Acid Changes in HIV-2 Protease Confer Class-Wide Sensitivity to Protease Inhibitors

Protease is essential for retroviral replication, and protease inhibitors (PI) are important for treating HIV infection. HIV-2 exhibits intrinsic resistance to most FDA-approved HIV-1 PI, retaining clinically useful susceptibility only to lopinavir, darunavir, and saquinavir. The mechanisms for this resistance are unclear; although HIV-1 and HIV-2 proteases share just 38 to 49% sequence identity, all critical structural features of proteases are conserved. Structural studies have implicated four amino acids in the ligand-binding pocket (positions 32, 47, 76, and 82). We constructed HIV-2ROD9 molecular clones encoding the corresponding wild-type HIV-1 amino acids (I32V, V47I, M76L, and I82V) either individually or together (clone PRΔ4) and compared the phenotypic sensitivities (50% effective concentration [EC50]) of mutant and wild-type viruses to nine FDA-approved PI. Single amino acid replacements I32V, V47I, and M76L increased the susceptibility of HIV-2 to multiple PI, but no single change conferred class-wide sensitivity. In contrast, clone PRΔ4 showed PI susceptibility equivalent to or greater than that of HIV-1 for all PI. We also compared crystallographic structures of wild-type HIV-1 and HIV-2 proteases complexed with amprenavir and darunavir to models of the PRΔ4 enzyme. These models suggest that the amprenavir sensitivity of PRΔ4 is attributable to stabilizing enzyme-inhibitor interactions in the P2 and P2' pockets of the protease dimer. Together, our results show that the combination of four amino acid changes in HIV-2 protease confer a pattern of PI susceptibility comparable to that of HIV-1, providing a structural rationale for intrinsic HIV-2 PI resistance and resolving long-standing questions regarding the determinants of differential PI susceptibility in HIV-1 and HIV-2.

IMPORTANCE

Proteases are essential for retroviral replication, and HIV-1 and HIV-2 proteases share a great deal of structural similarity. However, only three of nine FDA-approved HIV-1 protease inhibitors (PI) are active against HIV-2. The underlying reasons for intrinsic PI resistance in HIV-2 are not known. We examined the contributions of four amino acids in the ligand-binding pocket of the enzyme that differ between HIV-1 and HIV-2 by constructing HIV-2 clones encoding the corresponding HIV-1 amino acids and testing the PI susceptibilities of the resulting viruses. We found that the HIV-2 clone containing all four changes (PRΔ4) was as susceptible as HIV-1 to all nine PI. We also modeled the PRΔ4 enzyme structure and compared it to existing crystallographic structures of HIV-1 and HIV-2 proteases complexed with amprenavir and darunavir. Our findings demonstrate that four positions in the ligand-binding cleft of protease are the primary cause of HIV-2 PI resistance.

Elucidation of the Molecular Mechanism Driving Duplication of the HIV-1 PTAP Late Domain

HIV-1 uses cellular machinery to bud from infected cells. This cellular machinery is comprised of several multiprotein complexes known as endosomal sorting complexes required for transport (ESCRTs). A conserved late domain motif, Pro-Thr-Ala-Pro (PTAP), located in the p6 region of Gag (p6(Gag)), plays a central role in ESCRT recruitment to the site of virus budding. Previous studies have demonstrated that PTAP duplications are selected in HIV-1-infected patients during antiretroviral therapy; however, the consequences of these duplications for HIV-1 biology and drug resistance are unclear. To address these questions, we constructed viruses carrying a patient-derived PTAP duplication with and without drug resistance mutations in the viral protease. We evaluated the effect of the PTAP duplication on viral release efficiency, viral infectivity, replication capacity, drug susceptibility, and Gag processing. In the presence of protease inhibitors, we observed that the PTAP duplication in p6(Gag) significantly increased the infectivity and replication capacity of the virus compared to those of viruses bearing only resistance mutations in protease. Our biochemical analysis showed that the PTAP duplication, in combination with mutations in protease, enhances processing between the nucleocapsid and p6 domains of Gag, resulting in more complete Gag cleavage in the presence of protease inhibitors. These results demonstrate that duplication of the PTAP motif in p6(Gag) confers a selective advantage in viral replication by increasing Gag processing efficiency in the context of protease inhibitor treatment, thereby enhancing the drug resistance of the virus. These findings highlight the interconnected role of PTAP duplications and protease mutations in the development of resistance to antiretroviral therapy.

IMPORTANCE

Resistance to current drug therapy limits treatment options in many HIV-1-infected patients. Duplications in a Pro-Thr-Ala-Pro (PTAP) motif in the p6 domain of Gag are frequently observed in viruses derived from patients on protease inhibitor (PI) therapy. However, the reason that these duplications arise and their consequences for virus replication remain to be established. In this study, we examined the effect of PTAP duplication on PI resistance in the context of wild-type protease or protease bearing PI resistance mutations. We observe that PTAP duplication markedly enhances resistance to a panel of PIs. Biochemical analysis reveals that the PTAP duplication reverses a Gag processing defect imposed by the PI resistance mutations in the context of PI treatment. The results provide a long-sought explanation for why PTAP duplications arise in PI-treated patients.

HIV-1 protease substrate-groove: Role in substrate recognition and inhibitor resistance

A key target in the treatment of HIV-1/AIDS has been the viral protease. Here we first studied in silico the evolution of protease resistance. Primary active site resistance mutations were found to weaken interactions between protease and both inhibitor and substrate P4-P4' residues. We next studied the effects of secondary resistance mutations, often distant from the active site, on protease binding to inhibitors and substrates. Those secondary mutations contributed to the rise of multi-drug resistance while also enhancing viral replicative capacity. Here many secondary resistance mutations were found in the HIV-1 protease substrate-grooves, one on each face of the symmetrical protease dimer. The protease active site binds substrate P4-P4' residues, while the substrate-groove allows the protease to bind residues P12-P5/P5'-P12', for a total of twenty-four residues. The substrate-groove secondary resistance mutations were found to compensate for the loss of interactions between the inhibitor resistant protease active site and substrate P4-P4' residues, due to primary resistance mutations, by increasing interactions with substrate P12-P5/P5'-P12' residues. In vitro experiments demonstrated that a multi-drug resistant protease with substrate-groove resistance mutations was slower than wild-type protease in cleaving a peptide substrate, which did not allow for substrate-groove interactions, while it had similar activity as wild-type protease when using a Gag polyprotein in which cleavage-site P12-P5/P5'-P12' residues could be bound by the protease substrate-grooves. When the Gag MA/CA cleavage site P12-P5/P5'-P12' residues were mutated the multi-drug resistant protease cleaved the mutant Gag significantly slower, indicating the importance of the protease S-grooves in binding to substrate.

A Direct Interaction with RNA Dramatically Enhances the Catalytic Activity of the HIV-1 Protease In Vitro

Though the steps of human immunodeficiency virus type 1 (HIV-1) virion maturation are well documented, the mechanisms regulating the proteolysis of the Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) remain obscure. One proposed mechanism argues that the maturation intermediate p15NC must interact with RNA for efficient cleavage by the PR. We investigated this phenomenon and found that processing of multiple substrates by the HIV-1 PR was enhanced in the presence of RNA. The acceleration of proteolysis occurred independently from the substrate's ability to interact with nucleic acid, indicating that a direct interaction between substrate and RNA is not necessary for enhancement. Gel-shift assays demonstrated the HIV-1 PR is capable of interacting with nucleic acids, suggesting that RNA accelerates processing reactions by interacting with the PR rather than the substrate. All HIV-1 PRs examined have this ability; however, the HIV-2 PR does not interact with RNA and does not exhibit enhanced catalytic activity in the presence of RNA. No specific sequence or structure was required in the RNA for a productive interaction with the HIV-1 PR, which appears to be principally, though not exclusively, driven by electrostatic forces. For a peptide substrate, RNA increased the kinetic efficiency of the HIV-1 PR by an order of magnitude, affecting both turnover rate (k(cat)) and substrate affinity (K(m)). These results suggest that an allosteric binding site exists on the HIV-1 PR and that HIV-1 PR activity during maturation could be regulated in part by the juxtaposition of the enzyme with virion-packaged RNA.

Drug-resistant HIV-1 protease regains functional dynamics through cleavage site coevolution

Drug resistance is caused by mutations that change the balance of recognition favoring substrate cleavage over inhibitor binding. Here, a structural dynamics perspective of the regained wild-type functioning in mutant HIV-1 proteases with coevolution of the natural substrates is provided. The collective dynamics of mutant structures of the protease bound to p1-p6 and NC-p1 substrates are assessed using the Anisotropic Network Model (ANM). The drug-induced protease mutations perturb the mechanistically crucial hinge axes that involve key sites for substrate binding and dimerization and mainly coordinate the intrinsic dynamics. Yet with substrate coevolution, while the wild-type dynamic behavior is restored in both p1-p6 ((LP) (1'F)p1-p6D30N/N88D) and NC-p1 ((AP) (2) (V)NC-p1V82A) bound proteases, the dynamic behavior of the NC-p1 bound protease variants (NC-p1V82A and (AP) (2) (V)NC-p1V82A) rather resemble those of the proteases bound to the other substrates, which is consistent with experimental studies. The orientational variations of residue fluctuations along the hinge axes in mutant structures justify the existence of coevolution in p1-p6 and NC-p1 substrates, that is, the dynamic behavior of hinge residues should contribute to the interdependent nature of substrate recognition. Overall, this study aids in the understanding of the structural dynamics basis of drug resistance and evolutionary optimization in the HIV-1 protease system.

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.

State of the art prediction of HIV-1 protease cleavage sites

MOTIVATION

Understanding the substrate specificity of human immunodeficiency virus (HIV)-1 protease is important when designing effective HIV-1 protease inhibitors. Furthermore, characterizing and predicting the cleavage profile of HIV-1 protease is essential to generate and test hypotheses of how HIV-1 affects proteins of the human host. Currently available tools for predicting cleavage by HIV-1 protease can be improved.

RESULTS

The linear support vector machine with orthogonal encoding is shown to be the best predictor for HIV-1 protease cleavage. It is considerably better than current publicly available predictor services. It is also found that schemes using physicochemical properties do not improve over the standard orthogonal encoding scheme. Some issues with the currently available data are discussed.

AVAILABILITY AND IMPLEMENTATION

The datasets used, which are the most important part, are available at the UCI Machine Learning Repository. The tools used are all standard and easily available.

CONTACT

thorsteinn.rognvaldsson@hh.se.

A modular system to evaluate the efficacy of protease inhibitors against HIV-2

The human immunodeficiency virus (HIV) protease is a homodimeric aspartyl protease that is crucial for the viral life-cycle, cleaving proviral polyproteins, hence creating mature protein components that are required for the formation of an infectious virus. With diagnostic measures and clinically used protease inhibitors focusing on HIV-1, due to its higher virulence and prevalence, studies of the efficacy of those inhibitors on HIV-2 protease remain widely lacking. Utilizing a wild-type HIV-2 vector backbone and cloning techniques we have developed a cassette system where the efficacy of clinically used protease inhibitors can be studied for various serotypes of HIV-2 protease both in enzymatic and cell culture assays. In our experiments, optimization of the expression protocol led to a relatively stable enzyme, for cell culture assays, the efficiency of transfection and transduction capability of the modified vector was tested and was not found to differ from that of the wild-type, moreover, a 2^nd generation protease inhibitor was used to demonstrate the usefulness of the system. The combination of assays performed with our cassette system is expected to provide an accurate measure of the efficacy of currently used; as well as experimental protease inhibitors on HIV-2.

In silico prediction of mutant HIV-1 proteases cleaving a target sequence

HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636–1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.

Quantitative assessment of in vivo HIV protease activity using genetically engineered QD-based FRET probes

HIV protease plays a central role in its life cycle leading to release of functional viral particles. It has been successfully used as a therapeutic target to block HIV infection. Several protease inhibitors (PIs) are currently being employed as a part of anti-HIV therapy. However, the constant genetic drift in the virus leads to accumulation of mutations in both cleavage site and the protease, resulting in resistance and failure of therapy. We reported the use of a quantum dot (QD)-based protein probe for the in vivo monitoring of HIV-1 protease activity based on fluorescence resonance energy transfer. In the current study, we demonstrate the utility of this approach by quantifying the in vivo cleavage rates of three known protease and cleavage site mutations in the presence or absence of different PIs. The changes in IC50 values for the different PIs were similar to that observed in patients, validating our assay as a rapid platform for PI screening.

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.

Bioinformatic approaches for modeling the substrate specificity of HIV-1 protease: an overview

HIV-1 protease has a broad and complex substrate specificity, which hitherto has escaped a simple comprehensive definition. This, and the relatively high mutation rate of the retroviral protease, makes it challenging to design effective protease inhibitors. Several attempts have been made during the last two decades to elucidate the enigmatic cleavage specificity of HIV-1 protease and to predict cleavage of novel substrates using bioinformatic analysis methods. This review describes the methods that have been utilized to date to address this important problem and the results achieved. The data sets used are also reviewed and important aspects of these are highlighted.

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.

Gag-Pol processing during HIV-1 virion maturation: a systems biology approach

Proteolytic processing of Gag and Gag-Pol polyproteins by the viral protease (PR) is crucial for the production of infectious HIV-1, and inhibitors of the viral PR are an integral part of current antiretroviral therapy. The process has several layers of complexity (multiple cleavage sites and substrates; multiple enzyme forms; PR auto-processing), which calls for a systems level approach to identify key vulnerabilities and optimal treatment strategies. Here we present the first full reaction kinetics model of proteolytic processing by HIV-1 PR, taking into account all canonical cleavage sites within Gag and Gag-Pol, intermediate products and enzyme forms, enzyme dimerization, the initial auto-cleavage of full-length Gag-Pol as well as self-cleavage of PR. The model allows us to identify the rate limiting step of virion maturation and the parameters with the strongest effect on maturation kinetics. Using the modelling framework, we predict interactions and compensatory potential between individual cleavage rates and drugs, characterize the time course of the process, explain the steep dose response curves associated with PR inhibitors and gain new insights into drug action. While the results of the model are subject to limitations arising from the simplifying assumptions used and from the uncertainties in the parameter estimates, the developed framework provides an extendable open-access platform to incorporate new data and hypotheses in the future.

Impact of gag genetic determinants on virological outcome to boosted lopinavir-containing regimen in HIV-2-infected patients

OBJECTIVE

This study investigated the impact on virological outcome of the gag cleavage sites and the protease-coding region mutations in protease inhibitor-naive and protease inhibitor-experienced patients infected with HIV-2 receiving lopinavir (LPV) containing regimen.

METHODS

Baseline gag and protease-coding region were sequenced in 46 HIV-2 group A-infected patients receiving lopinavir. Virological response was defined as plasma viral load less than 100 copies/ml at month 3. Associations between virological response and frequencies of mutations in gag [matrix/capsid (CA), CA/p2, p2/nucleocapsid (NC), NC/p1, p1/p6] and gag-pol (NC/p6) cleavage site and protease-coding region, with respect to the HIV-2ROD strain, were tested using Fisher's exact test.

RESULTS

Virological response occurred in 14 of 17 (82%) protease inhibitor-naive and 17 of 29 (59%) protease inhibitor-experienced patients. Virological failure was associated with higher baseline viral load (median: 6765 versus 1098 copies/ml, P = 0.02). More protease-coding region mutations were observed in protease inhibitor-experienced compared with protease inhibitor-naive patients (median: 8 versus 5, P = 0.003). In protease inhibitor-naive patients, T435A (NC/p6), V447M (p1/p6), and Y14H (protease-coding region) were associated with virological failure (P = 0.011, P = 0.033, P = 0.022, respectively). T435A and V447M were associated with Y14H (P = 0.018, P = 0.039, respectively). In protease inhibitor-experienced patients, D427E (NC/p1) was associated with virological response (P = 0.014). A430V (NC/p1) and I82F (protease-coding region) were associated with virological failure (P = 0.046, P = 0.050, respectively). Mutations at position 430 were associated with a higher number of mutations in protease-coding region (median: 10 versus 7, P = 0.008).

CONCLUSION

We have demonstrated, for the first time, an association between gag, gag-pol cleavage site and protease-coding region mutations, with distinct profiles between protease inhibitor-naive and protease inhibitor-experienced patients. These mutations might impact the virological outcome of HIV-2-infected patients receiving LPV-containing regimen.

Structural, kinetic, and thermodynamic studies of specificity designed HIV-1 protease

HIV-1 protease recognizes and cleaves more than 12 different substrates leading to viral maturation. While these substrates share no conserved motif, they are specifically selected for and cleaved by protease during viral life cycle. Drug resistant mutations evolve within the protease that compromise inhibitor binding but allow the continued recognition of all these substrates. While the substrate envelope defines a general shape for substrate recognition, successfully predicting the determinants of substrate binding specificity would provide additional insights into the mechanism of altered molecular recognition in resistant proteases. We designed a variant of HIV protease with altered specificity using positive computational design methods and validated the design using X-ray crystallography and enzyme biochemistry. The engineered variant, Pr3 (A28S/D30F/G48R), was designed to preferentially bind to one out of three of HIV protease's natural substrates; RT-RH over p2-NC and CA-p2. In kinetic assays, RT-RH binding specificity for Pr3 increased threefold compared to the wild-type (WT), which was further confirmed by isothermal titration calorimetry. Crystal structures of WT protease and the designed variant in complex with RT-RH, CA-p2, and p2-NC were determined. Structural analysis of the designed complexes revealed that one of the engineered substitutions (G48R) potentially stabilized heterogeneous flap conformations, thereby facilitating alternate modes of substrate binding. Our results demonstrate that while substrate specificity could be engineered in HIV protease, the structural pliability of protease restricted the propagation of interactions as predicted. These results offer new insights into the plasticity and structural determinants of substrate binding specificity of the HIV-1 protease.

Context surrounding processing sites is crucial in determining cleavage rate of a subset of processing sites in HIV-1 Gag and Gag-Pro-Pol polyprotein precursors by viral protease

Processing of the human immunodeficiency virus type 1 (HIV-1) Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) is essential for the production of infectious particles. However, the determinants governing the rates of processing of these substrates are not clearly understood. We studied the effect of substrate context on processing by utilizing a novel protease assay in which a substrate containing HIV-1 matrix (MA) and the N-terminal domain of capsid (CA) is labeled with a FlAsH (fluorescein arsenical hairpin) reagent. When the seven cleavage sites within the Gag and Gag-Pro-Pol polyproteins were placed at the MA/CA site, the rates of cleavage changed dramatically compared with that of the cognate sites in the natural context reported previously. The rate of processing was affected the most for three sites: CA/spacer peptide 1 (SP1) (≈10-fold increase), SP1/nucleocapsid (NC) (≈10-30-fold decrease), and SP2/p6 (≈30-fold decrease). One of two multidrug-resistant (MDR) PR variants altered the pattern of processing rates significantly. Cleavage sites within the Pro-Pol region were cleaved in a context-independent manner, suggesting for these sites that the sequence itself was the determinant of rate. In addition, a chimera consisting of SP1/NC P4-P1 and MA/CA P1'-P4' residues (ATIM↓PIVQ) abolished processing by wild type and MDR proteases, and the reciprocal chimera consisting of MA/CA P4-P1 and SP1/NC P1'-4' (SQNY↓IQKG) was cleaved only by one of the MDR proteases. These results suggest that complex substrate interactions both beyond the active site of the enzyme and across the scissile bond contribute to defining the rate of processing by the HIV-1 PR.

Critical differences in HIV-1 and HIV-2 protease specificity for clinical inhibitors

Clinical inhibitor amprenavir (APV) is less effective on HIV-2 protease (PR₂) than on HIV-1 protease (PR₁). We solved the crystal structure of PR₂ with APV at 1.5 Å resolution to identify structural changes associated with the lowered inhibition. Furthermore, we analyzed the PR₁ mutant (PR(1M) ) with substitutions V32I, I47V, and V82I that mimic the inhibitor binding site of PR₂. PR(1M) more closely resembled PR₂ than PR₁ in catalytic efficiency on four substrate peptides and inhibition by APV, whereas few differences were seen for two other substrates and inhibition by saquinavir (SQV) and darunavir (DRV). High resolution crystal structures of PR(1M) with APV, DRV, and SQV were compared with available PR₁ and PR₂ complexes. Val/Ile32 and Ile/Val47 showed compensating interactions with SQV in PR(1M) and PR₁, however, Ile82 interacted with a second SQV bound in an extension of the active site cavity of PR(1M). Residues 32 and 82 maintained similar interactions with DRV and APV in all the enzymes, whereas Val47 and Ile47 had opposing effects in the two subunits. Significantly diminished interactions were seen for the aniline of APV bound in PR₁ (M) and PR₂ relative to the strong hydrogen bonds observed in PR₁, consistent with 15- and 19-fold weaker inhibition, respectively. Overall, PR(1M) partially replicates the specificity of PR₂ and gives insight into drug resistant mutations at residues 32, 47, and 82. Moreover, this analysis provides a structural explanation for the weaker antiviral effects of APV on HIV-2.

Dynamics of preferential substrate recognition in HIV-1 protease: redefining the substrate envelope

Human immunodeficiency virus type 1 (HIV-1) protease (PR) permits viral maturation by processing the gag and gag-pro-pol polyproteins. HIV-1 PR inhibitors (PIs) are used in combination antiviral therapy but the emergence of drug resistance has limited their efficacy. The rapid evolution of HIV-1 necessitates consideration of drug resistance in novel drug design. Drug-resistant HIV-1 PR variants no longer inhibited efficiently, continue to hydrolyze the natural viral substrates. Though highly diverse in sequence, the HIV-1 PR substrates bind in a conserved three-dimensional shape we termed the substrate envelope. Earlier, we showed that resistance mutations arise where PIs protrude beyond the substrate envelope, because these regions are crucial for drug binding but not for substrate recognition. We extend this model by considering the role of protein dynamics in the interaction of HIV-1 PR with its substrates. We simulated the molecular dynamics of seven PR-substrate complexes to estimate the conformational flexibility of the bound substrates. Interdependence of substrate-protease interactions might compensate for variations in cleavage-site sequences and explain how a diverse set of sequences are recognized as substrates by the same enzyme. This diversity might be essential for regulating sequential processing of substrates. We define a dynamic substrate envelope as a more accurate representation of PR-substrate interactions. This dynamic substrate envelope, described by a probability distribution function, is a powerful tool for drug design efforts targeting ensembles of resistant HIV-1 PR variants with the aim of developing drugs that are less susceptible to resistance.

Reaction kinetics of catalyzed competitive heteropolymer cleavage

A theoretical formulation for complete heteropolymer degradation is developed in terms of Michaelis-Menten reaction kinetics under the quasi-steady-state approximation. This allows the concentration of the entire intermediate decomposition cascade to be accounted for as well as each species of emerging final product. The formulation is implemented computationally and results in stable reaction kinetics across a range of orders of magnitude for K(M) and k(cat). The model is compared with experiment, specifically in vitro HIV-1 protease-catalyzed retroviral Gag-polyprotein processing. Using an experimentally determined cleavage-polypeptide parameter set, good qualitative agreement is reached with Gag degradation kinetics, given the difference in experimental conditions. A parameter search within 1 order of magnitude of variation of the experimental set results in the determination of an optimal parameter set in complete agreement with experiment which allows the time evolution of each individual as well as intermediate species in Gag to be accurately followed. Future investigations that determine the required enzymatic parameters to populate such a scheme will allow for the model to be refined in order to track the time for viral maturation and infectivity.

The evolution of catalytic function in the HIV-1 protease

The evolution of species is a complex phenomenon based on the optimization of a multidimensional function referred to as fitness. At the level of biomolecular evolution, the fitness function can be reduced to include physiochemical properties relevant to the biological function of a particular molecule. In this work, questions involving the physical-chemical mechanisms underlying the evolution of HIV-1 protease are addressed through molecular simulation and subsequent analysis of thermodynamic properties related to the activity of the enzyme. Specifically, the impact of 40 single amino acid mutations on the binding affinity toward the matrix/capsid (MA/CA) substrate and corresponding transition state intermediate has been characterized using a molecular mechanics Poisson-Boltzmann surface area approach. We demonstrate that this approach is capable of extracting statistically significant information relevant to experimentally determined catalytic activity. Further, no correlation was observed between the effect of mutations on substrate and transition state binding, suggesting independent evolutionary pathways toward optimizing substrate specificity and catalytic activity. In addition, a detailed analysis of calculated binding affinity data suggests that ground-state destabilization (reduced binding affinity for the substrate) could be a contributing factor in the evolutionary optimization of HIV-1 protease. A numerical model is developed to demonstrate that ground-state destabilization is a valid mechanism for activity optimization given the high concentrations of substrate experienced by the functional enzyme in vivo.

The substrate specificity of Metarhizium anisopliae and Bos taurus carboxypeptidases A: insights into their use as tools for the removal of affinity tags

Carboxypeptidases may serve as tools for removal of C-terminal affinity tags. In the present study, we describe the expression and purification of an A-type carboxypeptidase from the fungal pathogen Metarhizium anisopliae (MeCPA) that has been genetically engineered to facilitate the removal of polyhistidine tags from the C-termini of recombinant proteins. A complete, systematic analysis of the specificity of MeCPA in comparison with that of bovine carboxypeptidase A (BoCPA) was carried out. Our results indicate that the specificity of the two enzymes is similar but not identical. Histidine residues are removed more efficiently by MeCPA. The very inefficient digestion of peptides with C-terminal lysine or arginine residues, along with the complete inability of the enzyme to remove a C-terminal proline, suggests a strategy for designing C-terminal affinity tags that can be trimmed by MeCPA (or BoCPA) to produce a digestion product with a homogeneous endpoint.

A switch-on mechanism to activate maize ribosome-inactivating protein for targeting HIV-infected cells

Maize ribosome-inactivating protein (RIP) is a plant toxin that inactivates eukaryotic ribosomes by depurinating a specific adenine residue at the α-sarcin/ricin loop of 28S rRNA. Maize RIP is first produced as a proenzyme with a 25-amino acid internal inactivation region on the protein surface. During germination, proteolytic removal of this internal inactivation region generates the active heterodimeric maize RIP with full N-glycosidase activity. This naturally occurring switch-on mechanism provides an opportunity for targeting the cytotoxin to pathogen-infected cells. Here, we report the addition of HIV-1 protease recognition sequences to the internal inactivation region and the activation of the maize RIP variants by HIV-1 protease in vitro and in HIV-infected cells. Among the variants generated, two were cleaved efficiently by HIV-1 protease. The HIV-1 protease-activated variants showed enhanced N-glycosidase activity in vivo as compared to their un-activated counterparts. They also possessed potent inhibitory effect on p24 antigen production in human T cells infected by two HIV-1 strains. This switch-on strategy for activating the enzymatic activity of maize RIP in target cells provides a platform for combating pathogens with a specific protease.

Identification of structural mechanisms of HIV-1 protease specificity using computational peptide docking: implications for drug resistance

Drug-resistant mutations (DRMs) in HIV-1 protease are a major challenge to antiretroviral therapy. Protease-substrate interactions that are determined to be critical for native selectivity could serve as robust targets for drug design that are immune to DRMs. In order to identify the structural mechanisms of selectivity, we developed a peptide-docking algorithm to predict the atomic structure of protease-substrate complexes and applied it to a large and diverse set of cleavable and noncleavable peptides. Cleavable peptides showed significantly lower energies of interaction than noncleavable peptides with six protease active-site residues playing the most significant role in discrimination. Surprisingly, all six residues correspond to sequence positions associated with drug resistance mutations, demonstrating that the very residues that are responsible for native substrate specificity in HIV-1 protease are altered during its evolution to drug resistance, suggesting that drug resistance and substrate selectivity may share common mechanisms.

Revealing the dimer dissociation and existence of a folded monomer of the mature HIV-2 protease

Purification and in vitro protein-folding schemes were developed to produce monodisperse samples of the mature wild-type HIV-2 protease (PR2), enabling a comprehensive set of biochemical and biophysical studies to assess the dissociation of the dimeric protease. An E37K substitution in PR2 significantly retards autoproteolytic cleavage during expression. Furthermore, it permits convenient measurement of the dimer dissociation of PR2(E37K) (elevated K(d) approximately 20 nM) by enzyme kinetics. Differential scanning calorimetry reveals a T(m) of 60.5 for PR2 as compared with 65.7 degrees C for HIV-1 protease (PR1). Consistent with weaker binding of the clinical inhibitor darunavir (DRV) to PR2, the T(m) of PR2 increases by 14.8 degrees C in the presence of DRV as compared with 22.4 degrees C for PR1. Dimer interface mutations, such as a T26A substitution in the active site (PR2(T26A)) or a deletion of the C-terminal residues 96-99 (PR2(1-95)), drastically increase the K(d) (>10(5)-fold). PR2(T26A) and PR2(1-95) consist predominantly of folded monomers, as determined by nuclear magnetic resonance (NMR) and size-exclusion chromatography coupled with multiangle light scattering and refractive index measurements (SMR), whereas wild-type PR2 and its active-site mutant PR2(D25N) are folded dimers. Addition of twofold excess active-site inhibitor promotes dimerization of PR2(T26A) but not of PR2(1-95), indicating that subunit interactions involving the C-terminal residues are crucial for dimer formation. Use of SMR and NMR with PR2 facilitates probing for potential inhibitors that restrict protein folding and/or dimerization and, thus, may provide insights for the future design of inhibitors to circumvent drug resistance.

Comparative studies on retroviral proteases: substrate specificity

Exogenous retroviruses are subclassified into seven genera and include viruses that cause diseases in humans. The viral Gag and Gag-Pro-Pol polyproteins are processed by the retroviral protease in the last stage of replication and inhibitors of the HIV-1 protease are widely used in AIDS therapy. Resistant mutations occur in response to the drug therapy introducing residues that are frequently found in the equivalent position of other retroviral proteases. Therefore, besides helping to understand the general and specific features of these enzymes, comparative studies of retroviral proteases may help to understand the mutational capacity of the HIV-1 protease.

Specificity rule discovery in HIV-1 protease cleavage site analysis

Several machine learning algorithms have recently been applied to modeling the specificity of HIV-1 protease. The problem is challenging because of the three issues as follows: (1) datasets with high dimensionality and small number of samples could misguide classification modeling and its interpretation; (2) symbolic interpretation is desirable because it provides us insight to the specificity in the form of human-understandable rules, and thus helps us to design effective HIV inhibitors; (3) the interpretation should take into account complexity or dependency between positions in sequences. Therefore, it is necessary to investigate multivariate and feature-selective methods to model the specificity and to extract rules from the model. We have tested extensively various machine learning methods, and we have found that the combination of neural networks and decompositional approach can generate a set of effective rules. By validation to experimental results for the HIV-1 protease, the specificity rules outperform the ones generated by frequency-based, univariate or black-box methods.

Mutational analysis of the C-terminal gag cleavage sites in human immunodeficiency virus type 1

Human immunodeficiency virus type 1 (HIV-1) Gag is expressed as a polyprotein that is cleaved into six proteins by the viral protease in a maturation process that begins during assembly and budding. While processing of the N terminus of Gag is strictly required for virion maturation and infectivity, the necessity for the C-terminal cleavages of Gag is less well defined. To examine the importance of this process, we introduced a series of mutations into the C terminus of Gag that interrupted the cleavage sites that normally produce in the nucleocapsid (NC), spacer 2 (SP2), or p6(Gag) proteins. Protein analysis showed that all of the mutant constructs produced virions efficiently upon transfection of cells and appropriately processed Gag polyprotein at the nonmutated sites. Mutants that produced a p9(NC/SP2) protein exhibited only minor effects on HIV-1 infectivity and replication. In contrast, mutants that produced only the p8(SP2/p6) or p15(NC/SP2/p6) protein had severe defects in infectivity and replication. To identify the key defective step, we quantified reverse transcription and integration products isolated from infected cells by PCR. All mutants tested produced levels of reverse transcription products either similar to or only somewhat lower than that of wild type. In contrast, mutants that failed to cleave the SP2-p6(Gag) site produced drastically less provirus than the wild type. Together, our results show that processing of the SP2-p6(Gag) and not the NC-SP2 cleavage site is important for efficient viral DNA integration during infection in vitro. In turn, this finding suggests an important role for the p9(NC/SP2) species in some aspect of integration.

HIV-1 protease dimer interface mutations that compensate for viral reverse transcriptase instability in infectious virions

Mature enzymes encoded within the human immunodeficiency virus type 1 (HIV-1) genome (protease (PR), reverse transcriptase (RT) and integrase (IN)) derive from proteolytic processing of a large polyprotein (Gag-Pol). Gag-Pol processing is catalyzed by the viral PR, which is active as a homodimer. The HIV-1 RT functions as a heterodimer (p66/p51) composed of subunits of 560 and 440 amino acid residues, respectively. Both subunits have identical amino acid sequence, but p51 lacks 120 residues that are removed by the HIV-1 PR during viral maturation. While p66 is the catalytic subunit, p51 has a primarily structural role. Amino acid substitutions affecting the stability of p66/p51 (i.e. F130W) have a deleterious effect on viral fitness. Previously, we showed that the effects of F130W are mediated by p51 and can be compensated by mutation T58S. While studying the dynamics of emergence of the compensatory mutation, we observed that mutations in the viral PR-coding region were selected in HIV clones containing the RT substitution F130W, before the imposition of T58S/F130W mutations. The PR mutations identified (G94S and T96S) improved the replication capacity of the F130W mutant virus. By using a trans-complementation assay, we demonstrate that the loss of p66/p51 heterodimer stability caused by Trp130 can be attributed to an increased susceptibility of RT to viral PR degradation. Recombinant HIV-1 PRs bearing mutations G94S or T96S showed decreased dimer stability and reduced catalytic efficiency. These results were consistent with crystallographic data showing the location of both residues in the PR dimerization interface.

Computational proteomics analysis of HIV-1 protease interactome

HIV-1 protease is a small homodimeric enzyme that ensures maturation of HIV virions by cleaving the viral precursor Gag and Gag-Pol polyproteins into structural and functional elements. The cleavage sites in the viral polyproteins share neither sequence homology nor binding motif and the specificity of the HIV-1 protease is therefore only partially understood. Using an extensive data set collected from 16 years of HIV proteome research we have here created a general and predictive rule-based model for HIV-1 protease specificity based on rough sets. We demonstrate that HIV-1 protease specificity is much more complex than previously anticipated, which cannot be defined based solely on the amino acids at the substrate's scissile bond or by any other single substrate amino acid position only. Our results show that the combination of at least three particular amino acids is needed in the substrate for a cleavage event to occur. Only by combining and analyzing massive amounts of HIV proteome data it was possible to discover these novel and general patterns of physico-chemical substrate cleavage determinants. Our study is an example how computational biology methods can advance the understanding of the viral interactomes.