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

Unexpected Dynamic Binding May Rescue the Binding Affinity of Rivaroxaban in a Mutant of Coagulation Factor X

A novel coagulation factor X (FX) Tyr319Cys mutation (Y99C as chymotrypsin numbering) was identified in a patient with severe bleeding. Unlike the earlier reported Y99A mutant, this mutant can bind and cleave its specific chromogenetic substrate at a normal level, suggesting an intact binding pocket. Here, using molecular dynamics simulations and MM-PBSA calculations on a FX-rivaroxaban (RIV) complex, we confirmed a much stronger binding of RIV in Y99C than in Y99A on a molecular level, which is actually the average result of multiple binding poses in dynamics. Detailed structural analyses also indicated the moderate flexibility of the 99-loop and the importance of the flexible side chain of Trp215 in the different binding poses. This case again emphasizes that binding of ligands may not only be a dynamic process but also a dynamic state, which is often neglected in drug design and screening based on static X-ray structures. In addition, the computational results somewhat confirmed our hypothesis on the activated Tyr319Cys FX (Y99C FXa) with an impaired procoagulant function to bind inhibitors of FXa and to be developed into a potential reversal agent for novel oral anticoagulants (NOAC).

Molecular dynamic simulations to investigate the structural impact of known drug resistance mutations on HIV-1C Integrase-Dolutegravir binding

Resistance associated mutations (RAMs) threaten the long-term success of combination antiretroviral therapy (cART) outcomes for HIV-1 treatment. HIV-1 Integrase (IN) strand transfer inhibitors (INSTIs) have proven to be a viable option for highly specific HIV-1 therapy. The INSTI, Dolutegravir is recommended by the World Health Organization for use as first-line cART. This study aims to understand how RAMs affect the stability of IN, as well as the binding of the drug Dolutegravir to the catalytic pocket of the protein. A homology model of HIV-1 subtype C IN was successfully constructed and validated. The site directed mutator webserver was used to predict destabilizing and/or stabilizing effects of known RAMs while FoldX confirmed any changes in protein energy upon introduction of mutation. Also, interaction analysis was performed between neighbouring residues. Three mutations known to be associated with Raltegravir, Elvitegravir and Dolutegravir resistance were selected; E92Q, G140S and Y143R, for molecular dynamics simulations. The structural quality assessment indicated high reliability of the HIV-1C IN tetrameric structure, with more than 90% confidence in modelled regions. Change in free energy for the three mutants indicated different effects, while simulation analysis showed G140S to have the largest affect on protein stability and flexibility. This was further supported by weaker non-bonded pairwise interaction energy and binding free energy values between the drug DTG and E92Q, Y143R and G140S mutants suggesting reduced binding affinity, as indicated by interaction analysis in comparison to the WT. Our findings suggest the G140S mutant has the strongest effect on the HIV-1C IN protein structure and Dolutegravir binding. To the best of our knowledge, this is the first study that uses the consensus wild type HIV-1C IN sequence to build an accurate 3D model to understand the effect of three known mutations on DTG drug binding in a South Africa context.

Spectrometric and computational studies of the binding of HIV-1 integrase inhibitors to viral DNA extremities

Three integrase strand transfer inhibitors are in intensive clinical use, raltegravir (RAL), elvitegravir (EVG) and dolutegravir (DTG). The onset of integrase resistance mutations limits their therapeutic efficiency. As put forth earlier, the drug affinity for the intasome could be improved by targeting preferentially the retroviral nucleobases, which are little, if at all, mutation-prone. We report experimental results of anisotropy fluorescence titrations of viral DNA by these three drugs. These show the DTG > EVG > RAL ranking of their inhibitory activities of the intasome to correspond to that of their free energies of binding, ∆Gs, to retroviral DNA, and that such a ranking is only governed by the binding enthalpies, ∆H, the entropy undergoing marginal variations. We sought whether this ranking might be reproduced through quantum chemistry (QC) Density Functional Theory calculations of intermolecular interaction energies between simplified models consisting of sole halobenzene ring and the highly conserved retroviral nucleobases G4 and C16. These calculations showed that binding of EVG has a small preference over DTG, while RAL ranked third. This indicates that additional interactions of the diketoacid parts of the drugs with DNA could be necessary to further enable preferential binding of DTG. The corresponding ∆Etotvalues computed with a polarizable molecular mechanics/dynamics procedure, Sum of Interactions Between Fragments Ab initio computed (SIBFA), showed good correlations with this ∆E(QC) ranking. These validations are an important step toward the use of polarizable molecular dynamics simulations on DTG or EVG derivatives in their complexes with the complete intasome, an application now motivated and enabled by the advent of currently developed and improved massively parallel software.

An unexpected dynamic binding mode between coagulation factor X and Rivaroxaban reveals importance of flexibility in drug binding

