Next Article in Journal
Beyond Cholera: Characterization of zot-Encoding Filamentous Phages in the Marine Fish Pathogen Vibrio anguillarum
Next Article in Special Issue
Contribution of Human Retroviruses to Disease Development—A Focus on the HIV– and HERV–Cancer Relationships and Treatment Strategies
Previous Article in Journal
Histone Deacetylase 6 Knockout Mice Exhibit Higher Susceptibility to Influenza A Virus Infection
Previous Article in Special Issue
A High-Throughput HIV-1 Drug Screening Platform, Based on Lentiviral Vectors and Compatible with Biosafety Level-1
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 7
10.3390/v12070729
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Targeting HIV-1 RNase H: N’-(2-Hydroxy-benzylidene)-3,4,5-Trihydroxybenzoylhydrazone as Selective Inhibitor Active against NNRTIs-Resistant Variants

by
Angela Corona
1,*,
Ester Ballana
2,
Simona Distinto
1,
Dominga Rogolino
3,
Claudia Del Vecchio
4,
Mauro Carcelli
3,
Roger Badia
2,
Eva Riveira-Muñoz
2,
Francesca Esposito
1,
Cristina Parolin
4,
José A. Esté
2,5,
Nicole Grandi
1 and
Enzo Tramontano
1,6
1
Department of Life and Environmental Sciences University of Cagliari, Cittadella Universitaria di Monserrato, 09042 Monserrato, Italy
2
AIDS Research Institute—IrsiCaixa, 08916 Badalona, Spain
3
Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, 43124 Parma, Italy
4
Department of Molecular Medicine, University of Padova, 35122 Padova, Italy
5
CienciaTraducida, 08391 Barcelona, Spain
6
Istituto di Ricerca Genetica e Biomedica, Consiglio Nazionale delle Ricerche (CNR), 09042 Monserrato, Cagliari, Italy
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(7), 729; https://0-doi-org.brum.beds.ac.uk/10.3390/v12070729
Submission received: 30 April 2020 / Revised: 28 June 2020 / Accepted: 1 July 2020 / Published: 6 July 2020
(This article belongs to the Special Issue Antiretroviral Drug Development and HIV Cure Research)
Browse Figures
Versions Notes

Abstract

:
HIV-1 infection requires life-long treatment and with 2.1 million new infections/year, faces the challenge of an increased rate of transmitted drug-resistant mutations. Therefore, a constant and timely effort is needed to identify new HIV-1 inhibitors active against drug-resistant variants. The ribonuclease H (RNase H) activity of HIV-1 reverse transcriptase (RT) is a very promising target, but to date, still lacks an efficient inhibitor. Here, we characterize the mode of action of N’-(2-hydroxy-benzylidene)-3,4,5-trihydroxybenzoylhydrazone (compound 13), an N-acylhydrazone derivative that inhibited viral replication (EC50 = 10 µM), while retaining full potency against the NNRTI-resistant double mutant K103N-Y181C virus. Time-of-addition and biochemical assays showed that compound 13 targeted the reverse-transcription step in cell-based assays and inhibited the RT-associated RNase H function, being >20-fold less potent against the RT polymerase activity. Docking calculations revealed that compound 13 binds within the RNase H domain in a position different from other selective RNase H inhibitors; site-directed mutagenesis studies revealed interactions with conserved amino acid within the RNase H domain, suggesting that compound 13 can be taken as starting point to generate a new series of more potent RNase H selective inhibitors active against circulating drug-resistant variants.

1. Introduction

The number of people worldwide infected with HIV-1 in 2019 has been estimated to be around 37.9 million [1]. Optimum treatment regimens consist of a combination of two [2] or more antiretrovirals with a different mode of action, that allows the suppression of viremia to undetectable levels, giving a life-expectation comparable to uninfected people [3,4]. However, despite promising efforts, presently approved treatments do not allow immunization of people [5] or the eradication of the infection [6], requiring life-long treatment with optimal adherence. At the moment, only 64.6% of HIV-1 infected people are accessing antiretroviral therapy [1]. Lack of treatment adherence and suboptimal coverage of the infected population cause the emergence of drug-resistant variants (DRV) and treatment failure [7]. DRV transmission narrows the options for treatment and is a major concern that requires an optimized therapy [8,9,10]. Thus, there is a constant need to explore new drugs and new targets for salvage therapy.
Among the viral encoded enzymatic functions, the HIV-1 RT-associated ribonuclease H (RNase H) represents the sole target for which an effective antiviral drug has not yet been approved. HIV-1 RNase H function is essential for viral replication [11] and is characterized by a highly conserved catalytic domain, with Asp–Asp–Glu (DDE) motif, essential for coordination of two Mg2+ or Mn2+ ions and by a tertiary structure [12] conserved among Ortervirales [13,14,15,16,17].
Two decades of research has allowed the identification of many compounds able to block HIV-1 RNase H activity. These compounds may be grouped in two main families according to their binding mode [18,19,20,21]: (i) active site inhibitors, which chelate the two Mg2+ within the active site; and (ii) RNase H allosteric inhibitors. The most potent inhibitors belong to the first class [22]. Active-site inhibitors take advantage of the driving force represented by their chelating moiety to stably interact with the catalytic core domain [22,23]. Despite notable anti-HIV potency of active-site inhibitors, off-target toxicity specific to this mode of action [24] and an unclear competition with the substrate [25] have precluded their further development as therapeutic agents.
Attempts to target highly conserved residues in the RNase H domain have been successfully shown [26,27,28]. Given the structural similarities and overlaps with the HIV-1 IN catalytic core, many of the identified RNase H active-site inhibitors exhibit dual anti-RNase H and anti-integrase (IN) activity [20,29]. Nevertheless, compounds with partial or even total selectivity for HIV-1 RNase H function has been obtained while retained full potency of inhibition against multidrug-resistant RTs [10].
We previously investigated a series of N-acylhydrazone analogs to identify dual inhibitors of both HIV IN and RNase H [30]. Among them several compounds inhibited the RT-associated RNase H function in the low micromolar range, being inactive against HIV-1 IN. In the present work, we aim to dissect the mode of action of these inhibitors, laying the bases for further development as lead compounds. To this aim, two among previously reported compounds were chosen: compound 13 (Figure 1) as an example of active-site inhibitor selective for RNase H in enzymatic assays [30] and compound 21 (Figure 1), that inhibited also HIV-1 IN catalytic activity [30] as counterpart example of an RNase H/IN dual inhibitor.
Time-of-addition and antiviral assays showed that compound 13 blocks HIV-1 replication in cell culture by targeting the RNase H function during the reverse-transcription step and, importantly, retained full potency of inhibition against a non-nucleoside RT inhibitor (NNRTI) drug-resistant strain. The combined results of computational modeling and site-directed mutagenesis revealed that compound 13 interacts with the RNase H active site, coordinating the ionic cofactors, binding in an orientation different from the known inhibitor RDS1759 [26], and establishing additional interactions with residues of the RNase H domain. The analysis of the RNase H domain nucleotide sequences revealed the presence of highly conserved amino acid residues, also part of the RNase H primer grip motif interacting with compound 13 [31]. This demonstrated that amino acid residues involved in compound 13 interaction with RT are highly conserved and give insights for the optimization of the lead.

2. Materials and0 Methods

2.1. Biology

2.1.1. Cells Reagents and Viruses

Compounds 13 and 21 were synthesized as previously reported [30]. The inhibitors: azidothymidine (AZT); efavirenz (EFV); nevirapine (NVP) and raltegravir (RAL) were purchased form Sigma-Aldrich (St. Louis, MO, USA). RNase H inhibitor was kindly provided by Prof. Roberto Di Santo. MT4 (ECACC 08081402) and HEK 293 T cells (ATCC® CRL-3216TM) were cultured in RPMI 1640 medium and Dulbecco’s modified Eagle’s medium (DMEM) (Gibco, Waltham, MA, USA), respectively. Both media were supplemented with 10% fetal calf serum (FCS) (Gibco, Waltham, MA, USA) and 1% penicillin/streptomycin. Stocks of HIV-1 NL4-3 strain were prepared by transfecting 293 T cells with the HIV-1 pNL4-3. The VSV-pseudotyped, NL4-3–GFP expressing virus (VSV/HIV-1GFP) was generated as reported [32].

