Cyclic HIV protease inhibitors: synthesis, conformational analysis, P2/P2' structure-activity relationship, and molecular recognition of cyclic ureas

New developments in anti-HIV chemotherapy

Virtually all the compounds that are currently used, or under advanced clinical trial, for the treatment of HIV infections, belong to one of the following classes: (i) nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs): i.e. zidovudine, didanosine, zalcitabine, stavudine, lamivudine, abacavir, emtricitabine, tenofovir (PMPA) disoproxil fumarate; (ii) non-nucleoside reverse transcriptase inhibitors (NNRTIs): i.e. nevirapine, delavirdine, efavirenz, emivirine; and (iii) protease inhibitors (PIs): i.e. saquinavir, ritonavir, indinavir, nelfinavir and amprenavir. In addition, various other events in the HIV replicative cycle are potential targets for chemotherapeutic intervention: (i) viral adsorption, through binding to the viral envelope glycoprotein gp120; (ii) viral entry, through blockade of the viral coreceptors CXCR4 and CCR5; (iii) virus-cell fusion; (iv) viral assembly and disassembly; (v) proviral DNA integration; (vi) viral mRNA transcription. Also, new NRTIs, NNRTIs and PIs have been developed that possess respectively improved metabolic characteristics, or increased activity against NNRTI-resistant HIV strains or, as in the case of PIs, a different, non-peptidic scaffold. Given the multitude of molecular targets with which anti-HIV agents can interact, one should be cautious in extrapolating from cell-free enzymatic assays to the mode of action of these agents in intact cells.

Interpreting trends in the binding of cyclic ureas to HIV-1 protease

The design of new HIV protease inhibitors requires an improved understanding of the physical basis of inhibitor/protein binding. Here, the binding affinities of seven aliphatic cyclic ureas to HIV-1 protease are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed U-shaped trend of binding free energy as a function of aliphatic chain length. Interestingly, the decrease in affinity for the longest chains is attributable primarily to the energy cost of partly desolvating charged aspartic and arginine groups at the mouths of the active site. Even aliphatic chains too short to contact these charged groups directly are subject to considerable desolvation penalties. We are not aware of other systems where binding affinity trends have been attributed to long-ranged electrostatic desolvation of ionized groups. A generalized Born/surface area solvation model yields a much smaller change in desolvation energy with chain length and, therefore, does not reproduce the experimental binding affinity trends. This result suggests that the generalized Born model should be used with caution for complex, partly desolvated systems like protein binding sites. We also find that changing the assumed protonation state of the active site aspartyl dyad significantly affects the computed binding affinity trends. The protonation state of the aspartyl dyad in the presence of cyclic ureas is discussed in light of the observation that the monoprotonated state reproduces the experimental results best.

Inhibition of human endogenous retrovirus-K10 protease in cell-free and cell-based assays

A full-length and C-terminally truncated version of human endogenous retrovirus (HERV)-K10 protease were expressed in Escherichia coli and purified to homogeneity. Both versions of the protease efficiently processed HERV-K10 Gag polyprotein substrate. HERV-K10 Gag was also cleaved by human immunodeficiency virus, type 1 (HIV-1) protease, although at different sites. To identify compounds that could inhibit protein processing dependent on the HERV-K10 protease, a series of cyclic ureas that had previously been shown to inhibit HIV-1 protease was tested. Several symmetric bisamides acted as very potent inhibitors of both the truncated and full-length form of HERV-K10 protease, in subnanomolar or nanomolar range, respectively. One of the cyclic ureas, SD146, can inhibit the processing of in vitro translated HERV-K10 Gag polyprotein substrate by HERV-K10 protease. In addition, in virus-like particles isolated from the teratocarcinoma cell line NCCIT, there is significant accumulation of Gag and Gag-Pol precursors upon treatment with SD146, suggesting the compound efficiently blocks HERV-K Gag processing in cells. This is the first report of an inhibitor able to block cell-associated processing of Gag polypeptides of an endogenous retrovirus.

Novel compounds in preclinical/early clinical development for the treatment of HIV infections