Rivaroxaban (RIV) is a direct oral anticoagulant (DOAC) targeting activated coagulation factor X (FXa). An earlier study reported the F174A mutant of FXa resistant to a RIV-like inhibitor, Apixaban. In current study, the detailed molecular mechanism of the resistance has been explored by molecular dynamics simulations on the impaired interactions between RIV and FXa in the damaged S4 pocket of F174A mutant. Besides, an unexpected relative stable binding mode of S1'S1 was revealed, which required dynamic motions of Gln192 and Gln61 to allow the morpholinone moiety of RIV to shift into the S1' pocket and form strong interactions. These dynamic motions of RIV and critical residues might be important in drug design for direct inhibitors of coagulation factors.

Prediction of the binding mode and resistance profile for a dual-target pyrrolyl diketo acid scaffold against HIV-1 integrase and reverse-transcriptase-associated ribonuclease H

The rapid emergence of drug-resistant variants is one of the most common causes of highly active antiretroviral therapeutic (HAART) failure in patients infected with HIV-1. Compared with the existing HAART, the recently developed pyrrolyl diketo acid scaffold targeting both HIV-1 integrase (IN) and reverse transcriptase-associated ribonuclease H (RNase H) is an efficient approach to counteract the failure of anti-HIV treatment due to drug resistance. However, the binding mode and potential resistance profile of these inhibitors with important mechanistic principles remain poorly understood. To address this issue, an integrated computational method was employed to investigate the binding mode of inhibitor JMC6F with HIV-1 IN and RNase H. By using per-residue binding free energy decomposition analysis, the following residues: Asp64, Thr66, Leu68, Asp116, Tyr143, Gln148 and Glu152 in IN, Asp443, Glu478, Trp536, Lys541 and Asp549 in RNase H were identified as key residues for JMC6F binding. And then computational alanine scanning was carried to further verify the key residues. Moreover, the resistance profile of the currently known major mutations in HIV-1 IN and 2 mutations in RNase H against JMC6F was predicted by in silico mutagenesis studies. The results demonstrated that only three mutations in HIV-1 IN (Y143C, Q148R and N155H) and two mutations in HIV-1 RNase H (Y501R and Y501W) resulted in a reduction of JMC6F potency, thus indicating their potential role in providing resistance to JMC6F. These data provided important insights into the binding mode and resistance profile of the inhibitors with a pyrrolyl diketo acid scaffold in HIV-1 IN and RNase H, which would be helpful for the development of more effective dual HIV-1 IN and RNase H inhibitors.

Use multiscale simulation to explore the effects of the homodimerizations between different conformation states on the activation and allosteric pathway for the μ-opioid receptor

Recently, oligomers of G-protein coupled receptors (GPCRs) have been an important topic in the GPCR fields. However, knowledge about their structures and activation mechanisms is very limited due to the absence of crystal structures reported. In this work, we used multiscale simulations to study the effects of homodimerization between different conformation states on their activation, dynamic behaviors, and allosteric communication pathways for μ-OR. The results indicated that the dimerization of one inactive monomer with either one inactive monomer or one active one could enhance its constitutive activation. However, the conformation state of the other protomer (e.g., active or inactive) can influence the activated extent. The dimerization between the two inactive protomers leads to a negative cooperativity for their activation, which should contribute to the asymmetric activation of GPCR dimers observed in some experiments. On the other hand, for the active monomer, its dimerization with one inactive receptor could alleviate its deactivation, whereby negative and positive cooperativities can be observed between the two subunits of the dimer, depending on the different regions. Observations from protein structure network (PSN) analysis indicated that the dimerization of one inactive monomer with one active one would cause a significant drop in the number of main pathways from the ligand binding pocket to the G-protein coupled region for the inactive protomer, while the impact is minor for the active protomer. But, for the active monomer or the inactive one, its dimerization with one inactive monomer would significantly change the types of residues participating in the pathway with the highest frequency.

Exploring the mechanism of F282L mutation-caused constitutive activity of GPCR by a computational study

G-protein-coupled receptors (GPCRs) are important drug targets and generally activated by ligands. However, some experiments found that GPCRs also give rise to constitutive activity through some mutations (viz., CAM), which are usually associated with different kinds of diseases. However, the mechanisms of CAMs and their roles in interactions with drug-ligands are unclear in experiments. Herein, we used microsecond molecular dynamics simulations to study the effect of one important F282L mutation on β2AR in order to address the questions above. With the aid of principle component and correlation analysis, our results revealed that the F282L mutation could increase the instability of the overall structure, increase the dramatic fluctuations of NPxxY and extracellular loops, and decrease restraint of the helices through weakening interhelical H-bonding and correlations between residues, which could partly contribute to the constitutive activity reported by the experiments. The observations from the protein structure network (PSN) analysis indicate that the mutant exhibits less information flow than the wild β2AR and weakens the role of TM5 and TM6 in the signal transmission, but it enhances the impact of TM3 on the orthosteric pathway and TM4 on the allosteric one. In addition, the results from the virtual screening reveal that the mutant prefers to select agonists rather than antagonists, similar to the active state but opposite of the inactive state, further confirming that the F282L mutation advances the activation of β2AR. Our observations provide valuable information for understanding the mechanism of the mutation-caused constitutive activity of GPCR and related drug-design.

