Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

Quantum mechanical analysis of newly synthesized HIV-1 protease inhibitors: evaluation of wild-type and resistant strain binding interactions

Inhibition of HIV-1 protease is a cornerstone of antiretroviral therapy. However, the notorious ability of HIV-1 to develop resistance to protease inhibitors (PIs), particularly darunavir (DRV), poses a major challenge. Using quantum chemistry and computer simulations, this study aims to investigate the interactions between two novel PIs, GRL-004 and GRL-063, as well as a wild-type (WT) HIV strain and a DRV-resistant mutant strain. To do this, we used molecular docking, molecular dynamics simulations, and quantum mechanical calculations to check how well GRL-004 and GRL-063 bound to both WT and DRV-resistant proteases. The results show that GRL-004 and GRL-063 bind very well to ASP29 in the WT strain. ASP29 is an important amino acid in the HIV protease dimer. Remarkably, amino acids such as ILE50 in the WT strains showed substantial binding energies to both drugs. Quantum energy calculations showed a slight reduction in the energy affinity of the interaction between the MUT strain and the ligand GRL-063, compared to the WT strain. GRL-004 showed similar interaction energy for both strains, suggesting that it has greater plasticity than GRL-063 despite its lower interaction affinity. Furthermore, GLY49B demonstrated strong binding energies regardless of mutations. Other relevant residues with strong binding energies include GLY49B, PHE82A, PRO81A, ASP29A, ASP25A and ALA28B. This study improves our understanding of receptor-ligand dynamics and the adaptability of new protease inhibitors (PIs), which has profound implications for the innovation of future antiretroviral drugs.

Interplay between protease and reverse transcriptase dimerization in a model HIV-1 polyprotein

The Gag-Pol polyprotein in human immunodeficiency virus type I (HIV-1) encodes enzymes that are essential for virus replication: protease (PR), reverse transcriptase (RT), and integrase (IN). The mature forms of PR, RT and IN are homodimer, heterodimer and tetramer, respectively. The precise mechanism underlying the formation of dimer or tetramer is not yet understood. Here, to gain insight into the dimerization of PR and RT in the precursor, we prepared a model precursor, PR-RT, incorporating an inactivating mutation at the PR active site, D25A, and including two residues in the p6* region, fused to a SUMO-tag, at the N-terminus of the PR region. We also prepared two mutants of PR-RT containing a dimer dissociation mutation either in the PR region, PR(T26A)-RT, or in the RT region, PR-RT(W401A). Size exclusion chromatography showed both monomer and dimer fractions in PR-RT and PR(T26A)-RT, but only monomer in PR-RT(W401A). SEC experiments of PR-RT in the presence of protease inhibitor, darunavir, significantly enhanced the dimerization. Additionally, SEC results suggest an estimated PR-RT dimer dissociation constant that is higher than that of the mature RT heterodimer, p66/p51, but slightly lower than the premature RT homodimer, p66/p66. Reverse transcriptase assays and RT maturation assays were performed as tools to assess the effects of the PR dimer-interface on these functions. Our results consistently indicate that the RT dimer-interface plays a crucial role in the dimerization in PR-RT, whereas the PR dimer-interface has a lesser role.

Natural Polymorphisms D60E and I62V Stabilize a Closed Conformation in HIV-1 Protease in the Absence of an Inhibitor or Substrate

HIV infection remains a global health issue plagued by drug resistance and virological failure. Natural polymorphisms (NPs) contained within several African and Brazilian protease (PR) variants have been shown to induce a conformational landscape of more closed conformations compared to the sequence of subtype B prevalent in North America and Western Europe. Here we demonstrate through experimental pulsed EPR distance measurements and molecular dynamic (MD) simulations that the two common NPs D60E and I62V found within subtypes F and H can induce a closed conformation when introduced into HIV-1PR subtype B. Specifically, D60E alters the conformation in subtype B through the formation of a salt bridge with residue K43 contained within the nexus between the flap and hinge region of the HIV-1 PR fold. On the other hand, I62V modulates the packing of the hydrophobic cluster of the cantilever and fulcrum, also resulting in a more closed conformation.

Dynamics of Low-Level Viremia and Immune Activation after Switching to a Darunavir-Based Regimen

There is an ongoing debate regarding whether low-level viremia (LLV), in particular persistent LLV, during HIV treatment with optimal adherence originates from low-level viral replication, viral production, or both. We performed an observational study in 30 individuals with LLV who switched to a boosted darunavir (DRV)-based therapy. In-depth virological analyses were used to characterize the viral population and the (activity) of the viral reservoir. Immune activation was examined using cell-bound and soluble markers. The primary outcome was defined as the effect on HIV-RNA and was categorized by responders (<50 cp/mL) or non-responders (>50 cp/mL). At week 24, 53% of the individuals were considered responders, 40% non-responders, and 7% could not be assigned. Sequencing showed no evolution or selection of drug resistance in the non-responders. Production of defective virus with mutations in either the protease (D25N) or RT active site contributed to persistent LLV in two individuals. We show that in about half of the study participants, the switch to a DRV-based regimen resulted in a viral response indicative of ongoing low-level viral replication as the cause of LLV before the switch. Our data confirm that in clinical management, high genetic barrier drugs like DRV are a safe choice, irrespective of the source of LLV.

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

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

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

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

