Configurations of diastereomeric hydroxyethylene isosteres strongly affect biological activities of a series of specific inhibitors of human-immunodeficiency-virus proteinase

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

Background

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

Results

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

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

Conclusions

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

Ligand conformational and solvation/desolvation free energy in protein-ligand complex formation

In this study, an extensive sampling of the conformational space of nine HIV-1 protease inhibitors was performed to estimate the uncertainty with which a single-conformation scoring scheme approximates the ligand-protein binding free energy. The SMD implicit solvation/desolvation energy and gas-phase PM6-DH2 energy were calculated for a set of 1600 conformations of each ligand. The probability density functions of the energies were compared with the values obtained from the single-conformation approach and from a short ab initio molecular dynamics simulation. The relative uncertainty in the score within the set of nine inhibitors was calculated to be 3.5 kcal·mol(-1) and 2.7 kcal·mol(-1) for the single-conformation and short dynamics, respectively. These results, though limited to the consideration of flexible ligands, provide a valuable insight into the precision of rigid models in the current computer-aided drug design.

A reliable docking/scoring scheme based on the semiempirical quantum mechanical PM6-DH2 method accurately covering dispersion and H-bonding: HIV-1 protease with 22 ligands

In this study, we introduce a fast and reliable rescoring scheme for docked complexes based on a semiempirical quantum mechanical PM6-DH2 method. The method utilizes a PM6-based Hamiltonian with corrections for dispersion energy and hydrogen bonds. The total score is constructed as the sum of the PM6-DH2 interaction enthalpy, the empirical force field (AMBER) interaction entropy, and the sum of the deformation (PM6-DH2, SMD) and the desolvation (SMD) energies of the ligand. The main advantage of the procedure is the fact that we do not add any empirical parameter for either an individual component of the total score or an individual protein-ligand complex. This rescoring method is applied to a very challenging system, namely, the HIV-1 protease with a set of ligands. As opposed to the conventional DOCK procedure, the PM6-DH2 rescoring based on all of the terms distinguishes between binders and nonbinders and provides a reliable correlation of the theoretical and experimental binding free energies. Such a dramatic improvement, resulting from the PM6-DH2 rescoring of all the complexes, provides a valuable yet inexpensive tool for rational drug discovery and de novo ligand design.

Enzymatic and structural analysis of the I47A mutation contributing to the reduced susceptibility to HIV protease inhibitor lopinavir

Lopinavir (LPV) is a second-generation HIV protease inhibitor (PI) designed to overcome resistance development in patients undergoing long-term antiviral therapy. The mutation of isoleucine at position 47 of the HIV protease (PR) to alanine is associated with a high level of resistance to LPV. In this study, we show that recombinant PR containing a single I47A substitution has the inhibition constant (K(i) ) value for lopinavir by two orders of magnitude higher than for the wild-type PR. The addition of the I47A substitution to the background of a multiply mutated PR species from an AIDS patient showed a three-order-of-magnitude increase in K(i) in vitro relative to the patient PR without the I47A mutation. The crystal structure of I47A PR in complex with LPV showed the loss of van der Waals interactions in the S2/S2' subsites. This is caused by the loss of three side-chain methyl groups due to the I47A substitution and by structural changes in the A47 main chain that lead to structural changes in the flap antiparallel beta-strand. Furthermore, we analyzed possible interaction of the I47A mutation with secondary mutations V32I and I54V. We show that both mutations in combination with I47A synergistically increase the relative resistance to LPV in vitro. The crystal structure of the I47A/I54V PR double mutant in complex with LPV shows that the I54V mutation leads to a compaction of the flap, and molecular modeling suggests that the introduction of the I54V mutation indirectly affects the strain of the bound inhibitor in the PR binding cleft.

Molecular analysis of the HIV-1 resistance development: enzymatic activities, crystal structures, and thermodynamics of nelfinavir-resistant HIV protease mutants

