Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10?ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.     

用分子模拟方法研究HIV-1整合酶突变体的耐药性机理

二酮酸类化合物(DKAs)是目前最有前景的HIV-1整合酶(integrase, IN)抑制剂．为了解DKAs引起的多种耐药株共有的耐药性机理,选择3种S-1360引起的IN耐药突变体,用分子对接和分子动力学模拟,研究了野生型和突变型IN与S-1360的结合模式,基于该结合模式探讨了3种耐药突变体所共有的耐药性机理．结果表明：在突变体中,S-1360结合到耐药突变IN核心区中的位置靠近功能loop 3区却远离与 DNA结合的关键残基,结合位置不同导致S-1360的抑制作用部分丧失;残基138到166区域的柔性对IN发挥生物学功能很重要,S-1360能与DNA结合的关键残基N155及K159形成氢键,这2个氢键作用降低了该区域的柔性,突变体中无类似氢键,因而该区域柔性增高;在突变体中,S-1360的苯环远离病毒DNA结合区,不能阻止病毒DNA末端暴露给宿主DNA;T66I突变导致残基Ⅰ的长侧链占据IN的活性口袋,阻止抑制剂以与野生型中相同的方式结合到活性中心,这均是产生抗药性的重要原因．这些模拟结果与实验结果吻合,可为抗IN的抑制剂设计和改造提供帮助．     

Broad antiretroviral activity and resistance profile of the novel human immunodeficiency virus integrase inhibitor elvitegravir (JTK-303/GS-9137)

Integrase (IN), an essential enzyme of human immunodeficiency virus (HIV), is an attractive antiretroviral drug target. The antiviral activity and resistance profile in vitro of a novel IN inhibitor, elvitegravir (EVG) (also known as JTK-303/GS-9137), currently being developed for the treatment of HIV-1 infection are described. EVG blocked the integration of HIV-1 cDNA through the inhibition of DNA strand transfer. EVG inhibited the replication of HIV-1, including various subtypes and multiple-drug-resistant clinical isolates, and HIV-2 strains with a 50% effective concentration in the subnanomolar to nanomolar range. EVG-resistant variants were selected in two independent inductions, and a total of 8 amino acid substitutions in the catalytic core domain of IN were observed. Among the observed IN mutations, T66I and E92Q substitutions mainly contributed to EVG resistance. These two primary resistance mutations are located in the active site, and other secondary mutations identified are proximal to these primary mutations. The EVG-selected IN mutations, some of which represent novel IN inhibitor resistance mutations, conferred reduced susceptibility to other IN inhibitors, suggesting that a common mechanism is involved in resistance and potential cross-resistance. The replication capacity of EVG-resistant variants was significantly reduced relative to both wild-type virus and other IN inhibitor-resistant variants selected by L-870,810. EVG and L-870,810 both inhibited the replication of murine leukemia virus and simian immunodeficiency virus, suggesting that IN inhibitors bind to a conformationally conserved region of various retroviral IN enzymes and are an ideal drug for a range of retroviral infections.     

Identification of HIV-1 integrase inhibitors via three-dimensional database searching using ASV and HIV-1 integrases as targets

Integration of viral DNA into the host cell genome is a critical step in the life cycle of HIV. This essential reaction is catalyzed by integrase (IN) through two steps, 3'-processing and DNA strand transfer. Integrase is an attractive target for drug design because there is no known cellular analogue and integration is essential for successful replication of HIV. A computational three-dimensional (3-D) database search was used to identify novel HIV-1 integrase inhibitors. Starting from the previously identified Y3 (4-acetylamino-5-hydroxynaphthalene-2,7-disulfonic acid) binding site on the avian sarcoma virus integrase (ASV IN), a preliminary search of all compounds in the nonproprietary, open part of the National Cancer Institute 3-D database yielded a collection of 3100 compounds. A more rigorous scoring method was used to rescreen the 3100 compounds against both ASV IN and HIV-1 IN. Twenty-two of those compounds were selected for inhibition assays against HIV-1 IN. Thirteen of the 22 showed inhibitory activity against HIV-1 IN at concentrations less than 200 microM and three of them showed antiviral activities in HIV-1 infected CEM cells with effective concentrations (EC50) ranging from 0.8 to 200 microM. Analysis of the computer-generated binding modes of the active compounds to HIV-1 IN showed that simultaneous interaction with the Y3 site and the catalytic site is possible. In addition, interactions between the active compounds and the flexible loop involved in the binding of DNA by IN are indicated to occur. The structural details and the unique binding motif between the HIV-1 IN and its inhibitors identified in the present work may contribute to the future development of IN inhibitors.     