The structural, dynamic, and thermodynamic basis of darunavir resistance of a heavily mutated HIV-1 protease using molecular dynamics simulation

The human immunodeficiency virus type 1 protease (HIV-1 PR) is an important enzyme in the life cycle of the HIV virus. It cleaves inactive pre-proteins of the virus and changes them into active proteins. Darunavir (DRV) suppresses the wild-type HIV-1 PR (WT-Pr) activity but cannot inhibit some mutant resistant forms (MUT-Pr). Increasing knowledge about the resistance mechanism can be helpful for designing more effective inhibitors. In this study, the mechanism of resistance of a highly MUT-Pr strain against DRV was investigated. For this purpose, complexes of DRV with WT-Pr (WT-Pr-D) and MUT-Pr (MUT-Pr-D) were studied by all-atom molecular dynamics simulation in order to extract the dynamic and energetic properties. Our data revealed that mutations increased the flap-tip flexibility due to the reduction of the flap-flap hydrophobic interactions. So, the protease’s conformation changed from a closed state to a semi-open state that can facilitate the disjunction of DRV from the active site. On the other hand, energy analysis limited to the final basins of the energy landscape indicated that the entropy of binding of DRV to MUT-Pr was more favorable than that of WT-Pr. However, the enthalpy penalty overcomes it and makes binding more unfavorable relative to the WT-Pr. The unfavorable interaction of DRV with R8, I50, I84, D25′, and A28′ residues in MUT-Pr-D relative to WT-Pr-D is the reason for this enthalpy penalty. Thus, mutations drive resistance to DRV. The hydrogen bond analysis showed that compared with WT-Pr, the hydrogen bonds between DRV and the active-site residues of MUT-Pr were disrupted.

Non-active site mutations in the HIV protease: Diminished drug binding affinity is achieved through modulating the hydrophobic sliding mechanism

The global HIV/AIDS epidemic still currently affects approximately 38 million individuals globally. The protease enzyme of the human immunodeficiency virus is a major drug target in antiviral therapy, however, under the influence of reverse transcriptase and in the context of drug pressure, the rapid PR mutation rate contributes significantly to clinical failure. The set of cooperative non-active site mutations, I13V/I62V/V77I, have been associated with reduced inhibitor susceptibility and are the focus of the current study. When compared to the wild-type protease the mutant protease exhibited decreased binding affinities towards ATV and DRV by 64- and 12-fold, respectively, and decreased the overall favourable Gibbs free energy for ATV, DRV, RTV and SQV. Moreover, these mutations decreased the thermal stability of the protease when in complex with ATV and DRV by approximately 6.4 and 4.2 °C, respectively. The crystal structure of the mutant protease revealed that the location of these mutations and their effect on the hydrophobic sliding mechanism may be crucial in their role in resistance.

Co-evolution of drug resistance and broadened substrate recognition in HIV protease variants isolated from an Escherichia coli genetic selection system

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate L-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself.

DEER-PREdict: Software for efficient calculation of spin-labeling EPR and NMR data from conformational ensembles

Owing to their plasticity, intrinsically disordered and multidomain proteins require descriptions based on multiple conformations, thus calling for techniques and analysis tools that are capable of dealing with conformational ensembles rather than a single protein structure. Here, we introduce DEER-PREdict, a software program to predict Double Electron-Electron Resonance distance distributions as well as Paramagnetic Relaxation Enhancement rates from ensembles of protein conformations. DEER-PREdict uses an established rotamer library approach to describe the paramagnetic probes which are bound covalently to the protein.DEER-PREdict has been designed to operate efficiently on large conformational ensembles, such as those generated by molecular dynamics simulation, to facilitate the validation or refinement of molecular models as well as the interpretation of experimental data. The performance and accuracy of the software is demonstrated with experimentally characterized protein systems: HIV-1 protease, T4 Lysozyme and Acyl-CoA-binding protein. DEER-PREdict is open source (GPLv3) and available at github.com/KULL-Centre/DEERpredict and as a Python PyPI package pypi.org/project/DEERPREdict.

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

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

Probing the Interaction between HIV-1 Protease and the Homodimeric p66/p66' Reverse Transcriptase Precursor by Double Electron-Electron Resonance EPR Spectroscopy

Following excision from the Gag-Pol polyprotein, HIV-1 reverse transcriptase is released as an asymmetric homodimer comprising two p66 subunits that are structurally dissimilar but identical in amino acid sequence. Subsequent cleavage of the RNase H domain from only one of the subunits, denoted p66', results in the formation of the mature p66/p51 enzyme in which catalytic activity resides in the p66 subunit, and the p51 subunit (derived from p66') provides a supporting structural scaffold. Here, we probe the interaction of the p66/p66' asymmetric reverse transcriptase precursor with HIV-1 protease by pulsed Q-band double electron-electron resonance EPR spectroscopy to measure distances between nitroxide labels introduced at surface-engineered cysteine residues. The data suggest that the flexible, exposed linker between the RNaseH and connection domains in the open state of the p66' subunit binds to the active site of protease in a configuration that is similar to that of extended peptide substrates.

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

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