Human immunodeficiency virus (HIV) encodes an aspartic protease (PR) that cleaves viral polyproteins into mature proteins, thus leading to the formation of infectious particles. Protease inhibitors (PIs) are successful virostatics. However, their efficiency is compromised by antiviral resistance. In the PR sequence of viral variants resistant to the PI nelfinavir, the mutations D30N and L90M appear frequently. However, these two mutations are seldom found together in vivo, suggesting that there are two alternative evolutionary pathways leading to nelfinavir resistance. Here we analyze the proteolytic activities, X-ray structures, and thermodynamics of inhibitor binding to HIV-1 PRs harboring the D30N and L90M mutations alone and in combination with other compensatory mutations. Vitality values obtained for recombinant mutant proteases and selected PR inhibitors confirm the crucial role of mutations in positions 30 and 90 for nelfinavir resistance. The combination of the D30N and L90M mutations significantly increases the enzyme vitality in the presence of nelfinavir, without a dramatic decrease in the catalytic efficiency of the recombinant enzyme. Crystal structures, molecular dynamics simulations, and calorimetric data for four mutants (D30N, D30N/A71V, D30N/N88D, and D30N/L90M) were used to augment our kinetic data. Calorimetric analysis revealed that the entropic contribution to the mutant PR/nelfinavir interaction is less favorable than the entropic contribution to the binding of nelfinavir by wild-type PR. This finding is supported by the structural data and simulations; nelfinavir binds most strongly to the wild-type protease, which has the lowest number of protein-ligand hydrogen bonds and whose structure exhibits the greatest degree of fluctuation upon inhibitor binding.

HIV-1 Protease Mutations and Inhibitor Modifications Monitored on a Series of Complexes. Structural Basis for the Effect of the A71V Mutation on the Active Site

Two new X-ray structures of an HIV-1 protease mutant (A71V, V82T, I84V) in complex with inhibitors SE and SQ, pseudotetrapeptide inhibitors with an acyclic S-hydroxyethylamine isostere, were determined. Comparison of eight structures exploring the binding of four similar inhibitors--SE, SQ (S-hydroxyethylamine isostere), OE (ethyleneamine), and QF34 (hydroxyethylene)--to wild-type and A71V/V82T/I84V HIV-1 protease elucidates the principles of altered interaction with changing conditions. The A71V mutation, which is distant from the active site, causes changes in the structure of the enzyme detectable by the means of X-ray structure analysis, and a route of propagation of the effect toward the active site is proposed.

Role of hydroxyl group and R/S configuration of isostere in binding properties of HIV-1 protease inhibitors

The crystal structure of the complex between human immunodeficiency virus type 1 (HIV-1) protease and a peptidomimetic inhibitor of ethyleneamine type has been refined to R factor of 0.178 with diffraction limit 2.5 A. The peptidomimetic inhibitor Boc-Phe-Psi[CH2CH2NH]-Phe-Glu-Phe-NH2 (denoted here as OE) contains the ethyleneamine replacement of the scissile peptide bond. The inhibitor lacks the hydroxyl group which is believed to mimic tetrahedral transition state of proteolytic reaction and thus is suspected to be necessary for good properties of peptidomimetic HIV-1 protease inhibitors. Despite the missing hydroxyl group the inhibition constant of OE is 1.53 nm and it remains in the nanomolar range also towards several available mutants of HIV-1 protease. The inhibitor was found in the active site of protease in an extended conformation with a unique hydrogen bond pattern different from hydroxyethylene and hydroxyethylamine inhibitors. The isostere nitrogen forms a hydrogen bond to one catalytic aspartate only. The other aspartate forms two weak hydrogen bridges to the ethylene group of the isostere. A comparison with other inhibitors of this series containing isostere hydroxyl group in R or S configuration shows different ways of accommodation of inhibitor in the active site. Special attention is devoted to intermolecular contacts between neighbouring dimers responsible for mutual protein adhesion and for a special conformation of Met46 and Phe53 side chains not expected for free protein in water solution.

General Inverse Solid-Phase Synthesis Method for C-Terminally Modified Peptide Mimetics

Peptide mimetics are of considerable interest as bioactive agents and drugs. C-terminally modified peptide mimetics are of particular interest given the synthetic versatility of the carboxyl group and its derivatives. A general approach to C-terminally modified peptide mimetics, based on a urethane attachment strategy and amino acid t-butyl ester-based N-to-C peptide synthesis, is described. This approach is compatible with the reaction conditions generally employed for solution-phase peptide mimetic synthesis. To develop and demonstrate this approach, it was employed for the solid-phase synthesis of peptide trifluoromethyl ketones, peptide boronic acids, and peptide hydroxamic acids. The development of a versatile general approach to C-terminally modified peptides using readily available starting materials provides a basis for the combinatorial and parallel solid-phase synthesis of these peptide mimetic classes for bioactive agent screening and also provides a basis for the further development of solid-phase C-terminal functional group elaboration strategies.

