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.     

An immunoenzymatic solid-phase assay for quantitative determination of HIV-1 protease activity

A novel immunoenzymatic procedure for the quantitative determination of HIV protease activity is provided. An N-terminal biotinylated peptide (DU1) that comprises an HIV-1 protease (HIV-PR) cleavage sequence was bound to streptavidin-coated microtiter plates. The bound peptide can be quantified by an immunoenzymatic procedure (enzyme-linked immunosorbent assay, ELISA) that includes a monoclonal antibody (Mab 332) against the peptide (DU1) C-terminal. The incubation of the bound peptide with HIV-PR in solution resulted in a signal decrement, as the peptide was hydrolyzed and the released C-terminal segment washed away. An equation that relates the amount of added enzyme to the kinetics of the reaction was written in order to describe this heterogeneous enzyme-quasi-saturable system. This equation allows quantitative determination of protease activity, a feature widely underrated in previous similar assays. The assay also allows evaluation of the inhibitory activity of HIV-PR inhibitors. Due to the intrinsic advantages of the ELISA format, this method could be used in high-throughput screening of HIV protease inhibitors. The assay can be extended to other proteolytic enzymes.     

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.     

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.     

Synthesis and activity of N-acyl azacyclic urea HIV-1 protease inhibitors with high potency against multiple drug resistant viral strains

As part of our efforts to identify potent HIV-1 protease inhibitors that are active against resistant viral strains, structural modification of the azacyclic urea (I) was undertaken by incorporating acyl groups as P₁′ ligands. The extensive SAR study has yielded a series of N-acyl azacyclic ureas (II), which are highly potent against both wild-type and multiple PI-resistant viral strains.     

结合分子相似性、药效团和分子对接筛选新的HIV-1蛋白酶抑制剂

结合分子相似性、药效团和分子对接建立兼顾计算效率和预测准确度的HIV-1蛋白酶抑制剂筛选方法。首先通过对现有HIV-1蛋白酶抑制剂分子进行相似性分析,选取代表性的HIV-1蛋白酶抑制剂作为模板分子,构建和优化药效团模型,并从1万个化合物中优先筛选出500个化合物。而后采用分子对接方法进一步考察化合物与HIV-1蛋白酶结合情况,得到4个新的活性候选化合物,并进行其结合自由能计算和抗突变性分析。结果表明新候选化合物ST025723和HIV-1蛋白酶表现出较好的相互作用和抗突变性,具有深入研究的价值,同时也证明分子相似性、药效团和分子对接相结合能够快速有效地发现新颖活性候选化合物。     

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     

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.     

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.     

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.     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.     

Drug resistance in HIV-1 protease: Flexibility-assisted mechanism of compensatory mutations

The emergence of drug-resistant variants is a serious side effect associated with acquired immune deficiency syndrome therapies based on inhibition of human immunodeficiency virus type 1 protease (HIV-1 PR). In these variants, compensatory mutations, usually located far from the active site, are able to affect the enzymatic activity via molecular mechanisms that have been related to differences in the conformational flexibility, although the detailed mechanistic aspects have not been clarified so far. Here, we perform multinanosecond molecular dynamics simulations on L63P HIV-1 PR, corresponding to the wild type, and one of its most frequently occurring compensatory mutations, M46I, complexed with the substrate and an enzymatic intermediate. The quality of the calculations is established by comparison with the available nuclear magnetic resonance data. Our calculations indicate that the dynamical fluctuations of the mutated enzyme differ from those in the wild type. These differences in the dynamic properties of the adducts with the substrate and with the gem-diol intermediate might be directly related to variations in the enzymatic activity and therefore offer an explanation of the observed changes in catalytic rate between wild type and mutated enzyme. We anticipate that this "flexibility-assisted" mechanism might be effective in the vast majority of compensatory mutations, which do not change the electrostatic properties of the enzyme.     

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.     

