Molecular dynamics approaches estimate the binding energy of HIV-1 integrase inhibitors and correlate with in vitro activity

Exploring Natural Alkaloids from Brazilian Biodiversity as Potential Inhibitors of the Aedes aegypti Juvenile Hormone Enzyme: A Computational Approach for Vector Mosquito Control

This study explores the potential inhibitory activity of alkaloids, a class of natural compounds isolated from Brazilian biodiversity, against the mJHBP enzyme of the Aedes aegypti mosquito. This mosquito is a significant vector of diseases such as dengue, zika, and chikungunya. The interactions between the ligands and the enzyme at the molecular level were evaluated using computational techniques such as molecular docking, molecular dynamics (MD), and molecular mechanics with generalized Born surface area (MMGBSA) free energy calculation. The findings suggest that these compounds exhibit a high binding affinity with the enzyme, as confirmed by the binding free energies obtained in the simulation. Furthermore, the specific enzyme residues that contribute the most to the stability of the complex with the compounds were identified: specifically, Tyr33, Trp53, Tyr64, and Tyr129. Notably, Tyr129 residues were previously identified as crucial in the enzyme inhibition process. This observation underscores the significance of the research findings and the potential of the evaluated compounds as natural insecticides against Aedes aegypti mosquitoes. These results could stimulate the development of new vector control agents that are more efficient and environmentally friendly.

Analogue and structure based approaches for modelling HIV-1 integrase inhibitors

A set of 220 inhibitors belonging to different structure classes and having HIV-1 integrase activity were collected along with their experimental pIC50 values. Geometries of all the inhibitors were fully optimized using B3LYP/6-31 + G(d) level of theory. These ligands were docked against 4 different HIV-1 integrase receptors (PDB IDs: 4LH5, 5KRS, 3ZSQ and 3ZSV). 30 docked poses were generated for all 220 inhibitors and ligand interaction of the first docked pose and the docked pose with the highest score were analysed. Residue GLU170 of 4LH5 receptor shows the highest number of interactions followed by ALA169, GLN168, HIS171 and ASP167 residues. Hydrogen bonding and stacking are mainly responsible for the interactions of these inhibitors with the receptor. We performed Molecular Dynamics (MD) simulation to observe the root-mean-square deviation (RMSD), for measure the average change of displacement between the atoms for a particular frame with respect to a reference and The Root Mean Square Fluctuation (RMSF) for characterization of local changes along the protein chain of the docked complexes. Analogue based models were generated to predict the pIC50 values for integrase inhibitors using various types of descriptors such as constitutional, geometrical, topological, quantum chemical and docking based descriptors. The best models were selected on the basis of statistical parameters and were validated by training and test set division. A few new inhibitors were designed on the basis of structure activity relationship and their pIC50 values were predicted using the generated models. All the designed new inhibitors a very high potential and may be used as potent inhibitors of HIV integrase. These models may be useful for further design and development of new and potent HIV integrase inhibitors.Communicated by Ramaswamy H. Sarma.

In silico identification and study of potential anti-mosquito juvenile hormone binding protein (MJHBP) compounds as candidates for dengue virus - Vector insecticides

Dengue has become a huge global health burden. It is currently recognized as the most rapidly spreading mosquito-borne viral disease. Yet, there are currently no licensed vaccines or specific therapeutics to manage the virus, thus, scaling up vector control approaches is important in controlling this viral spread. This study aimed to identify and study in silico, potential anti-mosquito compounds targeting Juvenile hormone (JH) mediated pathways via the Mosquito Juvenile Hormone Binding Protein (MJHBP). The study was implemented using series of computational methods. The query compounds included pyrethroids and those derived from ZINC and ANPDB databases using a simple pharmacophore model in Molecular Operating Environment (MOE). Molecular docking of selected compounds' library was implemented in MOE. The resultant high-score compounds were further validated by molecular dynamics simulation via Maestro 12.3 module and the respective Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) binding energies computed. The study identified compounds-pyrethroids, natural and synthetic - with high docking energy scores (ranging from 10.91-12.34 kcal/mol). On further analysis of the high-ranking (in terms of docking scores) compounds using MD simulation, the compounds - Ekeberin D4, Maesanin, Silafluofen and ZINC16919139- revealed very low binding energies (-122.99, -72.91 -104.50 and,-74.94 kcal/mol respectively), fairly stable complex and interesting interaction with JH-binding site amino acid residues on MJHBP. Further studies can explore these compounds in vitro/in vivo in the search for more efficient mosquito vector control.

Non-Nucleoside Reverse Transcriptase Inhibitors Join Forces with Integrase Inhibitors to Combat HIV

