AutoDocking dinucleotides to the HIV-1 integrase core domain: exploring possible binding sites for viral and genomic DNA

A mini-review on integrase inhibitors: The cornerstone of next-generation HIV treatment

Integrase inhibitors represent one of the most remarkable and effective advances in the treatment of HIV-1 infection. Their lack of human cellular equivalence has established integrase as a unique and ideal target for HIV-1 treatment. Over the last two decades, a variety of drugs and small molecule inhibitors have been developed to control or treat HIV infection. Many of these FDA-approved drugs are considered first-line options for AIDS patients. Unfortunately, resistance to these drugs has dictated the development of novel and more efficacious antiretroviral drugs. In this review article, we illustrate the key classes of antiretroviral integrase inhibitors available. We provide a comprehensive analysis of recent advancements in the development of integrase inhibitors, focusing on novel compounds and their distinct mechanisms of action. Our literature review highlights emerging allosteric integrase inhibitors that offer improved efficacy, resistance profiles, and pharmacokinetics. By integrating these recent advancements and clinical insights, this review aims to provide a thorough and updated understanding of integrase inhibitors, emphasizing their evolving role in HIV treatment.

C-Terminal Domain of Integrase Binds between the Two Active Sites

HIV integrase (HIV-IN), one of three HIV enzymes, is a target for the treatment of AIDS, but the full biological assembly has been difficult to characterize, hampering inhibitor design. The recent crystallographic structures of integrase from prototype foamy virus (PFV-IN) with bound DNA were a breakthrough, revealing how viral DNA organizes two integrase dimers into a tetramer that has the two active sites appropriately spaced for insertion of the viral DNA into host DNA. The organization of domains within each PFV-IN protein chain, however, varies significantly from that found in HIV-IN structures. With the goal of identifying shared structural characteristics, the interactions among components of the PFV-IN and HIV-IN assemblies were investigated with the macromolecular docking program DOT. DOT performs an exhaustive, rigid-body search between two macromolecules. Computational docking reproduced the crystallographic interactions of the PFV-IN catalytic and N-terminal domains with viral DNA and found similar viral DNA interactions for HIV-IN. Computational docking did not reproduce the crystallographic interactions of the PFV-IN C-terminal domain (CTD). Instead, two symmetry-related positions were found for the PFV-IN CTD that indicate formation of a CTD dimer between the two active sites. Our predicted CTD dimer is consistent with cross-linking studies showing interactions of the CTD with viral DNA that appear to be blocked in the PFV-IN structures. The CTD dimer can insert two arginine-rich loops between the two bound vDNA molecules and the host DNA, a region that is unoccupied in the PFV-IN crystallographic structures. The positive potential from these two loops would alleviate the large negative potential created by the close proximity of two viral vDNA ends, helping to bring together the two active sites and assisting host DNA binding. This study demonstrates the ability of computational docking to evaluate complex crystallographic assemblies, identify interactions that are influenced by the crystal environment, and provide plausible alternatives.

Design, development and evaluation of novel dual PPARδ/PPARγ agonists

Type 2 diabetes is at epidemic proportions and thus development of novel pharmaceutical therapies for improving insulin sensitivity has become of paramount importance. The objectives of the current study were to develop novel dual PPARγ/δ agonists without the deleterious side effects associated with full PPARγ agonists. Docking simulations of 23 novel compounds within the ligand binding domain of PPARγ/δ were performed using AutoDock Vina which consistently reproduced experimental binding poses from known PPAR agonists. Comparisons were made and described with other docking programs AutoDock and Surflex-Dock (from SYBYL-X). Biological evaluation of compounds was accomplished by transcriptional promoter activity assays, quantitative PCR gene analysis for known PPARγ/δ targets as well as in vitro assays for lipid accumulation and mitochondrial biogenesis verses known PPAR agonists. We found one (compound 9) out of the 23 compounds evaluated, to be the most potent and selective dual PPARγ/δ agonist which did not display the deleterious side effects associated with full PPARγ agonists.

Mechanism of microRNA-target interaction: molecular dynamics simulations and thermodynamics analysis

