Computational studies of darunavir into HIV-1 protease and DMPC bilayer: necessary conditions for effective binding and the role of the flaps

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the main targets toward AIDS therapy. We have selected the potent drug darunavir and a weak inhibitor (fullerene analog) as HIV-1 PR substrates to compare protease's conformational features upon binding. Molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and quantum-mechanical (QM) calculations indicated the importance of the stability of HIV-1 PR flaps toward effective binding: a weak inhibitor may induce flexibility to the flaps, which convert between closed and semiopen states. A water molecule in the darunavir-HIV-1 PR complex bridged the two flap tips of the protease through hydrogen bonding (HB) interactions in a stable structure, a feature that was not observed for the fullerene-HIV-1 PR complex. Additionally, despite that van der Waals interactions and nonpolar contribution to solvation favored permanent fullerene entrapment into the cavity, these interactions alone were not sufficient for effective binding; enhanced electrostatic interactions as observed in the darunavir-complex were the crucial component of the binding energy. An alternative pathway to the usual way of a ligand to access the cavity was also observed for both compounds. Each ligand entered the binding cavity through an opening between the one flap of the protease and a neighboring loop. This suggested that access to the cavity is not necessarily regulated by flap opening. Darunavir exerts its biological action inside the cell, after crossing the membrane barrier. Thus, we also initiated a study on the interactions between darunavir and the DMPC bilayer to reveal that the drug was accommodated inside the bilayer in conformations that resembled its structure into HIV-1 PR, being stabilized via HBs with the lipids and water molecules.     

Binding free energy contributions of interfacial waters in HIV-1 protease/inhibitor complexes

Water molecules are commonly observed in crystal structures of protein-ligand complexes where they mediate protein-ligand binding. It is of considerable theoretical and practical importance to determine quantitatively the individual free energy contributions of these interfacial water molecules to protein-ligand binding and to elucidate factors that influence them. The double-decoupling free energy molecular dynamics simulation method has been used to calculate the binding free energy contribution for each of the four interfacial water molecules observed in the crystal structure of HIV-1 protease complexed with KNI-272, a potent inhibitor. While two of these water molecules contribute significantly to the binding free energy, the other two have close to zero contribution. It was further observed that the protonation states of two catalytic aspartate residues, Asp25 and Asp125, strongly influence the free energy contribution of a conserved water molecule Wat301 and that different inhibitors significantly influence the free energy contribution of Wat301. Our results have important implications on our understanding of the role of interfacial water molecules in protein-ligand binding and to structure-based drug design aimed at incorporating these interfacial water molecules into ligands.     

Binding of novel fullerene inhibitors to HIV-1 protease: insight through molecular dynamics and molecular mechanics Poisson–Boltzmann surface area calculations

The objectives of this study include the design of a series of novel fullerene-based inhibitors for HIV-1 protease (HIV-1 PR), by employing two strategies that can also be applied to the design of inhibitors for any other target. Additionally, the interactions which contribute to the observed exceptionally high binding free energies were analyzed. In particular, we investigated: (1) hydrogen bonding (H-bond) interactions between specific fullerene derivatives and the protease, (2) the regions of HIV-1 PR that play a significant role in binding, (3) protease changes upon binding and (4) various contributions to the binding free energy, in order to identify the most significant of them. This study has been performed by employing a docking technique, two 3D-QSAR models, molecular dynamics (MD) simulations and the molecular mechanics Poisson–Boltzmann surface area (MM–PBSA) method. Our computed binding free energies are in satisfactory agreement with the experimental results. The suitability of specific fullerene derivatives as drug candidates was further enhanced, after ADMET (absorption, distribution, metabolism, excretion and toxicity) properties have been estimated to be promising. The outcomes of this study revealed important protein–ligand interaction patterns that may lead towards the development of novel, potent HIV-1 PR inhibitors.     

