The structural stability of the HIV-1 protease

Premature Activation of the HIV-1 Protease Is Influenced by Polymorphisms in the Hinge Region

HIV-1 protease inhibitors are an essential component of antiretroviral therapy. However, drug resistance is a pervasive issue motivating a persistent search for novel therapies. Recent reports found that when protease activates within the host cell's cytosol, it facilitates the pyroptotic killing of infected cells. This has led to speculation that promoting protease activation, rather than inhibiting it, could help to eradicate infected cells and potentially cure HIV-1 infection. Here, we used a nanoscale flow cytometry-based assay to characterize protease resistance mutations and polymorphisms. We quantified protease activity, viral concentration, and premature protease activation and confirmed previous findings that major resistance mutations generally destabilize the protease structure. Intriguingly, we found evidence that common polymorphisms in the hinge domain of protease can influence its susceptibility to premature activation. This suggests that viral heterogeneity could pose a considerable challenge for therapeutic strategies aimed at inducing premature protease activation in the future.

Characterization of pH-induced conformational changes in recombinant DENV NS2B-NS3pro

The increasing frequency of Dengue is a cause of severe epidemics and therefore demands strategies for effective prevention, diagnosis, and treatment. DENV-protease is being investigated as a potential therapeutic target. However, due to the flat and highly charged active site of the DENV-protease, designing orthosteric medicines is very difficult. In this study, we have done a thorough analysis of pH-dependent conformational changes in recombinantly expressed DENV protease using various spectroscopic techniques. Our spectroscopic study of DENV protease (NS2B-NS3pro) at different pH conditions gives important insights into the dynamicity of structural conformation. At physiological pH, the DENV-protease exists in a random-coiled state. Lowering the pH promotes the formation of alpha-helical and beta-sheet structures i.e. gain of secondary structure as shown by Far-UV CD. The light scattering and Thioflavin T (ThT)-binding assay proved the aggregation-prone tendency of DENV-protease at pH 4.0. Further, the confocal microscopy image intensity showed the amorphous aggregate formation of DENV protease at pH 4.0. Thus, the DENV protease acquires different conformations with changes in pH conditions. Together, these results have the potential to facilitate the design of a conformation destabilizer-based therapeutic strategy for dengue fever.

Chemoselective Caging of Carboxyl Groups for On-Demand Protein Activation with Small Molecules

Tools for on-demand protein activation enable impactful gain-of-function studies in biological settings. Thus far, however, proteins have been chemically caged at primarily Lys, Tyr, and Sec, typically through the genetic encoding of unnatural amino acids. Herein, we report that the preferential reactivity of diazo compounds with protonated acids can be used to expand this toolbox to solvent-accessible carboxyl groups with an elevated pKa value. As a model protein, we employed lysozyme (Lyz), which has an active-site Glu35 residue with a pKa value of 6.2. A diazo compound with a bioorthogonal self-immolative handle esterified Glu35 selectively, inactivating Lyz. The hydrolytic activity of the caged Lyz on bacterial cell walls was restored with two small-molecule triggers. The decaging was more efficient by small molecules than by esterases. This simple chemical strategy was also applied to a hemeprotein and an aspartyl protease, setting the stage for broad applicability.

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.

Modulation of the monomer-dimer equilibrium and catalytic activity of SARS-CoV-2 main protease by a transition-state analog inhibitor

The role of dimer formation for the onset of catalytic activity of SARS-CoV-2 main protease (MPro^WT) was assessed using a predominantly monomeric mutant (MPro^M). Rates of MPro^WT and MPro^M catalyzed hydrolyses display substrate saturation kinetics and second-order dependency on the protein concentration. The addition of the prodrug GC376, an inhibitor of MPro^WT, to MPro^M leads to an increase in the dimer population and catalytic activity with increasing inhibitor concentration. The activity reaches a maximum corresponding to a dimer population in which one active site is occupied by the inhibitor and the other is available for catalytic activity. This phase is followed by a decrease in catalytic activity due to the inhibitor competing with the substrate. Detailed kinetics and equilibrium analyses are presented and a modified Michaelis-Menten equation accounts for the results. These observations provide conclusive evidence that dimer formation is coupled to catalytic activity represented by two equivalent active sites.