Efficiency of a second-generation HIV-1 protease inhibitor studied by molecular dynamics and absolute binding free energy calculations

A subnanomolar inhibitor of human immunodeficiency virus type 1 (HIV-1) protease, designated QF34, potently inhibits the wild-type and drug-resistant enzyme. To explain its broad activity, the binding of QF34 to the wild-type HIV-1 protease is investigated by molecular dynamics simulations and compared to the binding of two inhibitors that are used clinically, saquinavir (SQV) and indinavir (IDV). Analysis of the flexibility of protease residues and inhibitor segments in the complex reveals that segments of QF34 were more mobile during the dynamics studies than the segments of SQV and IDV. The dynamics of hydrogen bonding show that QF34 forms a larger number of stable hydrogen bonds than the two inhibitors that are used clinically. Absolute binding free energies were calculated with molecular mechanics-generalized Born surface area (MM-GBSA) methodology using three protocols. The most consistent results were obtained using the single-trajectory approach, due to cancellation of errors and inadequate sampling in the separate-trajectory protocols. For all three inhibitors, energy components in favor of binding include van der Waals and electrostatic terms, whereas polar solvation and entropy terms oppose binding. Decomposition of binding energies reveals that more protease residues contribute significantly to the binding of QF34 than to the binding of SQV and IDV. Moreover, contributions from protease main chains and side chains are balanced in the case of QF34 (52:48 ratio, respectively), whereas side chain contributions prevail in both SQV and IDV (main-chain:side-chain ratios of 41:59 and 45:55, respectively). The presented results help explain the ability of QF34 to inhibit multiple resistant mutants and should be considered in the design of broad-specificity second-generation HIV-1 protease inhibitors.

Kinetics of the dimerization of retroviral proteases: the "fireman's grip" and dimerization

All retroviral proteases belong to the family of aspartic proteases. They are active as homodimers, each unit contributing one catalytic aspartate to the active site dyad. An important feature of all aspartic proteases is a conserved complex scaffold of hydrogen bonds supporting the active site, called the "fireman's grip," which involves the hydroxyl groups of two threonine (serine) residues in the active site Asp-Thr(Ser)-Gly triplets. It was shown previously that the fireman's grip is indispensable for the dimer stability of HIV protease. The retroviral proteases harboring Ser in their active site triplet are less active and, under natural conditions, are expressed in higher enzyme/substrate ratio than those having Asp-Thr-Gly triplet. To analyze whether this observation can be attributed to the different influence of Thr or Ser on dimerization, we prepared two pairs of the wild-type and mutant proteases from HIV and myeloblastosis-associated virus harboring either Ser or Thr in their Asp-Thr(Ser)-Gly triplet. The equilibrium dimerization constants differed by an order of magnitude within the relevant pairs. The proteases with Thr in their active site triplets were found to be approximately 10 times more thermodynamically stable. The dimer association contributes to this difference more than does the dissociation. We propose that the fireman's grip might be important in the initial phases of dimer formation to help properly orientate the two subunits of a retroviral protease. The methyl group of threonine might contribute significantly to fixing such an intermediate conformation.

Characterisation of Mutated Proteinases Derived from HIV-Positive Patients: Enzyme Activity, Vitality and Inhibition

HIV protease (PR) specifically cleaves viral polyproteins to yield infectious progeny virus particles. Inactivation of PR leads to loss of virus infectivity and PR thus became an attractive pharmaceutic target. Indeed, seven protease inhibitors (PI) have been approved for clinical use to date. However, emerging resistant viral variants with reduced sensitivity to PIs become a major obstacle to successful control of viral replication. We have previously reported the design, testing and structural analysis of a pseudopeptide inhibitor, QF34, which efficiently inhibits a wide variety of PR variants. In a clinical study, we have monitored more than 100 HIV-positive patients in the Czech Republic undergoing highly active antiretroviral therapy including PI. In this paper we describe kinetic characterisation of two highly resistant PR species isolated from these patients. The mutated proteases accumulated as much as 14 amino acid exchanges and develop resistance to saquinavir, ritonavir, indinavir and nelfinavir with vitality value up to 150. Kinetic analyses revealed that second-generation PI lopinavir and QF34 retained their subnanomolar potency against both multidrug resistant PR variants. These results suggest a route to the design of PIs capable of inhibiting a variety of resistant PR mutants.

An ethylenamine inhibitor binds tightly to both wild type and mutant HIV-1 proteases. Structure and energy study