MicroRNAs (miRNAs) are endogenously produced ∼21-nt riboregulators that associate with Argonaute (Ago) proteins to direct mRNA cleavage or repress the translation of complementary RNAs. Capturing the molecular mechanisms of miRNA interacting with its target will not only reinforce the understanding of underlying RNA interference but also fuel the design of more effective small-interfering RNA strands. To address this, in the present work the RNA-bound (Ago-miRNA, Ago-miRNA-target) and RNA-free Ago forms were analyzed by performing both molecular dynamics simulations and thermodynamic analysis. Based on the principal component analysis results of the simulation trajectories as well as the correlation analysis in fluctuations of residues, we discover that: 1) three important (PAZ, Mid and PIWI) domains exist in Argonaute which define the global dynamics of the protein; 2) the interdomain correlated movements are so crucial for the interaction of Ago-RNAs that they not only facilitate the relaxation of the interactions between residues surrounding the RNA binding channel but also induce certain conformational changes; and 3) it is just these conformational changes that expand the cavity of the active site and open putative pathways for both the substrate uptake and product release. In addition, by thermodynamic analysis we also discover that for both the guide RNA 5′-end recognition and the facilitated site-specific cleavage of the target, the presence of two metal ions (of Mg²⁺) plays a predominant role, and this conclusion is consistent with the observed enzyme catalytic cleavage activity in the ternary complex (Ago-miRNA-mRNA). Our results find that it is the set of arginine amino acids concentrated in the nucleotide-binding channel in Ago, instead of the conventionally-deemed seed base-paring, that makes greater contributions in stabilizing the binding of the nucleic acids to Ago.

One of the biggest surprises at the beginning of the ‘post-genome era’ was the discovery of numerous genes encoding microRNAs. The number of microRNA genes is estimated to be nearly 1% of that of protein-coding genes, which were found in genomes of such diverse organisms as Caenorhabditis elegans, Drosophila melanogaster, Arabidopsis thaliana, and Homo sapiens. Their products, tiny RNAs (miRNAs and siRNAs), are thought to bind to Argonaute (Ago) proteins and form effector complexes to direct mRNA cleavage or repress translation of complementary RNAs, during development, organogenesis, and very likely during many other processes. The cellular interactions between the miRNAs and their target RNAs associating with Ago are only beginning to be revealed, and details of this interaction mechanism at molecular level are still poorly understood. In this article we propose the possible mechanisms of miRNA-target interaction with special emphasis on their structural dynamic and thermodynamic aspects. The results of our model suggest the chemical and physical factors and effects that may be responsible for the miRNA-Ago assembly and miRNA-target recognition.

Chemistry and structure-activity relationship of the styrylquinoline-type HIV integrase inhibitors

In spite of significant progress in anti-HIV-1 therapy, current antiviral chemo-therapy still suffers from deleterious side effects and emerging drug resistance. Therefore, the development of novel antiviral drugs remains a crucial issue for the fight against AIDS. HIV-1 integrase is a key enzyme in the replication cycle of the retrovirus since it catalyzes the integration of the reverse transcribed viral DNA into the chromosomal DNA. Efforts to develop anti-integrase drugs started during the early nineties, culminating with the recent approval of Raltegravir. The discovery and the development of the styrylquinoline inhibitor class was an important step in the overall process. In this review we have described the key synthetic issues and the structure-activity relationship of this family of integrase inhibitors. Crystallographic and docking studies that shed light on their mechanism of action are also examined.

Specificity of LTR DNA recognition by a peptide mimicking the HIV-1 integrase {alpha}4 helix

HIV-1 integrase integrates retroviral DNA through 3′-processing and strand transfer reactions in the presence of a divalent cation (Mg²⁺ or Mn²⁺). The α4 helix exposed at the catalytic core surface is essential to the specific recognition of viral DNA. To define group determinants of recognition, we used a model composed of a peptide analogue of the α4 helix, oligonucleotides mimicking processed and unprocessed U5 LTR end and 5 mM Mg²⁺. Circular dichroism, fluorescence and NMR experiments confirmed the implication of the α4 helix polar/charged face in specific and non-specific bindings to LTR ends. The specific binding requires unprocessed LTR ends—i.e. an unaltered 3′-processing site CA↓GT3′—and is reinforced by Mg²⁺ (K_d decreases from 2 to 0.8 nM). The latter likely interacts with the ApG and GpT3′ steps of the 3′-processing site. With deletion of GT3′, only persists non-specific binding (K_d of 100 μM). Proton chemical shift deviations showed that specific binding need conserved amino acids in the α4 helix and conserved nucleotide bases and backbone groups at LTR ends. We suggest a conserved recognition mechanism based on both direct and indirect readout and which is subject to evolutionary pressure.