Novel radial distribution function approach in the study of point mutations: the HIV-1 protease case study

Background: Mutations are one of the engines of evolution. Under constant stress pressure, mutations can lead to the emergence of unwanted, drug-resistant entities. Methodology: The radial distribution function weighted by the number of valence shell electrons is used to design quantitative structure-activity relationship (QSAR) model relating descriptors with the inhibition constant for a series of wild-type HIV-1 protease inhibitor complexes. The residuals of complexes with mutant HIV-1 protease were correlated with the energy of the highest occupied molecular orbitals of the residues introduced to enzyme via point mutations. Conclusion: Successful identification of residues Ile3, Asp25, Val32 and Ile50 as the one whose substitution influences the inhibition constant the most, demonstrates the potential of the proposed methodology for the study of the effects of point mutations.

The Influence of Hydrogen Atoms on the Performance of Radial Distribution Function-Based Descriptors in the Chemoinformatic Studies of HIV-1 Protease Complexes with Inhibitors

AIMS

The aim of this letter is to explore the influence of adding hydrogen atoms to the crystallographic structures of HIV-1 protease complexes with a series of inhibitors on the performance of radial distribution function based descriptors recently introduced in chemoinformatic studies.

BACKGROUND

Quite recently the successful application of molecular descriptors based on a radial distribution function to correlate it with biologically interesting properties of a ligand - enzyme complex was demonstrated. Except its predictive power, the analysis of atoms with dominant contributions to the RDFs can be used to identify relevant atoms and interactions. Since original paper was published on dataset consisting of the X-ray structures of complexes without hydrogen atoms, we wonder weather addition of light atoms can provide us new piece of information.

OBJECTIVE

The primarily objective is to create the model correlating the RDF based descriptors and physicochemical properties of the HIV-1 protease complexes with inhibitors with hydrogen atoms. Then, we will compare the performance of new model with previous one, where the hydrogen atoms were discarded. Information about interactions between the enzyme and the inhibitors will be extracted from the analysis of the RDF.

METHODS

The radial distribution function descriptor weighted by the number of valence shell electrons has proven to be sensitive to the changes in the structure of the enzyme and enzyme-ligand complexes. For each structure in our data set, RDF will be calculated and using multiple linear regression method the mathematical model will be designed correlating RDF based descriptors and the physicochemical properties. Statistical analysis of the atom's contribution to the total RDF will reveal relevant interactions.

RESULTS

The applicability of RDF based descriptor for the correlation with pKi and EC50 values is demonstrated, while simple models containing only two or three parameters are able to explain 78 and 86 % of the variance, respectively. The models with explicitly included hydrogens are of comparable quality with the previous models without hydrogens. The analysis of the atom's dominant contributions highlighted the importance of the hydroxyl groups of the inhibitor near the Asp25 and Asp25' residues when it is bounded to the protease.

CONCLUSION

Models based on the RDF weighted by the number of valence shell electrons for correlating small number of molecular descriptors and physicocehmical properties for structures with and without hydrogens are of comparable quality and both can be used for identification of relevant functional groups and interactions. Other: Our approach can be integrated to the next generation virtual screening methods, because is fast, reliable with high predictability potential.

Evidence of genomic information and structural restrictions of HIV-1 PR and RT gene regions from individuals experiencing antiretroviral virologic failure

OBJECTIVES

This study analyzed Protease-PR and Reverse Transcriptase-RT HIV-1 genomic information entropy metrics among patients under antiretroviral virologic failure, according to the numbers of virologic failures or resistance mutations.

METHODS

For this purpose, we used genomic sequences from PR and RT of HIV-1 from a cohort of chronic patients followed up at São Paulo Hospital.

RESULTS

Informational entropy proportionally increases with the number of antiretroviral virologic failures in PR and RT (p < .001). Affected regions of PR were related to catalytic and structural functions, such as Fulcrum (K20) Flap (M46) and Cantilever (A71). In RT, this occurred at Fingers (E44) and Palm (K219). Informational entropy increases according to the number of resistance mutations in PR and RT (p < .001). Higher PR entropy was proportional to the resistance mutation numbers in Fulcrum (L10), Active site (L24) Flap (M46), Cantilever (L63) and near Interface (L90). In RT, they related to regions responsible for protein stability such as Fingers (T39) and Palm (L100).

CONCLUSIONS

The antiretroviral selective pressure affects HIV genomic informational entropy at the PR and RT regions, leading to the emergence of more unstable virions. Mapping the three-dimensional structure in these HIV-1 proteins is relevant to designing new antiretroviral targeting resistant strains.