An X-ray structure (resolution 2.2 A) of mutant HIV-1 protease (A71V, V82T, I84V) complexed with a newly developed peptidomimetic inhibitor with an ethylenamine isostere Boc-Phe-Psi[CH(2)CH(2)NH]-Phe-Glu-Phe-NH(2), denoted as OE, is described and compared with the complex of wild-type HIV-1 protease with the same inhibitor (resolution 2.5 A). OE shows tight binding to the wild type (K(i) = 1.5 nM) as well as mutant (K(i) = 4.1 nM) protease. The hydrogen bonds formed, in the case of hydroxyethylamine inhibitors, by a hydroxyl group are, in the case of OE inhibitors, replaced by a bifurcated hydrogen bond from the isosteric NH group to both catalytic aspartates Asp 25 and Asp 125. The binding modes of OE inhibitor to the wild type and mutant protease are similar. However, in the mutant protease, weaker van der Waals interactions of the mutated residues Val 84 and Val 184 with OE were found. This lack of interaction energy is compensated by a new aromatic hydrogen bond between the phenyl ring of the inhibitor in position P1 and the mutated residue Thr 182. Energy analysis based on molecular mechanics has been performed to distinguish between the static and dynamic backgrounds of disorder observed at the mutation sites Thr 82, Val 84, Thr 182, and Val 184.

Unusual binding mode of an HIV-1 protease inhibitor explains its potency against multi-drug-resistant virus strains

Protease inhibitors (PIs) are an important class of drugs for the treatment of HIV infection. However, in the course of treatment, resistant viral variants with reduced sensitivity to PIs often emerge and become a major obstacle to successful control of viral load. On the basis of a compound equipotently inhibiting HIV-1 and 2 proteases (PR), we have designed a pseudopeptide inhibitor, QF34, that efficiently inhibits a wide variety of PR variants. In order to analyze the potency of the inhibitor, we constructed PR species harboring the typical (signature) mutations that confer resistance to commercially available PIs. Kinetic analyses showed that these mutated PRs were inhibited up to 1,000-fold less efficiently by the clinically approved PIs. In contrast, all PR species were effectively inhibited by QF34. In a clinical study, we have monitored 30 HIV-positive patients in the Czech Republic undergoing highly active antiretroviral therapy, and have identified highly PI resistant variants. Kinetic analyses revealed that QF34 retained its subnanomolar potency against multi-drug resistant PR variants. X-ray crystallographic analysis and molecular modeling experiments explained the wide specificity of QF34: this inhibitor binds to the PR in an unusual manner, thus avoiding contact sites that are mutated upon resistance development, and the unusual binding mode and consequently the binding energy is therefore preserved in the complex with a resistant variant. These results suggest a promising route for the design of second-generation PIs that are active against a variety of resistant PR variants.

Hydroxyethylamine isostere of an HIV-1 protease inhibitor prefers its amine to the hydroxy group in binding to catalytic aspartates. A synchrotron study of HIV-1 protease in complex with a peptidomimetic inhibitor

A complex structure of HIV-1 protease with a hydroxyethylamine-containing inhibitor Boc-Phe-Psi[(S)-CH(OH)CH2NH]-Phe-Gln-Phe-NH2 has been determined by X-ray diffraction to 1.8 A resolution. The inhibitor is bound in the active site of the protease dimer with its hydroxyethylamine isostere participating in hydrogen bonds to the catalytic aspartates 25 and 25' and glycine 27' of the active site triads via five hydrogen bonds. The isostere amine interactions with the catalytic aspartates result in a displacement of the isostere hydroxy group in comparison with the common position known for analogous hydroxyethylamine containing inhibitors. A comparison with another inhibitor of this series shows that the change of one atom of the P2' side chain (Glu/Gln) leads to an altered ability of creating hydrogen bonds to the active site and within the inhibitor molecule. The diffraction data collected at a synchrotron radiation source enabled a detailed analysis of the complex solvation and of alternative conformations of protein side chains.

The murine endogenous retrovirus MIA14 encodes an active aspartic proteinase that is functionally similar to proteinases from D-type retroviruses

