Active Site Of HIV-1 Protease Research Articles

Evaluation of darunavir-derived HIV-1 protease inhibitors incorporating P2′ amide-derivatives: Synthesis, biological evaluation and structural studies

We report here the synthesis and biological evaluation of darunavir derived HIV-1 protease inhibitors and their functional effect on enzyme inhibition and antiviral activity in MT-2 cell lines. The P2′ 4-amino functionality was modified to make a number of amide derivatives to interact with residues in the S2′ subsite of the HIV-1 protease active site. Several compounds exhibited picomolar enzyme inhibitory and low nanomolar antiviral activity. The X-ray crystal structure of the chloroacetate derivative bound to HIV-1 protease was determined. Interestingly, the active chloroacetate group converted to the acetate functionality during X-ray exposure. The structure revealed that the P2′ carboxamide functionality makes enhanced hydrogen bonding interactions with the backbone atoms in the S2′-subsite.

Beyond darunavir: recent development of next generation HIV-1 protease inhibitors to combat drug resistance.

We report our recent development of a conceptually new generation of exceptionally potent non-peptidic HIV-1 protease inhibitors that displayed excellent pharmacological and drug-resistance profiles. Our X-ray structural studies of darunavir and other designed inhibitors from our laboratories led us to create a variety of inhibitors incorporating fused ring polycyclic ethers and aromatic heterocycles to promote hydrogen bonding interactions with the backbone atoms of HIV-1 protease as well as van der Waals interactions with residues in the S2 and S2' subsites. We have also incorporated specific functionalities to enhance van der Waals interactions in the S1 and S1' subsites. The combined effects of these structural templates are critical to the inhibitors' exceptional potency and drug-like properties. We highlight here our molecular design strategies to promote backbone hydrogen bonding interactions to combat drug-resistance and specific design of polycyclic ether templates to mimic peptide-like bonds in the HIV-1 protease active site. Our medicinal chemistry and drug development efforts led to the development of new generation inhibitors significantly improved over darunavir and displaying unprecedented antiviral activity against multidrug-resistant HIV-1 variants.

HIV-1 Gag Non-Cleavage Site PI Resistance Mutations Stabilize Protease/Gag Substrate Complexes In Silico via a Substrate-Clamp

HIV-1 protease active site inhibitors are a key part of antiretroviral therapy, though resistance can evolve rendering therapy ineffective. Protease inhibitor resistance typically starts with primary mutations around the active site, which reduces inhibitor binding, protease affinity for substrate cleavage site residues P4-P4′, and viral replication. This is often followed by secondary mutations in the protease substrate-grooves which restore viral replication by increasing protease affinity for cleavage site residues P12-P5/P5′-P12′, while maintaining resistance. However, mutations in Gag alone can also result in resistance. The Gag resistance mutations can occur in cleavage sites (P12-P12′) to increase PR binding, as well as at non-cleavage sites. Here we show in silico that Gag non-cleavage site protease inhibitor resistance mutations can stabilize protease binding to Gag cleavage sites which contain structured subdomains on both sides: SP1/NC, SP2/p6, and MA/CA. The Gag non-cleavage site resistance mutations coordinated a network of H-bond interactions between the adjacent structured subdomains of the Gag substrates to form a substrate-clamp around the protease bound to cleavage site residues P12-P12′. The substrate-clamp likely slows protease disassociation from the substrate, restoring the cleavage rate in the presence of the inhibitor. Native Gag substrates can also form somewhat weaker substrate-clamps. This explains the 350-fold slower cleavage rate for the Gag CA/SP1 cleavage site in that the CA-SP1 substrate lacks structured subdomains on both sides of the cleavage site, and so cannot form a substrate-clamp around the PR.

Potent HIV-1 protease inhibitors incorporating squaramide-derived P2 ligands: Design, synthesis, and biological evaluation.

We describe the design, synthesis, and biological evaluation of novel HIV-1 protease inhibitors containing a squaramide-derived scaffold as the P2 ligand in combination with a (R)-hydroxyethylamine sulfonamide isostere. Inhibitor 3h with an N-methyl-3-(R)-aminotetrahydrofuranyl squaramide P2-ligand displayed an HIV-1 protease inhibitory Ki value of 0.51 nM. An energy minimized model of 3h revealed the major molecular interactions between HIV-1 protease active site and the tetrahydrofuranyl squaramide scaffold that may be responsible for its potent activity.

Design, synthesis, and X-ray studies of potent HIV-1 protease inhibitors incorporating aminothiochromane and aminotetrahydronaphthalene carboxamide derivatives as the P2 ligands