Analysis of binding parameters of HIV-1 integrase inhibitors: correlates of drug inhibition and resistance

This study undertook an exploratory data analysis of the binding parameters of HIV-1 integrase inhibitors. The study group involved inhibitors in preclinical development from the diketo acid, pyrroloquinoline and naphthyridine carboxamide families and the most advanced inhibitors Raltegravir and Elvitegravir. Distinct differences were observed in the energetics of binding between the studied classes of inhibitors that also correlated with drug resistant patterns. Quantitative-property-activity-relationships correlated experimental IC(50) values to the binding energy and the logarithm of the partition coefficient between n-octanol and water (clogP). The approach followed here serves as an improved basis for the development of 'second generation' integrase inhibitors.

Dinucleotides docking to scorpion polypeptide toxins: a molecular modeling method for protein functional site recognition

To understand the principles underlying protein folding, many molecular modeling methods are being developed for predicting functional positions. In this work, fully flexible dinucleotides d(pApA), d(pApC), d(pApG), d(pApT), d(pCpA), d(pCpC), d(pCpG), d(pCpT), d(pGpA), d(pGpC), d(pGpG), d(pGpT), d(pTpA), d(pTpC), d(pTpG), and d(pTpT) were first docked onto the surface of scorpion polypeptide toxins (LqhIT2, PDB ID:2I61) and homology modeled ANEPIII. Automated docking was able to identify sites on scorpion polypeptide toxins where favorable nucleotide interactions can occur, and those sites were in agreement with the mutation data of this protein published recently [I. Karbat, R. Kahn, L. Cohen, N. Ilan, N. Gilles, G. Corzo, O. Froy, M. Gur, G. Albrecht, S.H. Heinemann, D. Gordon, M. Gurevitz, The unique pharmacology of the scorpion alpha-like toxin Lqh3 is associated with its flexible C-tail, Febs J 274 (2007) 1918-1931]. Simulation results suggested that dinucleotides docking is a suitable molecular modeling method that could be developed for protein functional site recognition.

Identification of ssDNA aptamers specific for anti-neuroexcitation peptide III and molecular modeling studies: insights into structural interactions

Twelve ssDNA aptamers specific for a novel recombinant anti-neuroexcitation peptide (ANEPIII) were identified using the SELEX method from a 79-nucleotide ssDNA pool to purify ANEPIII in a more efficient way. To further understand the binding modes between ssDNA and ANEPIII, fully flexible dinucleotides were docked onto the homology-modeled ANEPIII structure. AutoDocking identified favorable binding sites on ANEPIII for nucleotides, which was valuable for designing more potent ligands.

Kinetic study of the HIV-1 DNA 3'-end processing

The 3'-processing of viral DNA extremities is the first step in the integration process catalysed by human immunodeficiency virus (HIV)-1 integrase (IN). This reaction is relatively inefficient and processed DNAs are usually detected in vitro under conditions of excess enzyme. Despite such experimental conditions, steady-state Michaelis-Menten formalism is often applied to calculate characteristic equilibrium/kinetic constants of IN. We found that the amount of processed product was not significantly affected under conditions of excess DNA substrate, indicating that IN has a limited turnover for DNA cleavage. Therefore, IN works principally in a single-turnover mode and is intrinsically very slow (single-turnover rate constant = 0.004 min(-1)), suggesting that IN activity is mainly limited at the chemistry step or at a stage that precedes chemistry. Moreover, fluorescence experiments showed that IN-DNA product complexes were very stable over the time-course of the reaction. Binding isotherms of IN to DNA substrate and product also indicate tight binding of IN to the reaction product. Therefore, the slow cleavage rate and limited product release prevent or greatly reduce subsequent turnover. Nevertheless, the time-course of product formation approximates to a straight line for 90 min (apparent initial velocity), but we show that this linear phase is due to the slow single-turnover rate constant and does not indicate steady-state multiple turnover. Finally, our data ruled out the possibility that there were large amounts of inactive proteins or dead-end complexes in the assay. Most of complexes initially formed were active although dramatically slow.

