Structural basis for coevolution of a human immunodeficiency virus type 1 nucleocapsid-p1 cleavage site with a V82A drug-resistant mutation in viral protease

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

HIV Protease Hinge Region Insertions at Codon 38 Affect Enzyme Kinetics, Conformational Stability and Dynamics

HIV-1 protease is essential for the production of mature, infectious virions and is a major target in antiretroviral therapy. We successfully purified a HIV-1 subtype C variant, L38↑N↑L- 4, containing an insertion of asparagine and leucine at position 38 without the four background mutations - K20R, E35D, R57K, V82I using a modified purification protocol. Isothermal titration calorimetry indicated that 50% of the variant protease sample was in the active conformation compared to 62% of the wild type protease. The secondary structure composition of the variant protease was unaffected by the double insertion. The specific activity and kcat values of the variant protease were approximately 50% lower than the wild type protease values. The variant protease also exhibited a 1.6-fold increase in kcat/KM when compared to the wild type protease. Differential scanning calorimetry showed a 5 °C increase in Tm of the variant protease, indicating the variant was more stable than the wild type. Molecular dynamics simulations indicated the variant was more stable and compact than the wild type protease. A 3-4% increase in the flexibility of the hinge regions of the variant protease was observed. In addition, increased flexibility of the flaps, cantilever and fulcrum regions of the variant protease B chain was observed. The variant protease sampled only the closed flap conformation indicating a potential mechanism for drug resistance. The present study highlights the direct impact of a double amino acid insertion in hinge region on enzyme kinetics, conformational stability and dynamics of an HIV-1 subtype C variant protease.

CHARMM-GUI Implicit Solvent Modeler for Various Generalized Born Models in Different Simulation Programs

Implicit solvent models are widely used because they are advantageous to speed up simulations by drastically decreasing the number of solvent degrees of freedom, which allows one to achieve long simulation time scales for large system sizes. CHARMM-GUI, a web-based platform, has been developed to support the setup of complex multicomponent molecular systems and prepare input files. This study describes an Implicit Solvent Modeler (ISM) in CHARMM-GUI for various generalized Born (GB) implicit solvent simulations in different molecular dynamics programs such as AMBER, CHARMM, GENESIS, NAMD, OpenMM, and Tinker. The GB models available in ISM include GB-HCT, GB-OBC, GB-neck, GBMV, and GBSW with the CHARMM and Amber force fields for protein, DNA, RNA, glycan, and ligand systems. Using the system and input files generated by ISM, implicit solvent simulations of protein, DNA, and RNA systems produce similar results for different simulation packages with the same input information. Protein-ligand systems are also considered to further validate the systems and input files generated by ISM. Simple ligand root-mean-square deviation (RMSD) and molecular mechanics generalized Born surface area (MM/GBSA) calculations show that the performance of implicit simulations is better than docking and can be used for early stage ligand screening. These reasonable results indicate that ISM is a useful and reliable tool to provide various implicit solvent simulation applications.

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.

Potential Resistance of SARS-CoV-2 Main Protease (Mpro) against Protease Inhibitors: Lessons Learned from HIV-1 Protease

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome 2 (SARS-CoV-2), has been one of the most devastating pandemics of recent times. The lack of potent novel antivirals had led to global health crises; however, emergence and approval of potent inhibitors of the viral main protease (Mpro), such as Pfizer's newly approved nirmatrelvir, offers hope not only in the therapeutic front but also in the context of prophylaxis against the infection. By their nature, RNA viruses including human immunodeficiency virus (HIV) have inherently high mutation rates, and lessons learnt from previous and currently ongoing pandemics have taught us that these viruses can easily escape selection pressure through mutation of vital target amino acid residues in monotherapeutic settings. In this paper, we review nirmatrelvir and its binding to SARS-CoV-2 Mpro and draw a comparison to inhibitors of HIV protease that were rendered obsolete by emergence of resistance mutations, emphasizing potential pitfalls in the design of inhibitors that may be of important relevance to the long-term use of novel inhibitors against SARS-CoV-2.

The HIV-1 Nucleocapsid Regulates Its Own Condensation by Phase-Separated Activity-Enhancing Sequestration of the Viral Protease during Maturation

A growing number of studies indicate that mRNAs and long ncRNAs can affect protein populations by assembling dynamic ribonucleoprotein (RNP) granules. These phase-separated molecular ‘sponges’, stabilized by quinary (transient and weak) interactions, control proteins involved in numerous biological functions. Retroviruses such as HIV-1 form by self-assembly when their genomic RNA (gRNA) traps Gag and GagPol polyprotein precursors. Infectivity requires extracellular budding of the particle followed by maturation, an ordered processing of ∼2400 Gag and ∼120 GagPol by the viral protease (PR). This leads to a condensed gRNA-NCp7 nucleocapsid and a CAp24-self-assembled capsid surrounding the RNP. The choreography by which all of these components dynamically interact during virus maturation is one of the missing milestones to fully depict the HIV life cycle. Here, we describe how HIV-1 has evolved a dynamic RNP granule with successive weak–strong–moderate quinary NC-gRNA networks during the sequential processing of the GagNC domain. We also reveal two palindromic RNA-binding triads on NC, KxxFxxQ and QxxFxxK, that provide quinary NC-gRNA interactions. Consequently, the nucleocapsid complex appears properly aggregated for capsid reassembly and reverse transcription, mandatory processes for viral infectivity. We show that PR is sequestered within this RNP and drives its maturation/condensation within minutes, this process being most effective at the end of budding. We anticipate such findings will stimulate further investigations of quinary interactions and emergent mechanisms in crowded environments throughout the wide and growing array of RNP granules.

