Constructing HIV-1 integrase tetramer and exploring influences of metal ions on forming integrase-DNA complex

HIV-1 integrase (IN) is essential for the replication of HIV-1 in human cells. At present, the complete structure of complex IN-DNA has not been resolved. In this paper, a HIV-1 IN tetramer model was built with homology modeling and molecular dynamics simulation approach, in which two Mg2+ ions were reasonably located in each catalytic core domain. Moreover, it was found that the AB and CD chains of HIV-1 IN tetramer were different in the structures and metal ions of HIV-1 IN tetramer might have great influences on DNA locating on IN. These findings may provide a more complete structural basis for guiding drug discovery and revealing integration mechanism.     

Structural and Dynamical Properties of a Full-length HIV-1 Integrase: Molecular Dynamics Simulations

Abstract

The structural and dynamical properties of the complete full-length structure of HIV-1 integrase were investigated using Molecular Dynamics approach. Simulations were carried out for the three systems, core domain only (CORE), full-length structure without (FULL) and with a Mg²⁺ (FULL+ION) in its active site, aimed to investigate the difference in the molecular properties of the full-length models due to their different construction procedures as well as the effects of the two ends, C- and N-terminal, on those properties in the core domain. The full-length structure was prepared from the two experimental structures of two-domain fragment. The following properties were observed to differ significantly from the previous reports: (i) relative topology formed by an angle between the three domains; (ii) the cavity size defined by the catalytic triad, Asp64, Asp116, and Glul52; (iii) distances and solvation of the Mg²⁺; and (iv) conformation of the catalytic residues. In addition, the presence of the two terminal domains decreases the mobility of the central core domain significantly.     

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

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.     

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整合酶核心区野生型和F185K突变型活性和溶解性的比较及分子动力学模拟分析

表达纯化了野生型(WT)及F185K突变型HIV-1整合酶核心区蛋白(IN_C),并对二者的溶解性和活性进行了比较．实验结果表明：F185K 突变后IN_C溶解性显著提高,活性有一定程度降低．对WT和F185K IN_C体系进行了1800 ps的分子动力学模拟．模拟结果表明：F185K IN_C功能loop区柔性和蛋白质整体运动性降低,使蛋白质活性降低,F185K突变后盐桥网络的变化驱动了IN_C局部构象改变,引起IN_C表面的部分疏水残基被包埋,亲水残基暴露,相对亲水溶剂可接近面积增大,同时,突变后IN_C与水之间形成氢键的数量增加,与水之间作用加强,以上变化使IN_C溶解性提高．分子动力学模拟与实验结果相吻合．为理解蛋白质溶解性和对蛋白质进行可溶性改造提供了一定的理论依据．     

Interchain hydrophobic clustering promotes rigidity in HIV-1 protease flap dynamics: new insights from Molecular Dynamics

The dynamics of HIV-1 protease (HIV-pr), a drug target for HIV infection, has been studied extensively by both computational and experimental methods. The flap dynamics of HIV-pr is considered to be more important for better ligand binding and enzymatic actions. Moreover, it has been demonstrated that the drug-induced mutations can change the flap dynamics of HIV-pr affecting the binding affinity of the ligands. Therefore, detailed understanding of flap dynamics is essential for designing better inhibitors. Previous computational investigations observed significant variation in the flap opening in nanosecond time scale indicating that the dynamics is highly sensitive to the simulation protocols. To understand the sensitivity of the flap dynamics on the force field and simulation protocol, molecular dynamics simulations of HIV-pr have been performed with two different AMBER force fields, ff99 and ff02. Two different trajectories (20?ns each) were obtained using the ff99 and ff02 force field. The results showed polarizable force field (ff02) make the flap tighter than the nonpolarizable force field (ff99). Some polar interactions and hydrogen bonds involving flap residues were found to be stronger with ff02 force field. The formation of interchain hydrophobic cluster (between flap tip of one chain and active site wall of another chain) was found to be dominant in the semi-open structures obtained from the simulations irrespective of the force field. It is proposed that an inhibitor, which will promote this interchain hydrophobic clustering, may make the flaps more rigid, and presumably the effect of mutation would be small on ligand binding.     

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.     

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.     

Crystal structure of the HIV-1 integrase core domain in complex with sucrose reveals details of an allosteric inhibitory binding site

HIV integrase (IN) is an essential enzyme in HIV replication and an important target for drug design. IN has been shown to interact with a number of cellular and viral proteins during the integration process. Disruption of these important interactions could provide a mechanism for allosteric inhibition of IN. We present the highest resolution crystal structure of the IN core domain to date. We also present a crystal structure of the IN core domain in complex with sucrose which is bound at the dimer interface in a region that has previously been reported to bind integrase inhibitors.

Structured summary

MINT-7713125: IN (uniprotkb:P04585) and IN (uniprotkb:P04585) bind (MI:0407) by X-ray crystallography (MI:0114)     

A novel short peptide is a specific inhibitor of the human immunodeficiency virus type 1 integrase