Target-Specific Prediction of Ligand Affinity with Structure-Based Interaction Fingerprints

Discovery and optimization of small molecule inhibitors as therapeutic drugs have immensely benefited from rational structure-based drug design. With recent advances in high-resolution structure determination, computational power, and machine learning methodology, it is becoming more tractable to elucidate the structural basis of drug potency. However, the applicability of machine learning models to drug design is limited by the interpretability of the resulting models in terms of feature importance. Here, we take advantage of the large number of available inhibitor-bound HIV-1 protease structures and associated potencies to evaluate inhibitor diversity and machine learning models to predict ligand affinity. First, using a hierarchical clustering approach, we grouped HIV-1 protease inhibitors and identified distinct core structures. Explicit features including protein-ligand interactions were extracted from high-resolution cocrystal structures as 3D-based fingerprints. We found that a gradient boosting machine learning model with this explicit feature attribution can predict binding affinity with high accuracy. Finally, Shapley values were derived to explain local feature importance. We found specific van der Waals (vdW) interactions of key protein residues are pivotal for the predicted potency. Protein-specific and interpretable prediction models can guide the optimization of many small molecule drugs for improved potency.

Halogen Bond Interactions of Novel HIV-1 Protease Inhibitors (PI) (GRL-001-15 and GRL-003-15) with the Flap of Protease Are Critical for Their Potent Activity against Wild-Type HIV-1 and Multi-PI-Resistant Variants

We generated two novel nonpeptidic HIV-1 protease inhibitors (PIs), GRL-001-15 and GRL-003-15, which contain unique crown-like tetrahydropyranofuran (Crn-THF) and P2'-cyclopropyl-aminobenzothiazole (Cp-Abt) moieties as P2 and P2' ligands, respectively. GRL-001-15 and GRL-003-15 have meta-monofluorophenyl and para-monofluorophenyl at the P1 site, respectively, exert highly potent activity against wild-type HIV-1 with 50% effective concentrations (EC50s) of 57 and 50 pM, respectively, and have favorable cytotoxicity profiles with 50% cytotoxic concentrations (CC50s) of 38 and 11 μM, respectively. The activity of GRL-001-15 against multi-PI-resistant HIV-1 variants was generally greater than that of GRL-003-15. The EC50 of GRL-001-15 against an HIV-1 variant that was highly resistant to multiple PIs, including darunavir (DRV) (HIV-1DRVRP30), was 0.17 nM, and that of GRL-003-15 was 3.3 nM, while DRV was much less active, with an EC50 of 216 nM. The emergence of HIV-1 variants resistant to GRL-001-15 and GRL-003-15 was significantly delayed compared to that of variants resistant to selected PIs, including DRV. Structural analyses of wild-type protease (PRWT) complexed with the novel PIs revealed that GRL-001-15's meta-fluorine atom forms halogen bond interactions (2.9 and 3.0 Å) with Gly49 and Ile50, respectively, of the protease flap region and with Pro81' (2.7 and 3.2 Å), which is located close to the protease active site, and that two fluorine atoms of GRL-142-13 form multiple halogen bond interactions with Gly49, Ile50, Pro81', Ile82', and Arg8'. In contrast, GRL-003-15 forms halogen bond interactions with Pro81' alone, suggesting that the reduced antiviral activity of GRL-003-15 is due to the loss of the interactions with the flap region.

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.

GRL-079, a Novel HIV-1 Protease Inhibitor, Is Extremely Potent against Multidrug-Resistant HIV-1 Variants and Has a High Genetic Barrier against the Emergence of Resistant Variants

We identified four novel nonpeptidic human immunodeficiency virus type 1 (HIV-1) protease inhibitors (PIs), GRL-078, -079, -077, and -058, containing an alkylamine at the C-5 position of P2 tetrahydropyrano-tetrahydrofuran (Tp-THF) and a P2' cyclopropyl (Cp) (or isopropyl)-aminobenzothiazole (Abt) moiety. Their 50% effective concentrations (EC₅₀s) were 2.5 to 30 nM against wild-type HIV-1_NL4-3, 0.3 to 6.7 nM against HIV-2_EHO, and 0.9 to 90 nM against laboratory-selected PI-resistant HIV-1 and clinical HIV-1 variants resistant to multiple FDA-approved PIs (HIV_MDR). GRL-078, -079, -077, and -058 also effectively blocked the replication of HIV-1 variants highly resistant to darunavir (DRV) (HIV_DRV^r_p51), with EC₅₀s of 38, 62, 61, and 90 nM, respectively, while four FDA-approved PIs examined (amprenavir, atazanavir, lopinavir [LPV], and DRV) had virtually no activity (EC₅₀s of >1,000 nM) against HIV_DRV^r_p51 Structurally, GRL-078, -079, and -058 form strong hydrogen bond interactions between Tp-THF modified at C-5 and Asp29/Asp30/Gly48 of wild-type protease, while the P2' Cp-Abt group forms strong hydrogen bonds with Asp30'. The Tp-THF and Cp-Abt moieties also have good nonpolar interactions with protease residues located in the flap region. For selection with LPV and DRV by use of a mixture of 11 HIV_MDR strains (HIV_11MIX), HIV_11MIX became highly resistant to LPV and DRV over 13 to 32 and 32 to 41 weeks, respectively. However, for selection with GRL-079 and GRL-058, HIV_11MIX failed to replicate at >0.08 μM and >0.2 μM, respectively. Thermal stability results supported the highly favorable anti-HIV-1 potency of GRL-079 as well as other PIs. The present data strongly suggest that the P2 Tp-THF group modified at C-5 and the P2' Abt group contribute to the potent anti-HIV-1 profiles of the four PIs against HIV-1_NL4-3 and a wide spectrum of HIV_MDR strains.