In the treatment of acquired immune deficiency syndrome (AIDS), the diarylpyrimidine (DAPY) analogs etravirine (ETR) and rilpivirine (RPV) have been widely effective against human immunodeficiency virus (HIV) variants that are resistant to other non-nucleoside reverse transcriptase inhibitors (NNRTIs). With non-inferior or improved efficacy, better safety profiles, and lower doses or pill burdens than other NNRTIs in the clinic, combination therapies including either of these two drugs have led to higher adherence than other NNRTI-containing treatments. In a separate development, HIV integrase strand transfer inhibitors (INSTIs) have shown efficacy in treating AIDS, including raltegravir (RAL), elvitegravir (EVG), cabotegravir (CAB), bictegravir (BIC), and dolutegravir (DTG). Of these, DTG and BIC perform better against a wide range of resistance mutations than other INSTIs. Nevertheless, drug-resistant combinations of mutations have begun to emerge against all DAPYs and INSTIs, attributable in part to non-adherence. New dual therapies that may promote better adherence combine ETR or RPV with an INSTI and have been safer and non-inferior to more traditional triple-drug treatments. Long-acting dual- and triple-therapies combining ETR or RPV with INSTIs are under study and may further improve adherence. Here, highly resistant emergent mutations and efficacy data on these novel treatments are reviewed. Overall, ETR or RPV, in combination with INSTIs, may be treatments of choice as long-term maintenance therapies that optimize efficacy, adherence, and safety.

Structural Insights on Retroviral DNA Integration: Learning from Foamy Viruses

Foamy viruses (FV) are retroviruses belonging to the Spumaretrovirinae subfamily. They are non-pathogenic viruses endemic in several mammalian hosts like non-human primates, felines, bovines, and equines. Retroviral DNA integration is a mandatory step and constitutes a prime target for antiretroviral therapy. This activity, conserved among retroviruses and long terminal repeat (LTR) retrotransposons, involves a viral nucleoprotein complex called intasome. In the last decade, a plethora of structural insights on retroviral DNA integration arose from the study of FV. Here, we review the biochemistry and the structural features of the FV integration apparatus and will also discuss the mechanism of action of strand transfer inhibitors.

An effective HIV-1 integrase inhibitor screening platform: Rationality validation of drug screening, conformational mobility and molecular recognition analysis for PFV integrase complex with viral DNA

As an important target for the development of novel anti-AIDS drugs, HIV-1 integrase (IN) has been widely concerned. However, the lack of a complete accurate crystal structure of HIV-1 IN greatly blocks the discovery of novel inhibitors. In this work, an effective HIV-1 IN inhibitor screening platform, namely PFV IN, was filtered from all species of INs. Next, the 40.8% similarity with HIV-1 IN, as well as the high efficiency of virtual screening and the good agreement between calculated binding free energies and experimental ones all proved PFV IN is a promising screening platform for HIV-1 IN inhibitors. Then, the molecular recognition mechanism of PFV IN by its substrate viral DNA and six naphthyridine derivatives (NRDs) inhibitors was investigated through molecular docking, molecular dynamics simulations and water-mediated interactions analyses. The functional partition of NRDs IN inhibitors could be divided into hydrophobic and hydrophilic ones, and the Mg²⁺ ions, water molecules and conserved DDE motif residues all interacted with the hydrophilic partition, while the bases in viral DNA and residues like Tyr212, Pro214 interacted with the hydrophobic one. Finally, the free energy landscape (FEL) and cluster analyses were performed to explore the molecular motion of PFV IN-DNA system. It is found that the association with NRDs inhibitors would obviously decrease the motion amplitude of PFV IN-DNA, which may be one of the most potential mechanisms of IN inhibitors. This work will provide a theoretical basis for the inhibitor design based on the structure of HIV-1 IN.

Cellulose Acetate Phthalate and Antiretroviral Nanoparticle Fabrications for HIV Pre-Exposure Prophylaxis

To adequately reduce new HIV infections, development of highly effective pre-exposure prophylaxis (PrEP) against HIV infection in women is necessary. Cellulose acetate phthalate (CAP) is a pH sensitive polymer with HIV-1 entry inhibitory properties. Dolutegravir (DTG) is an integrase strand transfer inhibitor with potent antiretroviral activity. DTG delivered in combination with CAP may significantly improve current PrEP against HIV. In the present study, the development of DTG-loaded CAP nanoparticles incorporated in thermosensitive (TMS) gel at vaginal pH 4.2 and seminal fluid pH 7.4 is presented as proof-of-concept for improved PrEP. Water–oil–in–water homogenization was used to fabricate DTG-loaded CAP nanoparticles (DTG–CAP–NPs). Size, polydispersity, and morphological analyses illustrate that DTG–CAP–NPs were smooth and spherical, ≤200 nm in size, and monodispersed with a polydispersity index PDI ≤ 0.2. The drug encapsulation (EE%) and release profile of DTG–CAP–NPs was determined by HPLC analysis. The EE% of DTG in DTG–CAP–NPs was evaluated to be ~70%. The thermal sensitivity of the TMS gel was optimized and the pH dependency was evaluated by rheological analysis. DTG release studies in TMS gel revealed that DTG–CAP–NPs were stable in TMS gel at pH 4.2 while DTG–CAP–NPs in TMS gel at pH 7.4 rapidly release DTG (≥80% release within 1 h). Cytotoxicity studies using vaginal cell lines revealed that DTG–CAP–NPs were relatively non-cytotoxic at concentration <1 µg/mL. Confocal microscopic studies illustrate that ≥98% cells retained DTG–CAP–NPs intracellularly over seven days. Antiretroviral drug loaded nanocellulose fabrications in TMS gel delivered intravaginally may enhance both microbicidal and antiretroviral drug efficacy and may present a novel option for female PrEP against HIV.