Improved Production of HIV-1 Subtype C Protease from Transgenic E. Coli

Background:

Human Immunodeficiency Virus 1 (HIV-1) subtype C is responsible for the majority of infections of patients in Southern Africa. The HIV protease is a primary target for the development of highly efficient anti-retroviral pharmaceuticals because of its pivotal role in the maturation of the virus in the host cell. For target validation of novel HIV protease inhibitors, there is a need for the availability of an abundance of this protease.

Objective:

This study reports an optimized method to produce HIV-1 protease derived from HIV-1 subtype C.

Methods:

It involves the use of a transgenic E. coli strain that overexpresses the native form of the enzyme via inclusion bodies. A stringent method for the isolation, purification, and renaturation resulted in the production of highly pure active HIV-1 protease. In order to facilitate an increase in protease yields, an optimized growth strategy was developed. In this regard, a chemically defined medium with lower glucose content and devoid of essential amino acids of the TCA cycle was used as an alternative to the widely used nutrient-rich Luria Bertani (LB) medium.

Results:

Results indicated an increase in protease yield up to twice the amount, thereby making this medium an attractive alternative for increasing biomass and HIV protease production for future research.

Conclusion:

An optimized method for HIV-1 protease derived from HIV-1 subtype C production using chemically defined media was established. This was achieved using a known method to isolate and purify the enzyme with the use of a specialized feeding strategy.

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.

Prediction of HIV drug resistance based on the 3D protein structure: Proposal of molecular field mapping

A method for predicting HIV drug resistance by using genotypes would greatly assist in selecting appropriate combinations of antiviral drugs. Models reported previously have had two major problems: lack of information on the 3D protein structure and processing of incomplete sequencing data in the modeling procedure. We propose obtaining the 3D structural information of viral proteins by using homology modeling and molecular field mapping, instead of just their primary amino acid sequences. The molecular field potential parameters reflect the physicochemical characteristics associated with the 3D structure of the proteins. We also introduce the Bayesian conditional mutual information theory to estimate the probabilities of occurrence of all possible protein candidates from an incomplete sequencing sample. This approach allows for the effective use of uncertain information for the modeling process. We applied these data analysis techniques to the HIV-1 protease inhibitor dataset and developed drug resistance prediction models with reasonable performance.

The denatured state of HIV-1 protease under native conditions

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high‐resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV‐1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant‐free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV‐1 protease under native conditions displays rich patterns of transient native and non‐native structures, which could be of relevance to its guidance through a complex folding process.

Experiments and simulation on ZIKV NS2B-NS3 protease reveal its complex folding

Zika virus has been identified in various body fluids such as semen, urine, saliva, cerebrospinal fluid, and vaginal secretion of an infected individual. The pH of these fluids varies from mildly acidic to mildly alkaline. So it is imperative to understand the impact of these conditions on viral protein functioning. We investigated the NS2B-NS3 protease stability and its activity in different denaturing environments. Finding indicates that NS2B-NS3 protease maintains stability at pH 4.8-8.7. Thus it suggests that the complex remains functionally active to hydrolyze the polyprotein within a diverse environmental condition such as variable pH. Despite a stable structure at a broad pH range, a change in environmental conditions dramatically influence its protease activity. Moreover, it is susceptible to structural transformation leading to increased β-strand or helix content in the presence of alcohol. This study may help further to understand the folding-function relationship of the general flaviviral protease complex.

A Single Amino Acid Substitution at the HIV-1 Protease Termini Dimer Interface Significantly Reduces Viral Particles Processing Efficiency

