Hybrid quantum mechanical/molecular mechanical molecular dynamics simulations of HIV-1 integrase/inhibitor complexes

Human immunodeficiency virus (HIV)-1 integrase (IN) is an attractive target for development of acquired immunodeficiency syndrome chemotherapy. In this study, conventional and coupled quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations of HIV-1 IN complexed with 5CITEP (IN-5CITEP) were carried out. In addition to differences in the bound position of 5CITEP, significant differences at the two levels of theory were observed in the metal coordination geometry and the areas involving residues 116-119 and 140-166. In the conventional MD simulation, the coordination of Mg(2+) was found to be a near-perfect octahedral geometry whereas a distorted octahedral complex was observed in QM/MM. All of the above reasons lead to a different pattern of protein-ligand salt link formation that was not observed in the classical MD simulation. Furthermore to provide a theoretical understanding of inhibition mechanisms of 5CITEP and its derivative (DKA), hybrid QM/MM MD simulations of the two complexes (IN-5CITEP and IN-DKA) have been performed. The results reveal that areas involving residues 60-68, 116-119, and 140-149 were substantially different among the two systems. The two systems show similar pattern of metal coordination geometry, i.e., a distorted octahedron. In IN-DKA, both OD1 and OD2 of Asp-64 coordinate the Mg(2+) in a monodentate fashion whereas only OD1 is chelated to the metal as observed in IN-5CITEP. The high potency of DKA as compared to 5CITEP is supported by a strong salt link formed between its carboxylate moiety and the ammonium group of Lys-159. Detailed comparisons between HIV-1 IN complexed with DKA and with 5CITEP provide information about ligand structure effects on protein-ligand interactions in particular with the Lys-159. This is useful for the design of new selective HIV-1 IN inhibitors.     

Calculation of binding energy using BLYP/MM for the HIV-1 integrase complexed with the S-1360 and two analogues

Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. S-1360 is a potent and selective inhibitor of HIV-1 IN. In this work, we have carried out molecular dynamics (MD) simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) approach, to determine the protein-ligand interaction energy for S-1360 and two analogues. Analysis of the MD trajectories reveals that the strongest protein-inhibitor interactions, observed in the three studied complexes, are established with Lys-159 residue and Mg(2+) cation. Calculations of binding energy using BLYP/MM level of theory reveal that there is a direct relationship between this theoretical computed property and the experimental determined anti-HIV activity.     

Understanding the effect of drug-resistant mutations of HIV-1 intasome on raltegravir action through molecular modeling study

Raltegravir is the first FDA-approved drug targeting the strand transfer step of HIV-1 integration. However, the rapid emergence of viral strains that are highly resistant to raltegravir has become a critical problem. Unfortunately, the detailed molecular mechanism of how HIV-1 integrase (IN) mutations actually confer drug resistance is not well understood. In the present study, starting from our previously constructed complex of HIV-1 IN and viral DNA, we employed molecular dynamics (MD) simulation and molecular mechanics generalized Born surface area (MM-GBSA) calculation, to uncover the molecular mechanism behind the resistant mechanism of HIV-1 IN to raltegravir. The values of the calculated binding free energy follow consistently the experimentally observed ranking of resistance levels. A detailed analysis of the results of MD simulation suggests that the Tyr143 located in the 140s loop (e.g., residues from Gly140 to Gly149) is a key anchoring residue that leads to stable raltegravir binding. The decrease in the interaction at this residue is one of the key reasons responsible for the resistance of HIV-1 IN to raltegravir. Additionally, the calculation results also proved that the 3' adenosine flip in different conformations in the wild-type and mutant HIV-1 IN-viral DNA complexes play an important role in raltegravir binding. Our results could provide a structural and energetic understanding of the raltegravir-resistant mechanism at the atomic level and provide some new clues on how to design new drugs that may circumvent the known resistance mutations.     

用分子模拟方法研究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的抑制剂设计和改造提供帮助．     

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.     

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.     

Mutations in the Human Immunodeficiency Virus Type 1 Integrase D,D(35)E Motif Do Not Eliminate Provirus Formation