Drug-resistant molecular mechanism of CRF01_AE HIV-1 protease due to V82F mutation

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the major targets of anti-AIDS drug discovery. The circulating recombinant form 01 A/E (CRF01_AE, abbreviated AE) subtype is one of the most common HIV-1 subtypes, which is infecting more humans and is expanding rapidly throughout the world. It is, therefore, necessary to develop inhibitors against subtype AE HIV-1 PR. In this work, we have performed computer simulation of subtype AE HIV-1 PR with the drugs lopinavir (LPV) and nelfinavir (NFV), and examined the mechanism of resistance of the V82F mutation of this protease against LPV both structurally and energetically. The V82F mutation at the active site results in a conformational change of 79′s loop region and displacement of LPV from its proper binding site, and these changes lead to rotation of the side-chains of residues D25 and I50′. Consequently, the conformation of the binding cavity is deformed asymmetrically and some interactions between PR and LPV are destroyed. Additionally, by comparing the interactive mechanisms of LPV and NFV with HIV-1 PR we discovered that the presence of a dodecahydroisoquinoline ring at the P1′ subsite, a [2-(2,6-dimethylphenoxy)acetyl]amino group at the P2′ subsite, and an N2 atom at the P2 subsite could improve the binding affinity of the drug with AE HIV-1 PR. These findings are helpful for promising drug design. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparison of azacyclic urea A-98881 as HIV-1 protease inhibitor with cage dimeric N-benzyl 4-(4-methoxyphenyl)-1,4- dihydropyridine as representative of a novel class of HIV-1 protease inhibitors: A molecular modeling study

The functional groups of cage dimeric N-alkyl substituted 3,5-bis(hydroxymethyl)-4-(4-methoxyphenyl)-1,4-dihydropyridines are similar to those of cyclic and azacyclic ureas that are potent inhibitors of HIV-1 protease of the dihydroxyethylene- and hydroxyethylene type, respectively. In the following study the conformity of common functional groups is investigated concerning their orientation in space as well as in the enzyme HIV-1 protease. Starting from X-ray crystal data of the centrosymmetric cage dimeric N-benzyl derivative with ester groups, the derivative with hydroxymethylene groups was built and a systematic conformational search was performed for the conformationally important torsion angles considering electrostatic and van der Waals interactions. From the huge number of conformations those comprising centrosymmetrical and C2-symmetrical energy minima were selected and minimized. The three remaining conformers were fitted to the azacyclic urea A-98881 selected from the HIV-1 protease enzyme- inhibitor complex using the centroids of the corresponding aromatic residues and additionally by the field fit option of the Advanced CoMFA module of SYBYL. Interestingly, the energetically most favourable one, which, additionally, possesses C2-symmetry like the active site cavity of HIV-1 protease, showed the best fit. Comparing the electrostatic potential (EP) of the latter with the EP of A-98881 the aromatic residues show excellent accordance. Slight differences in the extent of the EP were found in the areas of the hydroxymethylene groups of the cage dimer and the single hydroxy group as well as the urea carbonyl group of A- 98881, respectively. In order to compare the binding possibilities to the enzyme HIV-1 protease for the cage dimer and A-98881, their interaction fields with certain probes (CH3 for alkyl, NHamide, and carbonyl, O– of COO–), representing the decisive functional groups of the active site, have been calculated using GRID and projected into the enzyme placing the structures according to the position of A-98881 in the enzyme- inhibitor complex. The strongest calculated fields of the O– probe were found near Asp 25 for both structures. Another respective conformity consists in the overlap of the fields for the NHamide probe near Ile 50 and 50 for the investigated cage dimer and A-98881.     

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host–guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.     

Computational simulations of HIV-1 proteases--multi-drug resistance due to nonactive site mutation L90M

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the proteins that currently available anti-HIV-1 drugs target. Inhibitors of HIV-1 PR have become available, and they have lowered the rate of mortality from acquired immune deficiency syndrome (AIDS) in advanced countries. However, the rate of emergence of drug-resistant HIV-1 variants is quite high because of their short retroviral life cycle and their high mutation rate. Serious drug-resistant mutations against HIV-1 PR inhibitors (PIs) frequently appear at the active site of PR. Exceptionally, some other mutations such as L90M cause drug resistance, although these appear at nonactive sites. The mechanism of resistance due to nonactive site mutations is difficult to explain. In this study, we carried out computational simulations of L90M PR in complex with each of three kinds of inhibitors and one typical substrate, and we clarified the mechanism of resistance. The L90M mutation causes changes in interaction between the side chain atoms of the 90th residue and the main chain atoms of the 25th residue, and a slight dislocation of the 25th residue causes rotation of the side chain at the 84th residue. The rotation of the 84th residue leads to displacement of the inhibitor from the appropriate binding location, resulting in a collision with the flap or loop region. The difference in levels of resistance to the three inhibitors has been explained from energetic and structural viewpoints, which provides the suggestion for promising drugs keeping its efficacy even for the L90M mutant.     