Impact of Genetic Variations in HIV-1 Tat on LTR-Mediated Transcription via TAR RNA Interaction

HIV-1 evades host defense through mutations and recombination events, generating numerous variants in an infected patient. These variants with an undiminished virulence can multiply rapidly in order to progress to AIDS. One of the targets to intervene in HIV-1 replication is the trans-activator of transcription (Tat), a major regulatory protein that transactivates the long terminal repeat promoter through its interaction with trans-activation response (TAR) RNA. In this study, HIV-1 infected patients (n = 120) from North India revealed Ser46Phe (20%) and Ser61Arg (2%) mutations in the Tat variants with a strong interaction toward TAR leading to enhanced transactivation activities. Molecular dynamics simulation data verified that the variants with this mutation had a higher binding affinity for TAR than both the wild-type Tat and other variants that lacked Ser46Phe and Ser61Arg. Other mutations in Tat conferred varying affinities for TAR interaction leading to differential transactivation abilities. This is the first report from North India with a clinical validation of CD4 counts to demonstrate the influence of Tat genetic variations affecting the stability of Tat and its interaction with TAR. This study highlights the co-evolution pattern of Tat and predominant nucleotides for Tat activity, facilitating the identification of genetic determinants for the attenuation of viral gene expression.

Oligomeric Conjugated Polyelectrolytes Display Site-Preferential Binding to an MS2 Viral Capsid

Opportunistic bacteria and viruses are a worldwide health threat prompting the need to develop new targeting modalities. A class of novel synthetic poly(phenylene ethynylene) (PPE)-based oligomeric conjugated polyelectrolytes (OPEs) have demonstrated potent wide-spectrum biocidal activity. A subset of cationic OPEs display high antiviral activity against the MS2 bacteriophage. The oligomers have been found to inactivate the bacteriophage and perturb the morphology of the MS2 viral capsid. However, details of the initial binding and interactions between the OPEs and the viruses are not well understood. In this study, we use a multiscale computational approach, including random sampling, molecular dynamics, and electronic structure calculations, to gain an understanding of the molecular-level interactions of a series of OPEs that vary in length, charge, and functional groups with the MS2 capsid. Our results show that OPEs strongly bind to the MS2 capsid protein assembly with binding energies of up to -30 kcal/mol. Free-energy analysis shows that the binding is dominated by strong van der Waals interactions between the hydrophobic OPE backbone and the capsid surface and strong electrostatic free energy contributions between the OPE charged moieties and charged residues on the capsid surface. This knowledge provides molecular-level insight into how to tailor the OPEs to optimize viral capsid disruption and increase OPE efficacy to target amphiphilic protein coats of icosahedral-based viruses.

6,7-Dihydroxyisoindolin-1-one and 7,8-Dihydroxy-3,4-Dihydroisoquinolin- 1(2H)-one Based HIV-1 Integrase Inhibitors

Integrase (IN) is an essential viral enzyme required for HIV-1 replication, which has been targeted by anti-AIDS therapeutics. Integrase strand transfer inhibitors (INSTIs) represent a new class of antiretroviral agents developed for the treatment of HIV-1 infections. Important structural features that are shared by many INSTIs include a coplanar arrangement of three heteroatoms that chelate two catalytic Mg(2+) ions in the IN active site and a linked halophenyl ring that binds in the hydrophobic pocket formed by the complex of IN with viral DNA. We recently reported bicyclic 6,7-dihydroxyoxoisoindolin-1-one-based IN inhibitors. In the current study, we modified these inhibitors in three ways. First, we increased the spacer length between the metalchelating triad and the halophenyl group. Second, we replaced the indoline [5,6] bicycle with a fused dihydroxyisoquinolinones [6,6] ring system. Finally, we prepared bis-6,7-dihydroxyisoindolin-1-one-4-sulfonamides as dimeric HIV-1 IN inhibitors. These new analogues showed low micromolar inhibitory potency in in vitro HIV-1 integrase assays.

Selectivity for strand-transfer over 3'-processing and susceptibility to clinical resistance of HIV-1 integrase inhibitors are driven by key enzyme-DNA interactions in the active site