Murine intracisternal A-type particles (IAPs) are endogenous retroviruses showing sequence homologies to B/D- and avian C-type retroviruses and a gene expression strategy similar to that of D-type retroviruses. These viruses form immature particles in the endoplasmic reticulum and do not release extracellular virions, but are competent for retrotransposition within the virus-producing cell. It had been assumed that lack of polyprotein processing and maturation is due to a defect in the viral proteinase (PR), but recent experiments have shown that polyprotein processing occurs when assembly of the mouse IAP MIA14 is artificially directed to the plasma membrane. We have expressed and purified recombinant MIA14 PR and show that it undergoes N- and C-terminal autoprocessing at defined sites. Using peptide cleavage and inhibition assays and in vitro cleavage of recombinant HIV-1 and MIA14 Gag polyproteins, we show that MIA14 PR is a catalytically competent enzyme comparable in its efficiency to PRs from type D exogenous retroviruses. MIA14 PR is related to the PR of Mason-Pfizer monkey virus both functionally and with respect to its expression strategy, and is distinct from HIV-1 PR with respect to substrate specificity and catalytic efficiency. These findings reveal a functional and possibly evolutionary relationship between MIA14 and D-type retroviruses and imply that a functional PR may be relevant for intracellular retrotransposition even in the case of an endogenous retrovirus that does not produce extracellular virus.

Cell-based fluorescence assay for human immunodeficiency virus type 1 protease activity

The human immunodeficiency virus type 1 (HIV-1) protease is essential for production of infectious virus and is therefore a major target for the development of drugs against AIDS. Cellular proteins are also cleaved by the protease, which explains its cytotoxic activity and the consequent failure to establish convenient cell-based protease assays. We have exploited this toxicity to develop a new protease assay that relies on transient expression of an artificial protease precursor harboring the green fluorescent protein (GFP-PR). The precursor is activated in vivo by autocatalytic cleavage, resulting in rapid elimination of protease-expressing cells. Treatment with therapeutic doses of HIV-1 protease inhibitors results in a dose-dependent accumulation of the fluorescent precursor that can be easily detected and quantified by flow cytometric and fluorimetric assays. The precursor provides a convenient and noninfectious model for high-throughput screenings of substances that can interfere with the activity of the protease in living cells.

A picomolar inhibitor of resistant strains of human immunodeficiency virus protease identified by a combinatorial approach

In order to identify inhibitors of various drug-resistant forms of the human immunodeficiency virus protease (HIV PR), we have designed and synthesized pseudopeptide libraries with a general structure Z-mimetic-Aa1-Aa2-NH2. Five different chemistries for peptide bond replacement have been employed and the resulting five individual sublibraries tested with the HIV PR and its drug-resistant mutants. Each mutant contains amino acid substitutions that have previously been shown to be associated with resistance to protease inhibitors, including Ritonavir, Indinavir, and Saquinavir. We have mapped the subsite preferences of resistant HIV PR species with the aim of selecting a pluripotent pharmaceutical lead. All of the enzyme species in this study manifest clear preference for an L-Glu residue in the P2' position. Slight, but significant, differences in P3' subsite specificity among individual resistant PR species have been documented. We have identified three compounds, combining the most favorable features of the inhibitor array, that exhibit low-nanomolar or picomolar Ki values for all three mutant PR species tested.

Systematic mutational analysis of the active-site threonine of HIV-1 proteinase: rethinking the "fireman's grip" hypothesis

Aspartic proteinases share a conserved network of hydrogen bonds (termed "fireman's grip"), which involves the hydroxyl groups of two threonine residues in the active site Asp-Thr-Gly triplets (Thr26 in the case of human immunodeficiency virus type 1 (HIV-1) PR). In the case of retroviral proteinases (PRs), which are active as symmetrical homodimers, these interactions occur at the dimer interface. For a systematic analysis of the "fireman's grip," Thr26 of HIV-1 PR was changed to either Ser, Cys, or Ala. The variant enzymes were tested for cleavage of HIV-1 derived peptide and polyprotein substrates. PR(T26S) and PR(T26C) showed similar or slightly reduced activity compared to wild-type HIV-1 PR, indicating that the sulfhydryl group of cysteine can substitute for the hydroxyl of the conserved threonine in this position. PR(T26A), which lacks the "fireman's grip" interaction, was virtually inactive and was monomeric in solution at conditions where wild-type PR exhibited a monomer-dimer equilibrium. All three mutations had little effect when introduced into only one chain of a linked dimer of HIV-1 PR. In this case, even changing both Thr residues to Ala yielded residual activity suggesting that the "fireman's grip" is not essential for activity but contributes significantly to dimer formation. Taken together, these results indicate that the "fireman's grip" is crucial for stabilization of the retroviral PR dimer and for overall stability of the enzyme.