Elasticity-Associated Functionality and Inhibition of the HIV Protease

HIV protease plays a critical role in the life cycle of the virus through the generation of mature and infectious virions. Detailed knowledge of the structure of the enzyme and its substrate has led to the development of protease inhibitors. However, the development of resistance to all currently available protease inhibitors has contributed greatly to the decreased success of antiretroviral therapy. When therapy failure occurs, multiple mutations are found within the protease sequence starting with primary mutations, which directly impact inhibitor binding, which can also negatively impact viral fitness and replicative capacity by decreasing the binding affinity of the natural substrates to the protease. As such, secondary mutations which are located outside of the active site region accumulate to compensate for the recurrently deleterious effects of primary mutations. However, the resistance mechanism of these secondary mutations is not well understood, but what is known is that these secondary mutations contribute to resistance in one of two ways, either through increasing the energetic penalty associated with bringing the protease into the closed conformation, or, through decreasing the stability of the protein/drug complex in a manner that increases the dissociation rate of the drug, leading to diminished inhibition. As a result, the elasticity of the enzyme-substrate complex has been implicated in the successful recognition and catalysis of the substrates which may be inferred to suggest that the elasticity of the enzyme/drug complex plays a role in resistance. A realistic representation of the dynamic nature of the protease may provide a more powerful tool in structure-based drug design algorithms.

Incorporating the Coevolving Information of Substrates in Predicting HIV-1 Protease Cleavage Sites

Human immunodeficiency virus 1 (HIV-1) protease (PR) plays a crucial role in the maturation of the virus. The study of substrate specificity of HIV-1 PR as a new endeavor strives to increase our ability to understand how HIV-1 PR recognizes its various cleavage sites. To predict HIV-1 PR cleavage sites, most of the existing approaches have been developed solely based on the homogeneity of substrate sequence information with supervised classification techniques. Although efficient, these approaches are found to be restricted to the ability of explaining their results and probably provide few insights into the mechanisms by which HIV-1 PR cleaves the substrates in a site-specific manner. In this work, a coevolutionary pattern-based prediction model for HIV-1 PR cleavage sites, namely EvoCleave, is proposed by integrating the coevolving information obtained from substrate sequences with a linear SVM classifier. The experiment results showed that EvoCleave yielded a very promising performance in terms of ROC analysis and f-measure. We also prospectively assessed the biological significance of coevolutionary patterns by applying them to study three fundamental issues of HIV-1 PR cleavage site. The analysis results demonstrated that the coevolutionary patterns offered valuable insights into the understanding of substrate specificity of HIV-1 PR.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 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.

Gag-protease coevolution shapes the outcome of lopinavir-inclusive treatment regimens in chronically infected HIV-1 subtype C patients

MOTIVATION

Commonly, protease inhibitor failure is characterized by the development of multiple protease resistance mutations (PRMs). While the impact of PRMs on therapy failure are understood, the introduction of Gag mutations with protease remains largely unclear.

RESULTS

Here, we utilized phylogenetic analyses and Bayesian network learning as tools to understand Gag-protease coevolution and elucidate the pathways leading to Lopinavir failure in HIV-1 subtype C infected patients. Our analyses indicate that while PRMs coevolve in response to drug selection pressure within protease, the Gag mutations added to the existing network while specifically interacting with known Lopinavir failure PRMs. Additionally, the selection of mutations at specific positions in Gag-protease suggests that these coevolving mutational changes occurs to maintain structural integrity during Gag cleavage.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

Amino-acid inserts of HIV-1 capsid (CA) induce CA degradation and abrogate viral infectivity: Insights for the dynamics and mechanisms of HIV-1 CA decomposition

Accumulation of amino acid (AA) insertions/substitutions are observed in the Gag-protein of HIV-1 variants resistant to HIV-1 protease inhibitors. Here, we found that HIV-1 carrying AA insertions in capsid protein (CA) undergoes aberrant CA degradation. When we generated recombinant HIV-1s (rHIV-1s) containing 19-AAs in Gag, such insertions caused significant CA degradation, which initiated in CA's C-terminal. Such rHIV-1s had remarkable morphological abnormality, decreased infectivity, and no replicative ability, which correlated with levels of CA degradation. The CA degradation observed was energy-independent and had no association with cellular/viral proteolytic mechanisms, suggesting that the CA degradation occurs due to conformational/structural incompatibility caused by the 19-AA insertions. The incorporation of degradation-prone CA into the wild-type CA resulted in significant disruption of replication competence in "chimeric" virions. The data should allow better understanding of the dynamics and mechanisms of CA decomposition/degradation and retroviral uncoating, which may lead to new approach for antiretroviral modalities.

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.

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.

Selection analyses of paired HIV-1 gag and gp41 sequences obtained before and after antiretroviral therapy