2.1.2. Antiviral and Cytotoxicity Assay

The antiviral assay was performed as previously described [33,34]. Briefly, drug-mediated inhibition of virus-induced cytotoxicity was assayed in MT4cells at 120 h post infection. Triplicate wells of 96-well plates containing serial dilutions of drugs were seeded with 3 × 104 MT4 cells and infected with HIV-1 NL4.3 wt or with the HIV-1 NL4-3 K103N-Y181C strain at multiplicity of infection (MOI) of 0.003 or 100 culture infective dose (CCID50) calculated using the Reed and Muench method [35]. Five days post infection, 20 μL 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) per well were added. After 1 h of incubation at 37 °C in a humidified atmosphere with 5% CO2, 150 μL of supernatant were carefully removed and 100 μL of acidified Triton X-100 isopropanol solution were added to each well and allow the complete dissolution of the formazan crystals by placing the plates on a vibrating platform shaker for 10 min.
The amount of formazan produced was quantified spectrophotometrically at 550/620 nm. The 50% effective concentration (EC50) and 50% cytotoxic concentration (CC50) values were calculated for each compound by applying a nonlinear regression with log dose vs. normalized response to the data.

2.1.3. Time of Addition Assay

Time-of-addition assay was performed according to [36] using the VSV/HIV-1GFP virus. The MT4 cells/well were pelleted and resuspended with RPMI containing VSV/HIV-1GFP at a MOI of 0.2, centrifuged for 1 h at 1200× g at 15 °C and washed 2 times with cold PBS. A total of 1 × 104 infected cells were seeded (100 μL) in a 96-well plate round bottom, after addition of 50 μL of RPMI to the whole plate and 50 μL 4X-dilutions of drugs to the first column, representing time 0. Drugs were added at 11 time-points post infection (1 h, 1.5 h, 2 h, 2.5 h, 3 h, 4 h, 5 h, 6 h, 7 h, 8 h, 24 h). Drugs were used at the following final concentrations: AZT—3.7 μM; NVP—3.7 μM; EFV—0.3 μM; RAL—2.25 μM; compound 13—25 μM; compound 21—25 μM; compound 21—25 μM. Twenty-four hours post infection, cells were washed twice with PBS and fixed with 1% formaldehyde. GPF expression was quantified by fluorescence-activated cell sorting (FACS) with a LSRII flow cytometer (BD Biosciences, San José, CA, USA), recording FSC, SSC and FITC data. Data were analyzed with FlowJo software version10.6.1 (Tree Star Inc., Ashland, OR, USA).

2.2. Molecular Modeling

2.2.1. Ligand Preparation

The ligand was built within the Mestro platform, the geometry was optimized with MacroModel (Schrödinger, LLC, New York, NY, USA, 2020) [37] using the Merck molecular force fields (MMFFs) [38], the GB/SA solvation model [39] and the Polak–Ribière conjugate gradient (PRCG) method and a convergence threshold of 0.05 kJ/(molÅ).

2.2.2. Protein Preparation