Integrase strand transfer inhibitors (INSTIs) are highly effective against HIV infections. Co-crystal structures of the prototype foamy virus intasome have shown that all three FDA-approved drugs, raltegravir (RAL), elvitegravir and dolutegravir (DTG), act as interfacial inhibitors during the strand transfer (ST) integration step. However, these structures give only a partial sense for the limited inhibition of the 3′-processing reaction by INSTIs and how INSTIs can be modified to overcome drug resistance, notably against the G140S-Q148H double mutation. Based on biochemical experiments with modified oligonucleotides, we demonstrate that both the viral DNA +1 and −1 bases, which flank the 3′-processing site, play a critical role for 3′-processing efficiency and inhibition by RAL and DTG. In addition, the G140S-Q148H (SH) mutant integrase, which has a reduced 3′-processing activity, becomes more active and more resistant to inhibition of 3′-processing by RAL and DTG in the absence of the −1 and +1 bases. Molecular modeling of HIV-1 integrase, together with biochemical data, indicate that the conserved residue Q146 in the flexible loop of HIV-1 integrase is critical for productive viral DNA binding through specific contacts with the virus DNA ends in the 3′-processing and ST reactions. The potency of integrase inhibitors against 3′-processing and their ability to overcome resistance is discussed.

Effect of HIV-1 Subtype C integrase mutations implied using molecular modeling and docking data

The degree of sequence variation in HIV-1 integrase genes among infected patients and their impact on clinical response to Anti retroviral therapy (ART) is of interest. Therefore, we collected plasma samples from 161 HIV-1 infected individuals for subsequent integrase gene amplification (1087 bp). Thus, 102 complete integrase gene sequences identified as HIV-1 subtype-C was assembled. This sequence data was further used for sequence analysis and multiple sequence alignment (MSA) to assess position specific frequency of mutations within pol gene among infected individuals. We also used biophysical geometric optimization technique based molecular modeling and docking (Schrodinger suite) methods to infer differential function caused by position specific sequence mutations towards improved inhibitor selection. We thus identified accessory mutations (usually reduce susceptibility) leading to the resistance of some known integrase inhibitors in 14% of sequences in this data set. The Stanford HIV-1 drug resistance database provided complementary information on integrase resistance mutations to deduce molecular basis for such observation. Modeling and docking analysis show reduced binding by mutants for known compounds. The predicted binding values further reduced for models with combination of mutations among subtype C clinical strains. Thus, the molecular basis implied for the consequence of mutations in different variants of integrase genes of HIV-1 subtype C clinical strains from South India is reported. This data finds utility in the design, modification and development of a representative yet an improved inhibitor for HIV-1 integrase.

Comparison of Newly Assembled Full Length HIV-1 Integrase With Prototype Foamy Virus Integrase: Structure-Function Prospective

Background

Drug design against human immunodeficiency virus type 1 (HIV-1) integrase through its mechanistic study is of great interest in the area in biological research. The main obstacle in this area is the absence of the full-length crystal structure for HIV-1 integrase to be used as a model. A complete structure, similar to HIV-1 of a prototype foamy virus integrase in complex with DNA, including all conservative residues, is available and has been extensively used in recent investigations.

Objectives

The aim of this study was to determine whether the above model is precisely representative of HIV-1 integrase. This would critically determine the success of any designed drug using the model in deactivation of integrase and AIDS treatment.

Materials and Methods

Primarily, a new structure for HIV-1 was constructed, using a crystal structure of prototype foamy virus as the starting structure. The constructed structure of HIV-1 integrase was simultaneously simulated with a prototype foamy virus integrase on a separate occasion.

Results

Our results indicate that the HIV-1 system behaves differently from the prototype foamy virus in terms of folding, hydration, hydrophobicity of binding site and stability.

Conclusions

Based on our findings, we can conclude that HIV-1 integrase is vastly different from the prototype foamy virus integrase and does not resemble it, and the modeling output of the prototype foamy virus simulations could not be simply generalized to HIV-1 integrase. Therefore, our HIV-1 model seems to be more representative and more useful for future research.

Development of HIV-1 integrase inhibitors: recent molecular modeling perspectives

Of the three viral enzymes essential to HIV replication, HIV-1 integrase (IN) is gaining popularity as a target for the antiviral therapy of AIDS. Substantial work focusing on IN has been done over the past three decades, which has facilitated and led to the approval of three drugs. Here, we discuss in detail the development of IN inhibitors between January 2012 and May 2014, with a particular focus on molecular simulation. We highlight controversial aspects of computational drug design and refer to alternative practices where appropriate. The analysis of these computational approaches provides some useful clues to the possible future discovery of novel IN inhibitors.

In silico modification of oseltamivir as neuraminidase inhibitor of influenza A virus subtype H1N1