Binding aspects of baicalein to HIV-1 integrase

Human immunodeficiency virus type 1 (HIV-1) integrase is an essential enzyme in the life cycle of the virus. It is responsible for catalyzing the insertion of the viral genome into the host cell chromosome. This integrase is an attractive target for the design of a HIV antiviral drug, because integrase has no human counterpart. In order to know the interaction mode of HIV-1 integrase with its inhibitor, we investigated the effect of the inhibitor, baicalein, on the conformation of the HIV-1 integrase catalytic domain [IN-(50-212/F185K)] using fluorescence and circular dichroism (CD) spectroscopy. We found that baicalein binds to the hydrophobic region of the HIV-1 integrase catalytic core domain. This binding of baicalein induces the conformational change of the enzyme. We also found that the binding ratio of baicalein to the HIV-1 integrase catalytic domain is 2:1.     

Hydrophobic sliding: a possible mechanism for drug resistance in human immunodeficiency virus type 1 protease

Hydrophobic residues outside the active site of HIV-1 protease frequently mutate in patients undergoing protease inhibitor therapy; however, the mechanism by which these mutations confer drug resistance is not understood. From analysis of molecular dynamics simulations, 19 core hydrophobic residues appear to facilitate the conformational changes that occur in HIV-1 protease. The hydrophobic core residues slide by each other, exchanging one hydrophobic van der Waal contact for another, with little energy penalty, while maintaining many structurally important hydrogen bonds. Such hydrophobic sliding may represent a general mechanism by which proteins undergo conformational changes. Mutation of these residues in HIV-1 protease would alter the packing of the hydrophobic core, affecting the conformational flexibility of the protease. Therefore these residues impact the dynamic balance between processing substrates and binding inhibitors, and thus contribute to drug resistance.     

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.     

Comparison of multiple molecular dynamics trajectories calculated for the drug-resistant HIV-1 integrase T66I/M154I catalytic domain

HIV-1 integrase (IN) is an essential enzyme for the viral replication and an interesting target for the design of new pharmaceuticals for multidrug therapy of AIDS. Single and multiple mutations of IN at residues T66, S153, or M154 confer degrees of resistance to several inhibitors that prevent the enzyme from performing its normal strand transfer activity. Four different conformations of IN were chosen from a prior molecular dynamics (MD) simulation on the modeled IN T66I/M154I catalytic core domain as starting points for additional MD studies. The aim of this article is to understand the dynamic features that may play roles in the catalytic activity of the double mutant enzyme in the absence of any inhibitor. Moreover, we want to verify the influence of using different starting points on the MD trajectories and associated dynamical properties. By comparison of the trajectories obtained from these MD simulations we have demonstrated that the starting point does not affect the conformational space explored by this protein and that the time of the simulation is long enough to achieve convergence for this system.     

Active site binding modes of the beta-diketoacids: a multi-active site approach in HIV-1 integrase inhibitor design

Predicting a bioactive conformation of a ligand is of paramount importance in rational drug design. The task becomes very difficult when the receptor site possesses a region with unusual conformational flexibility. Significant conformational differences are present in the active site regions in the available crystal structures of the core domains of HIV-1 integrase (IN). Among all reported IN inhibitors, the β-diketoacid class of compounds has proved to be of most promise and indeed S-1360 was the first IN inhibitor to enter clinical studies. With an aim to predict the bioactive (active site bound) conformation of S-1360, we performed extensive docking studies using three different reported crystal structures where the active site or partial active site region was resolved. For comparison we extended our studies to include 5CITEP (the first compound cocrystallized with IN core domain) and a bis-diketoacid (BDKA). We found that the conformation of S-1360 when bound in one of the active sites matches that of the experimentally observed results of IN escape mutants resistant to S-1360. Therefore, we propose that this active site conformation is the biologically relevant conformation and can be used for the future structure-based drug design studies selectively targeting IN.     