Virtually all the compounds that are currently used, or under advanced clinical trial, for the treatment of HIV infections, belong to one of the following classes: (i) nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs), (ii) non-nucleoside reverse transcriptase inhibitors (NNRTIs) and (iii) protease inhibitors (PIs). In addition to the reverse transcriptase and protease step, various other events in the HIV replicative cycle are potential targets for chemotherapeutic intervention: (i) viral adsorption, through binding to the viral envelope glycoprotein gp120 (polysulphates, polysulphonates, polyoxometalates, zintevir, negatively charged albumins); (ii) viral entry, through blockade of the viral coreceptors CXCR4 and CCR5 [bicyclams (AMD3100), polyphemusins (T22), TAK-779]; (iii) virus-cell fusion, through binding to the viral glycoprotein gp41 [T-20 (DP-178), siamycins, betulinic acid derivatives]; (iv) viral assembly and disassembly, through NCp7 zinc finger-targeted agents [2,2'-dithiobisbenzamides (DIBAs), azadicarbonamide (ADA)]; (v) proviral DNA integration, through integrase inhibitors such as L-chicoric acid; (vi) viral mRNA transcription, through inhibitors of the transcription (transactivation) process (peptoid CGP64222, fluoroquinolone K-12, Streptomyces product EM2487). Also, in recent years new NRTIs, NNRTIs and PIs have been developed that possess, respectively, improved metabolic characteristics (i.e. phosphoramidate and cyclosaligenyl pronucleotides of d4T), or increased activity against NNRTI-resistant HIV strains, or, in the case of PIs, a different, non-peptidic scaffold. Given the multitude of molecular targets with which anti-HIV agents can interact, one should be cautious in extrapolating from cell-free enzymatic assays to the mode of action of these agents in intact cells. A number of compounds (i.e. zintevir and L-chicoric acid, on the one hand; and CGP64222 on the other hand) have recently been found to interact with virus-cell binding and viral entry in contrast to their proposed modes of action targeted at the integrase and transactivation process, respectively.

Molecular and biological constraints on ligand-binding affinity and specificity

Interactions between biological macromolecules have characteristic values of affinity and specificity that are set according to the biological function that is served by the interaction in the organism. Here we examine the molecular mechanisms that are used to achieve the required values of affinity and specificity in various biological systems.

Computational design of new cyclic urea inhibitors for improved binding of HIV-1 aspartic protease

We report in this paper the design, by means of computational techniques, of new cyclic urea inhibitors of the HIV aspartic protease. The relationship between the complexation energies of the enzyme with known inhibitors and the experimentally determined log K(i) have been studied and used to predict inhibition constants for new inhibitors.

7-Azabicycloheptane carboxylic acid: a proline replacement in a boroarginine thrombin inhibitor

[formula: see text] The synthesis of thrombin inhibitor 3, which incorporates conformationally constrained 7-azabicycloheptane carboxylic acid (1) as a proline replacement, is described. The inhibition constant (Ki(thrombin) = 2.9 nM) indicates that 1 is a reasonable replacement of proline in the formation of a beta-turn tripeptide mimetic.

Quantitative structure-activity relationship studies on cyclic cyanoguanidines acting as HIV-1 protease inhibitors

A quantitative structure-activity relationship study has been performed on some cyclic cyanoguanidines that inhibit the enzyme HIV-1 protease (HIV-1-PR) and exhibit antiviral potency, and the results have been compared with those of cyclic urea derivatives. Both the enzyme inhibition activity and antiviral potency in cyclic cyanoguanidines as well as in cyclic urea derivatives are found to be primarily governed by hydrophobic property of substituents attached to nitrogen (P2/P2') and further enhanced by OH or NH2 group, if any, present in the substituents. However, aromatic substituents are found to be unfavourable to both the activities of cyclic cyanoguanidines but not to any activity of cyclic urea derivatives. Cyclic urea derivatives are indicated to be more potent than cyclic cyanoguanidines. A model for the interaction of cyclic cyanoguanidines with the receptor is proposed.

Flap opening and dimer-interface flexibility in the free and inhibitor-bound HIV protease, and their implications for function

BACKGROUND

(1)H and (15)N transverse relaxation measurements on perdeuterated proteins are ideally suited for detecting backbone conformational fluctuations on the millisecond-microsecond timescale. The identification of conformational exchange on this timescale by measuring the relaxation of both (1)H and (15)N holds great promise for the elucidation of functionally relevant conformational changes in proteins.

RESULTS

We measured the transverse (1)H and (15)N relaxation rates of backbone amides of HIV-1 protease in its free and inhibitor-bound forms. An analysis of these rates, obtained as a function of the effective rotating frame field, provided information about the timescale of structural fluctuations in several regions of the protein. The flaps that cover the active site of the inhibitor-bound protein undergo significant changes of backbone (φ,psi) angles, on the 100 micros timescale, in the free protein. In addition, the intermonomer beta-sheet interface of the bound form, which from protease structure studies appears to be rigid, was found to fluctuate on the millisecond timescale.

CONCLUSIONS

We present a working model of the flap-opening mechanism in free HIV-1 protease which involves a transition from a semi-open to an open conformation that is facilitated by interaction of the Phe53 ring with the substrate. We also identify a surprising fluctuation of the beta-sheet intermonomer interface that suggests a structural requirement for maturation of the protease. Thus, slow conformational fluctuations identified by (1)H and (15)N transverse relaxation measurements can be related to the biological functions of proteins.

Probing the role of structural water in a duplex oligodeoxyribonucleotide containing a water-mimicking base analog

Molecular dynamics simulations were performed on models of the dodecamer DNA double-stranded segment, [d(CGCGAATTCGCG)](2), in which each of the adenine residues, individually or jointly, was replaced by the water-mimicking analog 2'-deoxy-7-(hydroxy-methyl)-7-deazaadenosine (hm(7)c(7)dA) [Rockhill, J.K., Wilson,S.R. and Gumport,R.I. (1996) J. Am. Chem. Soc.,118, 10065-10068]. The simulations, when compared with those of the dodecamer itself, show that incorporation of the analog affects neither the overall DNA structure nor its hydrogen-bonding and stacking interactions when it replaces a single individual base. Furthermore, the water molecules near the bases in the singly-substituted oligonucleotides are similarly unaffected. Double substitutions lead to differences in all the aforementioned parameters with respect to the reference sequence. The results suggest that the analog provides a good mimic of specific 'ordered' water molecules observed in contact with DNA itself and at the interface between protein and DNA in specific complexes.

The maximal affinity of ligands

We explore the question of what are the best ligands for macromolecular targets. A survey of experimental data on a large number of the strongest-binding ligands indicates that the free energy of binding increases with the number of nonhydrogen atoms with an initial slope of approximately -1.5 kcal/mol (1 cal = 4.18 J) per atom. For ligands that contain more than 15 nonhydrogen atoms, the free energy of binding increases very little with relative molecular mass. This nonlinearity is largely ascribed to nonthermodynamic factors. An analysis of the dominant interactions suggests that van der Waals interactions and hydrophobic effects provide a reasonable basis for understanding binding affinities across the entire set of ligands. Interesting outliers that bind unusually strongly on a per atom basis include metal ions, covalently attached ligands, and a few well known complexes such as biotin-avidin.

Structure-based design of diaminopyranosides as a novel inhibitor core unit of HIV proteases

Novel HIV PR inhibitors, which contain a diaminopyranoside moiety as an inhibitor core unit, were designed based on the 3D structures of complexes of HIV-1 PR with transition-state mimics. These compounds were examined for their ability to inhibit the hydrolytic activity of a recombinant HIV-1 PR.

Thermodynamic linkage between the binding of protons and inhibitors to HIV-1 protease

The aspartyl dyad of free HIV-1 protease has apparent pK(a)s of approximately 3 and approximately 6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson-Boltzmann model of electrostatics is used to compute the apparent pK(a)s of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns epsilon = 4 to the solute interior and uses a detailed model of ionization; the other uses epsilon = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pK(a)s of approximately 3 and approximately 6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pK(a)s in complex protein environments. Because binding of a cyclic urea inhibitor shifts the pK(a)s of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.

Tricyclic ureas: a new class of HIV-1 protease inhibitors

A new class of tricyclic ureas containing a conformationally constrained proline was designed with the aid of molecular modeling. Efficient stereoselective intermolecular pinacol coupling represented the highlight of the synthesis. These rigid cyclic ureas are active towards HIV-1 protease, with 9 being the most potent compound (Ki = 9 nM) despite interacting with only three side chain binding pockets of HIV protease.

Design and selection of DMP 850 and DMP 851: the next generation of cyclic urea HIV protease inhibitors

BACKGROUND

Recent clinical trials have demonstrated that HIV protease inhibitors are useful in the treatment of AIDS. It is necessary, however, to use HIV protease inhibitors in combination with other antiviral agents to inhibit the development of resistance. The daunting ability of the virus to rapidly generate resistant mutants suggests that there is an ongoing need for new HIV protease inhibitors with superior pharmacokinetic and efficacy profiles. In our attempts to design and select improved cyclic urea HIV protease inhibitors, we have simultaneously optimized potency, resistance profile, protein binding and oral bioavailability.

RESULTS

We have discovered that nonsymmetrical cyclic ureas containing a 3-aminoindazole P2 group are potent inhibitors of HIV protease with excellent oral bioavailability. Furthermore, the 3-aminoindazole group forms four hydrogen bonds with the enzyme and imparts a good resistance profile. The nonsymmetrical 3-aminoindazoles DMP 850 and DMP 851 were selected as our next generation of cyclic urea HIV protease inhibitors because they achieve 8 h trough blood levels in dog, with a 10 mg/kg dose, at or above the protein-binding-adjusted IC90 value for the worst single mutant--that containing the Ile84-->Val mutation.

CONCLUSIONS

In selecting our next generation of cyclic urea HIV protease inhibitors, we established a rigorous set of criteria designed to maximize chances for a sustained antiviral effect in HIV-infected individuals. As DMP 850 and DMP 851 provide plasma levels of free drug that are sufficient to inhibit wild-type HIV and several mutant forms of HIV, they could show improved ability to decrease viral load for clinically significant time periods. The ultimate success of DMP 850 and DMP 851 in clinical trials might depend on achieving or exceeding the oral bioavailability seen in dog.

Quantitative structure-activity relationship studies on cyclic urea-based HIV protease inhibitors

A quantitative structure-activity relationship (QSAR) study is described on some cyclic ureas that inhibit the enzyme HIV-1 protease (HIV-1-PR) and exhibit antiviral potency. Both the enzyme inhibition activity and the antiviral potency were found to be primarily governed by the hydrophobic property of the substituents at the nitrogens (N2/N2') of the urea. Adjacent to the nitrogens, the C1/C1'-substituents are, however, found to affect the activity (inhibition) by their molecular size. The essential binding of the ureas with the receptor is, however, through multiple hydrogen bonding, where the substituents, too, can participate in such binding if they are capable of doing so. A schematic diagram of the overall interaction of the inhibitors with the receptor is presented.

Inhibitors of HIV-1 protease: a major success of structure-assisted drug design

Retroviral protease (PR) from the human immunodeficiency virus type 1 (HIV-1) was identified over a decade ago as a potential target for structure-based drug design. This effort was very successful. Four drugs are already approved, and others are undergoing clinical trials. The techniques utilized in this remarkable example of structure-assisted drug design included crystallography, NMR, computational studies, and advanced chemical synthesis. The development of these drugs is discussed in detail. Other approaches to designing HIV-1 PR inhibitors, based on the concepts of symmetry and on the replacement of a water molecule that had been found tetrahedrally coordinated between the enzyme and the inhibitors, are also discussed. The emergence of drug-induced mutations of HIV-1 PR leads to rapid loss of potency of the existing drugs and to the need to continue the development process. The structural basis of drug resistance and the ways of overcoming this phenomenon are mentioned.

Cyclic HIV protease inhibitors: design and synthesis of orally bioavailable, pyrazole P2/P2' cyclic ureas with improved potency

Highly potent HIV-1 protease (HIVPR) inhibitors have been designed and synthesized by introducing bidentate hydrogen-bonding oxime and pyrazole groups at the meta-position of the phenyl ring on the P2/P2' substituents of cyclic ureas. Nonsymmetrical cyclic ureas incorporating 3(1H)-pyrazolylbenzyl as P2 and hydrophilic functionalities as P2' show potent protease inhibition and antiviral activities against HIV and have good oral bioavailabilities. The X-ray structure of HIVPR.10A complex confirms that the two pyrazole rings of 10A form bidentate hydrogen bonds with the side-chain oxygen (C=O) and backbone nitrogen (N-H) of Asp30/30' of HIVPR.

Molecular recognition of cyclic urea HIV-1 protease inhibitors

As long as the threat of human immunodeficiency virus (HIV) protease drug resistance still exists, there will be a need for more potent antiretroviral agents. We have therefore determined the crystal structures of HIV-1 protease in complex with six cyclic urea inhibitors: XK216, XK263, DMP323, DMP450, XV638, and SD146, in an attempt to identify 1) the key interactions responsible for their high potency and 2) new interactions that might improve their therapeutic benefit. The structures reveal that the preorganized, C2 symmetric scaffolds of the inhibitors are anchored in the active site of the protease by six hydrogen bonds and that their P1 and P2 substituents participate in extensive van der Waals interactions and hydrogen bonds. Because all of our inhibitors possess benzyl groups at P1 and P1', their relative binding affinities are modulated by the extent of their P2 interactions, e.g. XK216, the least potent inhibitor (Ki (inhibition constant) = 4.70 nM), possesses the smallest P2 and the lowest number of P2-S2 interactions; whereas SD146, the most potent inhibitor (Ki = 0.02 nM), contains a benzimidazolylbenzamide at P2 and participates in fourteen hydrogen bonds and approximately 200 van der Waals interactions. This analysis identifies the strongest interactions between the protease and the inhibitors, suggests ways to improve potency by building into the S2 subsite, and reveals how conformational changes and unique features of the viral protease increase the binding affinity of HIV protease inhibitors.

The synthesis of symmetrical and unsymmetrical P1/P1' cyclic ureas as HIV protease inhibitors

Cyclic urea SD146, a potent HIV protease inhibitor bearing a flat resistance profile, possessed poor solubility and bioavailability, which precluded further development of the compound. In an effort to improve upon the pharmacokinetic profile of the compound, several analogs modified at the P1/P1' residues were prepared and evaluated. Several of those compounds displayed significant improvement of physical properties.