This research focused on the modification of the functional groups of oseltamivir as neuraminidase inhibitor against influenza A virus subtype H1N1. Interactions of three of the best ligands were evaluated in the hydrated state using molecular dynamics simulation at two different temperatures. The docking result showed that AD3BF2D ligand (N-[(1S,6R)-5-amino-5-{[(2R,3S,4S)-3,4-dihydroxy-4-(hydroxymethyl) tetrahydrofuran-2-yl]oxy}-4-formylcyclohex-3-en-1-yl]acetamide-3-(1-ethylpropoxy)-1-cyclohexene-1-carboxylate) had better binding energy values than standard oseltamivir. AD3BF2D had several interactions, including hydrogen bonds, with the residues in the catalytic site of neuraminidase as identified by molecular dynamics simulation. The results showed that AD3BF2D ligand can be used as a good candidate for neuraminidase inhibitor to cope with influenza A virus subtype H1N1.

Development of a receptor model for efficient in silico screening of HIV-1 integrase inhibitors

Integrase (IN) is a key viral enzyme for the replication of the type-1 human immunodeficiency virus (HIV-1), and as such constitutes a relevant therapeutic target for the development of anti-HIV agents. However, the lack of crystallographic data of HIV IN complexed with the corresponding viral DNA has historically hindered the application of modern structure-based drug design techniques to the discovery of new potent IN inhibitors (INIs). Consequently, the development and validation of reliable HIV IN structural models that may be useful for the screening of large databases of chemical compounds is of particular interest. In this study, four HIV-1 IN homology models were evaluated respect to their capability to predict the inhibition potency of a training set comprising 36 previously reported INIs with IC50 values in the low nanomolar to the high micromolar range. Also, 9 inactive structurally related compounds were included in this training set. In addition, a crystallographic structure of the IN-DNA complex corresponding to the prototype foamy virus (PFV) was also evaluated as structural model for the screening of inhibitors. The applicability of high throughput screening techniques, such as blind and ligand-guided exhaustive rigid docking was assessed. The receptor models were also refined by molecular dynamics and clustering techniques to assess protein sidechain flexibility and solvent effect on inhibitor binding. Among the studied models, we conclude that the one derived from the X-ray structure of the PFV integrase exhibited the best performance to rank the potencies of the compounds in the training set, with the predictive power being further improved by explicitly modeling five water molecules within the catalytic side of IN. Also, accounting for protein sidechain flexibility enhanced the prediction of inhibition potencies among the studied compounds. Finally, an interaction fingerprint pattern was established for the fast identification of potent IN inhibitors. In conclusion, we report an exhaustively validated receptor model if IN that is useful for the efficient screening of large chemical compounds databases in the search of potent HIV-1 IN inhibitors.

Impact of resistance mutations on inhibitor binding to HIV-1 integrase

HIV-1 integrase (IN) is essential for HIV-1 replication, catalyzing two key reaction steps termed 3' processing and strand transfer. Therefore, IN has become an important target for antiviral drug discovery. However, mutants have emerged, such as E92Q/N155H and G140S/Q148H, which confer resistance to raltegravir (RAL), the first IN strand transfer inhibitor (INSTI) approved by the FDA, and to the recently approved elvitegravir (EVG). To gain insights into the molecular mechanisms of ligand binding and drug resistance, we performed molecular dynamics (MD) simulations of homology models of the HIV-1 IN and four relevant mutants complexed with viral DNA and RAL. The results show that the structure and dynamics of the 140s' loop, comprising residues 140 to 149, are strongly influenced by the IN mutations. In the simulation of the G140S/Q148H double mutant, we observe spontaneous dissociation of RAL from the active site, followed by an intrahelical swing-back of the 3'-OH group of nucleotide A17, consistent with the experimental observation that the G140S/Q148H mutant exhibits the highest resistance to RAL compared to other IN mutants. An important hydrogen bond between residues 145 and 148 is present in the wild-type IN but not in the G140S/Q148H mutant, accounting for the structural and dynamical differences of the 140s' loop and ultimately impairing RAL binding in the double mutant. End-point free energy calculations that broadly capture the experimentally known RAL binding profiles elucidate the contributions of the 140s' loop to RAL binding free energies and suggest possible approaches to overcoming drug resistance.

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.

Biochemical analysis of the role of G118R-linked dolutegravir drug resistance substitutions in HIV-1 integrase

