A quantitative structure-activity relationship study for structurally diverse HIV-1 protease inhibitors: Contribution of conformational flexibility to inhibitory activity

In this study, we investigated by linear regression model the SAR data of the 15 HIV-1 protease inhibitors possessing structurally diverse scaffolds. First, a regression model was developed only using the enzyme-inhibitor interaction energy as a term of the model, but did not provide a good correlation with the inhibitory activity (R2 = 0.580 and Q2 = 0.500). Then, we focused on the conformational flexibility of the inhibitors which may represent the diversity of the inhibitors, and added two conformational parameters into the model, respectively: the number of rotatable bonds of ligands (deltaSrot) and the distortion energy of ligands (deltaElig). The regression model by adding deltaElig successfully improved the quality of the model (R2 = 0.771 and Q2 = 0.713) while the model with deltaSrot was unsuccessful. The prediction for a training inhibitor by the deltaElig model also showed good agreement with experimental activity. These results suggest that the conformational flexibility of HIV-1 protease inhibitors directly contributes to the enzyme inhibition.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Prediction of HIV-1 protease inhibitory activity of 4-hydroxy-5,6-dihydropyran-2-ones: QSAR study

Inhibition of human immunodeficiency virus 1 (HIV-1) protease is an important strategy for the treatment of HIV and acquired immune deficiency syndrome (AIDS). Therefore, HIV-1 protease inhibitory activity of dihydropyranone derivatives has been analyzed with different physico-chemical parameters. In the present work, QSAR studies were performed on a series of 4-hydroxy-5,6-dihydropyran-2-ones to explore the physico-chemical parameters responsible for their HIV-1 protease inhibitory activity. Physico-chemical parameters were calculated using WIN CAChe 6.1. Stepwise multiple linear regression analysis was performed to derive QSAR models which were further evaluated for statistical significance and predictive power by internal and external validation. The selected best QSAR model was having correlation coefficient (R)?=?0.875 and cross-validated squared correlation coefficient (Q²)?=?0.707. The developed significant QSAR model indicates that hydrophobicity of whole molecule and the substituent present at sixth position of dihydropyranones play an important role in the HIV-1 protease inhibitory activities of 4-hydroxy-5,6-dihydropyran-2-ones.     

Isolation, structure, and HIV-1-integrase inhibitory activity of structurally diverse fungal metabolites

HIV-1 integrase is a critical enzyme for replication of HIV, and its inhibition is one of the most promising new drug strategies for anti-retroviral therapy, with potentially significant advantages over existing therapies. In this report, a series of HIV-1 inhibitors isolated from the organic extract of fermentations from terrestrial fungi is described. These fungal species, belonging to a variety of genera, were collected from throughout the world following the strict guidelines of Rio Convention on Biodiversity. The polyketide- and terpenoid-derived inhibitors are represented by two naphthoquinones, a biphenyl and two triphenyls, a benzophenone, four aromatics with or without catechol units, a linear aliphatic terpenoid, a diterpenoid, and a sesterterpenoid. These compounds inhibited the coupled and strand-transfer reaction of HIV-1 integrase with an IC₅₀ value of 0.5–120 µM. The bioassay-directed isolation, structure elucidation, and HIV-1 inhibitory activity of these compounds are described.Abbreviations CH ₃ CN Acetonitrile - TFA Trifluoroacetic acid - CD ₃ CN Deuterio acetonitrile - CD ₂ Cl ₂ Deuterio methylene chloride - MEK Methyl ethyl ketone - ZnSe Zinc selenide - UV Ultraviolet - IR Infrared - NMR Nuclear magnetic resonance - COSY Homonuclear H-H correlation spectroscopy - HMQC Heteronuclear multiple-quantum coherence - HMBC Heteronuclear multiple-bond correlation - EIMS Electron-impact mass spectrometry - HREIMS High-resolution electron impact mass spectrometry - ESIMS Electrospray ionization mass spectrometry - HRFABMS High-resolution fast atom bombardment mass spectrometry - MES 2-(N-morpholino)-ethane sulfonic acid     