The retroviral encoded protein integrase (IN) is required for the insertion of the human immunodeficiency virus type 1 (HIV-1) proviral DNA into the host genome. In spite of the crucial role played by IN in the retroviral life cycle, which makes this enzyme an attractive target for the development of new anti-AIDS agents, very few inhibitors have been described and none seems to have a potential use in anti-HIV therapy. To obtain potent and specific IN inhibitors, we used the two-hybrid system to isolate short peptides. Using HIV-1 IN as a bait and a yeast genomic library as the source of inhibitory peptides (prey), we isolated a 33-mer peptide (I33) that bound tightly to the enzyme. I33 inhibited both in vitro IN activities, i.e. 3' end processing and strand transfer. Further analysis led us to select a shorter peptide, EBR28, corresponding to the N-terminal region of I33. Truncated variants showed that EBR28 interacted with the catalytic domain of IN interfering with the binding of the DNA substrate. Alanine single substitution of each EBR28 residue (alanine scanning) allowed the identification of essential amino acids involved in the inhibition. The EBR28 NMR structure shows that this peptide adopts an alpha-helical conformation with amphipathic properties. Additionally, EBR28 showed a significant antiviral effect when assayed on HIV-1 infected human cells. Thus, this potentially important short lead peptide may not only be helpful to design new anti-HIV agents, but also could prove very useful in further studies of the structural and functional characteristics of HIV-1 IN.     

Molecular dynamics studies on HIV-1 protease: a comparison of the flap motions between wild type protease and the M46I/G51D double mutant

The emergence of drug-resistant mutants of HIV-1 is a tragic effect associated with conventional long-treatment therapies against acquired immunodeficiency syndrome. These mutations frequently involve the aspartic protease encoded by the virus; knowledge of the molecular mechanisms underlying the conformational changes of HIV-1 protease mutants may be useful in developing more effective and longer lasting treatment regimes. The flap regions of the protease are the target of a particular type of mutations occurring far from the active site. These mutations modify the affinity for both substrate and ligands, thus conferring resistance. In this work, molecular dynamics simulations were performed on a native wild type HIV-1 protease and on the drug-resistant M46I/G51D double mutant. The simulation was carried out for a time of 3.5 ns using the GROMOS96 force field, with implementation of the SPC216 explicit solvation model. The results show that the flaps may exist in an ensemble of conformations between a “closed” and an “open” conformation. The behaviour of the flap tips during simulations is different between the native enzyme and the mutant. The mutation pattern leads to stabilization of the flaps in a semi-open configuration.     

The strand transfer oligonucleotide inhibitors of HIV-integrase

Retroviral integrase participates in two catalytic reactions, which require interactions with the two ends of the viral DNA in the 3′processing reaction, and with a targeted host DNA in the strand transfer reaction. The 3′-hydroxyl group of 2′-deoxyadenosine resulting from the specific removing of GT dinucleotide from the viral DNA in the processing reaction provides the attachment site for the host DNA in a transesterification reaction. We synthesized oligonucleotides (ONs) of various lengths that mimic the processed HIV-1 U5 terminus of the proviral long terminal repeat (LTR) and are ended by 2′-deoxyadenosine containing a 3′-O-phosphonomethyl group. The duplex stability of phosphonomethyl ONs was increased by covalent linkage of the modified strand with its complementary strand by a triethylene glycol loop (TEG). Modified ONs containing up to 10 bases inhibited in vitro the strand transfer reaction catalyzed by HIV-1 integrase at nanomolar concentrations.     

Enzymatic capability of HIS-tagged HIV-1 integrase using oligonucleotide disintegration substrates

Disintegration, wherein a half-site integration substrate is resolved into separate viral and host DNA components via DNA strand transfer, is one of three well-established in vitro activities of HIV-1 integrase. The role of disintegration in the HIV-1 replicative cycle, however, remains a mystery. In this report, we describe the expression inEscherichia coli and purification of HIV-1 integrase as a fusion protein containing a 6×His tag at its amino terminus. Integrase resolved dumbbell and Y-substrates optimally at pH 6.8–7.2 in the presence of 2 mM MnCl₂. Substrate requirements for intramolecular disintegration included a 10 base pair viral U5 LTR arm and a CA dinucleotide located at the 3 end of the LTR. Disintegration was not sensitive to changes in the host DNA portion of the substrate. A dumbbell substrate with a 5 oligo-dA tail also underwent disintegration. The released LTR arm with an oligo-dA tail was utilized as a template primer by several DNA polymerases indicating that disintegration occurred via nucleophilic attack on the phosphodiester bond located immediately adjacent to the CA dinucleotide at the 3 end of the LTR. Coupled disintegration-DNA polymerase reactions provided a highly efficient and sensitive means of detecting disintegration activity. Integrase also catalyzed an apparently concerted disintegration-5-end joining reaction in which an LTR arm was transferred from one dumbbell substrate molecule to another.     

HIV-1整合酶四聚体结构模拟及其活性位点分析