Drug resistance mutations (DRMs) have been reported for all currently approved anti-HIV drugs, including the latest integrase strand transfer inhibitors (INSTIs). We previously used the new INSTI dolutegravir (DTG) to select a G118R integrase resistance substitution in tissue culture and also showed that secondary substitutions emerged at positions H51Y and E138K. Now, we have characterized the impact of the G118R substitution, alone or in combination with either H51Y or E138K, on 3' processing and integrase strand transfer activity. The results show that G118R primarily impacted the strand transfer step of integration by diminishing the ability of integrase-long terminal repeat (LTR) complexes to bind target DNA. The addition of H51Y and E138K to G118R partially restored strand transfer activity by modulating the formation of integrase-LTR complexes through increasing LTR DNA affinity and total DNA binding, respectively. This unique mechanism, in which one function of HIV integrase partially compensates for the defect in another function, has not been previously reported. The G118R substitution resulted in low-level resistance to DTG, raltegravir (RAL), and elvitegravir (EVG). The addition of either of H51Y or E138K to G118R did not enhance resistance to DTG, RAL, or EVG. Homology modeling provided insight into the mechanism of resistance conferred by G118R as well as the effects of H51Y or E138K on enzyme activity. The G118R substitution therefore represents a potential avenue for resistance to DTG, similar to that previously described for the R263K substitution. For both pathways, secondary substitutions can lead to either diminished integrase activity and/or increased INSTI susceptibility.

Computational design of a full-length model of HIV-1 integrase: modeling of new inhibitors and comparison of their calculated binding energies with those previously studied

A full-length model of integrase (IN) of the human immunodeficiency virus type 1 (HIV-1) was constructed based on the distinctly resolved X-ray crystal structures of its three domains, named N-terminal, catalytic core and C-terminal. Thirty-one already known inhibitors with varieties of structural differences as well as nine newly tested ones were docked into the catalytic core. The molecular dynamic (MD) and binding properties of these complexes were obtained by MD calculations. The binding energies calculated by molecular mechanic/Poisson Boltzmann solvation area were significantly correlationed with available IC50. Four inhibitors including two newly designed were also docked into the full-length model and their MD behaviors and binding properties were calculated. It was found that one of the newly designed compounds forms a better complex with HIV-1 IN compared to the rest including raltegravir. MD calculations were performed with AMBER suite of programs using ff99SB force field for the proteins and the general Amber force field for the ligands. In conclusion, the results have produced a promising standpoint not only in the construction of the full-length model but also in development of new drugs against it. However, the role of multimer formation and the involvement of DNAs, and their subsequent effect on the complexation and inhibition, are required to arrive at a conclusive decision.

Molecular dynamics simulation of complex Histones Deacetylase (HDAC) Class II Homo Sapiens with suberoylanilide hydroxamic acid (SAHA) and its derivatives as inhibitors of cervical cancer

Cervical cancer is second most common cancer in woman worldwide. Cervical cancer caused by human papillomavirus (HPV) oncogene. Inhibition of histone deacetylase (HDAC) activity has been known as a potential strategy for cancer therapy. SAHA is an HDAC inhibitor that has been used in cancer therapy but still has side effects. SAHA modification proposed to minimize side effects. Triazole attachment on the chain of SAHA has been known to enhance the inhibition ability of SAHA and less toxic. In this study, it will be carried out with molecular dynamic simulations of SAHA modifications consisting ligand 1a, 2a and, 2c to interact with six HDAC in hydrated conditions. To all six HDAC Class II, performed docking with SAHA and a modified inhibitor. The docking results were then carried out molecular dynamics simulations to determine the inhibitor affinities in hydrated conditions. The molecular dynamic simulations results show better affinities of ligand 2c with HDAC 4, 6, and 7 than SAHA itself, and good affinity was also shown by ligand 2a and 1c on HDAC 5 and 9. The results of this study can be a reference to obtain better inhibitors.

Raltegravir flexibility and its impact on recognition by the HIV-1 IN targets

HIV-1 IN is a pertinent target for the development of AIDS chemotherapy. The first IN-specific inhibitor approved for the treatment of HIV/AIDS, RAL, was designed to block the ST reaction. We characterized the structural and conformational features of RAL and its recognition by putative HIV-1 targets - the unbound IN, the vDNA, and the IN•vDNA complex - mimicking the IN states over the integration process. RAL binding to the targets was studied by performing an extensive sampling of the inhibitor conformational landscape and by using four different docking algorithms: Glide, Autodock, VINA, and SurFlex. The obtained data evidenced that: (i) a large binding pocket delineated by the active site and an extended loop in the unbound IN accommodates RAL in distinct conformational states all lacking specific interactions with the target; (ii) a well-defined cavity formed by the active site, the vDNA, and the shortened loop in the IN•vDNA complex provide a more optimized inhibitor binding site in which RAL chelates Mg(2+) cations; (iii) a specific recognition between RAL and the unpaired cytosine of the processed DNA is governed by a pair of strong H-bonds similar to those observed in DNA base pair G-C. The identified RAL pose at the cleaved vDNA shed light on a putative step of RAL inhibition mechanism. This modeling study indicates that the inhibition process may include as a first step RAL recognition by the processed vDNA bound to a transient intermediate IN state, and thus provides a potentially promising route to the design of IN inhibitors with improved affinity and selectivity.