The core domain of human immunodeficiency virus type 1 (HIV-1) integrase (IN) contains a D,D(35)E motif, named for the phylogenetically conserved glutamic acid and aspartic acid residues and the invariant 35 amino acid spacing between the second and third acidic residues. Each acidic residue of the D,D(35)E motif is independently essential for the 3′-processing and strand transfer activities of purified HIV-1 IN protein. Using a replication-defective viral genome with a hygromycin selectable marker, we recently reported that a mutation at any of the three residues of the D,D(35)E motif produces a 10³- to 10⁴-fold reduction in infectious titer compared with virus encoding wild-type IN (A. D. Leavitt et al., J. Virol. 70:721–728. 1996). The infectious titer, as measured by the number of hygromycin-resistant colonies formed following infection of cells in culture, was less than a few hundred colonies per μg of p24. To understand the mechanism by which the mutant virions conferred hygromycin resistance, we characterized the integrated viral DNA in cells infected with virus encoding mutations at each of the three residues of the D,D(35)E motif. We found the integrated viral DNA to be colinear with the incoming viral genome. DNA sequencing of the junctions between integrated viral DNA and host DNA showed that (i) the characteristic 5-bp direct repeat of host DNA flanking the HIV-1 provirus was not maintained, (ii) integration often produced a deletion of host DNA, (iii) integration sometimes occurred without the viral DNA first undergoing 3′-processing, (iv) integration sites showed a strong bias for a G residue immediately adjacent to the conserved viral CA dinucleotide, and (v) mutations at each of the residues of the D,D(35)E motif produced essentially identical phenotypes. We conclude that mutations at any of the three acidic residues of the conserved D,D(35)E motif so severely impair IN activity that most, if not all, integration events by virus encoding such mutations are not IN mediated. IN-independent provirus formation may have implications for anti-IN therapeutic agents that target the IN active site.     

Molecular dynamics studies of the wild-type and double mutant HIV-1 integrase complexed with the 5CITEP inhibitor: mechanism for inhibition and drug resistance

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the life cycle of the virus and is an attractive target for the development of new drugs useful in acquired immunodeficiency syndrome multidrug therapy. Starting from the crystal structure of the 5CITEP inhibitor bound to the active site in the catalytic domain of the HIV-1 IN, two different molecular dynamics simulations in water have been carried out. In the first simulation the wild-type IN was used, whereas in the second one the double mutation T66I/M154I, described to lead to drug resistance, was introduced in the protein. Compelling differences have been observed in these two structures during analyses of the molecular dynamics trajectories, particularly in the inhibitor binding modes and in the conformational flexibility of the loop (residues 138-149) located near the three catalytic residues in the active site (Asp(64), Asp(116), Glu(152)). Because the conformational flexibility of this region is important for efficient biological activity and its behavior is quite different in the two models, we suggest a hypothetical mechanism for the inhibition and drug resistance of HIV-1 IN. These results can be useful for the rational design of more potent and selective integrase inhibitors and may allow for the design of inhibitors that will be more robust against known resistance mutations.     

Binding mode prediction of biologically active compounds from plant Salvia Miltiorrhiza as integrase inhibitor

Integrase (IN), an essential enzyme for HIV-1 replication, has been targeted in antiretroviral drug therapy. The emergence of HIV-1 variants clinically resistant to antiretroviral agents has lead to the development of alternative IN inhibitors. In the present work, binding modes of a high potent IN inhibitor, M₅22 and M₅32, within the catalytic binding site of wild type (WT) IN were determined using molecular docking calculation. Both M₅22 and M₅32 displayed similar modes of binding within the IN putative binding pocket and exhibited favorable interactions with the catalytic Mg²⁺ ions, the nearby amino acids and viral DNA through metal-ligand chelation, hydrogen bonding and π-π stacking interactions. Furthermore, the modes of action of these two compounds against the mutated Y212R, N224H and S217H PFV IN were also predicted. Although the replacement of amino acid could somehow disturb inhibitor binding mode, almost key interactions which detected in the WT complexes were fairly conserved. Detailed information could highlight the application of M₅22 and M₅32 as candidate IN inhibitors for drug development against drug resistant strains.     

Insights into saquinavir resistance in the G48V HIV-1 protease: quantum calculations and molecular dynamic simulations

The spread of acquired immune deficiency syndrome has increasingly become a great concern owing largely to the failure of chemotherapies. The G48V is considered the key signature residue mutation of HIV-1 protease developing with saquinavir therapy. Molecular dynamics simulations of the wild-type and the G48V HIV-1 protease complexed with saquinavir were carried out to explore structure and interactions of the drug resistance. The molecular dynamics results combined with the quantum-based and molecular mechanics Poisson-Boltzmann surface area calculations indicated a monoprotonation took place on D25, one of the triad active site residues. The inhibitor binding of the triad residues and its interaction energy in the mutant were similar to those in the wild-type. The overall structure of both complexes is almost identical. However, the steric conflict of the substituted valine results in the conformational change of the P2 subsite and the disruption of hydrogen bonding between the -NH of the P2 subsite and the backbone -CO of the mutated residue. The magnitude of interaction energy changes was comparable to the experimental K(i) data. The designing for a new drug should consider a reduction of steric repulsion on P2 to enhance the activity toward this mutant strain.     