The dimeric form of HIV-1 protease (PR) is required for its full proteolytic activity. The stability of the dimer primarily depends on the termini interface, with N-terminal residues 1-4 of one monomer encountering C-terminal residues 96-99 of another. We made an alanine substitution for valine 3 (V3) or leucine 97 (L97) at the termini dimer interface and tested their proteolytic activity. We found that an alanine substitution for L97 (PRL97A) completely inhibited the proteolytic activity of the PR. However, an alanine substitution for V3 (PRV3A) partially impaired the proteolytic activity. We then introduced two forced-dimerization systems involving nucleocapsid (NC) replacement or the addition of 1-2 leucine zippers to determine whether the proteolytic activity of dimer-defective PRs could be restored. We found that two forced-dimerization systems compensated for the defect in PRV3A, but not in PRL97A. This implies that PRV3A and PRL97A potentially impair the PR via different mechanisms or cause defects in PR activity to different extents. These novel findings will likely serve as a foundation for developing new PR inhibitors for treating drug-resistant HIV-1 infections in the future.

A synergy of activity, stability, and inhibitor-interaction of HIV-1 protease mutants evolved under drug-pressure

A clinically-relevant, drug-resistant mutant of HIV-1 protease (PR), termed Flap+_(I54V) and containing L10I, G48V, I54V and V82A mutations, is known to produce significant changes in the entropy and enthalpy balance of drug-PR interactions, compared to wild-type PR. A similar mutant, Flap+_(I54A) , which evolves from Flap+_(I54V) and contains the single change at residue 54 relative to Flap+_(I54V) , does not. Yet, how Flap+_(I54A) behaves in solution is not known. To understand the molecular basis of V54A evolution, we compared nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy, isothermal titration calorimetry, and enzymatic assay data from four PR proteins: PR (pWT), Flap+_(I54V) , Flap+_(I54A) , and Flap+_(I54) , a control mutant that contains only L10I, G48V and V82A mutations. Our data consistently show that selection to the smaller side chain at residue 54, not only decreases inhibitor affinity, but also restores the catalytic activity.

Drug susceptibility and replication capacity of a rare HIV-1 subtype C protease hinge region variant

BACKGROUND

Protease inhibitors form the main component of second-line antiretroviral treatment in South Africa. Despite their efficacy, mutations arising within the HIV-1 gag and protease coding regions contribute to the development of resistance against this class of drug. In this paper we investigate a South African HIV-1 subtype C Gag-protease that contains a hinge region mutation and insertion (N37T↑V).

METHODS

In vitro single-cycle drug susceptibility and viral replication capacity assays were performed on W1201i, a wild-type reference isolate (MJ4) and a chimeric construct (MJ4GagN37T↑VPR). Additionally, enzyme assays were performed on the N37T↑V protease and a wild-type reference protease.

RESULTS

W1201i showed a small (threefold), but significant (P<0.0001) reduction in drug susceptibility to darunavir compared with MJ4. Substitution of W1201i-Gag with MJ4-Gag resulted in an additional small (twofold), but significant (P<0.01) reduction in susceptibility to lopinavir and atazanavir. The W1201i pseudovirus had a significantly (P<0.01) reduced replication capacity (16.4%) compared with the MJ4. However, this was dramatically increased to 164% (P<0.05) when W1201i-Gag was substituted with MJ4-Gag. Furthermore, the N37T↑V protease displayed reduced catalytic processing compared with the SK154 protease.

CONCLUSIONS

Collectively, these data suggest that the N37T↑V mutation and insertion increases viral infectivity and decreases drug susceptibility. These variations are classified as secondary mutations, and indirectly impact inhibitor binding, enzyme fitness and enzyme stability. Additionally, polymorphisms arising in Gag can modify the impact of protease with regards to viral replication and susceptibility to 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.

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

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