Active site binding modes of curcumin in HIV-1 protease and integrase

Structure models for the interaction of curcumin with HIV-1 integrase (IN) and protease (PR) were investigated using computational docking. Curcumin was found to bind preferentially in similar ways to the active sites of both IN and PR. For IN, the binding site is formed by residues Asp64, His67, Thr66, Glu92, Thr93, Asp116, Ser119, Asn120, and Lys159. Docked curcumin contacts the catalytic residues adjacent to Asp116 and Asp64, and near the divalent metal (Mg2+). In the PR docking, the curcumin structure fitted well to the active site, interacting with residues Asp25, Asp29, Asp30, Gly27', Asp29', and Asp30'. The results suggest that o-hydroxyl and/or keto-enol structures are important for both IN and PR inhibitory actions. The symmetrical structure of curcumin seems to play an important role for binding to the PR protein, whereas the keto-enol and only one side of the terminal o-hydroxyl showed tight binding to the IN active site.     

Large-scale conformational dynamics of the HIV-1 integrase core domain and its catalytic loop mutants

HIV-1 integrase is one of the three essential enzymes required for viral replication and has great potential as a novel target for anti-HIV drugs. Although tremendous efforts have been devoted to understanding this protein, the conformation of the catalytic core domain around the active site, particularly the catalytic loop overhanging the active site, is still not well characterized by experimental methods due to its high degree of flexibility. Recent studies have suggested that this conformational dynamics is directly correlated with enzymatic activity, but the details of this dynamics is not known. In this study, we conducted a series of extended-time molecular dynamics simulations and locally enhanced sampling simulations of the wild-type and three loop hinge mutants to investigate the conformational dynamics of the core domain. A combined total of >480 ns of simulation data was collected which allowed us to study the conformational changes that were not possible to observe in the previously reported short-time molecular dynamics simulations. Among the main findings are a major conformational change (>20 A) in the catalytic loop, which revealed a gatinglike dynamics, and a transient intraloop structure, which provided a rationale for the mutational effects of several residues on the loop including Q(148), P(145), and Y(143). Further, clustering analyses have identified seven major conformational states of the wild-type catalytic loop. Their implications for catalytic function and ligand interaction are discussed. The findings reported here provide a detailed view of the active site conformational dynamics and should be useful for structure-based inhibitor design for integrase.     

Development of an AlphaScreen-based HIV-1 integrase dimerization assay for discovery of novel allosteric inhibitors

In recent years, HIV-1 integrase (IN) has become an established target in the field of antiretroviral drug discovery. However, its sole clinically approved inhibitor, the integrase strand transfer inhibitor (INSTI) raltegravir, has a surprisingly low genetic barrier for resistance. Furthermore, the only two other integrase inhibitors currently in advanced clinical trials, elvitegravir and dolutegravir, share its mechanism of action and certain resistance pathways. To maintain a range of treatment options, drug discovery efforts are now turning toward allosteric IN inhibitors, which should be devoid of cross-resistance with INSTIs. As IN requires a precise and dynamic equilibrium between several oligomeric species for its activities, the modulation of this equilibrium presents an interesting allosteric target. We report on the development, characterization, and validation of an AlphaScreen-based assay for high-throughput screening for modulators of HIV-1 IN dimerization. Compounds identified as hits in this assay proved to act as allosteric IN inhibitors. Additionally, the assay offers a flexible platform to study IN dimerization.     

Identification of amino acids in HIV-1 and avian sarcoma virus integrase subsites required for specific recognition of the long terminal repeat Ends