Most HIV-1-infected individuals with virological failure on a pharmacologically-boosted protease inhibitor (PI) regimen do not develop PI-resistance protease mutations. One proposed explanation is that HIV-1 gag or gp41 cytoplasmic domain mutations might also reduce PI susceptibility. In a recent study of paired gag and gp41 sequences from individuals with virological failure on a PI regimen, we did not identify PI-selected mutations and concluded that if such mutations existed, larger numbers of paired sequences from multiple studies would be needed for their identification. In this study, we generated site-specific amino acid profiles using gag and gp41 published sequences from 5,338 and 4,242 ART-naïve individuals, respectively, to assist researchers identify unusual mutations arising during therapy and to provide scripts for performing established and novel maximal likelihood estimates of dN/dS substitution rates in paired sequences. The pipelines used to generate the curated sequences, amino acid profiles, and dN/dS analyses will facilitate the application of consistent methods to paired gag and gp41 sequence datasets and expedite the identification of potential sites under PI-selection pressure.

Evolution of gag and gp41 in Patients Receiving Ritonavir-Boosted Protease Inhibitors

Several groups have proposed that genotypic determinants in gag and the gp41 cytoplasmic domain (gp41-CD) reduce protease inhibitor (PI) susceptibility without PI-resistance mutations in protease. However, no gag and gp41-CD mutations definitively responsible for reduced PI susceptibility have been identified in individuals with virological failure (VF) while receiving a boosted PI (PI/r)-containing regimen. To identify gag and gp41 mutations under selective PI pressure, we sequenced gag and/or gp41 in 61 individuals with VF on a PI/r (n = 40) or NNRTI (n = 20) containing regimen. We quantified nonsynonymous and synonymous changes in both genes and identified sites exhibiting signal for directional or diversifying selection. We also used published gag and gp41 polymorphism data to highlight mutations displaying a high selection index, defined as changing from a conserved to an uncommon amino acid. Many amino acid mutations developed in gag and in gp41-CD in both the PI- and NNRTI-treated groups. However, in neither gene, were there discernable differences between the two groups in overall numbers of mutations, mutations displaying evidence of diversifying or directional selection, or mutations with a high selection index. If gag and/or gp41 encode PI-resistance mutations, they may not be confined to consistent mutations at a few sites.

Improving Viral Protease Inhibitors to Counter Drug Resistance

Drug resistance is a major problem in health care, undermining therapy outcomes and necessitating novel approaches to drug design. Extensive studies on resistance to viral protease inhibitors, particularly those of HIV-1 and hepatitis C virus (HCV) protease, revealed a plethora of information on the structural and molecular mechanisms underlying resistance. These insights led to several strategies to improve viral protease inhibitors to counter resistance, such as exploiting the essential biological function and leveraging evolutionary constraints. Incorporation of these strategies into structure-based drug design can minimize vulnerability to resistance, not only for viral proteases but for other quickly evolving drug targets as well, toward designing inhibitors one step ahead of evolution to counter resistance with more intelligent and rational design.

Modulation of HIV protease flexibility by the T80N mutation

The flexibility of HIV protease (HIVp) plays a critical role in enabling enzymatic activity and is required for substrate access to the active site. While the importance of flexibility in the flaps that cover the active site is well known, flexibility in other parts of the enzyme is also critical for function. One key region is a loop containing Thr 80, which forms the walls of the active site. Although not situated within the active site, amino acid Thr80 is absolutely conserved. The mutation T80N preserves the structure of the enzyme but catalytic activity is completely lost. To investigate the potential influence of the T80N mutation on HIVp flexibility, wide-angle X-ray scattering (WAXS) data was measured for a series of HIVp variants. Starting with a calculated WAXS pattern from a rigid atomic model, the modulations in the intensity distribution caused by structural fluctuations in the protein were predicted by simple analytic methods and compared with the experimental data. An analysis of T80N WAXS data shows that this variant is significantly more rigid than the WT across all length scales. The effects of this single point mutation extend throughout the protein, to alter the mobility of amino acids in the enzymatic core. These results support the contentions that significant protein flexibility extends throughout HIVp and is critical to catalytic function.

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.

Function-based mutation-resistant synthetic signaling device activated by HIV-1 proteolysis

The high mutation rate of the human immunodeficiency virus type 1 (HIV-1) virus is a major problem since it evades the function of antibodies and chemical inhibitors. Here, we demonstrate a viral detection strategy based on synthetic biology principles to detect a specific viral function rather than a particular viral protein. The resistance caused by mutations can be circumvented since the mutations that cause the loss of function also incapacitate the virus. Many pathogens encode proteases that are essential for their replication and that have a defined substrate specificity. A genetically encoded sensor composed of a fused membrane anchor, viral protease target site, and an orthogonal transcriptional activator was engineered into a human cell line. The HIV-1 protease released the transcriptional activator from the membrane, thereby inducing transcription of the selected genes. The device was still strongly activated by clinically relevant protease mutants that are resistant to protease inhibitors. In the future, a similar principle could be applied to detect also other pathogens and functions.

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.

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.

HIV-1 protease-substrate coevolution in nelfinavir resistance