NMR and MD studies combined to elucidate inhibitor and water interactions of HIV-1 protease and their modulations with resistance mutations

Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.

Exploring the flap dynamics of the South African HIV subtype C protease in presence of FDA-approved inhibitors: MD study

HIV-1 protease (HIV PR) is considered as one of the most attractive targets for the treatment of HIV and the impact of flap dynamics of HIV PR on the binding affinities of protease inhibitors (PIs) is a crucial ongoing research field. Recently, our research group evaluated the binding affinities of different FDA approved PIs against the South African HIV-1 subtype C (C-SA) protease (PR). The CSA-HIV PR displayed weaker binding affinity for most of the clinical PIs compared to HIV-1 B subtype for West and Central Europe, the Americas. In the current work, the flap dynamics of four different systems of HIV-1 C-SA PR complexed to FDA approved second generation PIs and its impact on binding was explored over the molecular dynamic trajectories. It was observed that the interactions of the selected drugs with the binding site residues of the protease may not be the major contributor for affinity towards PIs. Various post-MD analyses were performed, also entropic contributions, solvation free energies and hydrophobic core formation interactions were studied to assess how the flap dynamics of C-SA PR which is affected by such factors. From these contributions, large van der Waals interactions and low solvation free energies were found to be major factors for the higher activity of ATV against C-SA HIV PR. Furthermore, a comparatively stable hydrophobic core may be responsible for higher stability of the PR flaps of the ATV complex. The outcome of this study provides significant guidance to how the flap dynamics of C-SA PR is affected by various factors as a result of the binding affinity of various protease inhibitors. It will also assist with the design of potent inhibitors against C-SA HIV PR that apart from binding in the active site of PR can interacts with the flaps to prevent opening of the flaps resulting in inactivation of the protease.

Using the Relative Energy Gradient Method with Interacting Quantum Atoms to Determine the Reaction Mechanism and Catalytic Effects in the Peptide Hydrolysis in HIV-1 Protease

The reaction mechanism in an active site is of the utmost importance when trying to understand the role that an enzyme plays in biological processes. In a recently published paper [Theor. Chem. Acc. 2017, 136, 86], we formalised the Relative Energy Gradient (REG) method for automating an Interacting Quantum Atoms (IQA) analysis. Here, the REG method is utilised to determine the mechanism of peptide hydrolysis in the aspartic active site of the enzyme HIV-1 Protease. Using the REG method along with the IQA approach we determine the mechanism of peptide hydrolysis without employing any arbitrary parameters and with remarkable ease (albeit at large computational cost: the system contains 133 atoms, which means that there are 17 689 individual IQA terms to be calculated). When REG and IQA work together it is possible to determine a reaction mechanism at atomistic resolution from data directly derived from quantum calculations, without arbitrary parameters. Moreover, the mechanism determined by this novel method gives concrete insight into how the active site residues catalyse peptide hydrolysis.

Evolution of inhibitor-resistant natural mutant forms of HIV-1 protease probed by pre-steady state kinetic analysis

Pre-steady state kinetic analysis of mechanistic features of substrate binding and processing is crucial for insight into the evolution of inhibitor-resistant forms of HIV-1 protease. These data may provide a correct vector for rational drug design assuming possible intrinsic dynamic effects. These data should also give some clues to the molecular mechanism of protease action and resistance to inhibitors. Here we report pre-steady state kinetics of the interaction of wild type or mutant forms of HIV-1 protease with a FRET-labeled peptide. The three-stage "minimal" kinetic scheme with first and second reversible steps of substrate binding and with following irreversible peptide cleavage step adequately described experimental data. For the first time, a set of "elementary" kinetic parameters of wild type HIV-1 protease and its natural mutant inhibitor-resistant forms MDR-HM, ANAM-11 and prDRV4 were compared. Inhibitors of the first and second generation were used to estimate the inhibitory effects on HIV-1 protease activity. The resulting set of kinetic data supported that the mutant forms are kinetically unaffected by inhibitors of the first generation, proving their functional resistance to these compounds. The second generation inhibitor darunavir inhibited mutant forms MDR-HM and ANAM-11, but was ineffective against prDRV4. Our kinetic data revealed that these inhibitors induced different conformational changes in the enzyme and, thereby they have different mode of binding in the enzyme active site. These data confirmed hypothesis that the driving force of the inhibitor-resistance evolution is disruption of enzyme-inhibitor complex by changing of the contact network in the inhibitor binding site.