We describe the design, synthesis, and biological evaluation of a series of novel HIV-1 protease inhibitors with carboxamide derivatives as the P2 ligands. We have specifically designed aminothiochromane and aminotetrahydronaphthalene-based carboxamide ligands to promote hydrogen bonding and van der Waals interactions in the active site of HIV-1 protease. Inhibitors 4e and 4j have shown potent enzyme inhibitory and antiviral activity. High resolution X-ray crystal structures of 4d- and 4k-bound HIV-1 protease revealed molecular insights into the ligand-binding site interactions.

The ANANKE relative energy gradient (REG) method to automate IQA analysis over configurational change

Much chemical insight ultimately comes down to finding out which fragment of a total system behaves like the total system, in terms of an energy profile. A simple example is that of the water dimer, where this system is regarded as held together by a hydrogen bond. The hydrogen bond consists of two atoms (H···O), which energetically behave similarly to the total system (H2O)2. However, from a quantum mechanical point of view, each atom in the total system interacts with any other atom. Thus, the view that the hydrogen bond by itself governs the energetic stability of the water dimer needs rigorous justification. In this work, we propose a method that provides such a justification, in general, but only illustrated on the water dimer here. This method is based on the topological energy partitioning method called interacting quantum atoms (IQA). The method is implemented in the program ANANKE, which calculates correlations between the energy profile of the total system and those of subsystems (or fragments). ANANKE acts on the IQA energy contributions obtained for a sequence of full-system geometries controlled by a coordinate of interest (e.g. the O···H distance in the water dimer). Although applied only for the water dimer in this work, the method is general and able to explain the gauche effect, the torsional barrier in biphenyl, the arrow-pushing scheme of an enzymatic reaction (peptide hydrolysis in the HIV-1 Protease active site), and halogen-alkane nucleophilic substitution (SN2) reactions. Those applications will appear elsewhere as separate and elaborated case studies; here we focus on the details of the ANANKE method and its justification, using the water dimer as a concrete case.

Effects of drug-resistant mutations on the dynamic properties of HIV-1 protease and inhibition by Amprenavir and Darunavir.

Molecular dynamics simulations are performed to investigate the dynamic properties of wild-type HIV-1 protease and its two multi-drug-resistant variants (Flap + (L10I/G48V/I54V/V82A) and Act (V82T/I84V)) as well as their binding with APV and DRV inhibitors. The hydrophobic interactions between flap and 80 s (80’s) loop residues (mainly I50-I84’ and I50’-I84) play an important role in maintaining the closed conformation of HIV-1 protease. The double mutation in Act variant weakens the hydrophobic interactions, leading to the transition from closed to semi-open conformation of apo Act. APV or DRV binds with HIV-1 protease via both hydrophobic and hydrogen bonding interactions. The hydrophobic interactions from the inhibitor is aimed to the residues of I50 (I50’), I84 (I84’), and V82 (V82’) which create hydrophobic core clusters to further stabilize the closed conformation of flaps, and the hydrogen bonding interactions are mainly focused with the active site of HIV-1 protease. The combined change in the two kinds of protease-inhibitor interactions is correlated with the observed resistance mutations. The present study sheds light on the microscopic mechanism underlying the mutation effects on the dynamics of HIV-1 protease and the inhibition by APV and DRV, providing useful information to the design of more potent and effective HIV-1 protease inhibitors.

QSAR prediction of HIV-1 protease inhibitory activities using docking derived molecular descriptors.

In this study, application of a new hybrid docking-quantitative structure activity relationship (QSAR) methodology to model and predict the HIV-1 protease inhibitory activities of a series of newly synthesized chemicals is reported. This hybrid docking-QSAR approach can provide valuable information about the most important chemical and structural features of the ligands that affect their inhibitory activities. Docking studies were used to find the actual conformations of chemicals in active site of HIV-1 protease. Then the molecular descriptors were calculated from these conformations. Multiple linear regression (MLR) and least square support vector machine (LS-SVM) were used as QSAR models, respectively. The obtained results reveal that statistical parameters of the LS-SVM model are better than the MLR model, which indicate that there are some non-linear relations between selected molecular descriptors and anti-HIV activities of interested chemicals. The correlation coefficient (R), root mean square error (RMSE) and average absolute error (AAE) for LS-SVM are: R=0.988, RMSE=0.207 and AAE=0.145 for the training set, and R=0.965, RMSE=0.403 and AAE=0.338 for the test set. Leave one out cross validation test was used for assessment of the predictive power and validity of models which led to cross-validation correlation coefficient QUOTE of 0.864 and 0.850 and standardized predicted relative error sum of squares (SPRESS) of 0.553 and 0.581 for LS-SVM and MLR models, respectively.