A three-dimensional model of the human immunodeficiency virus type 1 integration complex

While the general features of HIV-1 integrase function are understood, there is still uncertainty about the composition of the integration complex and how integrase interacts with viral and host DNA. We propose an improved model of the integration complex based on current experimental evidence including a comparison with the homologous Tn5 transposase containing bound DNA and an analysis of DNA binding sites using Goodford's GRID. Our model comprises a pair of integrase dimers, two strands of DNA to represent the viral DNA ends and a strand of bent DNA representing the host chromosome. In our model, the terminal four base pairs of each of the viral DNA strands interact with the integrase dimer providing the active site, while bases one turn away interact with a flexible loop (residues 186-194) on the second integrase dimer. We propose that residues E152, Q148 and K156 are involved in the specific recognition of the conserved CA dinucleotide and that the active site mobile loop (residues 140-149) stabilises the integration complex by acting as a barrier to separate the two viral DNA ends. In addition, the residues responsible for DNA binding in our model show a high level of amino acid conservation.

Genetic analyses of DNA-binding mutants in the catalytic core domain of human immunodeficiency virus type 1 integrase

The catalytic core domain (CCD) of human immunodeficiency virus type 1 (HIV-1) integrase (IN) harbors the enzyme active site and binds viral and chromosomal DNA during integration. Thirty-five CCD mutant viruses were constructed, paying particular attention to conserved residues in the Phe(139)-Gln(146) flexible loop and abutting Ser(147)-Val(165) amphipathic alpha helix that were implicated from previous in vitro work as important for DNA binding. Defective viruses were typed as class I mutants (specifically blocked at integration) or pleiotropic class II mutants (additional particle assembly and/or reverse transcription defects). Whereas HIV-1(P145A) and HIV-1(Q146K) grew like the wild type, HIV-1(N144K) and HIV-1(Q148L) were class I mutants, reinforcing previous results that Gln-148 is important for DNA binding and uncovering for the first time an important role for Asn-144 in integration. HIV-1(Q62K), HIV-1(H67E), HIV-1(N120K), and HIV-1(N155K) were also class I mutants, supporting findings that Gln-62 and Asn-120 interact with viral and target DNA, respectively, and suggesting similar integration-specific roles for His-67 and Asn-155. Although results from complementation analyses established that IN functions as a multimer, the interplay between active-site and CCD DNA binding functions was unknown. By using Vpr-IN complementation, we determined that the CCD protomer that catalyzes integration also preferentially binds to viral and target DNA. We additionally characterized E138K as an intramolecular suppressor of Gln-62 mutant virus and IN. The results of these analyses highlight conserved CCD residues that are important for HIV-1 replication and integration and define the relationship between DNA binding and catalysis that occurs during integration in vivo.

Docking dinucleotides to HIV-1 integrase carboxyl-terminal domain to find possible DNA binding sites

HIV-1 integrase mediates a processed viral DNA into the host genome DNA. The binding modes between HIV-1 integrase and viral/host DNA have not been clarified until now. The C-terminal domain of integrase has been found to be a DNA-binding domain. In this work, we have explored the possible DNA binding sites of dimeric C-terminal domain by docking dinucleotides to HIV-1 integrase. The docking results suggest that two symmetrical DNA-binding sites are likely located on the outside surface of the dimeric C-terminal domain and not located in the groove. Those sites are in agreement with the experimental data.

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 beta-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.

HIV-1 integrase complexes with DNA dissociate in the presence of short oligonucleotides conjugated to acridine

The human immunodeficiency virus type 1 (HIV-1) integrase is an essential enzyme in the life cycle of the virus and is therefore an attractive target for the development of new antiviral drugs. Among them, inhibitors which are capable of targeting the preassembled integrase/DNA complex are of particular interest, because they could suppress integrase activity in the context of the HIV-1 preintegration complex. Here, we study the mechanism of action of 11-mer oligonucleotides, which are efficient inhibitors of the catalytic activity of integrase, provided that they are conjugated to a hydrophobic compound, acridine. To understand the mechanism of the conjugate inhibitory action, we used a steady-state fluorescence anisotropy assay, which allowed us to study the stability of the integrase/DNA complex in various conditions. We found that oligonucleotide-acridine conjugates induced the efficient dissociation of preassembled integrase/DNA complexes. The simultaneous presence of both acridine and an oligonucleotidic moiety is required for the inhibitory activity of conjugates. However, the dissociation effect is not dependent on the oligonucleotide sequence. Finally, our results suggest that the conjugates bind directly to integrase within its complex with DNA at a site different from the viral DNA binding site.