Structural Stability of the Coiled-Coil Domain of Tumor Susceptibility Gene (TSG)-101

The tumor susceptibility gene-101 coiled coil domain (TSG101cc) is an integral component of the endosomal maturation machinery and cytokinesis, and also interacts with several transcription factors. The TSG101cc has been crystallized as a homotetramer but is known to interact with two of its binding partners as a heterotrimer. To investigate this apparent discrepancy, we examined the solution thermodynamics of the TSG101cc. Here, we use circular dichroism, differential scanning calorimetry, analytical ultracentrifugation, fluorescence, and structural thermodynamic analysis to investigate the structural stability and the unfolding of the TSG101cc. We demonstrate that TSG101cc exists in solution primarily as a tetramer, which unfolds in a two-state manner. Surprisingly, no homodimeric or homotrimeric species were detected. Structural thermodynamic analysis of the homotetrameric structure and comparison with known oligomeric coiled-coils suggests that the TSG101cc homotetramer is comparatively unstable on a per residue basis. Furthermore, the homotrimeric coiled-coil is predicted to be much less stable than the functional heterotrimeric coiled-coil in the endosomal sorting complex required for transport 1 (ESCRT1). These results support a model whereby the tetramer–monomer equilibrium of TSG101 serves as the cellular reservoir of TSG101, which is effectively outcompeted when its binding partners are present and the heteroternary complex can form.

Development of peptide inhibitors of HIV transmission

Treatment of HIV has long faced the challenge of high mutation rates leading to rapid development of resistance, with ongoing need to develop new methods to effectively fight the infection. Traditionally, early HIV medications were designed to inhibit RNA replication and protein production through small molecular drugs. Peptide based therapeutics are a versatile, promising field in HIV therapy, which continues to develop as we expand our understanding of key protein-protein interactions that occur in HIV replication and infection. This review begins with an introduction to HIV, followed by the biological basis of disease, current clinical management of the disease, therapeutics on the market, and finally potential avenues for improved drug development.

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.

Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS

HIV-1 protease inhibitors continue to play an important role in the treatment of HIV/AIDS, transforming this deadly ailment into a more manageable chronic infection. Over the years, intensive research has led to a variety of approved protease inhibitors for the treatment of HIV/AIDS. In this review, we outline current drug design and medicinal chemistry efforts toward the development of next-generation protease inhibitors beyond the currently approved drugs.

HIV-1 Protease Dimerization Dynamics Reveals a Transient Druggable Binding Pocket at the Interface

The binding mechanism of HIV-1 protease monomers leading to the catalytically competent dimeric enzyme has been investigated by means of state-of-the-art atomistic simulations. The emerging picture allows a deeper understanding of experimental observations and reveals that water molecules trapped at the interface have an important role in slowing down the kinetics of the association process. Unexpectedly, a cryptic binding pocket is identified at the interface of the complex, corresponding to a partially bound dimer that lacks enzymatic function. The pocket has a transient nature with a lifetime longer than 1 μs, and it displays very favorable druggability features. Docking as well as MM-GBSA free-energy calculations further support the possibility to target the new binding site by means of inhibitors able to prevent the complete dimerization by capturing the inactive conformation. This discovery could open the way to the rational design of a new class of anti-HIV drugs.

Conformational stability and self-association equilibrium in biologics