A tetramer model for HIV-1 integrase (IN) with DNA representing 20 bp of the U3 and U5 long terminal repeats (LTR) termini was assembled using structural and biochemical data and molecular dynamics simulations. It predicted amino acid residues on the enzyme surface that can interact with the LTR termini. A separate structural alignment of HIV-1, simian sarcoma virus (SIV), and avian sarcoma virus (ASV) INs predicted which of these residues were unique. To determine whether these residues were responsible for specific recognition of the LTR termini, the amino acids from ASV IN were substituted into the structurally equivalent positions of HIV-1 IN, and the ability of the chimeras to 3 ' process U5 HIV-1 or ASV duplex oligos was determined. This analysis demonstrated that there are multiple amino acid contacts with the LTRs and that substitution of ASV IN amino acids at many of the analogous positions in HIV-1 IN conferred partial ability to cleave ASV substrates with a concomitant loss in the ability to cleave the homologous HIV-1 substrate. HIV-1 IN residues that changed specificity include Val(72), Ser(153), Lys(160)-Ile(161), Gly(163)-Val(165), and His(171)-Leu(172). Because a chimera that combines several of these substitutions showed a specificity of cleavage of the U5 ASV substrate closer to wild type ASV IN compared with chimeras with individual amino acid substitutions, it appears that the sum of the IN interactions with the LTRs determines the specificity. Finally, residues Ser(153) and Val(72) in HIV-1 IN are among those that change in enzymes that develop resistance to naphthyridine carboxamide- and diketo acid-related inhibitors in cells. Thus, amino acid residues involved in recognition of the LTRs are among these positions that change in development of drug resistance.     

Structural analysis of a mutant of the HIV-1 integrase zinc finger domain that forms a single conformation

HIV-1 integrase consists of three functional domains, an N-terminal zinc finger domain, a catalytic core domain and a C-terminal DNA binding domain. NMR analysis of an isolated N-terminal domain (IN(1-55)) has shown that IN(1-55) exists in two conformational states [E and D forms; Cai et al. (1997) Nat. Struct. Biol. 4, 567-577]. The two forms differ in the coordination of the zinc ion by two histidine residues. In the present study, structural analysis of a mutant of IN(1-55), Y15A, by NMR spectroscopy indicated that the mutant protein folds correctly but takes only the E form. Since the Y15A mutation abrogates the HIV-1 infectivity, Y15 might have some important role in the full-length integrase activity during the virus infection cycle. Our results suggest a possible role of Y15 in structural transition between the E and D forms of HIV-1 integrase to allow the optimal tetramerization.     

Insights into the interactions between HIV-1 integrase and human LEDGF/p75 by molecular dynamics simulation and free energy calculation

Human cellular protein LEDGF/p75 (lens epithelium-derived growth factor) is an important binding partner of human immunodeficiency virus type 1 (HIV-1) integrase (IN). Without LEDGF/p75, HIV-1 can not complete its life cycle. To study the detailed interactions between LEDGF/p75 and HIV-1 IN, and then obtain the hotspots at the binding interface, 13 ns molecular dynamics simulations were carried out here. One-hundred snapshots extracted from the last 4 ns trajectories were used for calculation of binding free energy and decomposition of the energy by residue. First, the structural changes and their dynamic interactions were investigated focused on the production stage. And then, the free energy was discussed. On the basis of the above results, it could be suggested that residues Gln168, Glu170, and Thr174 in chain A of IN, Thr125, and Trp131 in chain B of IN as well as Ile365, Asp366, Phe406, and Val408 in LEDGF/p75 were responsible for their binding. These results might be helpful for discovery and design of small molecules to interrupt the interaction between HIV-1 IN and LEDGF/p75.     