Dolutegravir Interactions with HIV-1 Integrase-DNA: Structural Rationale for Drug Resistance and Dissociation Kinetics

Signature HIV-1 integrase mutations associated with clinical raltegravir resistance involve 1 of 3 primary genetic pathways, Y143C/R, Q148H/K/R and N155H, the latter 2 of which confer cross-resistance to elvitegravir. In accord with clinical findings, in vitro drug resistance profiling studies with wild-type and site-directed integrase mutant viruses have shown significant fold increases in raltegravir and elvitegravir resistance for the specified viral mutants relative to wild-type HIV-1. Dolutegravir, in contrast, has demonstrated clinical efficacy in subjects failing raltegravir therapy due to integrase mutations at Y143, Q148 or N155, which is consistent with its distinct in vitro resistance profile as dolutegravir’s antiviral activity against these viral mutants is equivalent to its activity against wild-type HIV-1. Kinetic studies of inhibitor dissociation from wild-type and mutant integrase-viral DNA complexes have shown that dolutegravir also has a distinct off-rate profile with dissociative half-lives substantially longer than those of raltegravir and elvitegravir, suggesting that dolutegravir’s prolonged binding may be an important contributing factor to its distinct resistance profile. To provide a structural rationale for these observations, we constructed several molecular models of wild-type and clinically relevant mutant HIV-1 integrase enzymes in complex with viral DNA and dolutegravir, raltegravir or elvitegravir. Here, we discuss our structural models and the posited effects that the integrase mutations and the structural and electronic properties of the integrase inhibitors may have on the catalytic pocket and inhibitor binding and, consequently, on antiviral potency in vitro and in the clinic.     

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.     

Determining the Functions of HIV-1 Tat and a Second Magnesium Ion in the CDK9/Cyclin T1 Complex: A Molecular Dynamics Simulation Study

The current paradigm of cyclin-dependent kinase (CDK) regulation based on the well-established CDK2 has been recently expanded. The determination of CDK9 crystal structures suggests the requirement of an additional regulatory protein, such as human immunodeficiency virus type 1 (HIV-1) Tat, to exert its physiological functions. In most kinases, the exact number and roles of the cofactor metal ions remain unappreciated, and the repertoire has thus gained increasing attention recently. Here, molecular dynamics (MD) simulations were implemented on CDK9 to explore the functional roles of HIV-1 Tat and the second Mg²⁺ ion at site 1 (Mg₁ ²⁺). The simulations unveiled that binding of HIV-1 Tat to CDK9 not only stabilized hydrogen bonds (H-bonds) between ATP and hinge residues Asp104 and Cys106, as well as between ATP and invariant Lys48, but also facilitated the salt bridge network pertaining to the phosphorylated Thr186 at the activation loop. By contrast, these H-bonds cannot be formed in CDK9 owing to the absence of HIV-1 Tat. MD simulations further revealed that the Mg₁ ²⁺ ion, coupled with the Mg₂ ²⁺ ion, anchored to the triphosphate moiety of ATP in its catalytic competent conformation. This observation indicates the requirement of the Mg₁ ²⁺ ion for CDK9 to realize its function. Overall, the introduction of HIV-1 Tat and Mg₁ ²⁺ ion resulted in the active site architectural characteristics of phosphorylated CDK9. These data highlighted the functional roles of HIV-1 Tat and Mg₁ ²⁺ ion in the regulation of CDK9 activity, which contributes an important complementary understanding of CDK molecular underpinnings.     

Why is Substrate Peptide Binding Unsusceptible to Multidrug-Resistant Mutations in HIV-1 Protease? A Structural and Energetic Analysis

Understanding of the molecular mechanism and biological implication underlying the difference in binding of substrate peptides and small-molecule inhibitors to multidrug-resistant mutants of HIV-1 protease would help to develop new anti-HIV agents combating drug resistance. Here, an integration of rigorous quantum mechanics/molecular mechanics (QM/MM) analysis and empirical Poisson–Boltzmann/surface area (PB/SA) model is described to investigate the structural basis and energetic property of wild-type HIV-1 protease and its mutants in recognizing and binding with a wide variety of ligands, including the peptides derived from its cognate cleavage sites and the cleavage site variants as well as a number of FDA-approved protease inhibitors, attempting to explain why is substrate binding unsusceptible to most observed HIV-1 protease mutations. A preliminary test study demonstrates that the combined QM/MM–PB/SA scheme is able to effectively reproduce the relative ligand binding energy changes upon protease single- and double-mutations, albeit the absolute values appear to be different significantly between the calculated and experimental results. With the QM/MM–PB/SA calculations a complete mutation energy map of HIV-1 protease–ligand interactions is created, which unravels distinct affinity pictures of wild-type substrates, substrate variants and, particularly, the protease inhibitors bound to HIV-1 protease mutants, suggesting that, on the one hand, the evaluation pressure under anti-HIV chemotherapies addresses site-directed protease mutations that impair and undermine the intermolecular interactions specific to inhibitors but not substrates; on the other hand, co-evaluation of protease and its substrate peptides provides a more effective mechanism to avoid therapeutic surveillance. Further, nonbonded interaction analysis and computational alanine scanning reveal 12 key residues that is critical for substrate binding, from which the Asn25, Gly27, Ala28, Asp29 and Pro81 are identified that have not yet been found to cause drug resistance and hence would be the promising sites targeted by new protease inhibitors.     