Resistance to various human immunodeficiency virus type 1 (HIV-1) protease inhibitors (PIs) challenges the effectiveness of therapies in treating HIV-1-infected individuals and AIDS patients. The virus accumulates mutations within the protease (PR) that render the PIs less potent. Occasionally, Gag sequences also coevolve with mutations at PR cleavage sites contributing to drug resistance. In this study, we investigated the structural basis of coevolution of the p1-p6 cleavage site with the nelfinavir (NFV) resistance D30N/N88D protease mutations by determining crystal structures of wild-type and NFV-resistant HIV-1 protease in complex with p1-p6 substrate peptide variants with L449F and/or S451N. Alterations of residue 30's interaction with the substrate are compensated by the coevolving L449F and S451N cleavage site mutations. This interdependency in the PR-p1-p6 interactions enhances intermolecular contacts and reinforces the overall fit of the substrate within the substrate envelope, likely enabling coevolution to sustain substrate recognition and cleavage in the presence of PR resistance mutations.

IMPORTANCE

Resistance to human immunodeficiency virus type 1 (HIV-1) protease inhibitors challenges the effectiveness of therapies in treating HIV-1-infected individuals and AIDS patients. Mutations in HIV-1 protease selected under the pressure of protease inhibitors render the inhibitors less potent. Occasionally, Gag sequences also mutate and coevolve with protease, contributing to maintenance of viral fitness and to drug resistance. In this study, we investigated the structural basis of coevolution at the Gag p1-p6 cleavage site with the nelfinavir (NFV) resistance D30N/N88D protease mutations. Our structural analysis reveals the interdependency of protease-substrate interactions and how coevolution may restore substrate recognition and cleavage in the presence of protease drug resistance mutations.

HIV-1 assembly, budding, and maturation

A defining property of retroviruses is their ability to assemble into particles that can leave producer cells and spread infection to susceptible cells and hosts. Virion morphogenesis can be divided into three stages: assembly, wherein the virion is created and essential components are packaged; budding, wherein the virion crosses the plasma membrane and obtains its lipid envelope; and maturation, wherein the virion changes structure and becomes infectious. All of these stages are coordinated by the Gag polyprotein and its proteolytic maturation products, which function as the major structural proteins of the virus. Here, we review our current understanding of the mechanisms of HIV-1 assembly, budding, and maturation, starting with a general overview and then providing detailed descriptions of each of the different stages of virion morphogenesis.

Improved Generalized Born Solvent Model Parameters for Protein Simulations

The generalized Born (GB) model is one of the fastest implicit solvent models and it has become widely adopted for Molecular Dynamics (MD) simulations. This speed comes with tradeoffs, and many reports in the literature have pointed out weaknesses with GB models. Because the quality of a GB model is heavily affected by empirical parameters used in calculating solvation energy, in this work we have refit these parameters for GB-Neck, a recently developed GB model, in order to improve the accuracy of both the solvation energy and effective radii calculations. The data sets used for fitting are significantly larger than those used in the past. Comparing to other pairwise GB models like GB-OBC and the original GB-Neck, the new GB model (GB-Neck2) has better agreement to Poisson-Boltzmann (PB) in terms of reproducing solvation energies for a variety of systems ranging from peptides to proteins. Secondary structure preferences are also in much better agreement with those obtained from explicit solvent MD simulations. We also obtain near-quantitative reproduction of experimental structure and thermal stability profiles for several model peptides with varying secondary structure motifs. Extension to non-protein systems will be explored in the future.

Structural and thermodynamic basis of amprenavir/darunavir and atazanavir resistance in HIV-1 protease with mutations at residue 50

Drug resistance occurs through a series of subtle changes that maintain substrate recognition but no longer permit inhibitor binding. In HIV-1 protease, mutations at I50 are associated with such subtle changes that confer differential resistance to specific inhibitors. Residue I50 is located at the protease flap tips, closing the active site upon ligand binding. Under selective drug pressure, I50V/L substitutions emerge in patients, compromising drug susceptibility and leading to treatment failure. The I50V substitution is often associated with amprenavir (APV) and darunavir (DRV) resistance, while the I50L substitution is observed in patients failing atazanavir (ATV) therapy. To explain how APV, DRV, and ATV susceptibility are influenced by mutations at residue 50 in HIV-1 protease, structural and binding thermodynamics studies were carried out on I50V/L-substituted protease variants in the compensatory mutation A71V background. Reduced affinity to both I50V/A71V and I50L/A71V double mutants is largely due to decreased binding entropy, which is compensated for by enhanced enthalpy for ATV binding to I50V variants and APV binding to I50L variants, leading to hypersusceptibility in these two cases. Analysis of the crystal structures showed that the substitutions at residue 50 affect how APV, DRV, and ATV bind the protease with altered van der Waals interactions and that the selection of I50V versus I50L is greatly influenced by the chemical moieties at the P1 position for APV/DRV and the P2 position for ATV. Thus, the varied inhibitor susceptibilities of I50V/L protease variants are largely a direct consequence of the interdependent changes in protease inhibitor interactions.

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.

Mutations in multiple domains of Gag drive the emergence of in vitro resistance to the phosphonate-containing HIV-1 protease inhibitor GS-8374