Biologics exist in equilibrium between native, partially denatured, and denatured conformational states. The population of any of these states is dictated by their Gibbs energy and can be altered by changes in physical and solution conditions. Some conformations have a tendency to self-associate and aggregate, an undesirable phenomenon in protein therapeutics. Conformational equilibrium and self-association are linked thermodynamic functions. Given that any associative reaction is concentration dependent, conformational stability studies performed at different protein concentrations can provide early clues to future aggregation problems. This analysis can be applied to the selection of protein variants or the identification of better formulation solutions. In this review, we discuss three different aggregation situations and their manifestation in the observed conformational equilibrium of a protein.

Synthetic, structural mimetics of the β-hairpin flap of HIV-1 protease inhibit enzyme function

Small-molecule mimetics of the β-hairpin flap of HIV-1 protease (HIV-1 PR) were designed based on a 1,4-benzodiazepine scaffold as a strategy to interfere with the flap-flap protein-protein interaction, which functions as a gated mechanism to control access to the active site. Michaelis-Menten kinetics suggested our small-molecules are competitive inhibitors, which indicates the mode of inhibition is through binding the active site or sterically blocking access to the active site and preventing flap closure, as designed. More generally, a new bioactive scaffold for HIV-1PR inhibition has been discovered, with the most potent compound inhibiting the protease with a modest K(i) of 11 μM.

Thermodynamics and kinetics of inhibitor binding to human equilibrative nucleoside transporter subtype-1

Many nucleoside transport inhibitors are in clinical use as anti-cancer, vasodilator and cardioprotective drugs. However, little is known about the binding energetics of these inhibitors to nucleoside transporters (NTs) due to their low endogenous expression levels and difficulties in the biophysical characterization of purified protein with ligands. Here, we present kinetics and thermodynamic analyses of inhibitor binding to the human equilibrative nucleoside transporter-1 (hENT1), also known as SLC29A1. Using a radioligand binding assay, we obtained equilibrium binding and kinetic rate constants of well-known NT inhibitors--[(3)H]nitrobenzylmercaptopurine ribonucleoside ([(3)H]NBMPR), dilazep, and dipyridamole--and the native permeant, adenosine, to hENT1. We observed that the equilibrium binding affinities for all inhibitors decreased whereas, the kinetic rate constants increased with increasing temperature. Furthermore, we found that binding is enthalpy driven and thus, an exothermic reaction, implying that the transporter does not discriminate between its inhibitors and substrates thermodynamically. This predominantly enthalpy-driven binding by four chemically distinct ligands suggests that the transporter may not tolerate diversity in the type of interactions that lead to high affinity binding. Consistent with this, the measured activation energy of [(3)H]NBMPR association was relatively large (20 kcal mol(-1)) suggesting a conformational change upon inhibitor binding. For all three inhibitors the enthalpy (ΔH°) and entropy (ΔS°) contributions to the reaction energetics were determined by van't Hoff analysis to be roughly similar (25-75% ΔG°). Gains in enthalpy with increasing polar surface area of inhibitors suggest that the binding is favored by electrostatic or polar interactions between the ligands and the transporter.

Fluorogenic Assay for Inhibitors of HIV-1 Protease with Sub-picomolar Affinity

A fluorogenic substrate for HIV-1 protease was designed and used as the basis for a hypersensitive assay. The substrate exhibits a k_cat of 7.4 s⁻¹, K_M of 15 μM, and an increase in fluorescence intensity of 104-fold upon cleavage, thus providing sensitivity that is unmatched in a continuous assay of HIV-1 protease. These properties enabled the enzyme concentration in an activity assay to be reduced to 25 pM, which is close to the K_d value of the protease dimer. By fitting inhibition data to Morrison’s equation, K_i values of amprenavir, darunavir, and tipranavir were determined to be 135, 10, and 82 pM, respectively. This assay, which is capable of measuring K_i values as low as 0.25 pM, is well-suited for characterizing the next generation of HIV-1 protease inhibitors.