Study on the molecular mechanism of inhibiting HIV-1 integrase by EBR28 peptide via molecular modeling approach

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the HIV-1 lifecycle which aids the integration of viral DNA into the host chromosome. Recently synthesized 12-mer peptide EBR28, which can strongly bind to IN, is one of the most potential small peptide leading compounds inhibiting IN binding with viral DNA. However, the binding mode between EBR28 peptide with HIV-1 IN and the inhibition mechanism remain uncertain. In this paper, the binding modes of EBR28 with HIV-1 IN monomer core domain (IN(1)) and dimmer core domain (IN(2)) were investigated by using molecular docking and molecular dynamics (MD) simulation methods. The results indicated that EBR28 bound to the interfaces of the IN(1) and IN(2) systems mainly through the hydrophobic interactions with the beta3, alpha1 and alpha5 regions of the proteins. The binding free energies for IN(1) with a series of EBR28 mutated peptides were calculated with the MM/GBSA model, and the correlation between the calculated and experimental binding free energies is very good (r=0.88). Thus, the validity of the binding mode of IN(1) with EBR28 was confirmed. Based on the binding modes, the inhibition mechanism of EBR28 was explored by analyzing the essential dynamics (ED), energy decomposition and the mobility of EBR28 in the two docked complexes. The proposed inhibition mechanism is represented that EBR28 binds to the interface of IN(1) to form the IN(1)_EBR28 complex and preventes the formation of IN dimmer, finally leads to the partial loss of binding potency for IN with viral DNA. All of the above simulation results agree well with experimental data, which provide us with some helpful information for designing anti-HIV small peptide drugs.     

Substrate Recognition and Motion Mode Analyses of PFV Integrase in Complex with Viral DNA via Coarse-Grained Models

HIV-1 integrase (IN) is an important target in the development of drugs against the AIDS virus. Drug design based on the structure of IN was markedly hampered due to the lack of three-dimensional structure information of HIV-1 IN-viral DNA complex. The prototype foamy virus (PFV) IN has a highly functional and structural homology with HIV-1 IN. Recently, the X-ray crystal complex structure of PFV IN with its cognate viral DNA has been obtained. In this study, both Gaussian network model (GNM) and anisotropy network model (ANM) have been applied to comparatively investigate the motion modes of PFV DNA-free and DNA-bound IN. The results show that the motion mode of PFV IN has only a slight change after binding with DNA. The motion of this enzyme is in favor of association with DNA, and the binding ability is determined by its intrinsic structural topology. Molecular docking experiments were performed to gain the binding modes of a series of diketo acid (DKA) inhibitors with PFV IN obtained from ANM, from which the dependability of PFV IN-DNA used in the drug screen for strand transfer (ST) inhibitors was confirmed. It is also found that the functional groups of keto-enol, bis-diketo, tetrazole and azido play a key role in aiding the recognition of viral DNA, and thus finally increase the inhibition capability for the corresponding DKA inhibitor. Our study provides some theoretical information and helps to design anti-AIDS drug based on the structure of IN.     

Development of a high-throughput assay for the HIV-1 integrase disintegration reaction

Both HIV-1 integrase (IN) and the central catalytic domain of IN (IN-CCD) catalyze the disintegration reaction in vitro. In this study, IN and IN-CCD proteins were expressed and purified, and a high-throughput format enzyme-linked immunosorbent assay (ELISA) was developed for the disintegration reaction. IN exhibited a marked preference for Mn²⁺ over Mg²⁺ as the divalent cation cofactor in disintegration. Baicalein, a known IN inhibitor, was found to be an IN-CCD inhibitor. The assay is sensitive and specific for the study of disintegration reaction as well as for the in vitro identification of antiviral drugs targeting IN, especially targeting IN-CCD.