Structural evidence for effectiveness of darunavir and two related antiviral inhibitors against HIV-2 protease

No drug has been targeted specifically for HIV-2 (human immunodeficiency virus type 2) infection despite its increasing prevalence worldwide. The antiviral HIV-1 (human immunodeficiency virus type 1) protease (PR) inhibitor darunavir and the chemically related GRL98065 and GRL06579A were designed with the same chemical scaffold and different substituents at P2 and P2′ to optimize polar interactions for HIV-1 PR (PR1). These inhibitors are also effective antiviral agents for HIV-2-infected cells. Therefore, crystal structures of HIV-2 PR (PR2) complexes with the three inhibitors have been solved at 1.2-Å resolution to analyze the molecular basis for their antiviral potency. Unusually, the crystals were grown in imidazole and zinc acetate buffer, which formed interactions with the PR2 and the inhibitors. Overall, the structures were very similar to the corresponding inhibitor complexes of PR1 with an RMSD of 1.1 Å on main-chain atoms. Most hydrogen-bond and weaker C-H…O interactions with inhibitors were conserved in the PR2 and PR1 complexes, except for small changes in interactions with water or disordered side chains. Small differences were observed in the hydrophobic contacts for the darunavir complexes, in agreement with relative inhibition of the two PRs. These near-atomic-resolution crystal structures verify the inhibitor potency for PR1 and PR2 and will provide the basis for the development of antiviral inhibitors targeting PR2.     

A quantitative structure-activity relationship study on some HIV-1 protease inhibitors using molecular connectivity index