HIV-1复制需要HIV-1整合酶将其环状DNA整合进宿主DNA中,这其中包括2个重要反应,即“3′-加工”和“链转移”,两者均由HIV-1整合酶催化完成.阻断其中的任一反应,都能达到抑制HIV-1复制的目的.因此,了解HIV-1整合酶的完整结构和聚合状态,对深入探讨其作用机理及设计新型抑制剂具有重要的指导作用.然而,迄今为止仅有HIV-1整合酶单独结构域的晶体结构可供参考,而其全酶晶体结构尚未获得解析.本研究利用分子模拟技术,通过蛋白质 蛋白质/DNA分子对接、动力学模拟等方法,构建了全长整合酶四聚体的结构模型、HIV-1 DNA与整合酶复合物的结构模型,进一步从理论上证实HIV-1整合酶是以四聚体形态发挥催化作用,明确“3′-加工”和“链转移”在HIV-1整合酶上的催化位点.同时,通过与作用机理相似的细菌转座子Tn5转座酶等的结构比对,推测HIV-1整合酶的核心结构域中应有第2个Mg^２＋存在,其位置螯合于Asp64与Glu152之间.在HIV-1整合酶结构研究的基础上,有望进一步设计出新的抗艾滋病药物.     

Molecular dynamics of the full-length p53 monomer

The p53 protein is frequently mutated in a very large proportion of human tumors, where it seems to acquire gain-of-function activity that facilitates tumor onset and progression. A possible mechanism is the ability of mutant p53 proteins to physically interact with other proteins, including members of the same family, namely p63 and p73, inactivating their function. Assuming that this interaction might occurs at the level of the monomer, to investigate the molecular basis for this interaction, here, we sample the structural flexibility of the wild-type p53 monomeric protein. The results show a strong stability up to 850 ns in the DNA binding domain, with major flexibility in the N-terminal transactivations domains (TAD1 and TAD2) as well as in the C-terminal region (tetramerization domain). Several stable hydrogen bonds have been detected between N-terminal or C-terminal and DNA binding domain, and also between N-terminal and C-terminal. Essential dynamics analysis highlights strongly correlated movements involving TAD1 and the proline-rich region in the N-terminal domain, the tetramerization region in the C-terminal domain; Lys120 in the DNA binding region. The herein presented model is a starting point for further investigation of the whole protein tetramer as well as of its mutants.     

Molecular dynamics simulations applied to the study of subtypes of HIV-1 protease common to Brazil, Africa, and Asia

Africa accounts for the majority of HIV-1 infections worldwide caused mainly by the A and C viral subtypes rather than B subtype, which prevails in the United States and Western Europe. In Brazil, B subtype is the major subtype, but F, C, and A also circulate. These non-B subtypes present polymorphisms, and some of them occur at sites that have been associated with drug resistance, including the HIV-1 protease (PR), one important drug target. Here, we report a Molecular Dynamics study of the B and non-B PR complexed with the inhibitor ritonavir to delineate the behavior of each subtype. We compare root mean squared deviation, binding free energy by linear interaction energy approach, hydrogen bonds, and intermolecular contact surface area between inhibitor and PR. From our results, we can provide a basis to understand the molecular mechanism of drug resistance in non-B subtypes. In this sense, we found a decrease of approx 4 kcal/mol in ΔG of binding between B and non-B subtypes. This corresponds to the loss of one hydrogen bond, which is in agreement with our H-bond analysis. Previous experimental affinity studies reported analogous results with inhibition constant values for non-B PR.     

NMR structure of the Tn916 integrase-DNA complex

The integrase protein catalyzes the excision and integration of the Tn916 conjugative transposon, a promiscuous genetic element that spreads antibiotic resistance in pathogenic bacteria. The solution structure of the N-terminal domain of the Tn916 integrase protein bound to its DNA-binding site within the transposon arm has been determined. The structure reveals an interesting mode of DNA recognition, in which the face of a three-stranded antiparallel beta-sheet is positioned within the major groove. A comparison to the structure of the homing endonuclease I-Ppol-DNA complex suggests that the three-stranded sheet may represent a new DNA-binding motif whose residue composition and position within the major groove are varied to alter specificity. The structure also provides insights into the mechanism of conjugative transposition. The DNA in the complex is bent approximately 35 degrees and may, together with potential interactions between bound integrase proteins at directly repeated sites, significantly bend the arms of the transposon.     

Molecular dynamics studies of the conformation of sorbitol

Molecular dynamics simulations of a 3 molal aqueous solution of d-sorbitol (also called d-glucitol) have been performed at 300 K, as well as at two elevated temperatures to promote conformational transitions. In principle, sorbitol is more flexible than glucose since it does not contain a constraining ring. However, a conformational analysis revealed that the sorbitol chain remains extended in solution, in contrast to the bent conformation found experimentally in the crystalline form. While there are 243 staggered conformations of the backbone possible for this open-chain polyol, only a very limited number were found to be stable in the simulations. Although many conformers were briefly sampled, only eight were significantly populated in the simulation. The carbon backbones of all but two of these eight conformers were completely extended, unlike the bent crystal conformation. These extended conformers were stabilized by a quite persistent intramolecular hydrogen bond between the hydroxyl groups of carbon C-2 and C-4. The conformational populations were found to be in good agreement with the limited available NMR data except for the C-2–C-3 torsion (spanned by the O-2–O-4 hydrogen bond), where the NMR data support a more bent structure.