A homology model of HIV-1 integrase and analysis of mutations designed to test the model

Although there are structures of the different domains of human immunodeficiency virus type 1 (HIV-1) integrase (IN), there is no structure of the entire protein. The recently determined crystal structures of the prototype foamy virus (PFV) IN tetramer, in complexes with viral DNA, led to the generation of models of full-length HIV-1 IN. These models were generated, in part, by superimposing the structures of the domains of HIV-1 IN onto the structure of full-length PFV IN. We developed a model for HIV-1 IN-based solely on its sequence alignment with PFV IN-that differs in several ways from the previous models. Specifically, in our model, the junction between the catalytic core domain and C-terminal domain adopts a helix-loop-helix motif that is similar to the corresponding segment of PFV IN and differs from the crystal structures of these two HIV-1 IN domains. The alignment of residues in the C-terminal domain also differs from the previous models. Our model can be used to explain the phenotype of previously published HIV-1 IN mutants. We made additional mutants, and the behavior of these new mutants provides additional support for the model.

Activities, crystal structures, and molecular dynamics of dihydro-1H-isoindole derivatives, inhibitors of HIV-1 integrase

On the basis of a series of lactam and phthalimide derivatives that inhibit HIV-1 integrase, we developed a new molecule, XZ-259, with biochemical and antiviral activities comparable to raltegravir. We determined the crystal structures of XZ-259 and four other derivatives in complex with the prototype foamy virus intasome. The compounds bind at the integrase-Mg(2+)-DNA interface of the integrase active site. In biochemical and antiviral assays, XZ-259 inhibits raltegravir-resistant HIV-1 integrases harboring the Y143R mutation. Molecular modeling is also presented suggesting that XZ-259 can bind in the HIV-1 intasome with its dimethyl sulfonamide group adopting two opposite orientations. Molecular dynamics analyses of the HIV-1 intasome highlight the importance of the viral DNA in drug potency.

Structural-Functional Analysis of 2,1,3-Benzoxadiazoles and Their N-oxides As HIV-1 Integrase Inhibitors

Human immunodeficiency virus type 1 integrase is one of the most attractive targets for the development of anti-HIV-1 inhibitors. The capacity of a series of 2,1,3-benzoxadiazoles (benzofurazans) and their N-oxides (benzofuroxans) selected using the PASS software to inhibit the catalytic activity of HIV-1 integrase was studied in the present work. Only the nitro-derivatives of these compounds were found to display inhibitory activity. The study of the mechanism of inhibition by nitro-benzofurazans/benzofuroxans showed that they impede the substrate DNA binding at the integrase active site. These inhibitors were also active against integrase mutants resistant to raltegravir, which is the first HIV-1 integrase inhibitor approved for clinical use. The comparison of computer-aided estimations of the pharmacodynamic and pharmacokinetic properties of the compounds studied and raltegravir led us to conclude that these compounds show promise and need to be further studied as potential HIV-1 integrase inhibitors.

HIV integrase inhibitors: 20-year landmark and challenges

Since the discovery of HIV as the cause for AIDS 30 years ago, major progress has been made, including the discovery of drugs that now control the disease. Here, we review the integrase (IN) inhibitors from the discovery of the first compounds 20 years ago to the approval of two highly effective IN strand transfer inhibitors (INSTIs), raltegravir (Isentress) and elvitegravir (Stribild), and the promising clinical activity of dolutegravir. After summarizing the molecular mechanism of action of the INSTIs as interfacial inhibitors, we discuss the remaining challenges. Those include: overcoming resistance to clinical INSTIs, long-term safety of INSTIs, cost of therapy, place of the INSTIs in prophylactic treatments, and the development of new classes of inhibitors (the LEDGINs) targeting IN outside its catalytic site. We also discuss the role of chromatin and host DNA repair factor for the completion of integration.

From base pair to bedside: molecular simulation and the translation of genomics to personalized medicine

Despite the promises made that genomic sequencing would transform therapy by introducing a new era of personalized medicine, relatively few tangible breakthroughs have been made. This has led to the recognition that complex interactions at multiple spatial, temporal, and organizational levels may often combine to produce disease. Understanding this complexity requires that existing and future models are used and interpreted within a framework that incorporates knowledge derived from investigations at multiple levels of biological function. It also requires a computational infrastructure capable of dealing with the vast quantities of data generated by genomic approaches. In this review, we discuss the use of molecular modeling to generate quantitative and qualitative insights at the smallest scales of the systems biology hierarchy, how it can play an important role in the development of a systems understanding of disease and in the application of such knowledge to help discover new therapies and target existing ones on a personal level.

Peptides that inhibit HIV-1 integrase by blocking its protein-protein interactions