Transient HIV-1 Gag-protease interactions revealed by paramagnetic NMR suggest origins of compensatory drug resistance mutations

Cleavage of the group-specific antigen (Gag) polyprotein by HIV-1 protease represents the critical first step in the conversion of immature noninfectious viral particles to mature infectious virions. Selective pressure exerted by HIV-1 protease inhibitors, a mainstay of current anti-HIV-1 therapies, results in the accumulation of drug resistance mutations in both protease and Gag. Surprisingly, a large number of these mutations (known as secondary or compensatory mutations) occur outside the active site of protease or the cleavage sites of Gag (located within intrinsically disordered linkers connecting the globular domains of Gag to one another), suggesting that transient encounter complexes involving the globular domains of Gag may play a role in guiding and facilitating access of the protease to the Gag cleavage sites. Here, using large fragments of Gag, as well as catalytically inactive and active variants of protease, we probe the nature of such rare encounter complexes using intermolecular paramagnetic relaxation enhancement, a highly sensitive technique for detecting sparsely populated states. We show that Gag-protease encounter complexes are primarily mediated by interactions between protease and the globular domains of Gag and that the sites of transient interactions are correlated with surface exposed regions that exhibit a high propensity to mutate in the presence of HIV-1 protease inhibitors.

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

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

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

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

ARTSY-J: Convenient and precise measurement of (3)JHNHα couplings in medium-size proteins from TROSY-HSQC spectra

A new and convenient method, named ARTSY-J, is introduced that permits extraction of the (3)JHNHα couplings in proteins from the relative intensities in a pair of (15)N-(1)H TROSY-HSQC spectra. The pulse scheme includes (3)JHNHα dephasing of the narrower TROSY (1)H(N)-{(15)N} doublet component during a delay, integrated into the regular two-dimensional TROSY-HSQC pulse scheme, and compares the obtained intensity with a reference spectrum where (3)JHNHα dephasing is suppressed. The effect of passive (1)H(α) spin flips downscales the apparent (3)JHNHα coupling by a uniform factor that depends approximately linearly on both the duration of the (3)JHNHα dephasing delay and the (1)H-(1)H cross relaxation rate. Using such a correction factor, which accounts for the effects of both inhomogeneity of the radiofrequency field and (1)H(α) spin flips, agreement between prior and newly measured values for the small model protein GB3 is better than 0.3Hz. Measurement for the HIV-1 protease homodimer (22kDa) yields (3)JHNHα values that agree to better than 0.7Hz with predictions made on the basis of a previously parameterized Karplus equation. Although for Gly residues the two individual (3)JHNHα couplings cannot be extracted from a single set of ARTSY-J spectra, the measurement provides valuable ϕ angle information.

Improved darunavir genotypic mutation score predicting treatment response for patients infected with HIV-1 subtype B and non-subtype B receiving a salvage regimen

OBJECTIVES

The objective of this study was to improve the prediction of the impact of HIV-1 protease mutations in different viral subtypes on virological response to darunavir.

METHODS

Darunavir-containing treatment change episodes (TCEs) in patients previously failing PIs were selected from large European databases. HIV-1 subtype B-infected patients were used as the derivation dataset and HIV-1 non-B-infected patients were used as the validation dataset. The adjusted association of each mutation with week 8 HIV RNA change from baseline was analysed by linear regression. A prediction model was derived based on best subset least squares estimation with mutational weights corresponding to regression coefficients. Virological outcome prediction accuracy was compared with that from existing genotypic resistance interpretation systems (GISs) (ANRS 2013, Rega 9.1.0 and HIVdb 7.0).

RESULTS

TCEs were selected from 681 subtype B-infected and 199 non-B-infected adults. Accompanying drugs were NRTIs in 87%, NNRTIs in 27% and raltegravir or maraviroc or enfuvirtide in 53%. The prediction model included weighted protease mutations, HIV RNA, CD4 and activity of accompanying drugs. The model's association with week 8 HIV RNA change in the subtype B (derivation) set was R(2) = 0.47 [average squared error (ASE) = 0.67, P < 10(-6)]; in the non-B (validation) set, ASE was 0.91. Accuracy investigated by means of area under the receiver operating characteristic curves with a binary response (above the threshold value of HIV RNA reduction) showed that our final model outperformed models with existing interpretation systems in both training and validation sets.