Denatured state aggregation parameters derived from concentration dependence of protein stability

Protein aggregation is a major issue affecting the long-term stability of protein preparations. Proteins exist in equilibrium between the native and denatured or partially denatured conformations. Often denatured or partially denatured conformations are prone to aggregate because they expose to solvent the hydrophobic core of the protein. The aggregation of denatured protein gradually shifts the protein equilibrium toward increasing amounts of denatured and ultimately aggregated protein. Recognizing and quantitating the presence of denatured protein and its aggregation at the earliest possible time will bring enormous benefits to the identification and selection of optimal solvent conditions or the engineering of proteins with the best stability/aggregation profile. In this article, a new approach that allows simultaneous determination of structural stability and the amount of denatured and aggregated protein is presented. This approach is based on the analysis of the concentration dependence of the Gibbs energy (ΔG) of protein stability. It is shown that three important quantities can be evaluated simultaneously: (i) the population of denatured protein, (ii) the population of aggregated protein, and (iii) the fraction of denatured protein that is aggregated.

Switches of hydrogen bonds during ligand-protein association processes determine binding kinetics

Revealing the processes of ligand-protein associations deepens our understanding of molecular recognition and binding kinetics. Hydrogen bonds (H-bonds) play a crucial role in optimizing ligand-protein interactions and ligand specificity. In addition to the formation of stable H-bonds in the final bound state, the formation of transient H-bonds during binding processes contributes binding kinetics that define a ligand as a fast or slow binder, which also affects drug action. However, the effect of forming the transient H-bonds on the kinetic properties is little understood. Guided by results from coarse-grained Brownian dynamics simulations, we used classical molecular dynamics simulations in an implicit solvent model and accelerated molecular dynamics simulations in explicit waters to show that the position and distribution of the H-bond donor or acceptor of a drug result in switching intermolecular and intramolecular H-bond pairs during ligand recognition processes. We studied two major types of HIV-1 protease ligands: a fast binder, xk263, and a slow binder, ritonavir. The slow association rate in ritonavir can be attributed to increased flexibility of ritonavir, which yields multistep transitions and stepwise entering patterns and the formation and breaking of complex H-bond pairs during the binding process. This model suggests the importance of conversions of spatiotemporal H-bonds during the association of ligands and proteins, which helps in designing inhibitors with preferred binding kinetics.

Human immunodeficiency virus-1 (HIV-1)-mediated apoptosis: new therapeutic targets

HIV has posed a significant challenge due to the ability of the virus to both impair and evade the host’s immune system. One of the most important mechanisms it has employed to do so is the modulation of the host’s native apoptotic pathways and mechanisms. Viral proteins alter normal apoptotic signaling resulting in increased viral load and the formation of viral reservoirs which ultimately increase infectivity. Both the host’s pro- and anti-apoptotic responses are regulated by the interactions of viral proteins with cell surface receptors or apoptotic pathway components. This dynamic has led to the development of therapies aimed at altering the ability of the virus to modulate apoptotic pathways. These therapies are aimed at preventing or inhibiting viral infection, or treating viral associated pathologies. These drugs target both the viral proteins and the apoptotic pathways of the host. This review will examine the cell types targeted by HIV, the surface receptors exploited by the virus and the mechanisms whereby HIV encoded proteins influence the apoptotic pathways. The viral manipulation of the hosts’ cell type to evade the immune system, establish viral reservoirs and enhance viral proliferation will be reviewed. The pathologies associated with the ability of HIV to alter apoptotic signaling and the drugs and therapies currently under development that target the ability of apoptotic signaling within HIV infection will also be discussed.

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.

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.

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.

Introduction of polar groups on the naphthalene scaffold of molecular tongs inhibiting wild-type and mutated HIV-1 protease dimerization