Structure-based QSAR analysis of a set of 4-hydroxy-5,6-dihydropyrones as inhibitors of HIV-1 protease: an application of the receptor-dependent (RD) 4D-QSAR formalism

Receptor-dependent (RD) 4D-QSAR models were constructed for a set of 39 4-hydroxy-5,6-dihydropyrone analogue HIV-1 protease inhibitors. The receptor model used in this QSAR analysis was derived from the HIV-1 protease (PDB ID ) crystal structure. The bound ligand in the active site of the enzyme, also a 4-hydroxy-5,6-dihydropyrone analogue, was used as the reference ligand for docking the data set compounds. The optimized RD 4D-QSAR models are not only statistically significant (r(2) = 0.86, q(2) = 0.80 for four- and greater-term models) but also possess reasonable predictivity based on test set predictions. The proposed "active" conformations of the docked analogues in the active site of the enzyme are consistent in overall molecular shape with those suggested from crystallographic studies. Moreover, the RD 4D-QSAR models also "capture" the existence of specific induced-fit interactions between the enzyme active site and each specific inhibitor. Hydrophobic interactions, steric shape requirements, and hydrogen bonding of the 4-hydroxy-5,6-dihydropyrone analogues with the HIV-1 protease binding site model dominate the RD 4D-QSAR models in a manner again consistent with experimental conclusions. Some possible hypotheses for the development of new lead HIV-1 protease inhibitors can be inferred from the RD 4D-QSAR models.     

CONCEPTS: New dynamic algorithm for de novo drug suggestion

We describe a new method for de novo design of molecules that bind to protein active sites. The method, CONCEPTS (Creation of Novel Compounds by Evaluation of Particles at Target Sites), places a group of atom-like particles in the site. The particles are free to move within the site to improve binding to the protein. A key innovation of this technique is that covalent connections are made among the particles in a stochastic and dynamically reversible manner. These changes in the topology are either accepted or rejected depending on their ability to improve the total energy of the enzyme–inhibitor complex. The method is applied to two test systems: The FK506 binding protein (FKBP-12) and HIV-1 aspartyl protease. In both cases, we are able to predict, de novo, drugs that have striking similarities to known potent inhibitors and that can successfully be used to generate “hits” of the known inhibitors from a data base. © John Wiley & Sons, Inc.     

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25相似文献   

Prediction and analysis of binding affinities for chemically diverse HIV-1 PR inhibitors by the modified SAFE_p approach

One of the biggest challenges in the "in silico" screening of enzyme ligands is to have a protocol that could predict the ligand binding free energies. In our group we have developed a very simple screening function (referred to as solvent accessibility free energy of binding predictor, SAFE_p) which we have applied previously to the study of peptidic HIV-1 protease (HIV-1 PR) inhibitors and later to cyclic urea type HIV-1 PR inhibitors. In this work, we have extended the SAFE_p protocol to a chemically diverse set of HIV-1 PR inhibitors with binding constants that differ by several orders of magnitude. The resulting function is able to reproduce the ranking and in many cases the value of the inhibitor binding affinities for the HIV-1 PR, with accuracy comparable with that of costlier protocols. We also demonstrate that the binding pocket SAFE_p analysis can contribute to the understanding of the physical forces that participate in ligand binding. The analysis tools afforded by our protocol have allowed us to identify an induced fit phenomena mediated by the inhibitor and have demonstrated that larger fragments do not necessarily contribute the most to the binding free energy, an outcome partially brought about by the substantial role the desolvation penalty plays in the energetics of binding. Finally, we have revisited the effect of the Asp dyad protonation state on the predicted binding affinities.     

Conformational analysis of TMC114, a novel HIV-1 protease inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.     

Virtual screening for anti-HIV-1 RT and anti-HIV-1 PR inhibitors from the Thai medicinal plants database: a combined docking with neural networks approach

The virtual screening approach for docking small molecules into a known protein structure is a powerful tool for drug design. In this work, a combined docking and neural network approach, using a self-organizing map, has been developed and applied to screen anti-HIV-1 inhibitors for two targets, HIV-1 RT and HIV-1 PR, from active compounds available in the Thai Medicinal Plants Database. Based on nevirapine and calanolide A as reference structures in the HIV-1 RT binding site and XK-263 in the HIV-1 PR binding site, 2,684 compounds in the database were docked into the target enzymes. Self-organizing maps were then generated with respect to three types of pharmacophoric groups. The map of the reference structures were then superimposed on the feature maps of all screened compounds. Only the structures having similar features to the reference compounds were accepted. By using the SOMs, the number of candidates for HIV-1 RT was reduced to six and nine compounds consistent with nevirapine and calanolide A, respectively, as references. For the HIV-1 PR target, there are 135 screened compounds showed good agreement with the XK-263 feature map. These screened compounds will be further tested for their HIV-1 inhibitory affinities. The obtained results indicate that this combined method is clearly helpful to perform the successive screening and to reduce the analyzing step from AutoDock and scoring procedure.     