The QSAR and docking calculations of fullerene derivatives as HIV-1 protease inhibitors

The inhibition of HIV-1 protease is considered as one of the most important targets for drug design and the deactivation of HIV-1. In the present work, the fullerene surface (C60) is modified by adding oxygen atoms as well as hydroxymethylcarbonyl (HMC) groups to form 6 investigated fullerene derivative compounds. These compounds have one, two, three, four or five O atoms+HMC groups at different positions on phenyl ring. The effect of the repeating of these groups on the ability of suggested compounds to inhibit the HIV protease is studied by calculating both Quantitative Structure Activity Relationship (QSAR) properties and docking simulation. Based on the QSAR descriptors, the solubility and the hydrophilicity of studied fullerene derivatives increased with increasing the number of oxygen atoms+HMC groups in the compound. While docking calculations indicate that, the compound with two oxygen atoms+HMC groups could interact and binds with HIV-1 protease active site. This is could be attributed to the active site residues of HIV-1 protease are hydrophobic except the two aspartic acids. So that, the increase in the hydrophilicity and polarity of the compound is preventing and/or decreasing the hydrophobic interaction between the compound and HIV-1 protease active site.

Drug resistance conferred by mutations outside the active site through alterations in the dynamic and structural ensemble of HIV-1 protease.

HIV-1 protease inhibitors are part of the highly active antiretroviral therapy effectively used in the treatment of HIV infection and AIDS. Darunavir (DRV) is the most potent of these inhibitors, soliciting drug resistance only when a complex combination of mutations occur both inside and outside the protease active site. With few exceptions, the role of mutations outside the active site in conferring resistance remains largely elusive. Through a series of DRV–protease complex crystal structures, inhibition assays, and molecular dynamics simulations, we find that single and double site mutations outside the active site often associated with DRV resistance alter the structure and dynamic ensemble of HIV-1 protease active site. These alterations correlate with the observed inhibitor binding affinities for the mutants, and suggest a network hypothesis on how the effect of distal mutations are propagated to pivotal residues at the active site and may contribute to conferring drug resistance.

Quantized Water Access to the HIV-1 Protease Active Site as a Proposed Mechanism for Cooperative Mutations in Drug Affinity

The development of resistance to different drugs remains a major problem for a wide range of infections. In particular, combinations of specific mutations, which individually demonstrate no effect, exhibit significant cooperativity. Here we show that changes to the energy of ligand binding in different resistant HIV-1 proteases are correlated with the creation of water binding sites in the active site. This correlation is conserved across two drugs (ritonavir and lopinavir). We propose that individual mutations induce changes in flap packing that are insufficient to allow water binding but in combination allow access, leading to the observed cooperative resistance.

Higher Desolvation Energy Reduces Molecular Recognition in Multi-Drug Resistant HIV-1 Protease.

Designing HIV-1 protease inhibitors that overcome drug-resistance is still a challenging task. In this study, four clinical isolates of multi-drug resistant HIV-1 proteases that exhibit resistance to all the US FDA-approved HIV-1 protease inhibitors and also reduce the substrate recognition ability were examined. A multi-drug resistant HIV-1 protease isolate, MDR 769, was co-crystallized with the p2/NC substrate and the mutated CA/p2 substrate, CA/p2 P1’F. Both substrates display different levels of molecular recognition by the wild-type and multi-drug resistant HIV-1 protease. From the crystal structures, only limited differences can be identified between the wild-type and multi-drug resistant protease. Therefore, a wild-type HIV-1 protease and four multi-drug resistant HIV-1 proteases in complex with the two peptides were modeled based on the crystal structures and examined during a 10 ns-molecular dynamics simulation. The simulation results reveal that the multi-drug resistant HIV-1 proteases require higher desolvation energy to form complexes with the peptides. This result suggests that the desolvation of the HIV-1 protease active site is an important step of protease-ligand complex formation as well as drug resistance. Therefore, desolvation energy could be considered as a parameter in the evaluation of future HIV-1 protease inhibitor candidates.

Design, Synthesis, and X-ray Structure of Substituted Bis-tetrahydrofuran (Bis-THF)-Derived Potent HIV-1 Protease Inhibitors