A quantum mechanic/molecular mechanic study of the wild-type and N155S mutant HIV-1 integrase complexed with diketo acid

Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. Recently, mutation studies have been reported that have shown that a certain degree of viral resistance to diketo acids (DKAs) appears when some amino acid residues of the IN active site are mutated. Mutations represent a fascinating experimental challenge, and we invite theoretical simulations for the disclosure of still unexplored features of enzyme reactions. The aim of this work is to understand the molecular mechanisms of HIV-1 IN drug resistance, which will be useful for designing anti-HIV inhibitors with unique resistance profiles. In this study, we use molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy for wild-type and N155S mutant HIV-1 IN, both complexed with a DKA. This hybrid methodology has the advantage of the inclusion of quantum effects such as ligand polarization upon binding, which can be very important when highly polarizable groups are embedded in anisotropic environments, for example in metal-containing active sites. Furthermore, an energy terms decomposition analysis was performed to determine contributions of individual residues to the enzyme-inhibitor interactions. The results reveal that there is a strong interaction between the Lys-159, Lys-156, and Asn-155 residues and Mg²⁺ cation and the DKA inhibitor. Our calculations show that the binding energy is higher in wild-type than in the N155S mutant, in accordance with the experimental results. The role of the mutated residue has thus been checked as maintaining the structure of the ternary complex formed by the protein, the Mg²⁺ cation, and the inhibitor. These results might be useful to design compounds with more interesting anti-HIV-1 IN activity on the basis of its three-dimensional structure.     

Litchi chinensis-derived terpenoid as anti-HIV-1 protease agent: structural design from molecular dynamics simulations

The molecular structures of the binding between human immunodeficiency virus-1 protease (HIV-1PR) and various inhibitors including existing extensive natural products extracts have been investigated for anti-HIV drug development. In this study, the binding of HIV-1PR and a terpenoid from Litchi chinensis extracts (3-oxotrirucalla-7,24-dien-21-oic acid) was investigated in order to clarify the inhibition effectiveness of this compound. Molecular dynamics (MD) simulations of HIV-1PR complex with 3-oxotrirucalla-7,24-dien-21-oic acid were performed including water molecules. The MD simulation results indicated the formation of hydrogen bonds between the oxygen atoms of the inhibitor and the catalytic aspartates, which are commonly found in inhibitors–protease complexes. On the other hand, no hydrogen bonding of this particular inhibitor to the flap region was found. In addition, the radial distribution function of water oxygens around the catalytic carboxylate nitrogens of Asp29 and Asp30 suggests that at least one or two water molecules are in the active site region whereas direct interaction of the inhibitor was found for catalytic carboxylate oxygen of Asp25. The results of this simulation, in comparison with the structures of other HIV-PR inhibitor complexes, could lead to a better understanding of the activity of 3-oxotrirucalla-7,24-dien-21-oic acid.     

Mechanism of inhibition of HIV-1 integrase by G-tetrad-forming oligonucleotides in Vitro

The G-tetrad-forming oligonucleotides and have been identified as potent inhibitors of human immunodeficiency virus type 1 integrase (HIV-1 IN) activity (Rando, R. F., Ojwang, J., Elbaggari, A., Reyes, G. R., Tinder, R., McGrath, M. S., and Hogan, M. E. (1995) J. Biol. Chem. 270, 1754-1760; Mazumder, A., Neamati, N., Ojwang, J. O., Sunder, S., Rando, R. F., and Pommier, Y. (1996) Biochemistry 35, 13762-13771; Jing, N., and Hogan, M. E. (1998) J. Biol. Chem. 273, 34992-34999). To understand the inhibition of HIV-1 IN activity by the G-quartet inhibitors, we have designed the oligonucleotides and, composed of three and four G-quartets with stem lengths of 19 and 24 A, respectively. The fact that increasing the G-quartet stem length from 15 to 24 A kept inhibition of HIV-1 IN activity unchanged suggests that the binding interaction occurs between a GTGT loop domain of the G-quartet inhibitors and a catalytic site of HIV-1 IN, referred to as a face-to-face interaction. Docking the NMR structure of (Jing and Hogan (1998)) into the x-ray structure of the core domain of HIV-1 IN, HIV-1 IN-(51-209) (Maignan, S., Guilloteau, J.-P. , Qing, Z.-L., Clement-Mella, C., and Mikol, V. (1998) J. Mol. Biol. 282, 359-368), was performed using the GRAMM program. The statistical distributions of hydrogen bonding between HIV-1 IN and were obtained from the analyses of 1000 random docking structures. The docking results show a high probability of interaction between the GTGT loop residues of the G-quartet inhibitors and the catalytic site of HIV-1 IN, in agreement with the experimental observation.