Structure-based virtual screening toward the discovery of novel inhibitors for impeding the protein-protein interaction between HIV-1 integrase and human lens epithelium-derived growth factor (LEDGF/p75)

HIV-1 integrase is a unique promising component of the viral replication cycle, catalyzing the integration of reverse transcribed viral cDNA into the host cell genome. Generally, IN activity requires both viral as well as a cellular co-factor in the processing replication cycle. Among them, the human lens epithelium-derived growth factor (LEDGF/p75) represented as promising cellular co-factor which supports the viral replication by tethering IN to the chromatin. Due to its major importance in the early steps of HIV replication, the interaction between IN and LEDGF/p75 has become a pleasing target for anti-HIV drug discovery. The present study involves the finding of novel inhibitor based on the information of dimeric CCD of IN in complex with known inhibitor, which were carried out by applying a structure-based virtual screening concept with molecular docking. Additionally, Free binding energy, ADME properties, PAINS analysis, Density Functional Theory, and Enrichment Calculations were performed on selected compounds for getting a best lead molecule. On the basis of these analyses, the current study proposes top 3 compounds: Enamine-Z742267384, Maybridge-HTS02400, and Specs-AE-848/37125099 with acceptable pharmacological properties and enhanced binding affinity to inhibit the interaction between IN and LEDGF/p75. Furthermore, Simulation studies were carried out on these molecules to expose their dynamics behavior and stability. We expect that the findings obtained here could be future therapeutic agents and may provide an outline for the experimental studies to stimulate the innovative strategy for research community.     

Design and discovery of flavonoid-based HIV-1 integrase inhibitors targeting both the active site and the interaction with LEDGF/p75

HIV integrase (IN) is an essential enzyme for the viral replication. Currently, three IN inhibitors have been approved for treating HIV-1 infection. All three drugs selectively inhibit the strand transfer reaction by chelating a divalent metal ion in the enzyme active site. Flavonoids are a well-known class of natural products endowed with versatile biological activities. Their β-ketoenol or catechol structures can serve as a metal chelation motif and be exploited for the design of novel IN inhibitors. Using the metal chelation as a common pharmacophore, we introduced appropriate hydrophobic moieties into the flavonol core to design natural product-based novel IN inhibitors. We developed selective and efficient syntheses to generate a series of mono 3/5/7/3′/4′-substituted flavonoid derivatives. Most of these new compounds showed excellent HIV-1 IN inhibitory activity in enzyme-based assays and protected against HIV-1 infection in cell-based assays. The 7-morpholino substituted 7c showed effective antiviral activity (EC₅₀ = 0.826 μg/mL) and high therapeutic index (TI > 242). More significantly, these hydroxyflavones block the IN–LEDGF/p75 interaction with low- to sub-micromolar IC₅₀ values and represent a novel scaffold to design new generation of drugs simultaneously targeting the catalytic site as well as protein–protein interaction domains.     

Identification of small peptides inhibiting the integrase‐LEDGF/p75 interaction through targeting the cellular co‐factor

The integration of the viral DNA into the host genome is one of the essential steps in the HIV replication cycle. This process is mediated by the viral enzyme integrase (IN) and lens epithelium‐derived growth factor (LEDGF/p75). LEDGF/p75 has been identified as a crucial cellular co‐factor of integration that acts by tethering IN to the cellular chromatin. Recently, circular peptides were identified that bind to the C‐terminal domain of IN and disrupt the interaction with LEDGF/p75. Starting from the circular peptides, we identified a short peptidic sequence able to inhibit the LEDGF/p75‐IN interaction at low μM concentration through its binding to the IN binding site of LEDGF/p75. This discovery can lead to the synthesis of peptidomimetics with high anti‐HIV activity targeting the cellular co‐factor LEDGF/p75 and not the viral protein IN. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Identification of the LEDGF/p75 binding site in HIV-1 integrase

Lens epithelium-derived growth factor (LEDGF)/p75 is an important cellular co-factor for human immunodeficiency virus (HIV) replication. We originally identified LEDGF/p75 as a binding partner of integrase (IN) in human cells. The interaction has been mapped to the integrase-binding domain (IBD) of LEDGF/p75 located in the C-terminal part. We have subsequently shown that IN carrying the Q168A mutation remains enzymatically active but is impaired for interaction with LEDGF/p75. To map the integrase/LEDGF interface in more detail, we have now identified and characterized two regions within the enzyme involved in the interaction with LEDGF/p75. The first region centers around residues W131 and W132 while the second extends from I161 up to E170. For the different IN mutants the interaction with LEDGF/p75 and the enzymatic activities were determined. IN(W131A), IN(I161A), IN(R166A), IN(Q168A) and IN(E170A) are impaired for interaction with LEDGF/p75, but retain 3' processing and strand transfer activities. Due to impaired integration, an HIV-1 strain containing the W131A mutation in IN displays reduced replication capacity, whereas virus carrying IN(Q168A) is replication defective. Comparison of the wild-type IN-LEDGF/p75 co-crystal structure with that of the modelled structure of the IN(Q168A) and IN(W131A) mutant integrases corroborated our experimental data.     