Secondary structure propensity and chirality of the amyloidophilic peptide p5 and its analogues impacts ligand binding - In vitro characterization

BACKGROUND

Polybasic helical peptides, such as peptide p5, bind human amyloid extracts and synthetic amyloid fibrils. When radiolabeled, peptide p5 has been shown to specifically bind amyloid in vivo thereby allowing imaging of the disease. Structural requirements for heparin and amyloid binding have been studied using analogs of p5 that modify helicity and chirality.

METHODS

Peptide-ligand interactions were studied using CD spectroscopy and solution-phase binding assays with radiolabeled p5 analogs. The interaction of a subset of peptides was further studied by using molecular dynamics simulations.

RESULTS

Disruption of the peptide helical structure reduced peptide binding to heparin and human amyloid extracts. The all-D enantiomer and the β-sheet-structured peptide bound all substrates as well as, or better than, p5. The interaction of helical and β-sheet structured peptides with Aβ fibrils was modeled and shown to involve both ionic and non-ionic interactions.

CONCLUSIONS

The α-helical secondary structure of peptide p5 is important for heparin and amyloid binding; however, helicity is not an absolute requirement as evidenced by the superior reactivity of a β-sheet peptide. The differential binding of the peptides with heparin and amyloid fibrils suggests that these molecular interactions are different. The all-D enantiomer of p5 and the β-sheet peptide are candidates for amyloid targeting reagents in vivo.

GENERAL SIGNIFICANCE

Efficient binding of polybasic peptides with amyloid is dependent on the linearity of charge spacing in the context of an α-helical secondary structure. Peptides with an α-helix or β-sheet propensity and with similar alignment of basic residues is optimal.

New insights into the interaction between pyrrolyl diketoacids and HIV-1 integrase active site and comparison with RNase H

HIV-1 integrase (IN) inhibitors are one of the most recent innovations in the treatment of HIV infection. The selection of drug resistance viral strains is however a still open issue requiring constant efforts to identify new anti-HIV-1 drugs. Pyrrolyl diketo acid (DKA) derivatives inhibit HIV-1 replication by interacting with the Mg²⁺ cofactors within the HIV-1 IN active site or within the HIV-1 reverse-transcriptase associated ribonuclease H (RNase H) active site. While the interaction mode of pyrrolyl DKAs with the RNase H active site has been recently reported and substantiated by mutagenesis experiments, their interaction within the IN active site still lacks a detailed understanding. In this study, we investigated the binding mode of four pyrrolyl DKAs to the HIV-1 IN active site by molecular modeling coupled with site-directed mutagenesis studies showing that the DKA pyrrolyl scaffold primarily interacts with the IN amino residues P145, Q146 and Q148. Importantly, the tested DKAs demonstrated good effectiveness against HIV-1 Raltegravir resistant Y143A and N155H INs, thus showing an interaction pattern with relevant differences if compared with the first generation IN inhibitors. These data provide precious insights for the design of new HIV inhibitors active on clinically selected Raltegravir resistant variants. Furthermore, this study provides new structural information to modulate IN and RNase H inhibitory activities for development of dual-acting anti-HIV agents.

Preclinical Validation of the Heparin-Reactive Peptide p5+14 as a Molecular Imaging Agent for Visceral Amyloidosis

Amyloid is a complex pathologic matrix comprised principally of paracrystalline protein fibrils and heparan sulfate proteoglycans. Systemic amyloid diseases are rare, thus, routine diagnosis is often challenging. The glycosaminoglycans ubiquitously present in amyloid deposits are biochemically and electrochemically distinct from those found in the healthy tissues due to the high degree of sulfation. We have exploited this unique property and evaluated heparin-reactive peptides, such as p5+14, as novel agents for specifically targeting and imaging amyloid. Herein, we demonstrate that radiolabeled p5+14 effectively bound murine AA amyloid in vivo by using molecular imaging. Biotinylated peptide also reacted with the major forms of human amyloid in tissue sections as evidenced immunohistochemically. Furthermore, we have demonstrated that the peptide also binds synthetic amyloid fibrils that lack glycosaminoglycans implying that the dense anionic motif present on heparin is mimicked by the amyloid protein fibril itself. These biochemical and functional data support the translation of radiolabeled peptide p5+14 for the clinical imaging of amyloid in patients.