CONCLUSIONS

A model with a new darunavir-weighted mutation score outperformed existing GISs in both B and non-B subtypes in predicting virological response to darunavir.

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

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

Conformational variation of an extreme drug resistant mutant of HIV protease

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

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

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

Mapping of functional domains and characterization of the transcription factor Cph1 that mediate morphogenesis in Candida albicans

Cph1, a transcription factor of the Mitogen Activated Protein (MAP) kinase pathway, regulates morphogenesis in human fungal pathogen Candida albicans. Here, by following a systemic deletion approach, we have identified functional domains and motifs of Cph1 that are involved in transcription factor activity and cellular morphogenesis. We found that the N-terminal homeodomain is essential for the DNA binding activity; however, C-terminal domain and polyglutamine motif (PQ) are indispensable for the transcriptional activation function. Complementation analysis of the cph1Δ null mutant using various deletion derivatives revealed functional significance of the N- and C-terminal domains and PQ motif in filamentation process, chlamydospore formation and sensitivity to the cell wall interfering compounds. Genome-wide identification of the Cph1 binding site and quantitative RT-PCR transcript analysis in cph1Δ null mutant revealed that a number of genes which are associated with the filamentous growth, maintaining cell wall organization and mitochondrial function, and the genes of the pH response pathway are the transcriptional targets of Cph1. The data also suggest that Cph1 may function as a positive or negative regulator depending on the morphological state and physiological conditions. Moreover, differential expression of the upstream MAP kinase pathway genes in wild type and cph1Δ null mutant indicated the existence of a feedback regulation.

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

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

Structural basis and distal effects of Gag substrate coevolution in drug resistance to HIV-1 protease

Drug resistance mutations in response to HIV-1 protease inhibitors are selected not only in the drug target but elsewhere in the viral genome, especially at the protease cleavage sites in the precursor protein Gag. To understand the molecular basis of this protease-substrate coevolution, we solved the crystal structures of drug resistant I50V/A71V HIV-1 protease with p1-p6 substrates bearing coevolved mutations. Analyses of the protease-substrate interactions reveal that compensatory coevolved mutations in the substrate do not restore interactions lost due to protease mutations, but instead establish other interactions that are not restricted to the site of mutation. Mutation of a substrate residue has distal effects on other residues' interactions as well, including through the induction of a conformational change in the protease. Additionally, molecular dynamics simulations suggest that restoration of active site dynamics is an additional constraint in the selection of coevolved mutations. Hence, protease-substrate coevolution permits mutational, structural, and dynamic changes via molecular mechanisms that involve distal effects contributing to drug resistance.

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

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

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

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

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

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.

Understanding HIV-1 protease autoprocessing for novel therapeutic development

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

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

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

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

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

Elucidating a relationship between conformational sampling and drug resistance in HIV-1 protease

Enzyme targets in rapidly replicating systems, such as retroviruses, commonly respond to drug-selective pressure with mutations arising in the active site pocket that limit inhibitor effectiveness by introducing steric hindrance or by eliminating essential molecular interactions. However, these primary mutations are disposed to compromising pathogenic fitness. Emerging secondary mutations, which are often found outside of the binding cavity, may or can restore fitness while maintaining drug resistance. The accumulated drug pressure selected mutations could have an indirect effect in the development of resistance, such as altering protein flexibility or the dynamics of protein-ligand interactions. Here, we show that accumulation of mutations in a drug-resistant HIV-1 protease (HIV-1 PR) variant, D30N/M36I/A71V, changes the fractional occupancy of the equilibrium conformational sampling ensemble. Correlations are made among populations of the conformational states, namely, closed-like, semiopen, and open-like, with inhibition constants, as well as kinetic parameters. Mutations that stabilize a closed-like conformation correlate with enzymes of lowered activity and with higher affinity for inhibitors, which is corroborated by a further increase in the fractional occupancy of the closed state upon addition of inhibitor or substrate-mimic. Cross-resistance is found to correlate with combinations of mutations that increase the population of the open-like conformations at the expense of the closed-like state while retaining native-like occupancy of the semiopen population. These correlations suggest that at least three states are required in the conformational sampling model to establish the emergence of drug resistance in HIV-1 PR. More importantly, these results shed light on a possible mechanism whereby mutations combine to impart drug resistance while maintaining catalytic activity.

Inhibitor-induced conformational shifts and ligand-exchange dynamics for HIV-1 protease measured by pulsed EPR and NMR spectroscopy