A new series of naphthalene-based molecular tongs containing polar groups at the 3-position of the naphthalene scaffold was synthesized and its anti-dimerization activity was evaluated against HIV-1 protease.

Exploring key orientations at protein-protein interfaces with small molecule probes

Small molecule probes that selectively perturb protein-protein interactions (PPIs) are pivotal to biomedical science, but their discovery is challenging. We hypothesized that conformational resemblance of semirigid scaffolds expressing amino acid side-chains to PPI-interface regions could guide this process. Consequently, a data mining algorithm was developed to sample huge numbers of PPIs to find ones that match preferred conformers of a selected semirigid scaffold. Conformations of one such chemotype (1aaa; all methyl side-chains) matched several biomedically significant PPIs, including the dimerization interface of HIV-1 protease. On the basis of these observations, four molecules 1 with side-chains corresponding to the matching HIV-1 dimerization interface regions were prepared; all four inhibited HIV-1 protease via perturbation of dimerization. These data indicate this approach may inspire design of small molecule interface probes to perturb PPIs.

Energetics of oligomeric protein folding and association

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

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.

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.

Effects of HIV-1 protease on cellular functions and their potential applications in antiretroviral therapy

Human Immunodeficiency Virus Type 1 (HIV-1) protease inhibitors (PIs) are the most potent class of drugs in antiretroviral therapies. However, viral drug resistance to PIs could emerge rapidly thus reducing the effectiveness of those drugs. Of note, all current FDA-approved PIs are competitive inhibitors, i.e., inhibitors that compete with substrates for the active enzymatic site. This common inhibitory approach increases the likelihood of developing drug resistant HIV-1 strains that are resistant to many or all current PIs. Hence, new PIs that move away from the current target of the active enzymatic site are needed. Specifically, allosteric inhibitors, inhibitors that prohibit PR enzymatic activities through non-competitive binding to PR, should be sought. Another common feature of current PIs is they were all developed based on the structure-based design. Drugs derived from a structure-based strategy may generate target specific and potent inhibitors. However, this type of drug design can only target one site at a time and drugs discovered by this method are often associated with strong side effects such as cellular toxicity, limiting its number of target choices, efficacy, and applicability. In contrast, a cell-based system may provide a useful alternative strategy that can overcome many of the inherited shortcomings associated with structure-based drug designs. For example, allosteric PIs can be sought using a cell-based system without considering the site or mechanism of inhibition. In addition, a cell-based system can eliminate those PIs that have strong cytotoxic effect. Most importantly, a simple, economical, and easy-to-maintained eukaryotic cellular system such as yeast will allow us to search for potential PIs in a large-scaled high throughput screening (HTS) system, thus increasing the chances of success. Based on our many years of experience in using fission yeast as a model system to study HIV-1 Vpr, we propose the use of fission yeast as a possible surrogate system to study the effects of HIV-1 protease on cellular functions and to explore its utility as a HTS system to search for new PIs to battle HIV-1 resistant strains.

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.

Carbonylhydrazide-based molecular tongs inhibit wild-type and mutated HIV-1 protease dimerization

We have designed and synthesized new molecular tongs based on a rigid naphthalene scaffold and evaluated their antidimer activity on HIV-1 protease (PR). We inserted carbonylhydrazide and oligohydrazide (azatide) fragments into their peptidomimetic arms to reduce hydrophobicity and increase metabolic stability. These fragments are designed to disrupt the protein-protein interactions by reproducing the hydrogen bond pattern found in the antiparallel β-sheet formed between the N- and C-ends of the two monomers in the native PR. Kinetic analyses and fluorescent probe binding studies showed that several molecular tongs can inhibit PR dimerization. The best nonpeptidic molecular tongs to date were obtained with an inhibition constant K(id) of 50 nM for PR and 80 nM for the multimutated protease ANAM-11. The PR inhibition was selective, the aspartic proteases renin and pepsin were not inhibited.

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.