GS-8374 is a potent HIV protease inhibitor (PI) with a unique diethyl-phosphonate moiety. Due to a balanced contribution of enthalpic and entropic components to its interaction with the protease (PR) active site, the compound retains activity against HIV mutants with high-level multi-PI resistance. We report here the in vitro selection and characterization of HIV variants resistant to GS-8374. While highly resistant viruses with multiple mutations in PR were isolated in the presence of control PIs, an HIV variant displaying moderate (14-fold) resistance to GS-8374 was generated only after prolonged passaging for >300 days. The isolate showed low-level cross-resistance to darunavir, atazanavir, lopinavir, and saquinavir, but not other PIs, and contained a single R41K mutation in PR combined with multiple genotypic changes in the Gag matrix, capsid, nucleocapsid, and SP2 domains. Mutations also occurred in the transframe peptide and p6* domain of the Gag-Pol polyprotein. Analysis of recombinant HIV variants indicated that mutations in Gag, but not the R41K in PR, conferred reduced susceptibility to GS-8374. The Gag mutations acted in concert, since they did not affect susceptibility when introduced individually. Analysis of viral particles revealed that the mutations rendered Gag more susceptible to PR-mediated cleavage in the presence of GS-8374. In summary, the emergence of resistance to GS-8374 involved a combination of substrate mutations without typical resistance mutations in PR. These substrate changes were distributed throughout Gag and acted in an additive manner. Thus, they are classified as primary resistance mutations indicating a unique mechanism and pathway of resistance development for GS-8374.

The choreography of HIV-1 proteolytic processing and virion assembly

HIV-1 has been the target of intensive research at the molecular and biochemical levels for >25 years. Collectively, this work has led to a detailed understanding of viral replication and the development of 24 approved drugs that have five different targets on various viral proteins and one cellular target (CCR5). Although most drugs target viral enzymatic activities, our detailed knowledge of so much of the viral life cycle is leading us into other types of inhibitors that can block or disrupt protein-protein interactions. Viruses have compact genomes and employ a strategy of using a small number of proteins that can form repeating structures to enclose space (i.e. condensing the viral genome inside of a protein shell), thus minimizing the need for a large protein coding capacity. This creates a relatively small number of critical protein-protein interactions that are essential for viral replication. For HIV-1, the Gag protein has the role of a polyprotein precursor that contains all of the structural proteins of the virion: matrix, capsid, spacer peptide 1, nucleocapsid, spacer peptide 2, and p6 (which contains protein-binding domains that interact with host proteins during budding). Similarly, the Gag-Pro-Pol precursor encodes most of the Gag protein but now includes the viral enzymes: protease, reverse transcriptase (with its associated RNase H activity), and integrase. Gag and Gag-Pro-Pol are the substrates of the viral protease, which is responsible for cleaving these precursors into their mature and fully active forms (see Fig. 1A).

Extreme entropy-enthalpy compensation in a drug-resistant variant of HIV-1 protease

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.

Human Immunodeficiency Virus Gag and protease: partners in resistance

Human Immunodeficiency Virus (HIV) maturation plays an essential role in the viral life cycle by enabling the generation of mature infectious virus particles through proteolytic processing of the viral Gag and GagPol precursor proteins. An impaired polyprotein processing results in the production of non-infectious virus particles. Consequently, particle maturation is an excellent drug target as exemplified by inhibitors specifically targeting the viral protease (protease inhibitors; PIs) and the experimental class of maturation inhibitors that target the precursor Gag and GagPol polyproteins. Considering the different target sites of the two drug classes, direct cross-resistance may seem unlikely. However, coevolution of protease and its substrate Gag during PI exposure has been observed both in vivo and in vitro. This review addresses in detail all mutations in Gag that are selected under PI pressure. We evaluate how polymorphisms and mutations in Gag affect PI therapy, an aspect of PI resistance that is currently not included in standard genotypic PI resistance testing. In addition, we consider the consequences of Gag mutations for the development and positioning of future maturation inhibitors.

Protease-Mediated Maturation of HIV: Inhibitors of Protease and the Maturation Process

Protease-mediated maturation of HIV-1 virus particles is essential for virus infectivity. Maturation occurs concomitant with immature virus particle release and is mediated by the viral protease (PR), which sequentially cleaves the Gag and Gag-Pol polyproteins into mature protein domains. Maturation triggers a second assembly event that generates a condensed conical capsid core. The capsid core organizes the viral RNA genome and viral proteins to facilitate viral replication in the next round of infection. The fundamental role of proteolytic maturation in the generation of mature infectious particles has made it an attractive target for therapeutic intervention. Development of small molecules that target the PR active site has been highly successful and nine protease inhibitors (PIs) have been approved for clinical use. This paper provides an overview of their development and clinical use together with a discussion of problems associated with drug resistance. The second-half of the paper discusses a novel class of antiretroviral drug termed maturation inhibitors, which target cleavage sites in Gag not PR itself. The paper focuses on bevirimat (BVM) the first-in-class maturation inhibitor: its mechanism of action and the implications of naturally occurring polymorphisms that confer reduced susceptibility to BVM in phase II clinical trials.

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.

Modulation of HIV-1 Gag NC/p1 cleavage efficiency affects protease inhibitor resistance and viral replicative capacity

Background