Starting crystal coordinates of the complex RT–RNase H inhibitor was downloaded from the Protein Data Bank (http://www.rcsb.org/) pdb accession code 3LP2 [40]. Afterward, the protein was prepared using the Schrödinger protein preparation wizard [41]: hydrogen atoms were added to the system. Partial atomic charges were assigned according to the optimized potential for liquid simulations (OPLS-2005) force field [42]. A restrained minimization was performed to optimize hydrogen atoms and to remove eventually any high-energy contacts or distorted bonds, angles and dihedrals with to converge atoms to RMSD 0.3.

2.2.3. Docking and Post-Docking

Compound 13 was docked into wt RTs through QM-polarized ligand docking protocol applying default settings [43]. In order to better take into account the induced fit phenomena, the most energy favored generated complexes were fully optimized with the OPLS force field in GB/SA implicit water [39], setting a 10,000-step-interaction analysis with Polak–Ribière conjugate gradient (PRCG) method and a convergence criterion of 0.1 kJ/(molÅ). The resulting complexes were considered for the binding modes graphical analysis with Pymol (Schrödinger, LLC, New York, NY, USA, 2020) [44] and Mestro (Schrödinger, LLC, New York, NY, USA, 2020) [41].

2.3. Pharmacophore Model Generation

All ligands reported in previous work [30] were modelled by means of LigandScout software (ver. 4.4) (Inte: Ligand, Vienna, Austria) [45,46]. Conformational analysis was performed before the alignment process applying best setting (i.e., max 200 conformations). Then compounds were aligned and the ligand-based share pharmacophore was automatically generated considering both scoring functions: pharmacophore fit and atom overlap.

2.4. Molecular Biology

2.4.1. Expression and Purification of Recombinant HIV-1 RTs Group M Subtype B wt and Mutants

The p6HRT-prot vector was kindly provided by Dr. Stuart Le Grice Laboratory. Heterodimeric RT was expressed and purified essentially as described [47]. Briefly, Escherichia coli strain M15 containing the p6HRT-prot vector were grown up to an OD600 of 0.7 and induced with isopropyl β-D-1-thiogalactopyranoside (IPTG) (Sigma-Aldrich, St. Louis, MO, USA) 1.7 mM for 4 h at 37 °C. Cells were pelleted and lysed with 50-mM sodium phosphate pH 7.8, 0.5-mg/mL lysozyme, sonicated and centrifuged at 30,000× g for 1 h. The proteins were expressed and purified by a combination of affinity and ion-exchange chromatography. First the supernatant was loaded into a Ni2+ Sepharose (GE Healthcare, Chicago, IL, USA) with loading buffer (50-mM sodium phosphate pH 7.8, 0.3-M NaCl, 10% glycerol, 10-mM imidazole) and washed thoroughly with the same buffer at 80-mM imidazole. RTs were eluted with a gradient of up to 0.5-M imidazole. Fractions were collected; protein purity was checked by SDS-PAGE and found to be higher than 90%. RT containing fractions were pooled and diluted 1:1 with dilute buffer (50-mM sodium phosphate pH 7.0, 10% glycerol) then loaded into a Hi-trap Heparin HP (GE Healthcare, Chicago, IL, USA) preequilibrated with 10-column volume of loading buffer 2 (50-mM sodium phosphate pH 7.0, 10% glycerol, 150-mM NaCl). The column was then washed with loading buffer 2 and RT was gradient-eluted with elute buffer 2 (50-mM sodium phosphate pH 7.0, 10% glycerol, 150-mM NaCl). Fractions were collected, the protein was dialyzed and stored in buffer containing 50-mM Tris HCl pH 7.0, 25-mM NaCl, 1-mM EDTA, 50% glycerol. Catalytic activities and protein concentration were determined. Enzyme-containing fractions were pooled, and aliquots were stored at −80 °C.

2.4.2. Site-Directed Mutagenesis

Amino acid substitutions were introduced into the p66 HIV-1 RT subunit coded in a p6HRT-prot plasmid using the QuikChange protocol (Agilent Technologies, Inc., Santa Clara, CA, USA).

2.4.3. HIV-1 DNA Polymerase-Independent RNase H Activity Determination

The RT-associated RNase H activity of wt and mutated HIV RTs was measured as described [48], in 100 µL reaction volume containing 50 mM Tris HCl pH 7.8, 6-mM MgCl2, 1-mM dithiothreitol (DTT), 80-mM KCl, hybrid RNA/DNA (5′-GTTTTCTTTTCCCCCCTGAC-3′-fluorescein, 5′-CAAAAGAAAAGGGGGGACUG-3′-dabcyl, from Metabion, Planegg, Germany) and different amounts of enzymes according to a linear range of dose–response curve (60 ng of wt RT; 60 ng R448A RT; 160 ng R557A RT; 600 ng Q475A RT). The reaction mixture was incubated for 1 h at 37 °C, the reaction was stopped by addition of EDTA and products were measured with a Victor 3 (PerkinElmer, Waltham, MA, USA) at 490/528 nm. Mean ± standard deviation of 50% inhibitory concentration (IC50) values was determined and p-values were calculated between IC50 value against the wt and IC50 value against the mutants by paired, two-tailed t-tests using GraphPad Prism 6.01 software (GraphPad Software, Inc.; San Diego, CA, USA). Figures were made with GraphPad Prism 6 version 6.01.

2.4.4. HIV-1 RNA-Dependent DNA Polymerase (RDDP) Activity Determination

The HIV-1 RT-associated RDDP activity was measured as described [49], in 25 µL volume containing 60-mM Tris-HCl pH 8.1, 8-mM MgCl2, 60-mM KCl, 13-mM DTT, poly(A)-oligo(dT), 100-µM dTTP and 6 ng wt RT. After enzyme addition, the reaction mixture was incubated for 30 min at 37 °C and enzymatic reaction was stopped by addition of EDTA. Reaction products were detected by PicoGreen addition and measured with a Victor 3 (PerkinElmer, Waltham, MA, USA) at 502/523 nm. Figures were made with GraphPad Prism 6 version 6.01.

2.5. Genetic Analysis

A nucleotide sequence collection including a total of 15.865 RT isolates was downloaded from the HIV-1 Single Genome Sequence (SGS) Database of the Stanford University HIV database [50,51] (https://hivdb.stanford.edu/project/sgs/). RT nucleotide sequences were aligned using Geneious software [52] implemented with MAFFT aligner [53], and RH nucleotide portion was extracted and translated to obtain the deriving RH proteins. The RT isolates devoid of the RNase H portions were excluded, as well as those including only a partial RNase H domain or showing internal frameshifts and stop codons, ending up with a total of 5798 aligned RH proteins. The latter were analyzed with respect to HIV-1 strain NL4.3 RT–RNase H, used as reference protein, indicating for each amino acid residue the percentage of conservation. A logo representing the degree of conservation of each amino acidic position was also produced from the RNase H domain protein alignment with the tool WebLogo 3 (http://weblogo.threeplusone.com/create.cgi).

3. Results

To determine the mode of action of N-acylhydrazone analogs that displayed selective inhibition of HIV-1 RT RNase H associated function and to acquire information useful to optimize the scaffold, a series of cell -based and enzymatic assays were performed on compound 13.

3.1. Antiviral Activity of N-acylhydrazone Analogs against wt and NNRTI-Resistant HIV-1 Strains

Antiviral activity and cytotoxicity of compounds 13 and 21 were determined at 5 days on MT4 cells against HIV-1 NL4.3 wt strain and HIV-1 NL4.3 carrying the K103N-Y181C double mutation, commonly selected in patients and that confers resistance to non-nucleoside RT inhibitors (NNRTIs). The nucleoside inhibitor AZT, the NNRTIs EFV and NVP, the IN inhibitor RAL and the RNase H inhibitor RDS1759 were used as positive controls (Table 1). Both N-acylhydrazone derivatives inhibited viral replication retaining similar potency of inhibition against the NNRTI-resistant double mutant NL4-3 K103N-Y181C virus, similarly to AZT and RAL and the RNase H inhibitor RDS1759, while the NNRTIs lost their antiviral activity. Unfortunately, the cytotoxicity showed by the compounds determined a selectivity index below 10.

3.2. Analysis of Compound 13 Target of Inhibition by Time-of-Addition and Comparative Inhibition of HIV-1 RT Associated Activities

To confirm that compound 13 was selectively targeting the RT-associated RNase H function during HIV-1 replication, we performed a time-of-addition of compound 13 at 25 μM using the RNase H/IN dual inhibitor compound 21, the nucleoside inhibitor AZT, the NNRTIs EFV and NVP and the IN inhibitor RAL as internal controls (Figure 2A). The two N-acylhydrazones showed a very different behavior: compound 13 lost activity 6 h after infection, in a time–window between AZT and EFV/NVP, in agreement with the hypothesis of targeting the retrotranscription step of the replicative viral cycle consistently with the inactivity against HIV-1 IN previously reported [30] and suggesting a mode of action different from nucleoside analogs and NNRTIs. Differently, compound 21 retained its potency over time, together with the IN inhibitor RAL affecting also the integration step, accordingly with the previously reported inhibition of IN catalytic activity [30].
Once shown that compound 13 is indeed selectively inhibiting the viral genome reverse transcription, since several RNase H inhibitors were demonstrated to be able to inhibit both RT-associated functions [20,54], we asked whether compound 13 could affect also the RT-associated RNA-dependent DNA polymerase (RDDP) activity. Hence, a comparative inhibition of both HIV-1 RT associated activities was performed (Figure 2B). To overcome any possibility of overestimating the potency of inhibition, the experiments were performed in conditions of competition between the compounds and the substrate, since both inhibitor and substrate were added in solution before the enzyme. Results showed that compound 13 inhibited RT-associated RNase H activity with an IC50 value of 3.6 μM, while only a residual inhibition of the RDDP activity was observed (IC50 value of 81 μM), with no inhibitory effect at 25 μM, clearly indicating that the inhibition observed at 25 μM in the cell-based time-of-addition experiment and the inhibition of viral replication at 5 days (EC50 10.1 µM) was specifically due to its inhibition of the HIV-1 RT-associated RNase H function.

3.3. Analysis of the Binding of Compound 13 within the HIV-1 RT RNase H Domain

In order to have insights on the binding mode of compound 13 within the HIV-1 RT RNase H domain, docking studies were performed.
Docking calculations showed that compound 13 binds within the RNase H domain in two possible orientations (Figure 3). In the first one (Figure 3A,B), compound 13 coordinates the two Mn2+ ions with the carbonyl moiety of the acylhydrazone group, establishing additional interactions with residues Ser553, Asp498, Asn474, Gln475, Glu478. In the second binding mode (Figure 3C,D) compound 13 coordinates the two Mn2+ cofactors with the pyrogallol moiety serving as a chelating motif, while the remaining atoms of the molecule accommodate between alpha-helix E (in yellow) and beta-sheet 1 (in green), interacting with charged and polar residues Asp443, Glu449, Glu478, Arg557, Asn474, Arg557, establishing further hydrogen bonds with the side chain of Gly444, Ala446 and Arg448 (Figure 4B). This last orientation of the chelating portion is different from the one reported for the selective RNase H inhibitor RDS1759, but similar to the chelating portion of other compounds as beta-tuhjaplicinol (pdb 3I1G) [55] and 2-hydroxyisoquinoline-1,3-dione (pdb code 5UV5) [56], while the binding mode exhibited by the tail portion of compound 13 significantly differs from all of them.

3.4. Site-Directed Mutagenesis Studies Detailing the Binding Mode of Compound 13

To confirm the model suggested by the docking calculation—and to determine which of the two binding orientations was the most probable—site-directed mutagenesis studies were performed on amino acid residues involved in critical interactions with either compound 13 or RDS1759 hypothesized to bind to RT differently from compound 13. Among the suggested residues, Arg448, Gln475 and Arg557 were chosen as the most significant to assess the binding mode. In particular, Arg448 plays a discriminant role since it is involved in the interaction between compound 13 and the RNase H domain in the sole second orientation, while both Arg448 and Arg557 would discriminate the differences with the known inhibitor RDS1759 that does not interact with that portion of the HIV-1 RT RNase H domain [26]. Hence, the single point mutations R448A, Q475A and R557A were independently introduced in the HIV-1 pol and mutant R448A-, Q475A- and R557A-RT were expressed and purified. Compound 13 and RDS1759 inhibitory profile of the HIV-1 RT RNase H function of the enzymes were determined in parallel with the wt (Figure 4). Results showed that compound 13 inhibited the wt enzyme with an IC50 value of 6.9 μM, with a consistent decrease in potency against all the mutated enzymes (Figure 4A). The IC50 value increases of 5.2-fold against Q475A (p-value = 0.0003); 9.2-fold against R448A (p-value <0.0001) and more than 16.2-fold against R557A (p-value < 0.0001). Differently, RDS1759 potency of inhibition was not significantly affected against RTs carrying the R448A or R557A mutations as compared to the wt RT (p-values equal to 0.5545 and 0.9682, respectively) while it significantly lost its potency when tested on the Q475A RT (p-value = 0.0008). Overall, the data confirm the binding pose number 2 suggested by docking calculation.

3.5. Genetic Analysis

To establish whether the residues involved in the interaction between the HIV-1 RNase H domain and compound 13 were conserved, we evaluated the existing variability of HIV-1 RT RNase H domain among circulating HIV-1 variants, using a nucleotide collection including a total of 15.865 RT isolates from the HIV-1 Single Genome Sequence (SGS) Database of the Stanford University HIV database [51]. After alignment of the sequences, we retained for analysis only the ones containing a complete HIV-1 RT RNase H domain, for a total number of 5798 isolates. The latter were analyzed using as reference protein HIV-1 strain NL4.3 RNase H, indicating for each amino acid residue the percentage of conservation (Table S1, Supplementary Materials). A logo representing the degree of conservation of each amino acid position was produced (Figure 5A). Results showed that the HIV-1 RT RNase H domain has a very conserved structure (Figure 5A, PDB code 3K2P), with a number of patterns of residues that display a variability lower than 0.5% (amino acid residues: 427–430; 438–445; 472–475; 485–490; 497–501; 535–542; 549–557) with Arg557 and Gln475 showing a 0% variability and Asn474 and Gln444 showing a 0.1% variability. Interestingly, residues Arg448 and Glu449 also displayed a good degree of conservation (96.70% and 95.40%, respectively), with Arg449 being replaced by Lys in the 3.3% of the cases and Glu449 being replaced by Asp in 4.6% of the cases, suggesting strong preservation of their structural functions and of the spatial characteristic needed for binding of compound 13.

3.6. Pharmacophore Model

Considering the biologic data previously reported for compounds 1–23 (Table S2) [30] and aiming at rationalizing the activity of this class of HIV-1 RNase H, we carried out pharmacophore model generation.
The share ligand-based pharmacophore model reported several common features required for inhibitory activity (Figure 6): the hydroxylic groups are key features and their presence in position 3,4,5 of phenyl ring (e.g., compounds 13 and 21, Figure 6C) is important. In fact, (i) the lack of this portion (e.g., compound 1), (ii) the hydroxylic group in other positions (compound 5) or (iii) the substitution with methoxy groups (e.g., compound 19) led to the loss of activity. Moreover, the presence in the linker of a moiety able to accept and donate hydrogen bonds is relevant, as it is a second aromatic group. These additional features indicate the additional interactions established by the “tail” of the compound while accommodated between alpha helix E and beta sheet 1, that stabilize the complex with the enzyme and are important to the binding, besides the driving force represented by metal binding. While this second aromatic group is a key factor for activity, there are no indications on its possible substitutions. Worth to note, most of the active compounds have hydroxylic group in position 2 (e.g., compounds 13, 14, 15, 16, 17, 18, 20, 21, 22), but others do not have it and their activity is still good (compounds 11, 12, 23), hence suggesting that this group is not essential for interaction with the enzyme. Furthermore, the excluded volume coat (in gray, Figure 6B), highlights the importance of steric occupancy and ideally defines the pocket. Docking experiments (Figure 3) also show that, while pyrogallol moiety drives the binding in the RNase H domain by chelating the ions, the tail is accommodated between alpha helix E and beta sheet 1. Hence the tail cannot be too bulky.

4. Discussion

Among the many molecules identified as HIV-1 RT-associated RNase H inhibitors the most potent inhibitors act by coordinating the ionic cofactors within the RNase H active sites [20,23]. Among them, a number of compounds were deeply investigated for their activity on RTs and viruses carrying mutations conferring resistance to approved inhibitors. Results showed that they retained full potency of inhibition against many drug-resistant viruses [10,57,58]. The development of these inhibitors, however, was hampered by the steep substrate barrier [25] or by loss of potency in cell-based assays [22].
Here, we analyzed the mode of action of an active-site RNase H inhibitor, compound 13, also analyzing the conserved patterns of HIV-1 RT RNase H domain suitable to be targeted to improve its potency.
Cell-based assays showed that compound 13 inhibits the HIV-1 replication with an EC50 value of 10.1 μM, in a five-day method with a full replicating virus in MT-4 cells, comparable with other RNase H inhibitors recently reported [22]. Of note, the inhibition of viral replication well reflects compound 13 efficacy in biochemical assays (HIV-1 RNase H IC50 = 6.8 μM).
We determined that reverse transcription is the actual target for compound 13, consistently with the inactivity against HIV-1 IN previously reported [30] and we showed that its activity cannot be correlated with the poor inhibition of HIV-1 RDDP activity (HIV-1 RDDP IC50 = 81 μM), possibly caused by some long-range effects already hypothesized for inhibitors binding within the RNase H domain [58]. Consistently, compound 13 retained full potency of inhibition against the NNRTI drug-resistant HIV-1 strain K103N-Y181C.
Despite a relatively low SI, this information encouraged the study of compound 13 binding mode by docking calculations and its validation by site-directed mutagenesis studies. Compound 13, contains a 3,4,5-trihydroxybenzohydrazone scaffold that shares remarkable similarities with other N-acylhydrazones previously reported as the first allosteric inhibitors of HIV-1 RT-associated RNase H function [59,60]. Our binding mode suggested that compound 13 binds within the RNase H active site with the pyrogallol group is serving as the chelating portion, similarly to what is hypothesized for other pyrogallol-derivatives RNase H inhibitors [61]. Compound 13 displayed an orientation different from the one reported for others active site RNase H inhibitors, whose binding mode involved mostly the highly conserved residue His539 [23,58,62] and from the diketo acid inhibitor RDS1759. Similar to RDS1759, compound 13 shows interactions with residues Asn474 and Gln475 involved in the RNase H primer grip motif [26,31] but, besides, it also interacts with Gly444, Arg448, Glu449 and Arg557 amino acid residues, some of them previously investigated and shown to be involved in the binding of RNase H active-site inhibitors [26,27,28].To prove compound 13 binding mode and its most critical interactions with RT, site-directed mutagenesis experiments were conducted, generating RTs carrying R448A, Q475A or R557A amino acid substitutions. All the three substitutions impair the catalytic efficiency of HIV-1 RT-associated RNase H activity by affecting both kcat and Km values leading to a reduced enzymatic efficiency [26,31].
The inhibitory activity of compound 13 and RDS1759 was assessed in conditions of competition with the substrate, proving that the compound potency was not affected by substrate displacement, as suggested elsewhere for other classes of compounds [25,63]. The significant loss of potency of compound 13 against all the generated mutant RTs proved the critical role of those amino acid residues for its binding to the RNase H domain and implies a binding orientation different from RDS1759, imidazolidinedione derivatives and beta-thujaplicinol [26,27] that were reported to be fully active against both R448A and R557A HIV-1 RTs.
Moreover, to provide useful information not only for the optimization of these molecules, but, in a broader context, for the design of new inhibitors, we performed the analysis of the conservation of the residues of the HIV-1 RT RNase H domain among the circulating variants, using the HIV-1 SGS Database of the Stanford University HIV database [50,51] that contains sequences of isolated from naïve and drug-experienced samples. The analysis revealed a high degree of conservation of the whole HIV-1 RT RNase H domain, particularly regarding the active site catalytic core. Interestingly, the analysis also revealed that more peripheric residues, such as Arg448 and Arg557, involved in the binding of compound 13, retain a good degree of conservation among naïve patients and patients treated with different classes of RT inhibitors, also carrying drug-resistant mutations [50,51], hence supporting the development of new RNase H inhibitors that would be potentially endowed with high drug-resistant barrier. This could be one of the main reason to justify the development of compounds against a virus with a wide therapeutic armamentarium already available such as HIV, because it may envision to a combination treatment including a lower number of drugs.
Finally, to summarize the structural features required for this new class of active-site HIV-1 RT RNase H inhibitors a pharmacophore model was obtained. The model highlights the chelating moiety already known to be the main feature for other classes of RNase H active site inhibitors such hydroxypyrimidine [20], beta-thujaplicinol derivatives [22,55], N-hydroximides [64]. Besides that, the pharmacophore also gave relevant knowledge on the importance and required spatial orientation of aromatic and hydrazone moieties. Of note, the model shows that the coat of the excluded volumes that delimitate the space is also relevant to the design of new molecules.

5. Conclusions

We characterized a new promising N-acylhydrazone derivative, compound 13, that selectively inhibits the HIV-1 RT-associated RNase H function during the HIV-1 replication, interacting with highly conserved residues within the RNase H domain of HIV-1 RT and that retains activity against an NNRTI-resistant strain. The binding mode showed different orientations form others previously reported suggesting new possibilities for targeting the highly conserved HIV-1 RT RNase H domain with drugs with potential high-genetic barrier. We also determined a pharmacophore model to understand the important features for the activity and can be used to drive the optimization of the series of molecule analyzed. Furthermore, it can be applied to select new scaffolds in commercial or other available databases. Overall, we have good bases to generate a new series of more potent, HIV-1 RT RNase H selective inhibitors, active against circulating drug-resistant variants, improving their binding and their safety profile.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/7/729/s1. Table S1: Percentage of conservation of the amino acid residue in the corresponding position, Table S2: Structures and anti-HIV-1 RT RNase H activities of compounds 1–23.

Author Contributions

Conceptualization, A.C. and E.T.; methodology, A.C., S.D., R.B., E.R.-M., F.E., J.A.E.; modeling, S.D.; investigation, A.C., S.D., N.G., R.B., E.B., C.D.V.; resources, E.T., J.A.E., M.C., D.R.; data curation, A.C., E.B., N.G., S.D.; writing—original draft preparation, A.C.; writing—review and editing, A.C., J.A.E., N.G., E.T., S.D., C.P.; visualization, A.C., N.G.; supervision, E.T., C.P., J.A.E.; project administration, E.T.; funding acquisition, E.T., J.A.E. All authors have read and agreed to the published version of the manuscript.

Funding

A.C. research was funded by the Sardinian Regional Government grant LR07/17 (F76C18000800002). A.C. acknowledges the support of the ISAR Chu Family Foundation Scholarships for Early Career Women in Science. The study was funded in part by project BFU2015-63800R (JAE) from the Spanish Ministerio de Ciencia e Innovación. D.R. research was supported by the national grant “Progetti di ricerca di Rilevante Interesse Nazionale” (PRIN 2017BMK8JR).

Acknowledgments

Authors thank Roberto di Santo for providing the internal control RDS1759. D.R. and M.C. thank the “Centro Interdipartimentale Misure Giuseppe Casnati” of the University of Parma for facilities.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. UNAIDS 2018 Global HIV Statistics. Available online: https://www.unaids.org/sites/default/files/media_asset/UNAIDS_FactSheet_en.pdf (accessed on 17 February 2020).
  2. Cahn, P.; Madero, J.S.; Arribas, J.R.; Antinori, A.; Ortiz, R.; Clarke, A.E.; Hung, C.C.; Rockstroh, J.K.; Girard, P.M.; Sievers, J.; et al. Durable Efficacy of Dolutegravir Plus Lamivudine in Antiretroviral Treatment-Naive Adults With HIV-1 Infection: 96-Week Results From the GEMINI-1 and GEMINI-2 Randomized Clinical Trials. J. Acquir. Immune Defic. Syndr. 2020, 83, 310–318. [Google Scholar] [CrossRef] [Green Version]
  3. Fauci, A.S.; Folkers, G.K. Toward an AIDS-free generation. JAMA J. Am. Med. Assoc. 2012, 308, 343–344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Deeks, S.G.; Lewin, S.R.; Havlir, D.V. The end of AIDS: HIV infection as a chronic disease. Lancet 2013, 382, 1525–1533. [Google Scholar] [CrossRef] [Green Version]
  5. Cohen, J. Combo of two HIV vaccines fails its big test. Science 2020, 367, 611–612. [Google Scholar] [CrossRef] [PubMed]
  6. Henderson, L.J.; Reoma, L.B.; Kovacs, J.A.; Nath, A. Advances toward Curing HIV-1 Infection in Tissue Reservoirs. J. Virol. 2019, 94, 1–21. [Google Scholar] [CrossRef]
  7. Stella-Ascariz, N.; Arribas, J.R.; Paredes, R.; Li, J.Z. The Role of HIV-1 Drug-Resistant Minority Variants in Treatment Failure. J. Infect. Dis. 2017, 216, S847–S850. [Google Scholar] [CrossRef]
  8. Gupta, R.K.; Gregson, J.; Parkin, N.; Haile-Selassie, H.; Tanuri, A.; Andrade Forero, L.; Kaleebu, P.; Watera, C.; Aghokeng, A.; Mutenda, N.; et al. HIV-1 drug resistance before initiation or re-initiation of first-line antiretroviral therapy in low-income and middle-income countries: A systematic review and meta-regression analysis. Lancet Infect. Dis. 2017, 18, 346–355. [Google Scholar] [CrossRef] [Green Version]
  9. Gregson, J.; Tang, M.; Ndembi, N.; Hamers, R.L.; Marconi, V.C.; Brooks, K.; Theys, K.; Arruda, M.; Garcia, F.; Monge, S.; et al. Global epidemiology of drug resistance after failure of WHO recommended first-line regimens for adult HIV-1 infection: A multicentre retrospective cohort study. Lancet Infect. Dis. 2016, 16, 565–575. [Google Scholar] [CrossRef] [Green Version]
  10. Schneider, A.; Corona, A.; Spöring, I.; Jordan, M.; Buchholz, B.; Maccioni, E.; Di Santo, R.; Bodem, J.; Tramontano, E.; Wöhrl, B.M. Biochemical characterization of a multi-drug resistant HIV-1 subtype AG reverse transcriptase: Antagonism of AZT discrimination and excision pathways and sensitivity to RNase H inhibitors. Nucleic Acids Res. 2016, 44, 2310–2322. [Google Scholar] [CrossRef] [Green Version]
  11. Schatz, O.; Cromme, F.V.; Naas, T.; Lindemann, D. Inactivation of the RNase H domain of HIV-1 reverse transcriptase blocks viral infectivity. In Gene Regulation and AIDS; Gulf Publishing Company: Houston, TX, USA, 1990; pp. 293–404. [Google Scholar]
  12. Davies, J.F.; Hostomska, Z.; Hostomsky, Z.; Jordan, S.R.; Matthews, D.A. Crystal structure of the ribonuclease H domain of HIV-1 reverse transcriptase. Science 1991, 252, 88–95. [Google Scholar] [CrossRef]
  13. Krupovic, M.; Blomberg, J.; Coffin, J.M.; Dasgupta, I.; Fan, H.; Geering, A.D.; Gifford, R.; Harrach, B.; Hull, R.; Johnson, W.; et al. Ortervirales: New Virus Order Unifying Five Families of Reverse-Transcribing Viruses. J. Virol. 2018, 92, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Tadokoro, T.; Kanaya, S. Ribonuclease H: Molecular diversities, substrate binding domains, and catalytic mechanism of the prokaryotic enzymes. FEBS J. 2009, 276, 1482–1493. [Google Scholar] [CrossRef] [PubMed]
  15. Engelman, A.; Cherepanov, P. Retroviral integrase structure and dna recombination mechanism. Microbiol. Spectr. 2014, 2, 1–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hyjek, M.; Figiel, M.; Nowotny, M. RNases H: Structure and mechanism. DNA Repair. 2019, 84, 102672. [Google Scholar] [CrossRef] [PubMed]
  17. Corona, A.; Schneider, A.; Schweimer, K.; Rösch, P.; Wöhrl, B.M.B.M.; Tramontano, E. Inhibition of foamy virus reverse transcriptase by human immunodeficiency virus type 1 ribonuclease H inhibitors. Antimicrob. Agents Chemother. 2014, 58, 4086–4093. [Google Scholar] [CrossRef] [Green Version]
  18. Corona, A.; Masaoka, T.; Tocco, G.; Tramontano, E.; Le Grice, S.F. Active site and allosteric inhibitors of the ribonuclease H activity of HIV reverse transcriptase. Futur. Med. Chem. 2013, 5, 2127–2139. [Google Scholar] [CrossRef]
  19. Wang, X.; Gao, P.; Menendez-Arias, L.; Liu, X.; Zhan, P.; Menéndez-Arias, L.; Liu, X.; Zhan, P. Update on Recent Developments in Small Molecular HIV-1 RNase H Inhibitors (2013–2016): Opportunities and Challenges. Curr. Med. Chem. 2018, 25, 1682–1702. [Google Scholar] [CrossRef]
  20. Tramontano, E.; Corona, A.; Menéndez-Arias, L. Ribonuclease H, an unexploited target for antiviral intervention against HIV and hepatitis B virus. Antivir. Res. 2019, 171, 104613. [Google Scholar] [CrossRef]
  21. Tramontano, E. HIV-1 RNase H: Recent progress in an exciting, yet little explored, drug target. Mini Rev. Med. Chem. 2006, 6, 727–737. [Google Scholar] [CrossRef]
  22. Wang, L.; Tang, J.; Huber, A.D.; Casey, M.C.; Kirby, K.A.; Wilson, D.J.; Kankanala, J.; Xie, J.; Parniak, M.A.; Sarafianos, S.G.; et al. 6-Arylthio-3-hydroxypyrimidine-2,4-diones potently inhibited HIV reverse transcriptase-associated RNase H with antiviral activity. Eur. J. Med. Chem. 2018, 156, 652–665. [Google Scholar] [CrossRef]
  23. Xi, Z.; Wang, Z.; Sarafianos, S.G.; Myshakina, N.S.; Ishima, R. Determinants of active-site inhibitor interaction with HIV-1 RNase H. ACS Infect. Dis. 2019, 5, 1963–1974. [Google Scholar] [CrossRef] [PubMed]
  24. Corona, A.; Esposito, F.; Tramontano, E. Can the ever-promising target HIV reverse transcriptase-associated RNase H become a success story for drug development? Future Virol. 2014, 9, 445–448. [Google Scholar] [CrossRef]
  25. Wang, L.; Sarafianos, S.G.; Wang, Z. Cutting into the Substrate Dominance: Pharmacophore and Structure-Based Approaches toward Inhibiting Human Immunodeficiency Virus Reverse Transcriptase-Associated Ribonuclease H. Acc. Chem. Res. 2020, 53, 218–230. [Google Scholar] [CrossRef] [PubMed]
  26. Corona, A.; Di Leva, F.S.; Thierry, S.; Pescatori, L.; Cuzzucoli Crucitti, G.; Subra, F.; Delelis, O.; Esposito, F.; Rigogliuso, G.; Costi, R.; et al. Identification of Highly Conserved Residues Involved in Inhibition of HIV-1 RNase H Function by Diketo Acid Derivatives. Antimicrob. Agents Chemother. 2014, 58, 6101–6110. [Google Scholar] [CrossRef] [Green Version]
  27. Poongavanam, V.; Corona, A.; Steinmann, C.; Scipione, L.; Grandi, N.; Pandolfi, F.; Di Santo, R.; Costi, R.; Esposito, F.; Tramontano, E.; et al. Structure-guided approach identifies a novel class of HIV-1 ribonuclease H inhibitors: Binding mode insights through magnesium complexation and site-directed mutagenesis studies. Medchemcomm 2018, 9, 562–575. [Google Scholar] [CrossRef] [Green Version]
  28. Messore, A.; Corona, A.; Madia, V.N.; Saccoliti, F.; Tudino, V.; De Leo, A.; Scipione, L.; De Vita, D.; Amendola, G.; Di Maro, S.; et al. Pyrrolyl Pyrazoles as Non-Diketo Acid Inhibitors of the HIV-1 Ribonuclease H Function of Reverse Transcriptase. ACS Med. Chem. Lett. 2020, 11, 798–805. [Google Scholar] [CrossRef]
  29. Corona, A.; Di Leva, F.S.; Rigogliuso, G.; Pescatori, L.; Noemi, V.; Subra, F.; Delelis, O.; Esposito, F.; Cadeddu, M.; Costi, R.; et al. New insights into the interaction between pyrrolyl diketoacids and HIV-1 integrase active site and comparison with RNase H. Antivir. Res. 2016, 134, 236–243. [Google Scholar] [CrossRef]
  30. Carcelli, M.; Rogolino, D.; Gatti, A.; Pala, N.; Corona, A.; Caredda, A.; Tramontano, E.; Pannecouque, C.; Naesens, L.; Esposito, F. Chelation motifs affecting metal-dependent viral enzymes: N-acylhydrazone ligands as dual target inhibitors of HIV-1 integrase and reverse transcriptase ribonuclease Hdomain. Front. Microbiol. 2017, 8, 440. [Google Scholar] [CrossRef]
  31. Rausch, J.W.; Lener, D.; Miller, J.T.; Julias, J.G.; Hughes, S.H.; Le Grice, S.F.J. Altering the RNase H primer grip of human immunodeficiency virus reverse transcriptase modifies cleavage specificity. Biochemistry 2002, 41, 4856–4865. [Google Scholar] [CrossRef]
  32. Clouser, C.L.; Patterson, S.E.; Mansky, L.M. Exploiting Drug Repositioning for Discovery of a Novel HIV Combination Therapy. J. Virol. 2010, 84, 9301–9309. [Google Scholar] [CrossRef] [Green Version]
  33. Pannecouque, C.; Daelemans, D.; De Clercq, E. Tetrazolium-based colorimetric assay for the detection of HIV replication inhibitors: Revisited 20 years later. Nat. Protoc. 2008, 3, 427–434. [Google Scholar] [CrossRef] [PubMed]
  34. Badia, R.; Grau, J.; Riveira-Muñoz, E.; Ballana, E.; Giannini, G.; Esté, J.A. The thioacetate-ω(γ-lactam carboxamide) HDAC inhibitor ST7612AA1 as HIV-1 latency reactivation agent. Antivir. Res. 2015, 123, 62–69. [Google Scholar] [CrossRef]
  35. Reed, L.J.; Muench, H. A simple method of estimating fifty per cent endpoints. Am. J. Epidemiol. 1938, 27, 493–497. [Google Scholar] [CrossRef]
  36. Daelemans, D.; Pauwels, R.; De Clercq, E.; Pannecouque, C. A time-of-drug addition approach to target identification of antiviral compounds. Nat. Protoc. 2011, 6, 925–933. [Google Scholar] [CrossRef] [PubMed]
  37. Mohamadi, F.; Richards, N.G.J.; Guida, W.C.; Liskamp, R.; Lipton, M.; Caufield, C.; Chang, G.; Hendrickson, T.; Still, W.C.; Giuda, W.C.; et al. Macromodel—An integrated software system for modeling organic and bioorganic molecules using molecular mechanics. J. Comput. Chem. 1990, 11, 440–467. [Google Scholar] [CrossRef]
  38. Halgren, T.A. Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular interactions. J. Comput. Chem. 1996, 17, 520–552. [Google Scholar] [CrossRef]
  39. Still, W.C.C.; Tempczyk, A.; Hawley, R.C.C.; Hendrickson, T. Semianalytical treatment of solvation for molecular mechanics and dynamics. J. Am. Chem. Soc. 1990, 112, 6127–6129. [Google Scholar] [CrossRef]
  40. Su, H.-P.; Yan, Y.; Prasad, G.S.; Smith, R.F.; Daniels, C.L.; Abeywickrema, P.D.; Reid, J.C.; Loughran, H.M.; Kornienko, M.; Sharma, S.; et al. Structural basis for the inhibition of RNase H activity of HIV-1 reverse transcriptase by RNase H active site-directed inhibitors. J. Virol. 2010, 84, 7625–7633. [Google Scholar] [CrossRef] [Green Version]
  41. Schrödinger LLC. Maestro GUI; Schrödinger LLC: New York, NY, USA, 2013. [Google Scholar]
  42. Jorgensen, W.L.; Jorgensen, W.L.; Maxwell, D.S.; Maxwell, D.S.; Tirado-Rives, J.; Tirado-Rives, J. Development and Testing of the OLPS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  43. Schrödinger LLC. QMPolarized Protocol; Schrödinger LLC: New York, NY, USA.
  44. Schrödinger LLC. The PyMOL Molecular Graphics System, Version 1.7.2.1 2014; Schrödinger LLC: New York, NY, USA.
  45. Wolber, G.; Langer, T. LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters. J. Chem. Inf. Model. 2005, 45, 160–169. [Google Scholar] [CrossRef]
  46. Wolber, G.; Dornhofer, A.A.; Langer, T. Efficient Overlay of Small Organic Molecules Using 3D Pharmacophores. J. Comput. Aided. Mol. Des. 2006, 20, 773–788. [Google Scholar] [CrossRef]
  47. Corona, A.; Meleddu, R.; Esposito, F.; Distinto, S.; Bianco, G.; Masaoka, T.; Maccioni, E.; Menéndez-Arias, L.; Alcaro, S.; Le Grice, S.F.J.; et al. Ribonuclease H/DNA Polymerase HIV-1 Reverse Transcriptase Dual Inhibitor: Mechanistic Studies on the Allosteric Mode of Action of Isatin-Based Compound RMNC6. PLoS ONE 2016, 11, e0147225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Pala, N.; Esposito, F.; Rogolino, D.; Carcelli, M.; Sanna, V.; Palomba, M.; Naesens, L.; Corona, A.; Grandi, N.; Tramontano, E.; et al. Inhibitory effect of 2,3,5,6-tetrafluoro-4-[4-(Aryl)-1H-1,2,3-triazol-1-yl]benzenesulfonamide derivatives on HIV reverse transcriptase associated RNase H activities. Int. J. Mol. Sci. 2016, 17, 1371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Costa, G.; Rocca, R.; Corona, A.; Grandi, N.; Moraca, F.; Romeo, I.; Talarico, C.; Giovanna, M.; Alessandra, F.; Ortuso, F.; et al. Novel natural non-nucleoside inhibitors of HIV-1 reverse transcriptase identified by shape- and structure-based virtual screening techniques. Eur. J. Med. Chem. 2019, 161, 1–10. [Google Scholar] [CrossRef]
  50. Rhee, S.-Y.; Kantor, R.; Katzenstein, D.A.; Camacho, R.; Morris, L.; Sirivichayakul, S.; Jorgensen, L.; Brigido, L.F.; Schapiro, J.M.; Shafer, R.W.; et al. HIV-1 pol mutation frequency by subtype and treatment experience: Extension of the HIVseq program to seven non-B subtypes. AIDS 2006, 20, 643–651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Rhee, S.-Y.; Gonzales, M.J.; Kantor, R.; Betts, B.J.; Ravela, J.; Shafer, R.W. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003, 31, 298–303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef]
  53. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  54. Distinto, S.; Maccioni, E.; Meleddu, R.; Corona, A.; Alcaro, S.; Tramontano, E. Molecular Aspects of the RT/drug Interactions. Perspective of Dual Inhibitors. Curr. Pharm. Des. 2013, 19, 1850–1859. [Google Scholar] [CrossRef]
  55. Himmel, D.M.; Maegley, K.A.; Pauly, T.A.; Bauman, J.D.; Das, K.; Dharia, C.; Clark, A.D.; Ryan, K.; Hickey, M.J.; Love, R.A.; et al. Structure of HIV-1 reverse transcriptase with the inhibitor beta-Thujaplicinol bound at the RNase H active site. Structure 2009, 17, 1625–1635. [Google Scholar] [CrossRef] [Green Version]
  56. Kirby, K.A.; Myshakina, N.A.; Christen, M.T.; Chen, Y.L.; Schmidt, H.A.; Huber, A.D.; Xi, Z.; Kim, S.; Rao, R.K.; Kramer, S.T.; et al. A 2-hydroxyisoquinoline-1,3-dione active-site RNase H inhibitor binds in multiple modes to HIV-1 reverse transcriptase. Antimicrob. Agents Chemother. 2017, 61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Tramontano, E.; Esposito, F.; Badas, R.; Di Santo, R.; Costi, R.; La Colla, P. 6-[1-(4-Fluorophenyl)methyl-1H-pyrrol-2-yl)]-2,4-dioxo-5-hexenoic acid ethyl ester a novel diketo acid derivative which selectively inhibits the HIV-1 viral replication in cell culture and the ribonuclease H activity in vitro. Antivir. Res. 2005, 65, 117–124. [Google Scholar] [CrossRef] [PubMed]
  58. Boyer, P.L.; Smith, S.J.; Zhao, X.Z.; Das, K.; Gruber, K.; Arnold, E.; Burke, T.R.; Hughes, S.H. Developing and Evaluating Inhibitors against the RNase H Active Site of HIV-1 Reverse Transcriptase. J. Virol. 2018, 92, 1–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Arion, D.; Sluis-Cremer, N.; Min, K.-L.L.; Abram, M.E.; Fletcher, R.S.; Parniak, M.A. Mutational analysis of Tyr-501 of HIV-1 reverse transcriptase. Effects on ribonuclease H activity and inhibition of this activity by N-acylhydrazones. J. Biol. Chem. 2002, 277, 1370–1374. [Google Scholar] [CrossRef] [Green Version]
  60. Sluis-Cremer, N. Destabilization of the HIV-1 Reverse Transcriptase Dimer upon Interaction with N-Acyl Hydrazone Inhibitors. Mol. Pharmacol. 2002, 62, 398–405. [Google Scholar] [CrossRef] [Green Version]
  61. Gao, P.; Wang, X.; Sun, L.; Cheng, X.; Poongavanam, V.; Kongsted, J.; Álvarez, M.; Luczkowiak, J.; Pannecouque, C.; De Clercq, E.; et al. Design, synthesis, and biologic evaluation of novel galloyl derivatives as HIV-1 RNase H inhibitors. Chem. Biol. Drug Des. 2019, 93, 582–589. [Google Scholar] [CrossRef] [Green Version]
  62. Wang, L.; Tang, J.; Huber, A.D.; Casey, M.C.; Kirby, K.A.; Wilson, D.J.; Kankanala, J.; Parniak, M.A.; Sarafianos, S.G.; Wang, Z. 6-Biphenylmethyl-3-hydroxypyrimidine-2,4-diones potently and selectively inhibited HIV reverse transcriptase-associated RNase H. Eur. J. Med. Chem. 2018, 156, 680–691. [Google Scholar] [CrossRef]
  63. Corona, A.; Desantis, J.; Massari, S.; Distinto, S.; Masaoka, T.; Sabatini, S.; Esposito, F.; Manfroni, G.; Maccioni, E.; Cecchetti, V.; et al. Studies on Cycloheptathiophene-3-carboxamide Derivatives as Allosteric HIV-1 Ribonuclease H Inhibitors. ChemMedChem 2016, 11, 1709–1720. [Google Scholar] [CrossRef]
  64. Klumpp, K. Two-metal ion mechanism of RNA cleavage by HIV RNase H and mechanism-based design of selective HIV RNase H inhibitors. Nucleic Acids Res. 2003, 31, 6852–6859. [Google Scholar] [CrossRef] [Green Version]