We investigated substituted bis-THF-derived HIV-1 protease inhibitors in order to enhance ligand-binding site interactions in the HIV-1 protease active site. In this context, we have carried out convenient syntheses of optically active bis-THF and C4-substituted bis-THF ligands using a [2,3]-sigmatropic rearrangement as the key step. The synthesis provided convenient access to a number of substituted bis-THF derivatives. Incorporation of these ligands led to a series of potent HIV-1 protease inhibitors. Inhibitor 23c turned out to be the most potent (K(i) = 2.9 pM; IC(50) = 2.4 nM) among the inhibitors. An X-ray structure of 23c-bound HIV-1 protease showed extensive interactions of the inhibitor with the protease active site, including a unique water-mediated hydrogen bond to the Gly-48 amide NH in the S2 site.

Design, Biologic Evaluation, and SAR of Novel Pseudo‐peptide Incorporating Benzheterocycles as HIV‐1 Protease Inhibitors

A series of novel HIV-1 protease inhibitors based on the (hydroxyethylamino)-sulfonamide isostere incorporating substituted phenyls and benzheterocycle derivatives bearing rich hydrogen bonding acceptors as P(2) ligands were synthesized. Prolonged chain linking the benzhereocycle to the carbonyl group resulted in partial loss of binding affinities. Introduction of a small alkyl substituent with appropriate size to the -CH2- of P(1)-P(2) linkage as a side chain resulted in improved inhibitory potency, and in this study, isopropyl was the best side chain. Replacement of the isobutyl substituent at P(1)'group with phenyl substituent decreased the inhibitory potency. One of the most potent inhibitor, compound 23 showing high affinity to HIV-1 protease with an IC(50) value of 5 nM, also exhibited good anti-SIV activity (EC(50) = 0.8 microM) with low toxicity (TC(50) > 100 microM). The flexible docking of inhibitor 23 to HIV-1 protease active site rationalized the interactions with protease.

Decomposing the Energetic Impact of Drug Resistant Mutations in HIV-1 Protease on Binding DRV

Darunavir (DRV) is a high affinity (4.5×10(-12) M, ΔG = -15.2 kcal/mol) HIV-1 protease inhibitor. Two drug-resistant protease variants FLAP+ (L10I, G48V, I54V, V82A) and ACT (V82T, I84V) decrease the binding affinity with DRV by 1.0 kcal/mol and 1.6 kcal/mol respectively. In this study the absolute and relative binding free energies of DRV with wild-type protease, FLAP+ and ACT were calculated with MM-PB/GBSA and thermodynamic integration methods, respectively. Free energy decomposition elucidated that the mutations conferred resistance by distorting the active site of HIV-1 protease so that the residues that lost binding free energy were not limited to the sites of mutation. Specifically the bis-tetrahydrofuranylurethane moiety of DRV maintained interactions with the FLAP+ and ACT variants, whereas the 4 - amino phenyl group lost more binding free energy with the protease in the FLAP+ and ACT complexes than in the wild-type protease which could account for the majority of the loss in binding free energy. This suggested that replacement of the 4 - amino phenyl group might generate new inhibitors less susceptible to the drug resistant mutations.

Substrate Recognition in HIV-1 Protease: A Computational Study

HIV-1 protease is a crucial enzyme for the life cycle of the human immunodeficiency virus, the retrovirus that triggers AIDS. It is well documented that HIV-1 protease mediates the cleavage of Gag, Gag-Pol, and Nef precursor polyproteins and is highly selective concerning the set of 12 different amino acid sequences that cleaves. However, the governing principles and physical parameters, which determine substrate recognition and specificity, remain poorly understood despite the many speculative proposals that abound in the literature. In fact, it has been difficult so far to circumvent the fact that protease's substrates share little sequence identity and lack an obvious consensus binding motif. We have used microsecond time scale MD simulations to quantitatively show that some sequences of the polyprotein Gag-Pol that are not cleaved (nonsubstrates) have in fact a higher affinity to the active site of HIV-1 protease than a substrate; i.e., recognition is not governed by affinity to the active site. On the basis of a detailed analysis of the results and experimental data, we propose that the recognition is based on the geometric specificity of PR:Gag and PR:Gag-Pol multiprotein complex, that selects which residues lie in the specific position that makes them accessible to the active site for cleavage.

Structural and electronic properties of new fullerene derivatives and their possible application as HIV-1 protease inhibitors