Double electron-electron resonance (DEER) spectroscopy was utilized to investigate shifts in conformational sampling induced by nine FDA-approved protease inhibitors (PIs) and a nonhydrolyzable substrate mimic for human immunodeficiency virus type 1 protease (HIV-1 PR) subtype B, subtype C, and CRF_01 A/E. The ligand-bound subtype C protease has broader DEER distance profiles, but trends for inhibitor-induced conformational shifts are comparable to those previously reported for subtype B. Ritonavir, one of the strong-binding inhibitors for subtypes B and C, induces less of the closed conformation in CRF_01 A/E. (1)H-(15)N heteronuclear single-quantum coherence (HSQC) spectra were acquired for each protease construct titrated with the same set of inhibitors. NMR (1)H-(15)N HSQC titration data show that inhibitor residence time in the protein binding pocket, inferred from resonance exchange broadening, shifting or splitting correlates with the degree of ligand-induced flap closure measured by DEER spectroscopy. These parallel results show that the ligand-induced conformational shifts resulting from protein-ligand interactions characterized by DEER spectroscopy of HIV-1 PR obtained at the cryogenic temperature are consistent with more physiological solution protein-ligand interactions observed by solution NMR spectroscopy.

Thermodynamics of strongly allosteric inhibition: a model study of HIV-1 protease

Protein inhibitors that shift the thermodynamic equilibrium towards a denatured state escape, in general, the straightforward framework of competitive or allosteric inhibitors. The equilibrium properties of peptides which compete with the folding, or more precisely destabilize the native state, of the human immunodeficiency virus (HIV)-1 protease monomer are studied within a structure-based model. The effect of peptides that disrupt the hydrophobic core of the protein can still be summarized in terms of an inhibition constant, which depends on the thermal stability of the protein. The state of the protein denatured by such a peptide is more structured than its intrinsic denatured state, but displays the same degree of compactness. Peptides that target less buried regions of the protein are less efficient and display a more complex thermodynamics that cannot be captured in a simple way.

Correlating conformational shift induction with altered inhibitor potency in a multidrug resistant HIV-1 protease variant

Inhibitor-induced conformational ensemble shifts in a multidrug resistant HIV-1 protease variant, MDR769, are characterized by site-directed spin labeling double electron-electron resonance spectroscopy. For MDR769 compared to the native enzyme, changes in inhibitor IC(50) values are related to a parameter defined as |ΔC|, which is the relative change in the inhibitor-induced shift to the closed state. Specifically, a linear correlation is found between |ΔC| and the magnitude of the change in IC(50), provided that inhibitor binding is not too weak. Moreover, inhibitors that exhibit MDR769 resistance no longer induce a strong shift to a closed conformational ensemble as seen previously in the native enzyme.

Mechanism of dissociative inhibition of HIV protease and its autoprocessing from a precursor

Dimerization is indispensible for release of the human immunodeficiency virus protease (PR) from its precursor (Gag-Pol) and ensuing mature-like catalytic activity that is crucial for virus maturation. We show that a single-chain Fv fragment (scFv) of a previously reported monoclonal antibody (mAb1696), which recognizes the N-terminus of PR, dissociates a dimeric mature D25N PR mutant with an enhanced dimer dissociation constant (K(d)) in the sub-micromolar range to form predominantly a monomer-scFv complex at a 1:1 ratio, along with small (5-10%) amounts of a dimer-scFv complex. Enzyme kinetics indicate a mixed mechanism of inhibition of the wild-type PR, which exhibits a K(d)<10nM, with effects both on K(m) and k(cat) at an scFv-to-PR ratio of 10:1. ScFv binds to the N-terminal peptide P(1)QITLW(6) of PR and to PR monomers with dissociation constants of ≤30 nM and ~100 nM, respectively. Consistent with an ~400-fold increase in the dissociation of the antibody (K(Ab)) on even addition of an acetyl group to P(1) of the peptide, the antibody fails to inhibit N-terminal autoprocessing of the PR from a model precursor (at ~5 μM). However, subsequent to this cleavage, it sequesters the PR, thus blocking autoprocessing at its C-terminus. A second monoclonal antibody [PRM1 (human monoclonal antibody to PR)], which recognizes part of the flap region (residues 41-47) of the mature PR and its precursor, does not inhibit autoprocessing and ensuing catalytic activity. However, its failure to recognize drug-resistant clinical mutants of PR may be beneficial to monitor the selection of mutations in this region under drug pressure.

Single point mutation induced alterations in the equilibrium structural transitions on the folding landscape of HIV-1 protease

Equilibrium folding-unfolding transitions are hard to study in HIV-1 protease (PR) because of its autolytic properties. Further, the protease exhibits many tolerant point mutations some of which also impart drug resistance to the protein. It is conceivable that the mutations affect protein's function by altering its folding characteristics; these would clearly depend on the nature of the mutations themselves. In this background, we report here NMR studies on the effects of D25 N mutation, which removes one negative charge from the protein at the active site, on the equilibrium folding behaviour of PR starting from its acetic acid denatured state. It is observed that in PRD25N two slowly exchanging conformations are present at the N-terminal. One of them is similar to that of PR. Though the conformational and dynamics preferences of PR and PRD25N are fairly similar in 9 M acetic acid, they seem to undergo different folding transitions when acetic acid concentration is reduced. The differences are seen in the active site, in the flap, and in the hinge of the flap regions. The present study suggests that such differences, though different in detail, would occur for other mutations as well, and also for different initial denatured states. These would have significant regulatory implications for the efficacy of protease function.