Viruses 12 00729 g001 550
Figure 1. Structure of N-acylhydrazone analogs.
Figure 1. Structure of N-acylhydrazone analogs.
Viruses 12 00729 g001
Viruses 12 00729 g002 550
Figure 2. Analysis of compound 13 target of inhibition. (A) Time-of-addition. Drugs were used at the following final concentrations: AZT 3.7 μM; NVP 3.7 μM; EFV 0.3 μM; RAL 2 25 μM; compound 13 25 μM; compound 21 25 μM. Twenty-four hours post infection GFP expression was quantified by FACS. A representative assay of two is shown; (B) comparative inhibition of HIV-1 RT-associated activities. RT-associated RNase H activity black line, RDDP blue line. Assays were run in duplicate; the graph represents values the mean and standard deviation of two independent experiments.
Figure 2. Analysis of compound 13 target of inhibition. (A) Time-of-addition. Drugs were used at the following final concentrations: AZT 3.7 μM; NVP 3.7 μM; EFV 0.3 μM; RAL 2 25 μM; compound 13 25 μM; compound 21 25 μM. Twenty-four hours post infection GFP expression was quantified by FACS. A representative assay of two is shown; (B) comparative inhibition of HIV-1 RT-associated activities. RT-associated RNase H activity black line, RDDP blue line. Assays were run in duplicate; the graph represents values the mean and standard deviation of two independent experiments.
Viruses 12 00729 g002
Viruses 12 00729 g003 550
Figure 3. Putative binding mode of compound 13 within the HIV-1 RT RNase H domain. (A–C) Compound 13 in cyan depicted in the wt-RT: p66 in green, p51 in orange; (B–D) 2D representation of compound 13 in complex with wt RT with indication of binding pocket interacting residues.
Figure 3. Putative binding mode of compound 13 within the HIV-1 RT RNase H domain. (A–C) Compound 13 in cyan depicted in the wt-RT: p66 in green, p51 in orange; (B–D) 2D representation of compound 13 in complex with wt RT with indication of binding pocket interacting residues.
Viruses 12 00729 g003
Viruses 12 00729 g004 550
Figure 4. Inhibitory profile of compound 13 against the HIV-1 RT RNase H function. (A) Compound 13 inhibition of HIV-1 RT RNase H activities; (B) RDS1759 inhibition of HIV-1 RT RNase H activities. Data represent the mean and SD of two experiments.
Figure 4. Inhibitory profile of compound 13 against the HIV-1 RT RNase H function. (A) Compound 13 inhibition of HIV-1 RT RNase H activities; (B) RDS1759 inhibition of HIV-1 RT RNase H activities. Data represent the mean and SD of two experiments.
Viruses 12 00729 g004
Viruses 12 00729 g005 550
Figure 5. Conservation profile of HIV-1 RT RNase H domain. (A) Logo of degree of conservation of amino acid residues: blue hydrophobic lateral chains, green neutral lateral chains, black hydrophilic lateral chains; (B) color scheme of the percentage of variability of HIV-1 RT RNase H domain (PDB code 3K2P). The color scale of the percentage of variability goes form red (low percentage of variation) to green (high percentage of variation).
Figure 5. Conservation profile of HIV-1 RT RNase H domain. (A) Logo of degree of conservation of amino acid residues: blue hydrophobic lateral chains, green neutral lateral chains, black hydrophilic lateral chains; (B) color scheme of the percentage of variability of HIV-1 RT RNase H domain (PDB code 3K2P). The color scale of the percentage of variability goes form red (low percentage of variation) to green (high percentage of variation).
Viruses 12 00729 g005
Viruses 12 00729 g006 550
Figure 6. (A) Share ligand-based pharmacophore model; (B) Share-based pharmacophore model with excluded volumes; (C) pharmacophore matching of compound 13 and 21. Red sphere, H bond acceptor; green sphere and arrow, H bond donor; blue circle, aromatic feature.
Figure 6. (A) Share ligand-based pharmacophore model; (B) Share-based pharmacophore model with excluded volumes; (C) pharmacophore matching of compound 13 and 21. Red sphere, H bond acceptor; green sphere and arrow, H bond donor; blue circle, aromatic feature.
Viruses 12 00729 g006
Table
Table 1. Antiviral activity of compounds against wt NL4-3 and NNRTI-resistant HIV-1 strains.
Table 1. Antiviral activity of compounds against wt NL4-3 and NNRTI-resistant HIV-1 strains.
HIV-1
wt
MT
EC50 1 (µM)		HIV-1
K103N-Y181C
MT4
EC50 1 (µM)		MT4
CC50 2
(µM)		SI 3
1310.1 ± 4.712.4 ± 0.861.5 ± 2.26.1
215.0 ± 1.48.7 ± 0.131.8 ± 0.46.4
RDS17598.2 ± 0.94.6 ± 1.139 ± 4.44.7
AZT0.007 ± 0.0050.0031 ± 0.0007>3.7>500
EFV0.050 ± 0.012>1.58>0.31>6.25
NVP0.67 ± 0.22>3.7>3.7>5.5
RAL0.0022 ± 0.00010.004 ± 0.001>2.25>1987
1 EC50: effective concentration 50 or concentration of compound required to inhibit 50% of HIV-induced cell death, evaluated with the MTT method in MT4 cells; 2 CC50: cytotoxic concentration 50 or concentration of compound required to induce 50% of death of non-infected cells, evaluated with the MTT method in MT4 cells; 3 SI: selectivity index, indicates the ratio between the CC50 and the EC50; AZT: azidothymidine, EFV: efavirenz, NVP: nevirapine, RAL: raltegravir; Data represent mean and SD of three independent experiments.