Mutations in the substrate of HIV-1 protease, especially changes in the NC/p1 cleavage site, can directly contribute to protease inhibitor (PI) resistance and also compensate for defects in viral replicative capacity (RC) due to a drug resistant protease. These NC/p1 changes are known to enhance processing of the Gag protein. To investigate the capacity of HIV-1 to modulate Gag cleavage and its consequences for PI resistance and RC, we performed a detailed enzymatic and virological analysis using a set of PI resistant NC/p1 variants (HXB2^431V, HXB2^436E+437T, HXB2^437T and HXB2^437V).

Results

Here, we demonstrate that single NC/p1 mutants, which displayed only a slight increase in PI resistance did not show an obvious change in RC. In contrast, the double NC/p1 mutant, which displayed a clear increase in processing efficiency and PI resistance, demonstrated a clear reduction in RC. Cleavage analysis showed that a tridecameric NC/p1 peptide representing the double NC/p1 mutant was cleaved in two specific ways instead of one.

The observed decrease in RC for the double NC/p1 mutant (HXB2^436E+437T) could (partially) be restored by either reversion of the 436E change or by acquisition of additional changes in the NC/p1 cleavage site at codon 435 or 438 as was revealed during in vitro evolution experiments. These changes not only restored RC but also reduced PI resistance levels. Furthermore these changes normalized Gag processing efficiency and obstructed the novel secondary cleavage site observed for the double NC/p1 mutant.

Conclusions

The results of this study clearly demonstrate that HIV-1 can modulate Gag processing and thereby PI resistance. Distinct increases in Gag cleavage and PI resistance result in a reduced RC that can only be restored by amino acid changes in NC/p1 which reduce Gag processing to an optimal rate.

HIV-1 Protease and Substrate Coevolution Validates the Substrate Envelope As the Substrate Recognition Pattern

Drug resistance of HIV-1 protease alters the balance in the molecular recognition events in favor of substrate processing versus inhibitor binding. To develop robust inhibitors targeting ensembles of drug-resistant variants, the code of this balance needs to be cracked. For this purpose, the principles governing the substrate recognition are required to be revealed. Previous crystallographic studies on the WT protease-substrate complexes showed that the substrates have a conserved consensus volume in the protease active site despite their low sequence homology. This consensus volume is termed as the substrate envelope. The substrate envelope was recently reevaluated by taking the substrate dynamics into account, and the dynamic substrate envelope was reported to better define the substrate specificity for HIV-1 protease. Drug resistance occurs mostly through mutations in the protease, occasionally accompanied by cleavage site mutations. In this study, three coevolved protease-substrate complexes (^AP2VNC-p1_V82A, ^LP1'Fp1-p6_D30N/N88D, and ^SP3'Np1-p6_D30N/N88D) were investigated for structural and dynamic properties by molecular modeling and dynamics simulations. The results show the substrate envelope is preserved by these cleavage site mutations in the presence of drug-resistance mutations in the protease, if not enhanced. This study on the conformational and mutational ensembles of protease-substrate complexes validates the substrate envelope as the substrate recognition motif for HIV-1 protease. The substrate envelope hypothesis allows for the elucidation of possible drug resistance mutation patterns in the polyprotein cleavage sites.

Mutational patterns in the frameshift-regulating site of HIV-1 selected by protease inhibitors

Sustained suppression of viral replication in HIV-1 infected patients is especially hampered by the emergence of HIV-1 drug resistance. The mechanisms of drug resistance mainly involve mutations directly altering the interaction of viral enzymes and inhibitors. However, protease inhibitors do not only select for mutations in the protease but also for mutations in the precursor Gag and Pol proteins. In this study, we analysed the frameshift-regulating site of HIV-1 subtype B isolates, which also encodes for Gag and Pol proteins, classified as either treatment-naïve (TN) or protease inhibitor resistant (PI-R). HIV-1 Gag cleavage site mutations (G435E, K436N, I437V, L449F/V) especially correlated with protease inhibitor resistance mutations, but also Pol cleavage site mutations (D05G, D05S) could be assigned to specific protease resistance profiles. Additionally, two Gag non-cleavage site mutations (S440F, H441P) were observed more often in HIV-1 isolates carrying protease resistance mutations. However, in dual luciferase assays, the frameshift efficiencies of specific clones did not reveal any effect from these mutations. Nevertheless, two patterns of mutations modestly increased the frameshift rates in vitro, but were not specifically accumulating in PI-resistant HIV-1 isolates. In summary, HIV-1 Gag cleavage site mutations were dominantly selected in PI-resistant HIV-1 isolates but also Pol cleavage site mutations influenced resistance profiles in the protease. Additionally, Gag non-cleavage site mutations accumulated in PI-resistant HIV-1 isolates, but were not related to an increased frameshift efficiency.

Characterization of HIV-1 from patients with virological failure to a boosted protease inhibitor regimen