Density functional theory (DFT) calculations have been carried out at the hybrid Becke 3-Lee–Yang–Parr; B3LYP/3-21G** level of theory to study two series of hydroxy-chalca-acetic acid-(4-pyrrolidin-1-yl-phenyl) ester [C 60-C 2H 4N-(4-XCOCH 2OH)C 6H 4] and hydroxy-chalcoacetic acid-[2-(2-hydroxy-acetylchalcanyl)-4-pyrrolidin-1-yl-phenyl] ester[C 60-C 2H 4N-(3,4-XCOCH 2OH)C 6H 4]. The X atom is O, S or Se for the two series. The vibrational spectra, physical, chemical, thermodynamics and Quantitative Structure Activity Relationship (QSAR) properties of the studied molecules are calculated and discussed. We have evaluated these molecules as HIV-1 protease inhibitors based on the hydrogenation interaction between the hydroxymethylcarbonyl (HMC) groups and the two aspartic acid of the HIV-1 protease active site. Results show that some of the investigated fullerene-based derivatives can be considered promising as HIV-1 protease inhibitors.

Catalytic Water Co-Existing with a Product Peptide in the Active Site of HIV-1 Protease Revealed by X-Ray Structure Analysis

BackgroundIt is known that HIV-1 protease is an important target for design of antiviral compounds in the treatment of Acquired Immuno Deficiency Syndrome (AIDS). In this context, understanding the catalytic mechanism of the enzyme is of crucial importance as transition state structure directs inhibitor design. Most mechanistic proposals invoke nucleophilic attack on the scissile peptide bond by a water molecule. But such a water molecule coexisting with any ligand in the active site has not been found so far in the crystal structures.Principal FindingsWe report here the first observation of the coexistence in the active site, of a water molecule WAT1, along with the carboxyl terminal product (Q product) peptide. The product peptide has been generated in situ through cleavage of the full-length substrate. The N-terminal product (P product) has diffused out and is replaced by a set of water molecules while the Q product is still held in the active site through hydrogen bonds. The position of WAT1, which hydrogen bonds to both the catalytic aspartates, is different from when there is no substrate bound in the active site. We propose WAT1 to be the position from where catalytic water attacks the scissile peptide bond. Comparison of structures of HIV-1 protease complexed with the same oligopeptide substrate, but at pH 2.0 and at pH 7.0 shows interesting changes in the conformation and hydrogen bonding interactions from the catalytic aspartates.Conclusions/SignificanceThe structure is suggestive of the repositioning, during substrate binding, of the catalytic water for activation and subsequent nucleophilic attack. The structure could be a snap shot of the enzyme active site primed for the next round of catalysis. This structure further suggests that to achieve the goal of designing inhibitors mimicking the transition-state, the hydrogen-bonding pattern between WAT1 and the enzyme should be replicated.

Inhibition of the Dimerization and Active Site of HIV-1 Protease by Secondary Metabolites from the Vietnamese Mushroom Ganoderma colossum

A new farnesyl hydroquinone, ganomycin I (1), was isolated along with ganomycin B (2) from the chloroform extract of the fruiting bodies of the Vietnamese mushroom Ganoderma colossum. These compounds inhibited HIV-1 protease with IC50 values of 7.5 and 1.0 microg/mL, respectively. Kinetic studies using Zhang-Poorman and Lineweaver plots revealed that compound 2 competitively inhibited the active site of the enzyme, whereas the tetracyclic triterpene schisanlactone A, previously isolated from the same fungus, was a dimerization inhibitor, with an IC50 value of 5.0 microg/mL. The previous findings were also confirmed by the virtual docking of both compounds with HIV-1 protease crystal structure.

Assessment of QM/MM Scoring Functions for Molecular Docking to HIV-1 Protease

We explore the ability of four quantum mechanical (QM)/molecular mechanical (MM) models to accurately identify the native pose of six HIV-1 protease inhibitors and compare them with the AMBER force field and ChemScore and GoldScore scoring functions. Three QM/MM scoring functions treated the ligand at the HF/6-31G*, AM1d, and PM3 levels; the fourth QM/MM function modeled the ligand and active site at the PM3-D level. For the discrimination of native from non-native poses, solvent-corrected HF/6-31G*:AMBER and AMBER functions exhibited the best overall performance. While the electrostatic component of the MM and QM/MM functions appears important for discriminating the native pose of the ligand, the polarization contribution in the QM/MM functions was relatively insensitive to a ligand's binding mode and, for one ligand, actually hindered discrimination. The inclusion of a desolvation penalty, here using a generalized Born solvent model, improved discrimination for the MM and QM/MM methods. There appeared to be no advantage to binding mode prediction by incorporating active site polarization at the PM3-D level. Finally, we found that choice of the protonation state of the aspartyl dyad in the HIV-1 protease active site influenced the ability of scoring methods to determine the native binding pose.