How Does Darunavir Prevent HIV-1 Protease Dimerization?

The drug Darunavir (DRV) is a potent inhibitor of HIV-1 protease (PR), a homodimeric essential enzyme of the AIDS virus. Recent experimental data suggest that DRV is able to prevent dimerization of HIV-1 PR, which, together with its high affinity for the mature enzyme, has been linked to the high genetic barrier to the development of viral resistance. The mechanism of dimerization inhibition and the binding mode(s) of DRV to monomeric HIV-1 PR are unknown. Here, multiple molecular dynamics simulations with explicit solvent (for a total of 11 μs with the CHARMM force field and 1 μs with the Amber force field) show that the monomer of HIV-1 PR is structurally stable and reveal a major binding mode of DRV. This binding mode is stabilized by favorable interactions between the apolar groups of DRV and the hydrophobic residues Ile32, Ile47, Ile50, Ile54, Pro79, Val82, and Ile84. The binding mode to monomeric HIV-1 PR identified by molecular dynamics is different from the two binding modes observed in the crystal structure of the complex with dimeric HIV-1 PR. As an example, there are no interactions between DRV and the catalytic Asp25 in the binding mode to monomeric HIV-1 PR revelead by the simulations. In contrast, the simulations show extensive and stable interactions between DRV and the flap (residues 46-55), which are likely to sterically hinder the formation of the flap interface as observed in the dimeric structure. Which of the two mechanisms of inhibition (dimerization inhibition by association with the flap or binding to the active site of the mature enzyme) dominates might depend on the HIV-1 PR mutations, and it is likely that dimerization inhibition is predominant for multiple mutations at the active site in multidrug resistant strains.

HIV-1 protease with 20 mutations exhibits extreme resistance to clinical inhibitors through coordinated structural rearrangements

The escape mutant of HIV-1 protease (PR) containing 20 mutations (PR20) undergoes efficient polyprotein processing even in the presence of clinical protease inhibitors (PIs). PR20 shows >3 orders of magnitude decreased affinity for PIs darunavir (DRV) and saquinavir (SQV) relative to PR. Crystal structures of PR20 crystallized with yttrium, substrate analogue p2-NC, DRV, and SQV reveal three distinct conformations of the flexible flaps and diminished interactions with inhibitors through the combination of multiple mutations. PR20 with yttrium at the active site exhibits widely separated flaps lacking the usual intersubunit contacts seen in other inhibitor-free dimers. Mutations of residues 35-37 in the hinge loop eliminate interactions and perturb the flap conformation. Crystals of PR20/p2-NC contain one uninhibited dimer with one very open flap and one closed flap and a second inhibitor-bound dimer in the closed form showing six fewer hydrogen bonds with the substrate analogue relative to wild-type PR. PR20 complexes with PIs exhibit expanded S2/S2' pockets and fewer PI interactions arising from coordinated effects of mutations throughout the structure, in agreement with the strikingly reduced affinity. In particular, insertion of the large aromatic side chains of L10F and L33F alters intersubunit interactions and widens the PI binding site through a network of hydrophobic contacts. The two very open conformations of PR20 as well as the expanded binding site of the inhibitor-bound closed form suggest possible approaches for modifying inhibitors to target extreme drug-resistant HIV.

Terminal interface conformations modulate dimer stability prior to amino terminal autoprocessing of HIV-1 protease

The HIV-1 protease (PR) mediates its own release (autoprocessing) from the polyprotein precursor, Gag-Pol, flanked by the transframe region (TFR) and reverse transcriptase at its N- and C-termini, respectively. Autoprocessing at the N-terminus of PR mediates stable dimer formation essential for catalytic activity, leading to the formation of infectious virus. An antiparallel β-sheet interface formed by the four N- and C-terminal residues of each subunit is important for dimer stability. Here, we present the first high-resolution crystal structures of model protease precursor-clinical inhibitor (PI darunavir or saquinavir) complexes, revealing varying conformations of the N-terminal flanking (S(-4)FNF(-1)) and interface residues (P(1)QIT(4)). A 180° rotation of the T(4)-L(5) peptide bond is accompanied by a new Q(2)-L(5) hydrogen bond and complete disengagement of PQIT from the β-sheet dimer interface, which may be a feature for intramolecular autoprocessing. This result is consistent with drastically lower thermal stability by 14-20 °C of PI complexes of precursors and the mature PR lacking its PQIT residues (by 18.3 °C). Similar to the TFR-PR precursor, this deletion also results in a darunavir dissociation constant (2 × 10(4))-fold higher and a markedly increased dimer dissociation constant relative to the mature PR. The terminal β-sheet perturbations of the dimeric structure likely account for the drastically poorer inhibition of autoprocessing of TFR-PR relative to the mature PR, even though significant differences in active site-PI interactions in these structures were not observed. The novel conformations of the dimer interface may be exploited to target selectively the protease precursor prior to its N-terminal cleavage.