The use of highly active antiretroviral treatment (HAART) regimens with unboosted protease inhibitors (PIs) has resulted in a high level of virological failure primarily due to the development of resistant virus. Current boosted PI regimens combine successfully low-dose ritonavir (r) with a second PI. The aim of the study was to estimate the proportion of patients, in a population based setting, who develop virological failure on a PI/r regimen. Through The Danish HIV Cohort Study 1,007 patients who received PI/r based treatment between 1995 and 2008 were identified. Twenty-three (2.3%) experienced virological failure, of whom 19 (83%) started PI/r treatment before 2001. Patients from Copenhagen (n=19) were selected to study the development of protease (PR) and gag cleavage site (CS) mutations during PI/r treatment and PI plasma levels at the time of virological failure. Three patients (16%) developed major PI resistance mutations. Mutations in the p7/p1 and p1/p6 gag CS only developed in patients with major or minor mutations in PR. Drug concentrations were low or undetectable in 10 out of the 19 patients. In total PR resistance mutations and low drug levels could account for 12 (63%) of the failure cases. In conclusion, virological failure to PI/r is a low and decreasing problem primarily caused by low plasma drug levels and to a lesser extent major PR mutations. Gag CS mutations did not contribute significantly to resistance development and virological failure.

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.

Nine crystal structures determine the substrate envelope of the MDR HIV-1 protease

Under drug selection pressure, emerging mutations render HIV-1 protease drug resistant, leading to the therapy failure in anti-HIV treatment. It is known that nine substrate cleavage site peptides bind to wild type (WT) HIV-1 protease in a conserved pattern. However, how the multidrug-resistant (MDR) HIV-1 protease binds to the substrate cleavage site peptides is yet to be determined. MDR769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90) was selected for present study to understand the binding to its natural substrates. MDR769 HIV-1 protease was co-crystallized with nine substrate cleavage site hepta-peptides. Crystallographic studies show that MDR769 HIV-1 protease has an expanded substrate envelope with wide open flaps. Furthermore, ligand binding energy calculations indicate weaker binding in MDR769 HIV-1 protease-substrate complexes. These results help in designing the next generation of HIV-1 protease inhibitors by targeting the MDR HIV-1 protease.

Clinical management of HIV drug resistance

Combination antiretroviral therapy for HIV-1 infection has resulted in profound reductions in viremia and is associated with marked improvements in morbidity and mortality. Therapy is not curative, however, and prolonged therapy is complicated by drug toxicity and the emergence of drug resistance. Management of clinical drug resistance requires in depth evaluation, and includes extensive history, physical examination and laboratory studies. Appropriate use of resistance testing provides valuable information useful in constructing regimens for treatment-experienced individuals with viremia during therapy. This review outlines the emergence of drug resistance in vivo, and describes clinical evaluation and therapeutic options of the individual with rebound viremia during therapy.

Within-host co-evolution of Gag P453L and protease D30N/N88D demonstrates virological advantage in a highly protease inhibitor-exposed HIV-1 case

To better understand the mechanism of HIV group-specific antigen (Gag) and protease (PR) co-evolution in drug-resistance acquisition, we analyzed a drug-resistance case by both bioinformatics and virological methods. We especially considered the quality of sequence data and analytical accuracy by introducing single-genome sequencing (SGS) and Spidermonkey/Bayesian graphical models (BGM) analysis, respectively. We analyzed 129 HIV-1 Gag-PR linkage sequences obtained from 8 time points, and the resulting sequences were applied to the Spidermonkey co-evolution analysis program, which identified ten mutation pairs as significantly co-evolving. Among these, we focused on associations between Gag-P453L, the P5' position of the p1/p6 cleavage-site mutation, and PR-D30N/N88D nelfinavir-resistant mutations, and attempted to clarify their virological significance in vitro by constructing recombinant clones. The results showed that P453L(Gag) has the potential to improve replication capacity and the Gag processing efficiency of viruses with D30N(PR)/N88D(PR) but has little effect on nelfinavir susceptibility. Homology modeling analysis suggested that hydrogen bonds between the 30th PR residue and the R452Gag are disturbed by the D30N(PR) mutation, but the impaired interaction is compensated by P453L(Gag) generating new hydrophobic interactions. Furthermore, database analysis indicated that the P453L(Gag)/D30N(PR)/N88D(PR) association was not specific only to our clinical case, but was common among AIDS patients.

Positive impact of HIV-1 gag cleavage site mutations on the virological response to darunavir boosted with ritonavir

We assessed the roles of baseline gag and gag-pol cleavage site mutations (CSM) on the virological outcome of a darunavir-based regimen in highly antiretroviral-experienced patients. We showed the association, in multivariate analysis, between the A431V gag CSM and the virological response, defined as a reduction in plasma HIV-1 RNA to <50 copies/ml at month 3 (P = 0.028). Our results suggest that a specific gag CSM might have a role on protease inhibitor susceptibility in an inhibitor-specific manner.

Explicit solvent dynamics and energetics of HIV-1 protease flap opening and closing