On the Function of the 14 Å Long Internal Cavity of Histone Deacetylase-Like Protein: Implications for the Design of Histone Deacetylase Inhibitors

Histone deacetylases (HDACs) play an important role in gene transcription. Inhibitors of HDACs induce cell differentiation and suppress cell proliferation in tumor cells. AutoDock calculations of known and novel HDAC inhibitors as well as of several probe molecules to histone deacetylase-like protein (HDLP), using a modified scoring function for metalloproteins, demonstrate excellent agreement (R = 0.92) between experimental and computed binding constants. Analysis of the docked structures allows a determination of the different binding motifs in known inhibitors. Such calculations are a useful tool for the prediction of binding constants for new HDAC inhibitors. Exploration of the 14 A long internal cavity adjacent to the active site by docking of small molecular probes suggest that it plays a crucial role by accepting the cleaved acetate and releasing it at the far side of the cavity. The importance of the findings for the design of new inhibitors is discussed.

Prediction of HIV-1 Integrase/Viral DNA Interactions in the Catalytic Domain by Fast Molecular Docking

This study details the separate analyses of binding specificity of HIV-1 integrase (IN) and viral B-DNA forms through ligand-receptor docking studies by means of a fast molecular docking method. The application of solvated electrostatics with the University of Houston Brownian Dynamics Program (UHBD) and configurational sampling by the Daughter of Turnip (DOT) docking program resulted in the computation of energies of more than 113 billion configurations for each ligand-receptor docking study, a procedure considered computationally intractable a few years ago. A specific binding pattern of viral DNA to the IN catalytic domain region has been predicted as a result of these calculations. In a representative docked configuration, we observe the 3'-hydroxyl of the conserved deoxyadenosine to be close to one of the two divalent metal ions that are necessary for catalysis. A superimposition of our energy-minimized docked complex on representative structures from a molecular dynamics (MD) simulation of a crystallographically resolved IN/inhibitor complex revealed an overlap of viral DNA with the inhibitor, indicating that the bound inhibitor might operate by blocking substrate binding. The DOT docking calculation also identified a second, adjacent DNA-binding site, which we believe is the nonspecific host DNA binding site. The binding pattern predicted by DOT complements previous electrostatics, MD simulation, photo-cross-linking, and mutagenesis studies and also provides a further refinement of the IN/viral DNA binding interaction as a basis for new structure-based design efforts.

Mechanism of HIV-1 Integrase Inhibition by Styrylquinoline Derivatives in Vitro

Styrylquinoline derivatives (SQ) efficiently inhibit the 3'-processing activity of integrase (IN) with IC50 values of between 0.5 and 5 microM. We studied the mechanism of action of these compounds in vitro. First, we used steady-state fluorescence anisotropy to assay the effects of the SQ derivatives on the formation of IN-viral DNA complexes independently of the catalytic process. The IC50 values obtained in activity and DNA-binding tests were similar, suggesting that the inhibition of 3'-processing can be fully explained by the prevention of IN-DNA recognition. SQ compounds act in a competitive manner, with Ki values of between 400 and 900 nM. In contrast, SQs did not inhibit 3'-processing when IN-DNA complexes were preassembled. Computational docking followed or not by molecular dynamics using the catalytic core of HIV-1 IN suggested a competitive inhibition mechanism, which is consistent with our previous data obtained with the corresponding Rous sarcoma virus domain. Second, we used preassembled IN-preprocessed DNA complexes to assay the potency of SQs against the strand transfer reaction, independently of 3'-processing. Inhibition occurred even if the efficiency was decreased by about 5- to 10-fold. Our results suggest that two inhibitor-binding modes exist: the first one prevents the binding of the viral DNA and then the two subsequent reactions (i.e., 3'-processing and strand transfer), whereas the second one prevents the binding of target DNA, thus inhibiting strand transfer. SQ derivatives have a higher affinity for the first site, in contrast to that observed for the diketo acids, which preferentially bind to the second one.