Previous Article | Next Article ![]()
Journal of Virology, December 2007, p. 13845-13851, Vol. 81, No. 24
0022-538X/07/$08.00+0 doi:10.1128/JVI.01184-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Tibotec BVBA, Mechelen, Belgium
Received 31 May 2007/ Accepted 30 September 2007
|
|
|---|
|
|
|---|
Currently, there are nine HIV-1 PIs in clinical use (amprenavir [APV], atazanavir [ATV], darunavir [DRV], indinavir [IDV], lopinavir [LPV], nelfinavir [NFV], ritonavir [RTV], saquinavir [SQV], and tipranavir [TPV]), and additional inhibitors are undergoing clinical trials (19). Despite the success of PIs in HAART, the emergence of multidrug-resistant (MDR) viruses has resulted in a high incidence of cross-resistance between PIs, which limits their sequential use, thereby reducing the treatment options available for HIV patients (4, 26). Thus, there is an ongoing clinical requirement for new antiviral agents with broad-spectrum activity against PI-resistant mutants for the long-term management of HIV infection. DRV, previously known as TMC114, is a very effective PI (3, 28) with high in vitro and in vivo potency against wild-type (WT) and MDR HIV-1. Moreover, DRV appears to have a very high genetic barrier to the development of resistance (2).
Although drug resistance is a complex phenomenon, the most important effect at the molecular level is the introduction of mutations within the protease gene, resulting in an enzyme with significantly decreased binding affinity for the inhibitor, with respect to the protease substrates (11). An inhibitor with a good resistance profile would possess high binding affinity to WT protease which would not be significantly reduced by mutations associated with drug resistance. Inhibition studies and in vitro resistance selection experiments are insufficient for delineating the interactions between HIV-1 protease and its inhibitors, in regard to resistance. Essential information on drug-target interactions can be derived from thermodynamic and kinetic binding studies.
By assessing the binding thermodynamics of DRV, using isothermal titration calorimetry, it has been shown that DRV binds tightly to WT HIV-1 protease with a favorable enthalpy due to the strong interactions between the bis-tetrahydrofuranyl urethane moiety of DRV and the main-chain atoms of Asp29 and Asp30 in the protease-active site (12, 13). This high affinity together with the close fit of DRV within the substrate envelope could explain DRV's high potency and much lower level of resistance selection.
Biosensor-based studies, or surface plasmon resonance (SPR) technology, are an efficient and sensitive methodology used to analyze the binding kinetics of an interaction (21, 25). This type of assay allows the measurement of an association rate (kon) next to a dissociation rate (koff), which are indicators for, respectively, the degree of recognition and the binding stability of the ligand-analyte interaction. The ratio of the two is a measure of the affinity or binding strength of the analyte for the ligand (KD = koff/kon), meaning that a high affinity can be caused by a high association or a slow dissociation rate (7).
To further establish the molecular basis for DRV's high antiviral potency and broad-spectrum activity, we measured the binding affinities of DRV and other currently approved PIs to WT and MDR proteases by using SPR. The MDR proteases were derived from clinical isolates containing 10 to 14 International AIDS Society-U.S.A. (IAS-USA) PI resistance-associated mutations (RAMs) (6) and included mutants with a decreased susceptibility to DRV. The MDR proteases contained 1 to 5 of the 11 mutations (V11I, V32I, L33F, I47V, I50V, I54L/M, G73S, L76V, I84V, and L89V) identified as those associated with a diminished DRV/ritonavir virologic response and present mostly with a high number of other PI RAMs (S. De Meyer, T. Vangeneugden, E. Lefebvre, H. Azijn, I. De Baere, B. Van Baelen, and M.-P. de Béthune, phenotypic and genotypic determinants of resistance to TMC114: pooled analysis of POWER 1, 2, and 3, presented at the XVth International HIV Drug Resistance Workshop, Sitges, Spain, 13 to 17 June 2006).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Genotypic profile of the MDR proteases derived from clinical isolatesa
|
The PIs ATV, APV, IDV, LPV, NFV, RTV, and SQV were purified from commercial formulations. DRV and TPV were produced as previously reported (28, 29). The purity of all PIs was >95% as determined by liquid chromatography-mass spectrometry.
Viruses. Recombinant HIV-1 isolates derived from clinical samples were constructed as previously described by cotransfection of MT4 cells with sample-derived viral protease and reverse transcriptase coding sequences and an HIV-1 HXB2-derived proviral clone deleted in the protease and reverse transcriptase coding region (5).
Antiviral assay. The antiviral activities of PIs against the laboratory-adapted WT (HIV-1-LAI) isolate and clinical sample-derived recombinant viruses were determined with a cell-based HIV-1 replication assay. MT4 cells equipped with long-terminal repeat-enhanced green fluorescent protein (EGFP) (150,000 cells/ml) were infected with HIV-1 (at a multiplicity of infection of 0.01) in the presence or absence of different inhibitor concentrations. Upon infection by HIV-1, expression of the viral Tat product increased transcription from the HIV-1 LTR promoter, leading to high-level expression of the EGFP reporter. After 3 days of incubation, the amount of HIV-1 replication was quantified by measuring the EGFP fluorescence and expressed as the 50% effective concentration (EC50) value or as the (n-fold) change in susceptibility by dividing the EC50 value for the virus tested by the EC50 value for the WT virus tested in parallel.
Interaction studies. Measurements of interactions between WT or mutant proteases and inhibitors were performed using a Biacore S51 instrument (Biacore AB, Uppsala, Sweden), essentially as described previously (17). Proteases were immobilized using amine coupling to CM5 sensor chips (Biacore AB) with an extra cross-linking step of 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC) and N-hydroxysuccinimide (NHS) for 7 min (WT) or 2 min (MDR) to stabilize the surface. Protease enzymes were dissolved at a concentration range of 2 to 10 µM in acetate buffer at pH 4.5 and immobilized at a flow rate of 10 µl/min and a contact time that varied between 10 and 20 min. In Biacore instruments, the amount of immobilized protein is measured in resonance units (RU), with one RU representing about 1 pg of protein bound per square millimeter of flow cell surface. In our system, we routinely reached a level of 1,500 to 3,000 RU of immobilized protein. Interaction studies were all performed at 20°C with a flow rate of 90 µl/min to minimize mass transport effects. Inhibitors were serially diluted in phosphate-buffered saline containing 5% dimethyl sulfoxide (DMSO) and 3 mM EDTA and supplemented with 0.005% or 0.05% P20, a nonionic surfactant, for studies with the WT and mutant proteases, respectively. The higher P20 level for mutant protease binding studies was needed to obtain binding levels that were high enough to study the weaker interactions of the mutant proteases. Increasing the P20 concentration from 0.005% to 0.05% did not influence PI binding characteristics when applied to WT protease binding studies (data not shown). The inhibitor concentration range was chosen based on the antiviral activity (EC50) of the inhibitor on the corresponding MDR virus. The compounds were injected for a period of 215 s and dissociated for a period of 1,000 s. A duplicate sample and a zero concentration sample were used as a positive and negative control, respectively. A DMSO concentration series was included to eliminate the contribution of a bulk solution to the measured response. Protease immobilization levels of 1,500 to 3,000 RU resulted in binding responses of 5 to 20 RU at the highest compound concentrations. The sensor surface was regenerated in between sample injections by injecting 100% ethylene glycol, and the flow system was washed with 50% DMSO. Buffer injections were conducted to detect possible carryover between samples.
Data analysis. As it was not possible to reach an equilibrium state in most of the experiments, all analyses were based on the evaluation of binding kinetics (8). Data were analyzed using simultaneous nonlinear regression analysis (global fitting) with Biacore S51 or Biacore T100 BiaEval evaluation software (Biacore AB).
The individual rate constants kon and koff and a derived affinity constant, KD = koff/kon, were determined by a kinetic evaluation of the sensorgrams fitting a 1:1 model or a heterogeneous ligand model. Experimentally, compound concentrations exceeding the KD value need to be used to allow a reliable analysis. For WT protease, where the KD value is low, experiments were done with correspondingly low concentrations of compound, and binding effects nonspecific to the protease target were undetectable. Since the evaluated inhibitors were expected to have only one binding site on the protease, the 1:1 Langmuir binding model was used (15). This model, assuming a single type of target binding site, was found to accurately fit the binding. For mutant proteases, the 1:1 model mismatched kinetic data due to additional binding components, particularly at high compound concentrations. Experiments using immobilized HIV-1 gp120 protein (data not shown) revealed that the additional binding components were not specific to the protease target but reflected the propensity of the tested PIs to bind nonspecifically to proteins. Therefore, the heterogeneous ligand model, which is essentially the sum of two independent 1:1 binding models (two kon, koff pairs), was used (10, 20). This model is a strong simplification, as nonspecific binding involves many different locations with many different binding constants, resulting in a binding level that does not saturate with increasing compound concentrations. However, the existing software did not allow more accurate modeling. Of the two sets of kinetic parameters generated, only one was considered relevant to compound activity. Components that occurred at the highest concentrations, that possessed unstable binding characteristics (high koff rates), and that did not saturate as a function of increasing concentration were categorized as nonspecific binding, but coverage of the 0.1 to 10x KD range was assured. The applied fit models accounted for limited mass transport effects.
To prevent suboptimal fitting, we required all important features of the data to be represented by the model and verified that no component of the model was dependent on minor features that could be due to experimental deviations. Every analysis was performed on at least two independent experiments, and it was verified that the fit values for kon, koff, or KD were reproduced within 1 order of magnitude.
Detection limits. Instrument limitations resulted in detection limits of 104 to 5 x 106 M–1s–1 for the kon value and 10–1 to 10–5 s–1 for the koff value (R. Karlsson, personal communication; 23). Values outside these detection limits (Table 2 and Table 3, bold type) should be interpreted as censored values, indicating kinetics that were slower or faster than the instrument was able to measure. The term kinetic limitation is used in Table 3 if the kon and koff values were so far outside the detection range that even an approximate determination was not possible. The fast kinetics of PI binding on MDR mutant proteases was the main reason for performing our measurements at the relatively low temperature of 20°C. Although a higher temperature would allow a more accurate determination of the DRV koff value for WT protease (koff value above the lower detection limit), it would compromise the measurements of mutants (koff value beyond the upper detection limit).
|
View this table: [in a new window] |
TABLE 2. Kinetic parameters for the binding of PIs to WT HIV-1 proteasea
|
|
View this table: [in a new window] |
TABLE 3. Kinetic parameters of the interaction between MDR HIV-1 protease (PR) mutants and PIsa
|
|
|
|---|
![]() View larger version (33K): [in a new window] |
FIG. 1. Overlay of the experimental sensorgrams of the interaction between WT protease and its inhibitors at a concentration of 40 nM (ATV, IDV, LPV, NFV, RTV, SQV, and TPV) or 50 nM (APV and DRV). The sensorgrams were normalized to a maximum response of 1. For clarity, sensorgrams are spread over two graphs, with DRV included in both graphs as a reference point.
|
2 x 106 M–1s–1). The dissociation of DRV from WT protease was extremely slow, up to 1,000 times slower than the dissociation for the other PIs tested, including APV. In fact, the koff value of the DRV-WT interaction was too slow to be reliably measured (koff of <10–5 s–1). This interaction profile of DRV resulted in a very tight binding to WT protease, with an affinity that was more than 2 orders of magnitude higher (KD < 10–12 M) than that for the other PIs. The extremely high stability of the binding of DRV to WT protease is also illustrated in the interaction kinetic diagram (Fig. 2, detection limits of the assay determined by biophysical and instrument limitations are shaded) (R. Karlsson, personal communication; 23), in which the diagonal lines represent different combinations of kon and koff values, yielding identical affinities. The binding affinity of TPV to WT protease was in the same range as that of APV, ATV, and LPV, whereas its in vitro antiviral activity was at least 10-fold lower. The kinetic parameters determined were in line with previously published results, as indicated in the KD reference value column in Table 2 (27).
![]() View larger version (24K): [in a new window] |
FIG. 2. WT protease interaction kinetic map for all currently used PIs shown as kon and koff values and the combinations of kon and koff that result in the same KD values (diagonal lines). The area beyond the detection limits due to biophysical and instrument limitations is shaded.
|
Kinetics of inhibitor binding to MDR protease mutants. The binding kinetics of the clinical PIs APV, ATV, DRV, LPV, and TPV with five MDR protease mutants were studied. These protease mutants were derived from clinical isolates with a decreased susceptibility to multiple HIV-1 PIs and contained 10 to 14 PI RAMs, based on the IAS-USA fall 2005 guidelines (6). The genotypes of the mutants are shown in Table 1.
In general, all PIs tested displayed a weaker binding affinity to the MDR proteases than to the WT protease. For DRV, the decrease in binding affinity to MDR proteases compared to that of WT protease was caused mainly by a faster dissociation rate, as the drop in association rate was less than 10-fold (Table 3 and Fig. 3) in all experiments. The association rates of DRV with the MDR proteases were higher than or similar to those of the other tested PIs, except for the binding of LPV to the mutant E (LPV, kon of
3.2 x 106 M–1s–1, versus DRV, kon of
5.6 x 105 M–1s–1). The dissociation rates of DRV from the MDR proteases were more than 1,000-fold faster than that from the WT protease but remained in most cases slower than or similar to those of the other tested PIs. The dissociation rates of ATV and TPV were slower than that of DRV for mutants B and C and mutants A and B, respectively, resulting in a higher binding affinity. The binding affinity of DRV was proportional to the number of mutations associated with a diminished susceptibility to DRV (De Meyer et al.). The mutants A, B, and C, containing four or five of such mutations, displayed the lowest binding affinity and hence a large (n-fold) change in antiviral activity. An accurate kinetic characterization of the binding of APV and LPV to mutants A and B was not possible because the interactions were too unstable to be measured, with very high association and dissociation rates (kinetic limitations). ATV and LPV showed a high decrease in binding affinity for mutant D, caused by a change in both the association and dissociation rates. Also for APV and TPV, there was a drop in both the association and the dissociation rates with mutant E, resulting in a decreased binding strength. The binding strength of DRV for mutant D and E proteases was at least 100-fold higher than for the other PIs, including APV, although it was lower than that for the WT protease. The association rates of TPV to the mutants A, B, and E were approximately 10-fold slower than that of DRV. For mutants C and D, extremely low binding levels (mutant D) and an association and/or dissociation rate outside the detection limits (mutants C and D) prevented accurate determination of TPV's kinetic parameters and were scored as not determinable (Table 3). Attempts to increase the immobilization level of mutant D and hence the binding level of TPV were unsuccessful (data not shown).
![]() View larger version (25K): [in a new window] |
FIG. 3. Interaction kinetic map of binding studies between APV, ATV, DRV, LPV, and TPV on WT and MDR proteases. The area beyond the detection limits due to biophysical and instrument limitations is shaded.
|
![]() View larger version (19K): [in a new window] |
FIG. 4. Relationship between binding affinity (KD) and antiviral activity (EC50) of PIs to HIV-1 MDR proteases.
|
|
|
|---|
With the WT protease, the binding affinity of DRV was more than 2 orders of magnitude higher than that of other approved PIs. This was due mainly to the extremely slow dissociation of DRV from WT protease, although its association also compared favorably to that of all tested PIs, with the exception of that of APV, which showed a similar association rate. However, while the association rates of DRV and its structural analogue APV were similar, the more than 100-fold difference in dissociation rates resulted in a much higher binding affinity of DRV to WT protease than with APV. Also, the thermodynamic binding characteristics determined with isothermal titration calorimetry reflect the high binding affinity of DRV to WT protease (12). DRV has been shown to have a significantly higher binding affinity than narrow-spectrum PIs, determined by both favorable enthalpy and entropy changes. In thermodynamic studies, the binding of DRV to WT protease was also approximately 2 orders of magnitude tighter than that of APV, with an extremely favorable enthalpy of binding for DRV compared to that of APV.
The koff value is believed to be a key indicator of in vivo duration of efficacy, more significant than the apparent affinity of the drug for its target, per se (1). The dissociative t1/2 value, which can be derived from the dissociation rate, is a direct measure of the interaction time between a drug and its target. DRV's extremely high t1/2 value with WT protease (>240 h) suggests that, in vivo, the drug remains bound and thus active throughout its elimination from plasma (terminal elimination half-life in plasma of
15 h when coadministered with RTV [V. Sekar, S. Spinosa-Guzman, E. Lefebvre, and R. Hoetelmans, clinical pharmacology of TMC114: a new HIV protease inhibitor, presented at the 16th International AIDS Conference, Toronto, Canada, 13 to 18 August 2006]). Thus, the duration of DRV action in vivo may be extended in comparison with that of the other PIs as ATV, LPV, SQV and TPV have a t1/2 between 1.2 and 1.8 h, while the interaction time of APV, IDV, and RTV with WT protease was limited to approximately 15 min.
PI resistance resulting from multiple protease mutations could arise as a result of decreases in association rates as well as increases in dissociation rates. In our experiments, the decrease in PI binding affinity with the studied MDR mutants was in general related to a faster dissociation rate. For DRV, a dissociation rate more than 1,000-fold faster than that of WT protease was the main cause of lower affinity as the decrease in association rate was, in all experiments, less than 10-fold. For APV, ATV, LPV, and TPV, the decrease in association rate was, for some mutants, also a contributing factor to the decreased affinity compared to that of WT protease. The association rates of DRV with the MDR proteases remained higher than or in the same range as those of the other PIs. Also the dissociation rates of DRV were slower than or similar to that of the other PIs, except for those mutants containing four to five protease mutations associated with a diminished susceptibility to DRV. For these mutants, ATV and/or TPV showed a slower dissociation rate and hence a higher binding affinity than that of DRV. The binding affinity of DRV to MDR proteases was proportional to the number of DRV RAMs. As such, decreased binding affinity to MDR proteases explains the molecular mechanism of resistance to DRV and the other PIs. DRV displayed a higher binding affinity to the MDR proteases than the chemically related APV. The binding of APV to the MDR mutants with five DRV-associated mutations was even too unstable to be measured. Although the antiviral activity of TPV was significantly lower than that of APV, ATV, and LPV, the binding affinity of TPV to WT protease was in the same range as that of these other PIs. This discrepancy between binding affinity and antiviral activity could be due to a limited permeability of TPV. The association rates of TPV with both WT and MDR proteases were also systematically 10-fold lower than that of DRV.
A more than 1,000-fold decrease in affinity for DRV compared to that of WT protease as observed for mutants A to D was needed to decrease its corresponding antiviral effect. For APV, ATV, LPV, and TPV, a decrease in binding affinity correlated from the onset with a drop in antiviral activity. The large drop in binding affinity of DRV that is sustained without compromising the corresponding antiviral activity provides a mechanistic background for the higher genetic barrier of DRV to development of resistance. The exceptionally high affinity of DRV for WT protease is a prerequisite for making this lag phase in affinity-activity correlation possible.
It has been suggested that DRV may have two binding sites on WT protease and that this would contribute to the high genetic barrier to resistance of DRV (14). Two binding sites with distinct binding kinetics would result from deviations from the 1:1 model on condition that the KD values for both binding interactions were within our studied concentration range. However, the 1:1 model accurately fitted the binding of all PIs studied, including DRV, to WT protease, indicating that any binding to alternate sites must be at least an order of magnitude weaker. A second binding component did appear in experiments with MDR proteases if higher compound concentrations had to be used, but the characteristics of this component indicated that it corresponded to nonspecific binding and not to binding to a protease target site.
The high impact of the dissociation rate of DRV on its affinity value confirms the relevance of the koff value for evaluation of the efficacy of a compound. As DRV already starts with an extremely high affinity to WT protease, it can sustain a substantial decrease in its binding affinity on MDR mutants without the loss of antiviral activity, allowing it to remain effective, as was shown in the POWER studies (9). The very high binding affinity of DRV to WT protease necessitates the simultaneous occurrence of multiple mutations before the virus becomes resistant to DRV. Further kinetic characterization will explore the influence of individual mutations on DRV binding to HIV-1 protease.
Published ahead of print on 10 October 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»