An accurate description of the conformational dynamics of the β-hairpin flaps of HIV-1 protease is of central importance in elucidating the functional recognition of the enzyme by ligands. Using all-atom molecular dynamics simulations in explicit solvent, with a total of 461 trajectories of ∼50 ns each, we report the closed, semiopen, open, and wide-open flap conformation of the free wild-type protease. The free energy of flap opening and closing from the semiopen state is 0.9 ± 0.2 and 2.4 ± 0.4 kcal/mol, respectively. The mean relaxation time of opening is ∼8 ns, in good agreement with NMR data. The explicit solvent simulations quantitatively confirm the hypothesis that the semiopen state is the dominant population in the free protease whilst fast flap tip fluctuations lead frequently to an open state. More pronounced flap rearrangements lead to a rare wide-open state with the catalytic site completely exposed to the solvent. The structures of the different flap conformations provided herein are of general interest for improved drug design of HIV-1 protease, in particular, the wide-open conformation could be favored by the large Gag and GagPol polyprotein chains. Strategies that take into account multiple flap-gating mechanisms may lead to more effective inhibitors.

Three residues in HIV-1 matrix contribute to protease inhibitor susceptibility and replication capacity

Other than cleavage site mutations, there is little data on specific positions within Gag that impact on HIV protease inhibitor susceptibility. We have recently shown that non-cleavage site mutations in gag, particularly within matrix protein can restore replication capacity and further reduce protease inhibitor drug susceptibility when coexpressed with a drug-resistant (mutant) protease. The matrix protein of this patient-derived virus was studied in order to identify specific changes responsible for this phenotype. Three amino acid changes in matrix (R76K, Y79F, and T81A) had an impact on replication capacity as well as drug susceptibility. Introduction of these three changes into wild-type (WT) matrix resulted in an increase in the replication capacity of the protease mutant virus to a level similar to that achieved by all the changes within the mutant matrix and part of the capsid protein. Pairs of changes to wild-type matrix led to an increased replication capacity of the protease mutant (although less than with all three changes). Having only these three changes to matrix in a wild-type virus (with wild-type protease) resulted in a 5- to 7-fold change in protease inhibitor 50% effective concentration (EC₅₀). Individual changes did not have as great an effect on replication capacity or drug susceptibility, demonstrating an interaction between these positions, also confirmed by sequence covariation analysis. Molecular modeling predicts that each of the three mutations would result in a loss of hydrogen bonds within α-helix-4 of matrix, leading to the hypothesis that more flexibility within this region or altered matrix structure would account for our findings.

Molecular Basis for Drug Resistance in HIV-1 Protease

HIV-1 protease is one of the major antiviral targets in the treatment of patients infected with HIV-1. The nine FDA approved HIV-1 protease inhibitors were developed with extensive use of structure-based drug design, thus the atomic details of how the inhibitors bind are well characterized. From this structural understanding the molecular basis for drug resistance in HIV-1 protease can be elucidated. Selected mutations in response to therapy and diversity between clades in HIV-1 protease have altered the shape of the active site, potentially altered the dynamics and even altered the sequence of the cleavage sites in the Gag polyprotein. All of these interdependent changes act in synergy to confer drug resistance while simultaneously maintaining the fitness of the virus. New strategies, such as incorporation of the substrate envelope constraint to design robust inhibitors that incorporate details of HIV-1 protease’s function and decrease the probability of drug resistance, are necessary to continue to effectively target this key protein in HIV-1 life cycle.

Rationale for more diverse inhibitors in competition with substrates in HIV-1 protease

The structural fluctuations of HIV-1 protease in interaction with its substrates versus inhibitors were analyzed using the anisotropic network model. The directions of fluctuations in the most cooperative functional modes differ mainly around the dynamically key regions, i.e., the hinge axes, which appear to be more flexible in substrate complexes. The flexibility of HIV-1 protease is likely optimized for the substrates' turnover, resulting in substrate complexes being dynamic. In contrast, in an inhibitor complex, the inhibitor should bind and lock down to inactivate the active site. Protease and ligands are not independent. Substrates are also more flexible than inhibitors and have the potential to meet the dynamic distributions that are inherent in the protease. This may suggest a rationale and guidelines for designing inhibitors that can better fit the ensemble of binding sites that are dynamically accessible to the protease.

The effect of clade-specific sequence polymorphisms on HIV-1 protease activity and inhibitor resistance pathways

The majority of HIV-1 infections around the world result from non-B clade HIV-1 strains. The CRF01_AE (AE) strain is seen principally in Southeast Asia. AE protease differs by approximately 10% in amino acid sequence from clade B protease and carries several naturally occurring polymorphisms that are associated with drug resistance in clade B. AE protease has been observed to develop resistance through a nonactive-site N88S mutation in response to nelfinavir (NFV) therapy, whereas clade B protease develops both the active-site mutation D30N and the nonactive-site mutation N88D. Structural and biochemical studies were carried out with wild-type and NFV-resistant clade B and AE protease variants. The relationship between clade-specific sequence variations and pathways to inhibitor resistance was also assessed. AE protease has a lower catalytic turnover rate than clade B protease, and it also has weaker affinity for both NFV and darunavir (DRV). This weaker affinity may lead to the nonactive-site N88S variant in AE, which exhibits significantly decreased affinity for both NFV and DRV. The D30N/N88D mutations in clade B resulted in a significant loss of affinity for NFV and, to a lesser extent, for DRV. A comparison of crystal structures of AE protease shows significant structural rearrangement in the flap hinge region compared with those of clade B protease and suggests insights into the alternative pathways to NFV resistance. In combination, our studies show that sequence polymorphisms within clades can alter protease activity and inhibitor binding and are capable of altering the pathway to inhibitor resistance.