Small-molecule dimerization inhibitors of wild-type and mutant HIV protease: a focused library approach

We demonstrate that a focused library based on truncated, cross-linked interfacial peptides of HIV-1 protease produces effective dimerization inhibitors of the enzyme. By combining individual changes of the library into a single compound, we obtained a significantly more potent agent and found that an additive increase in inhibitor efficacy was obtained. The good activity of library members against an active-site drug-resistant protease mutant bodes well for dimerization inhibition as a complementary method to targeting the active site.     

Does a diol cyclic urea inhibitor of HIV-1 protease bind tighter than its corresponding alcohol form? A study by free energy perturbation and continuum electrostatics calculations

The cyclic urea inhibitors of HIV-1 protease generally have two hydroxyl groups on the seven-membered ring. In this study, free energy perturbation and continuum electrostatic calculations were used to study the contributions of the two hydroxyl groups to the binding affinity and solubility of a cyclic urea inhibitor DMP323. The results indicated that the inhibitor with one hydroxyl group has better binding affinity and solubility than the inhibitor with two hydroxyl groups. Therefore, removal of one hydroxyl group from DMP323 may help to improve the properties of DMP323. This is also likely to be true for other cyclic urea inhibitors. The study also illustrated the difficulty in accurate modeling of the binding affinities of HIV-1 protease inhibitors, which involves many possible protonation states of the two catalytic aspartic acids in the active site of the enzyme.     

Enhancing protein backbone binding--a fruitful concept for combating drug-resistant HIV

The evolution of drug resistance is one of the most fundamental problems in medicine. In HIV/AIDS, the rapid emergence of drug-resistant HIV-1 variants is a major obstacle to current treatments. HIV-1 protease inhibitors are essential components of present antiretroviral therapies. However, with these protease inhibitors, resistance occurs through viral mutations that alter inhibitor binding, resulting in a loss of efficacy. This loss of potency has raised serious questions with regard to effective long-term antiretroviral therapy for HIV/AIDS. In this context, our research has focused on designing inhibitors that form extensive hydrogen-bonding interactions with the enzyme's backbone in the active site. In doing so, we limit the protease's ability to acquire drug resistance as the geometry of the catalytic site must be conserved to maintain functionality. In this Review, we examine the underlying principles of enzyme structure that support our backbone-binding concept as an effective means to combat drug resistance and highlight their application in our recent work on antiviral HIV-1 protease inhibitors.     

Macrocyclic Protease Inhibitors with Reduced Peptide Character

There is a real need for simple structures that define a β‐strand conformation, a secondary structure that is central to peptide–protein interactions. For example, protease substrates and inhibitors almost universally adopt this geometry on active site binding. A planar pyrrole is used to replace two amino acids of a peptide backbone to generate a simple macrocycle that retains the required geometry for active site binding. The resulting β‐strand templates have reduced peptide character and provide potent protease inhibitors with the attachment of an appropriate amino aldehyde to the C‐terminus. Picomolar inhibitors of cathepsin L and S are reported and the mode of binding of one example to the model protease chymotrypsin is defined by X‐ray crystallography.     

Rational design and synthesis of a novel class of active site-targeted HIV protease inhibitors containing a hydroxymethylcarbonyl isostere. Use of phenylnorstatine or allophenylnorstatine as a transition-state mimic.

A novel class of HIV-1 protease inhibitors containing a hydroxymethylcarbonyl (HMC) isostere were designed from the substrate transition state and synthesized. Phenylnorstatine [Pns; (2R,3S)-3-amino-2-hydroxy-4-phenylbutyric acid] and the 2S diastereomer, (2S,3S)-3-amino-2-hydroxy-4-phenylbutyric acid, named allophenylnorstatine (Apns) were effective transition-state mimics, and incorporation of Pns-Pro or Apns-Pro at the P1-P1' site gave potent and specific HIV-1 protease inhibitors. In the inhibitory assays, the chemically synthesized [Ala67,95] HIV-1 protease was used.