Inhibitory profile of a LEDGF/p75 peptide against HIV-1 integrase: insight into integrase-DNA complex formation and catalysis

A lens epithelium-derived growth factor (LEDGF)/p75 peptide was evaluated for human immunodeficiency virus type 1 integrase (IN) inhibitory activity. The LEDGF/p75 peptide modestly inhibited IN catalysis and was dependent on IN-DNA assembly. The peptide was also effective at disrupting LEDGF/p75-IN complex formation. We next investigated the activity of the LEDGF/p75 peptide on IN mutant proteins that are unable to catalyze the DNA strand transfer reaction. The LEDGF/p75 peptide displayed an increased potency on these IN proteins, from 5-fold to greater than 10-fold, indicating the IN multimeric state greatly influences the peptide inhibitory effects. These results shed light on IN-DNA multimeric formation, and how this process influences the LEDGF/p75-IN interaction.     

Effective prediction model and determination of binding residues influential for inhibitors targeting HIV-1 integrase-LEDGF/p75 interface by employing solvent accessible surface area energy as key determinant

Abstract

Development of a highly accurate prediction model for protein–ligand inhibition has been a major challenge in drug discovery. Herein, we describe a novel predictive model for the inhibition of HIV-1 integrase (IN)-LEDGF/p75 protein-protein interaction. The model was constructed using energy parameters approximated from molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations. Chemometric analysis using partial least squares (PLS) regression revealed that solvent accessible surface area energy (ΔG_SASA) is the major determinant parameter contributing greatly to the prediction accuracy. PLS prediction model on the ΔG_SASA values collected from 41 complexes yielded a strong correlation between the predicted and the actual inhibitory activities (R² = 0.9666, RMSEC of pIC₅₀ values = 0.0890). Additionally, for the test set of 14 complexes, the model performed satisfactorily with very low pIC₅₀ errors (Q² = 0.5168, RMSEP = 0.3325). A strong correlation between the buried surface areas on the IN protein, when bound with IN-LEDGF/p75 inhibitors, and the respective ΔG_SASA values was also obtained. Furthermore, the current method could identify ‘hot spots’of amino acid residues highly influential to the inhibitory activity prediction. This could present fruitful implications in binding site determination and future inhibitor developments targeting protein-protein interactions.

Communicated by Ramaswamy H. Sarma     

D77, one benzoic acid derivative, functions as a novel anti-HIV-1 inhibitor targeting the interaction between integrase and cellular LEDGF/p75

Integration of viral-DNA into host chromosome mediated by the viral protein HIV-1 integrase (IN) is an essential step in the HIV-1 life cycle. In this process, Lens epithelium-derived growth factor (LEDGF/p75) is discovered to function as a cellular co-factor for integration. Since LEDGF/p75 plays an important role in HIV integration, disruption of the LEDGF/p75 interaction with IN has provided a special interest for anti-HIV agent discovery. In this work, we reported that a benzoic acid derivative, 4-[(5-bromo-4-{[2,4-dioxo-3-(2-oxo-2-phenylethyl)-1,3-thiazolidin-5-ylidene]methyl}-2-ethoxyphenoxy)methyl]benzoic acid (D77) could potently inhibit the IN-LEDGF/p75 interaction and affect the HIV-1 IN nuclear distribution thus exhibiting antiretroviral activity. Molecular docking with site-directed mutagenesis analysis and surface plasmon resonance (SPR) binding assays has clarified possible binding mode of D77 against HIV-1 integrase. As the firstly discovered small molecular compound targeting HIV-1 integrase interaction with LEDGF/p75, D77 might supply useful structural information for further anti-HIV agent discovery.     

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.     

Design,synthesis and biological evaluation of imidazole and oxazole fragments as HIV-1 integrase-LEDGF/p75 disruptors and inhibitors of microbial pathogens