A quantitative structure-activity relationship (QSAR) study has been made on two different series of tetrahydropyrimidinones acting as HIV-1 protease inhibitors. A structural parameter, the first order valence molecular connectivity index ((1)chi(v)), has been used to account for the variation in the activity. The protease inhibition activity as well as the antiviral potency of the compounds are found to be significantly correlated with (1)chi(v) of P(2)/P(2') substituents attached to the two nitrogens N1 and N3, suggesting that substituents containing less electronegative and more saturated atoms, meaning thereby the less polar or more hydrophobic substituents, will be more advantageous. Further, if P(2) and P(2') are dissimilar, the former is found to be more effective than the latter. This difference is attributed to a conformational change in the enzyme that may be more favorable to P(2) binding than to P(2') binding.     

Purification of HIV-1 wild-type protease and characterization of proteolytically inactive HIV-1 protease mutants by pepstatin A affinity chromatography

Recombinant wild-type protease of human immunodeficiency virus, type [(HIV-1) expressed in E. coli was purified by pepstatin A affinity chromatography. An 88-fold purification was achieved giving a protease preparation with a specific enzymatic activity of approximately 3700 pmol/min/μg. Two proteolytically inactive HIV-1 mutant proteases (Arg-87 → Lys; Asn-88 → Glu) were found to bind to pepstatin A agarose, and they were purified as the wild-type protease. A third mutant protease (Arg-87 → Glu) was apparently unable to bind to pepstatin A under similar conditions. Binding to pepstatin A indicates the binding ability of the substrate binding site and the ability to form dimers. These features may be used to purify and to characterize other mutated HIV-1 proteases.     

High-performance liquid chromatographic assay of a new HIV-1 protease inhibitor, LB71350, in the plasma of dogs

A reliable reversed-phase high-performance liquid chromatographic method has been developed for the determination of LB71350 in the plasma of dogs. The analyte was deproteinized with 1.5 volumes of methanol and 0.5 volumes of 10% zinc sulfate, and the supernatant was injected into a 5-μm Capcell Pak C₁₈ column (150×4.6 mm I.D.). The mobile phase was a stepwise gradient mixture of acetonitrile and 0.2% triethylamine–HCl with a flow-rate of 1 ml/min and detection at UV 245 nm. The proportion of acetonitrile was kept at 52% for the first 6 min, increased to 100% for the next 0.5 min, kept at 100% for the next 2 min, decreased to 52% for the next 0.5 min, and finally kept at 52% for the next 7 min. The retention time of LB71350 was 6.9 min. The calibration was linear over the concentration range of 0.1–100 mg/l for dog plasma (r>0.997) and the limit of quantitation was 0.1 mg/l using 0.1 ml plasma. The quality control samples were reproducible with acceptable accuracy and precision at 0.1, 1, 10 and 100 mg/l concentrations. The within-day recovery (n=5) was 90.2–93.9%, the between-day recovery (n=5) was 89.5–93.5%, and the absolute between-day recovery (n=5) was 77–81%. The within-day precision (n=5) and between-day precision (n=5) were 2.59–5.82% and 3.17–4.55%, respectively. No interferences from endogenous substances were observed. Taken together, the above HPLC assay method by deproteinization and UV detection was suitable for the determination of LB71350 in the preclinical pharmacokinetics.     

Quantitative structure-activity relationship by CoMFA for cyclic urea and nonpeptide-cyclic cyanoguanidine derivatives on wild type and mutant HIV-1 protease

Abstact 3D-QSAR studies using the Comparative Molecular Field Analysis (CoMFA) methodology were conducted to predict the inhibition constants, K_i, and the inhibitor concentrations, IC₉₀ of 127 symmetrical and unsymmetrical cyclic urea and cyclic cyanoguanidine derivatives containing different substituent groups such as: benzyl, isopropyl, 4-hydroxybenzyl, ketone, oxime, pyrazole, imidazole, triazole and having anti-HIV-1 protease activities. A significant cross-validated correlation coefficient (q²) of 0.63 and a fitted correlation coefficient r² of 0.70 were obtained, indicating that the models can predict the anti-protease activity from poorly to highly active compounds reliably. The best predictions were obtained for: XV643 (predicted log 1/K_i=9.86), a 3,5-dimethoxy-benzyl cyclic urea derivate (molec60, predicted log 1/K_i=8.57) and a benzyl cyclic urea derivate (molec 61, predicted log 1/IC₉₀=6.87). Using the CoMFA method, we also predicted the biological activity of 14 cyclic urea derivatives that inhibit the HIV-1 protease mutants V82A, V82I and V82F. The predicted biological activities of the: (i) XNO63 (inhibitory activity on the mutant HIV-1 PR V82A), (ii) SB570 (inhibiting the mutant HIV-1 PR V82I) and also (iii) XV652 (during the interaction with the mutant HIV-1 PR V82F) were in good agreement with the experimental values.Figure Stereoview of the contour plots of the CoMFA steric and electrostatic fields. The favorable (indicated by blue polyhedra) and unfavorable (represented by red polyhedra) electrostatic areas and also the favorable (shown by green polyhedra) and unfavorable (shown by yellow polyhedra) steric areas formed around the most active molecule, 6a.     

Nitrogen-containing polyhydroxylated aromatics as HIV-1 integrase inhibitors: synthesis,structure-activity relationship analysis,and biological activity

Four series of forty-five nitrogen-containing polyhydroxylated aromatics based on caffeic acid phenethyl ester were designed and synthesized as HIV-1 integrase (IN) inhibitors. Most of these compounds inhibited IN catalytic activities in low micromolar range. Among these new analogues, compounds 9e and 9f were the most potent IN inhibitors with IC₅₀ value of 0.7 μM against strand transfer reaction. Their key structure-activity relationships were also discussed.     

The design,synthesis and structure-activity relationships associated with C28 amine-based betulinic acid derivatives as inhibitors of HIV-1 maturation

The design and synthesis of a series of C28 amine-based betulinic acid derivatives as HIV-1 maturation inhibitors is described. This series represents a continuation of efforts following on from previous studies of C-3 benzoic acid-substituted betulinic acid derivatives as HIV-1 maturation inhibitors (MIs) that were explored in the context of C-28 amide substituents. Compared to the C-28 amide series, the C-28 amine derivatives exhibited further improvements in HIV-1 inhibitory activity toward polymorphisms in the Gag polyprotein as well as improved activity in the presence of human serum. However, plasma exposure of basic amines following oral administration to rats was generally low, leading to a focus on moderating the basicity of the amine moiety distal from the triterpene core. The thiomorpholine dioxide (TMD) 20 emerged from this study as a compound with the optimal antiviral activity and an acceptable pharmacokinetic profile in the C-28 amine series. Compared to the C-28 amide 3, 20 offers a 2- to 4-fold improvement in potency towards the screening viruses, exhibits low shifts in the EC₅₀ values toward the V370A and ΔV370 viruses in the presence of human serum or human serum albumin, and demonstrates improved potency towards the polymorphic T371A and V362I virus variants.     

Analysis of the structure of chemically synthesized HIV-1 protease complexed with a hexapeptide inhibitor. Part I: Crystallographic refinement of 2 Å data

The structure of a complex between a hexapeptide-based inhibitor, MVT-101, and the chemically synthesized (Aba 67,95,167,195; Aba: l-α-amino-n-butyric acid) protease from the human immunodeficiency virus (HIV-1), reported previously at 2.3 Å has now been refined to a crystallographic R factor of 15.4% at 2.0 Å resolution. Root mean square deviations from ideality are 0.18 Å for bond lengths and 2.4° for the angles. The inhibitor can be fitted to the difference electron density map in two alternative orientations. Drastic differences are observed for positions and interactions at P3/S3 and P3′/S3′ subsites of the two orientations due to different crystallographic environments. © 1997 Wiley-Liss, Inc.     

Ultra-high resolution crystal structure of HIV-1 protease mutant reveals two binding sites for clinical inhibitor TMC114

TMC114 (darunavir) is a promising clinical inhibitor of HIV-1 protease (PR) for treatment of drug resistant HIV/AIDS. We report the ultra-high 0.84 A resolution crystal structure of the TMC114 complex with PR containing the drug-resistant mutation V32I (PR(V32I)), and the 1.22 A resolution structure of a complex with PR(M46L). These structures show TMC114 bound at two distinct sites, one in the active-site cavity and the second on the surface of one of the flexible flaps in the PR dimer. Remarkably, TMC114 binds at these two sites simultaneously in two diastereomers related by inversion of the sulfonamide nitrogen. Moreover, the flap site is shaped to accommodate the diastereomer with the S-enantiomeric nitrogen rather than the one with the R-enantiomeric nitrogen. The existence of the second binding site and two diastereomers suggest a mechanism for the high effectiveness of TMC114 on drug-resistant HIV and the potential design of new inhibitors.     

Biosensor analysis of the interaction of potential dimerization inhibitors with HIV-1 protease

Inhibitors of protein-protein interactions are currently considered as perspective prototypes of a new generation of drugs. The most attractive targets for such inhibitors are the oligomeric enzymes which active sites are formed by amino acid residues from different subunits. HIV-1 protease (HIVp), which is active only as a homodimer form, is the classic example of such enzymes. We have developed a new approach for experimental screening of HIVp dimerization inhibitors. It is based on an original biosensor test-system for differential analysis of interaction of tested substances with HIVp dimers and monomers. Using this test-system we have analyzed the most perspective candidate substances predicted by the method of virtual screening, and also some derivatives of glycyrrhizin, triterpenic and steroid glycosides. In the results of this study we have found one compound, which preferentially interacts with HIVp monomers and inhibits in vitro activity of this enzyme with the IC₅₀ value of about 10^?6 M.     

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.     

Structural insights into the South African HIV-1 subtype C protease: impact of hinge region dynamics and flap flexibility in drug resistance

The HIV protease plays a major role in the life cycle of the virus and has long been a target in antiviral therapy. Resistance of HIV protease to protease inhibitors (PIs) is problematic for the effective treatment of HIV infection. The South African HIV-1 subtype C protease (C-SA PR), which contains eight polymorphisms relative to the consensus HIV-1 subtype B protease, was expressed in Escherichia coli, purified, and crystallized. The crystal structure of the C-SA PR was resolved at 2.7?Å, which is the first crystal structure of a HIV-1 subtype C protease that predominates in Africa. Structural analyses of the C-SA PR in comparison to HIV-1 subtype B proteases indicated that polymorphisms at position 36 of the homodimeric HIV-1 protease may impact on the stability of the hinge region of the protease, and hence the dynamics of the flap region. Molecular dynamics simulations showed that the flap region of the C-SA PR displays a wider range of movements over time as compared to the subtype B proteases. Reduced stability in the hinge region resulting from the absent E35-R57 salt bridge in the C-SA PR, most likely contributes to the increased flexibility of the flaps which may be associated with reduced susceptibility to PIs.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:36     

Design,synthesis and biological evaluation of novel HIV-1 protease inhibitors with pentacyclic triterpenoids as P2-ligands

The design, synthesis and SAR study of a new series of HIV-1 protease inhibitors with pentacyclic triterpenoids as P2 ligands and phenylsulfonamide as P2′ ligands were discussed. These compounds exhibited micromolar inhibitory potency, among which compound T1c displayed HIV-1 protease inhibition with IC₅₀ values of 0.12?μM, which was 67 times the inhibitory activity of its raw material Ursolic acid (8.0?μM).     

The higher barrier of darunavir and tipranavir resistance for HIV-1 protease

Darunavir and tipranavir are two inhibitors that are active against multi-drug resistant (MDR) HIV-1 protease variants. In this study, the invitro inhibitory efficacy was tested against a MDR HIV-1 protease variant, MDR 769 82T, containing the drug resistance mutations of 46L/54V/82T/84V/90M. Crystallographic and enzymatic studies were performed to examine the mechanism of resistance and the relative maintenance of potency. The key findings are as follows: (i) The MDR protease exhibits decreased susceptibility to all nine HIV-1 protease inhibitors approved by the US Food and Drug Administration (FDA), among which darunavir and tipranavir are the most potent; (ii) the threonine 82 mutation on the protease greatly enhances drug resistance by altering the hydrophobicity of the binding pocket; (iii) darunavir or tipranavir binding facilitates closure of the wide-open flaps of the MDR protease; and (iv) the remaining potency of tipranavir may be preserved by stabilizing the flaps in the inhibitor-protease complex while darunavir maintains its potency by preserving protein main chain hydrogen bonds with the flexible P2 group. These results could provide new insights into drug design strategies to overcome multi-drug resistance of HIV-1 protease variants.     

1H NMR conformational study on n-terminal nonapeptide sequences of HIV-1 Tat protein: a contribution to structure–activity relationships

On the basis of our recent results, the N-terminal sequence of HIV-1 Tat protein as a natural competitive inhibitor of dipeptidyl peptidase IV (DP IV) is supposed to interact directly with the active site of DP IV hence mediating its immunosuppressive effects via specific DP IV interactions. Of special interest is the finding that amino acid substitutions of the Tat(1–9) peptide (MDPVDPNIE) in position 5 with S-isoleucine and in position 6 with S-leucine led to peptides with strongly reduced inhibitory activity suggesting differences in the solution conformation of the three analogues. Therefore, ¹H NMR techniques in conjunction with molecular modelling have been used here to determine the solution structure of Tat(1–9), I⁵-Tat(1–9) and L⁶-Tat(1–9) and to examine the influence of amino acid exchanges on structural features of these peptides. The defined structures revealed differences in the conformations what might be the reason for different interactions of these Tat(1–9) analogues with certain amino acids of the active site of DP IV. © 1998 European Peptide Society and John Wiley & Sons, Ltd.