HIV-1 integrase (IN) is one of the key enzymes in the viral replication cycle. It mediates the integration of viral cDNA into the host cell genome. IN activity requires interactions with several viral and cellular proteins, as well as IN oligomerization. Inhibition of IN is an important target for the development of anti-HIV therapies, but there is currently only one anti-HIV drug used in the clinic that targets IN. Several other small-molecule anti-IN drug leads are either undergoing clinical trials or in earlier stages of development. These molecules specifically inhibit one of the IN-mediated reactions necessary for successful integration. However, small-molecule inhibitors of protein-protein interactions are difficult to develop. In this review, we focus on peptides that inhibit IN. Peptides have advantages over small-molecule inhibitors of protein-protein interactions: they can mimic the structures of the binding domains within proteins, and are large enough to competitively inhibit protein-protein interactions. The development of peptides that bind IN and inhibit its protein-protein interactions will increase our understanding of the IN mode of action, and lead to the development of new drug leads, such as small molecules derived from these peptides, for better anti-HIV therapy.

Unprocessed viral DNA could be the primary target of the HIV-1 integrase inhibitor raltegravir

Integration of HIV DNA into host chromosome requires a 3'-processing (3'-P) and a strand transfer (ST) reactions catalyzed by virus integrase (IN). Raltegravir (RAL), commonly used in AIDS therapy, belongs to the family of IN ST inhibitors (INSTIs) acting on IN-viral DNA complexes (intasomes). However, studies show that RAL fails to bind IN alone, but nothing has been reported on the behaviour of RAL toward free viral DNA. Here, we assessed whether free viral DNA could be a primary target for RAL, assuming that the DNA molecule is a receptor for a huge number of pharmacological agents. Optical spectroscopy, molecular dynamics and free energy calculations, showed that RAL is a tight binder of both processed and unprocessed LTR (long terminal repeat) ends. Complex formation involved mainly van der Waals forces and was enthalpy driven. Dissociation constants (Kds) revealed that RAL affinity for unbound LTRs was stronger than for bound LTRs. Moreover, Kd value for binding of RAL to LTRs and IC50 value (half concentration for inhibition) were in same range, suggesting that RAL binding to DNA and ST inhibition are correlated events. Accommodation of RAL into terminal base-pairs of unprocessed LTR is facilitated by an extensive end fraying that lowers the RAL binding energy barrier. The RAL binding entails a weak damping of fraying and correlatively of 3'-P inhibition. Noteworthy, present calculated RAL structures bound to free viral DNA resemble those found in RAL-intasome crystals, especially concerning the contacts between the fluorobenzyl group and the conserved 5'C(4)pA(3)3' step. We propose that RAL inhibits IN, in binding first unprocessed DNA. Similarly to anticancer drug poisons acting on topoisomerases, its interaction with DNA does not alter the cut, but blocks the subsequent joining reaction. We also speculate that INSTIs having viral DNA rather IN as main target could induce less resistance.

The elvitegravir Quad pill: the first once-daily dual-target anti-HIV tablet

Anti-HIV combination therapies in a single formulation currently target only HIV-1 reverse transcriptase via two different mechanisms of action by associating a nucleoside and a non-nucleoside reverse transcriptase inhibitor. These combination therapies are therefore referred to as multi-class combination products. The elvitegravir Quad pill (Gilead Sciences), when approved by the Food and Drug Administration for the treatment of HIV/AIDS, will become the first once-daily dual-target anti-HIV tablet. This "4 in 1" tablet targets HIV-1 integrase by elvitegravir boosted by the pharmaco-enhancer cobicistat and HIV-1 reverse transcriptase by the two nucleoside reverse transcriptase inhibitors emtricitabine + tenofovir disoproxil fumarate. A second pill referred to as the 572-Trii pill (Shionogi-ViiV Healthcare, LLC), also based on the dual inhibition of integrase and reverse transcriptase, is currently in late-phase clinical trials. The availability of these novel once-daily anti-HIV tablets will improve treatment adherence and offer new perspective for patient failing existing antiviral regimens.

Novel therapeutic strategies targeting HIV integrase

Integration of the viral genome into host cell chromatin is a pivotal and unique step in the replication cycle of retroviruses, including HIV. Inhibiting HIV replication by specifically blocking the viral integrase enzyme that mediates this step is an obvious and attractive therapeutic strategy. After concerted efforts, the first viable integrase inhibitors were developed in the early 2000s, ultimately leading to the clinical licensure of the first integrase strand transfer inhibitor, raltegravir. Similarly structured compounds and derivative second generation integrase strand transfer inhibitors, such as elvitegravir and dolutegravir, are now in various stages of clinical development. Furthermore, other mechanisms aimed at the inhibition of viral integration are being explored in numerous preclinical studies, which include inhibition of 3' processing and chromatin targeting. The development of new clinically useful compounds will be aided by the characterization of the retroviral intasome crystal structure. This review considers the history of the clinical development of HIV integrase inhibitors, the development of antiviral drug resistance and the need for new antiviral compounds.