We describe here the synthesis of libraries of novel 1-subtituted-5-aryl-1H-imidazole, 5-aryl-4-tosyl-4,5-dihydro-1,3-oxazole and 5-aryl-1,3-oxazole fragments via microwave (MW)-assisted cycloaddition of para-toluenesulfonylmethyl isocyanide (TosMIC) to imines and aldehydes. The compounds obtained were biologically evaluated in an AlphaScreen HIV-1 IN-LEDGF/p75 inhibition assay with six imidazole-based compounds (16c, 16f, 17c, 17f, 20a and 20d) displaying more than 50% inhibition at 10 µM, with IC₅₀ values ranging from 7.0 to 30.4 µM. Additionally the hypothesis model developed predicts all active scaffolds except 20d to occupy similar areas as the N-heterocyclic (A) moiety and two aromatic rings (B and C) of previously identified inhibitor 5. These results indicate that the identified compounds represent a viable starting point for their use as templates in the design of next generation inhibitors targeting the HIV-1 IN and LEDGF/p75 protein-protein interaction. In addition, the in vitro antimicrobial properties of these fragments were tested by minimum inhibitory concentration (MIC) assays showing that compound 16f exhibited a MIC value of 15.6 μg/ml against S. aureus, while 17f displayed a similar MIC value against B. cereus, suggesting that these compounds could be further developed to specifically target those microbial pathogens.     

Fragment hopping approach directed at design of HIV IN-LEDGF/p75 interaction inhibitors

We recently identified a series of indole derivatives as active inhibitors of IN-LEDGF/p75 interaction through structure-based pharmacophore models generated from the crystal structure of dimeric catalytic core domain (CCD) of HIV-1 IN in complex with the LEDGF integrase binding domain (IBD). In this paper we used the fragment hopping approach to design small molecules able to prevent the IN-LEDGF/p75 interaction. By means of the proposed approach, we designed novel non-peptidyl compounds that mimic the biological function of some IBD residues and in particular the LEDGF hot spot residues Ile365 and Asp366. The biological results confirmed the importance of several structural requirements for the inhibitory effects of this class of compounds.     

Structural basis for HIV‐1 DNA integration in the human genome,role of the LEDGF/P75 cofactor

Integration of the human immunodeficiency virus (HIV‐1) cDNA into the human genome is catalysed by integrase. Several studies have shown the importance of the interaction of cellular cofactors with integrase for viral integration and infectivity. In this study, we produced a stable and functional complex between the wild‐type full‐length integrase (IN) and the cellular cofactor LEDGF/p75 that shows enhanced in vitro integration activity compared with the integrase alone. Mass spectrometry analysis and the fitting of known atomic structures in cryo negatively stain electron microscopy (EM) maps revealed that the functional unit comprises two asymmetric integrase dimers and two LEDGF/p75 molecules. In the presence of DNA, EM revealed the DNA‐binding sites and indicated that, in each asymmetric dimer, one integrase molecule performs the catalytic reaction, whereas the other one positions the viral DNA in the active site of the opposite dimer. The positions of the target and viral DNAs for the 3′ processing and integration reaction shed light on the integration mechanism, a process with wide implications for the understanding of viral‐induced pathologies.     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease

Two different structures of ligand-free HIV protease have been determined by X-ray crystallography. These structures differ in the position of two 12 residue, β-hairpin regions (or “flaps”) which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand-bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand-free enzyme, it is found that the two structures are nearly equally stable, with the ligand-bound-type structure being less stable, consistent with X-ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the β-sheet structure of the conformationally flexible, glycine-rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed. Proteins 32:7–16, 1998. © 1998 Wiley-Liss, Inc.     

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation

Abstract

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q² = 0.751, NOC = 6, r² = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q² = 0.744, NOC = 6, r² = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset.

Communicated by Ramaswamy H. Sarma     

Insight into the intrinsic flexibility of the SL1 stem-loop from genomic RNA of HIV-1 as probed by molecular dynamics simulation

The SL1 stem-loop is the dimerization initiation site for linking the two copies of the RNA forming the HIV-1 genome. The 26 nucleotides stem contains a defect consisting on a highly conserved G-rich 1-3 asymmetrical internal loop, which is a major site for nucleocapsid protein binding. Several NMR attempts were undertaken to determine the internal loop structure in the SL1 monomer. However, the RNA constructs used in the different studies were largely mutated, in particular with replacement of the nine nucleotides apical loop by a tetraloop, and divergent results were obtained ranging from a rigid structure to a particularly large flexibility. To investigate the reasons for such discrepancies, we used molecular dynamics simulation of the SL1 monomer to probe the effect of mutations on the conformational stability of the internal loop and of the whole stem. It is found that in the wild-type sequence, the internal loop displays conformational variability originating mainly from the nine nucleotides apical loop flexibility that causes large conformational fluctuations (without changing the average structure) in the 7 bp duplex linking the apical and internal loops. The large amplitude atomic motions in the duplex are transmitted to the internal loop in which they induce conformational changes characterized by a labile hydrogen bond network such as G5 successively H-bonded to A29 and G30. On the contrary, with a four nucleotides apical loop, conformational fluctuations in the duplex are reduced by a factor of 2 and are not sufficiently energizing for promoting changes in the internal loop structure at the time scale of the simulations.