Studies on flexibility and binding affinity of Asp25 of HIV-1 protease mutants

We have investigated and highlighted the behavior of binding residue, Asp25 by computational analysis, which play an important role in understanding docking process with drug molecule, Ritonavir (Norvir^®) and the flexibility nature of the Human Immunodeficiency Virus-1 (HIV-1) protease enzyme. It is well known that Ritonavir is a potent and a selective HIV-1 protease inhibitor. Molecular dockings were performed in order to gain insights regarding the binding mode of this inhibitor. In our analysis, we observed Ritonavir had different rank orders of scores against different mutant of this enzyme. Asp25 of the enzyme was found to be the active site for all the mutants. The results clearly suggest that Ritonavir is not able to appropriately bind at the active site of each HIV-1 protease mutant due to RMSD difference of the amino acid (Asp) at the position 25 of all mutants. These findings support the concept that 3D space of active site is a qualitative assessment for binding affinity of inhibitor with an enzyme. The investigation on the flexibility nature of Asp25 by normal mode analysis, show that binding residue posses less flexibility due to its solvation potential. The overall analysis of our study brings clarity to the binding behavior with respect to the different mutants with Ritonavir on the basis RMSD and also on the flexible nature of HIV-1 protease enzyme with respect to Asp25 position.     

Solvent accessibility as a predictive tool for the free energy of inhibitor binding to the HIV-1 protease.

We have developed a simple approach for the evaluation of the free energies of inhibitor binding to the protease of the human immunodeficiency virus (HIV-1 PR). Our algorithm is based on the observation that most groups that line the binding pockets of this enzyme are hydrophobic in nature. Based on this fact, we have likened the binding of an inhibitor to this enzyme to its transfer from water to a medium of lower polarity. The resulting expression produced values for the free energy of binding of inhibitors to the HIV-1 PR that are in good agreement with experimental values. The additive nature of this approach has enabled us to partition the free energy of binding into the contributions of single fragments. The resulting analysis clearly indicates the existence of a ranking in the participation of the enzyme's subsites in binding. Although all the enzyme's pockets contribute to binding, the ones that bind the P2-P'2 span of the inhibitor are in general the most critical for high inhibitor potency. Moreover, our method has allowed us to determine the nature of the functional groups that fit into given enzyme binding pockets. Perusal of the energy contributions of single side chains has shown that a large number of hydrophobic and aromatic groups located in the central portion of the HIV-1 PR inhibitors present optimal binding. All of these observations are in agreement with experimental evidence, providing a validation for the physical relevancy of our model.     

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.     

Structure-activity relationship study of tetrapeptide inhibitors of the Vascular Endothelial Growth Factor A binding to Neuropilin-1

Neuropilin-1 is considered as one of the key receptors responsible for signaling pathways involved in pathological angiogenesis necessary for tumor progression, therefore targeting of VEGF₁₆₅ binding to NRP-1 could be a relevant strategy for antiangiogenic treatment. It was shown before that the VEGF₁₆₅/NRP-1 interaction can be inhibited by short tetrapeptides with K/RXXR sequence.Here, we present a structure–activity relationship study of the systematic optimization of amino acid residues in positions 1–3 in the above tetrapeptides. All the 13 synthesized analogs possessed C-terminal arginine that is a necessary element for interaction with NRP-1. The obtained results of the inhibitory activity and modeling by molecular dynamics indicate that simultaneous interactions of the basic amino acid residues in position 1 and 4 (Arg) with Neuropilin-1 are crucial and their cooperation strongly affects the inhibitory activity. In addition, the binding strength is modulated by the flexibility of the peptide backbone (in the central part of the peptide), and the nature of the side chain of the amino acids at the second or third position. A dramatic decrease in the activity to the receptor is observed in flexible derivatives that are missing proline residues. The results described in this paper should prove useful for future studies aimed at establishing the best pharmacophore for inhibitors of VEGF₁₆₅ binding to NRP-1.