Share and Cite

MDPI and ACS Style

Corona, A.; Ballana, E.; Distinto, S.; Rogolino, D.; Del Vecchio, C.; Carcelli, M.; Badia, R.; Riveira-Muñoz, E.; Esposito, F.; Parolin, C.; et al. Targeting HIV-1 RNase H: N’-(2-Hydroxy-benzylidene)-3,4,5-Trihydroxybenzoylhydrazone as Selective Inhibitor Active against NNRTIs-Resistant Variants. Viruses 2020, 12, 729. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070729

AMA Style

Corona A, Ballana E, Distinto S, Rogolino D, Del Vecchio C, Carcelli M, Badia R, Riveira-Muñoz E, Esposito F, Parolin C, et al. Targeting HIV-1 RNase H: N’-(2-Hydroxy-benzylidene)-3,4,5-Trihydroxybenzoylhydrazone as Selective Inhibitor Active against NNRTIs-Resistant Variants. Viruses. 2020; 12(7):729. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070729

Chicago/Turabian Style

Corona, Angela, Ester Ballana, Simona Distinto, Dominga Rogolino, Claudia Del Vecchio, Mauro Carcelli, Roger Badia, Eva Riveira-Muñoz, Francesca Esposito, Cristina Parolin, and et al. 2020. "Targeting HIV-1 RNase H: N’-(2-Hydroxy-benzylidene)-3,4,5-Trihydroxybenzoylhydrazone as Selective Inhibitor Active against NNRTIs-Resistant Variants" Viruses 12, no. 7: 729. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070729

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop