A conserved hydrogen-bonding network of P2 bis-tetrahydrofuran-containing HIV-1 protease inhibitors (PIs) with a protease active-site amino acid backbone aids in their activity against PI-resistant HIV

Beyond darunavir: recent development of next generation HIV-1 protease inhibitors to combat drug resistance

We report our recent development of a conceptually new generation of exceptionally potent non-peptidic HIV-1 protease inhibitors that displayed excellent pharmacological and drug-resistance profiles. Our X-ray structural studies of darunavir and other designed inhibitors from our laboratories led us to create a variety of inhibitors incorporating fused ring polycyclic ethers and aromatic heterocycles to promote hydrogen bonding interactions with the backbone atoms of HIV-1 protease as well as van der Waals interactions with residues in the S2 and S2' subsites. We have also incorporated specific functionalities to enhance van der Waals interactions in the S1 and S1' subsites. The combined effects of these structural templates are critical to the inhibitors' exceptional potency and drug-like properties. We highlight here our molecular design strategies to promote backbone hydrogen bonding interactions to combat drug-resistance and specific design of polycyclic ether templates to mimic peptide-like bonds in the HIV-1 protease active site. Our medicinal chemistry and drug development efforts led to the development of new generation inhibitors significantly improved over darunavir and displaying unprecedented antiviral activity against multidrug-resistant HIV-1 variants.

N-Sulfonyl dipeptide nitriles as inhibitors of human cathepsin S: In silico design, synthesis and biochemical characterization

A library of cathepsin S inhibitors of the dipeptide nitrile chemotype, bearing a bioisosteric sulfonamide moiety, was synthesized. Kinetic investigations were performed at four human cysteine proteases, i.e. cathepsins S, B, K and L. Compound 12 with a terminal 3-biphenyl sulfonamide substituent was the most potent (K_i = 4.02 nM; selectivity ratio cathepsin S/K = 5.8; S/L = 67) and 24 with a 4'-fluoro-4-biphenyl sulfonamide substituent the most selective cathepsin S inhibitor (K_i = 35.5 nM; selectivity ratio cathepsin S/K = 57; S/L = 31). In silico design and biochemical evaluation emphasized the impact of the sulfonamide linkage on selectivity and a possible switch of P2 and P3 substituents with respect to the occupation of the corresponding binding sites of cathepsin S.

Single atom changes in newly synthesized HIV protease inhibitors reveal structural basis for extreme affinity, high genetic barrier, and adaptation to the HIV protease plasticity

HIV-1 protease inhibitors (PIs), such as darunavir (DRV), are the key component of antiretroviral therapy. However, HIV-1 often acquires resistance to PIs. Here, seven novel PIs were synthesized, by introducing single atom changes such as an exchange of a sulfur to an oxygen, scission of a single bond in P2′-cyclopropylaminobenzothiazole (or -oxazole), and/or P1-benzene ring with fluorine scan of mono- or bis-fluorine atoms around DRV’s scaffold. X-ray structural analyses of the PIs complexed with wild-type Protease (PR_WT) and highly-multi-PI-resistance-associated PR_DRV^R_P51 revealed that the PIs better adapt to structural plasticity in PR with resistance-associated amino acid substitutions by formation of optimal sulfur bond and adaptation of cyclopropyl ring in the S2′-subsite. Furthermore, these PIs displayed increased cell permeability and extreme anti-HIV-1 potency compared to DRV. Our work provides the basis for developing novel PIs with high potency against PI-resistant HIV-1 variants with a high genetic barrier.

Structural Adaptation of Darunavir Analogues against Primary Mutations in HIV-1 Protease

HIV-1 protease is one of the prime targets of agents used in antiretroviral therapy against HIV. However, under selective pressure of protease inhibitors, primary mutations at the active site weaken inhibitor binding to confer resistance. Darunavir (DRV) is the most potent HIV-1 protease inhibitor in clinic; resistance is limited, as DRV fits well within the substrate envelope. Nevertheless, resistance is observed due to hydrophobic changes at residues including I50, V82, and I84 that line the S1/S1' pocket within the active site. Through enzyme inhibition assays and a series of 12 crystal structures, we interrogated susceptibility of DRV and two potent analogues to primary S1' mutations. The analogues had modifications at the hydrophobic P1' moiety compared to DRV to better occupy the unexploited space in the S1' pocket where the primary mutations were located. Considerable losses of potency were observed against protease variants with I84V and I50V mutations for all three inhibitors. The crystal structures revealed an unexpected conformational change in the flap region of I50V protease bound to the analogue with the largest P1' moiety, indicating interdependency between the S1' subsite and the flap region. Collective analysis of protease-inhibitor interactions in the crystal structures using principle component analysis was able to distinguish inhibitor identity and relative potency solely based on van der Waals contacts. Our results reveal the complexity of the interplay between inhibitor P1' moiety and S1' mutations and validate principle component analyses as a useful tool for distinguishing resistance and inhibitor potency.

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

Design, Synthesis, Biological Evaluation, and X-ray Studies of HIV-1 Protease Inhibitors with Modified P2' Ligands of Darunavir

The structure-based design, synthesis, and biological evaluation of a series of nonpeptidic HIV-1 protease inhibitors with rationally designed P2' ligands are described. The inhibitors are designed to enhance backbone binding interactions, particularly at the S2' subsite. Synthesis of inhibitors was carried out efficiently. The stereochemistry of alcohol functionalities of the P2' ligands was set by asymmetric reduction of the corresponding ketone using (R,R)- or (S,S)-Noyori catalysts. A number of inhibitors displayed very potent enzyme inhibitory and antiviral activity. Inhibitors 3g and 3h showed enzyme K_i values of 27.9 and 49.7 pm and antiviral activity of 6.2 and 3.9 nm, respectively. These inhibitors also remained quite potent against darunavir-resistant HIV-1 variants. An X-ray structure of inhibitor 3g in complex with HIV-1 protease revealed key interactions in the S2' subsite.

GRL-09510, a Unique P2-Crown-Tetrahydrofuranylurethane -Containing HIV-1 Protease Inhibitor, Maintains Its Favorable Antiviral Activity against Highly-Drug-Resistant HIV-1 Variants in vitro

We report that GRL-09510, a novel HIV-1 protease inhibitor (PI) containing a newly-generated P2-crown-tetrahydrofuranylurethane (Crwn-THF), a P2′-methoxybenzene, and a sulfonamide isostere, is highly active against laboratory and primary clinical HIV-1 isolates (EC₅₀: 0.0014–0.0028 μM) with minimal cytotoxicity (CC₅₀: 39.0 μM). Similarly, GRL-09510 efficiently blocked the replication of HIV-1_NL4-3 variants, which were capable of propagating at high-concentrations of atazanavir, lopinavir, and amprenavir (APV). GRL-09510 was also potent against multi-drug-resistant clinical HIV-1 variants and HIV-2_ROD. Under the selection condition, where HIV-1_NL4-3 rapidly acquired significant resistance to APV, an integrase inhibitor raltegravir, and a GRL-09510 congener (GRL-09610), no variants highly resistant against GRL-09510 emerged over long-term in vitro passage of the virus. Crystallographic analysis demonstrated that the Crwn-THF moiety of GRL-09510 forms strong hydrogen-bond-interactions with HIV-1 protease (PR) active-site amino acids and is bulkier with a larger contact surface, making greater van der Waals contacts with PR than the bis-THF moiety of darunavir. The present data demonstrate that GRL-09510 has favorable features for treating patients infected with wild-type and/or multi-drug-resistant HIV-1 variants, that the newly generated P2-Crwn-THF moiety confers highly desirable anti-HIV-1 potency. The use of the novel Crwn-THF moiety sheds lights in the design of novel PIs.

Investigation of the binding free energies of FDA approved drugs against subtype B and C-SA HIV PR: ONIOM approach

Human immune virus subtype C is the most widely spread HIV subtype in Sub-Sahara Africa and South Africa. A profound structural insight on finding potential lead compounds is therefore necessary for drug discovery. The focus of this study is to rationalize the nine Food and Drugs Administration (FDA) HIV antiviral drugs complexed to subtype B and C-SA PR using ONIOM approach. To achieve this, an integrated two-layered ONIOM model was used to optimize the geometrics of the FDA approved HIV-1 PR inhibitors for subtype B. In our hybrid ONIOM model, the HIV-1 PR inhibitors as well as the ASP 25/25' catalytic active residues were treated at high level quantum mechanics (QM) theory using B3LYP/6-31G(d), and the remaining HIV PR residues were considered using the AMBER force field. The experimental binding energies of the PR inhibitors were compared to the ONIOM calculated results. The theoretical binding free energies (?G_bind) for subtype B follow a similar trend to the experimental results, with one exemption. The computational model was less suitable for C-SA PR. Analysis of the results provided valuable information about the shortcomings of this approach. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide much improved binding energies for complex enzyme drug interactions.

Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution

Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.

Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS

HIV-1 protease inhibitors continue to play an important role in the treatment of HIV/AIDS, transforming this deadly ailment into a more manageable chronic infection. Over the years, intensive research has led to a variety of approved protease inhibitors for the treatment of HIV/AIDS. In this review, we outline current drug design and medicinal chemistry efforts toward the development of next-generation protease inhibitors beyond the currently approved drugs.

Highly resistant HIV-1 proteases and strategies for their inhibition

The virally encoded protease is an important drug target for AIDS therapy. Despite the potency of the current drugs, infections with resistant viral strains limit the long-term effectiveness of therapy. Highly resistant variants of HIV protease from clinical isolates have different combinations of about 20 mutations and several orders of magnitude worse binding affinity for clinical inhibitors. Strategies are being explored to inhibit these highly resistant mutants. The existing inhibitors can be modified by introducing groups with the potential to form new interactions with conserved protease residues, and the flexible flaps. Alternative strategies are discussed, including designing inhibitors to bind to the open conformation of the protease dimer, and inhibition of the protease-catalyzed processing of the Gag-Pol precursor.

C-5-Modified Tetrahydropyrano-Tetrahydofuran-Derived Protease Inhibitors (PIs) Exert Potent Inhibition of the Replication of HIV-1 Variants Highly Resistant to Various PIs, including Darunavir

We identified three nonpeptidic HIV-1 protease inhibitors (PIs), GRL-015, -085, and -097, containing tetrahydropyrano-tetrahydrofuran (Tp-THF) with a C-5 hydroxyl. The three compounds were potent against a wild-type laboratory HIV-1 strain (HIV-1(WT)), with 50% effective concentrations (EC50s) of 3.0 to 49 nM, and exhibited minimal cytotoxicity, with 50% cytotoxic concentrations (CC50) for GRL-015, -085, and -097 of 80, >100, and >100 μM, respectively. All the three compounds potently inhibited the replication of highly PI-resistant HIV-1 variants selected with each of the currently available PIs and recombinant clinical HIV-1 isolates obtained from patients harboring multidrug-resistant HIV-1 variants (HIVMDR). Importantly, darunavir (DRV) was >1,000 times less active against a highly DRV-resistant HIV-1 variant (HIV-1DRV(R) P51); the three compounds remained active against HIV-1DRV(R) P51 with only a 6.8- to 68-fold reduction. Moreover, the emergence of HIV-1 variants resistant to the three compounds was considerably delayed compared to the case of DRV. In particular, HIV-1 variants resistant to GRL-085 and -097 did not emerge even when two different highly DRV-resistant HIV-1 variants were used as a starting population. In the structural analyses, Tp-THF of GRL-015, -085, and -097 showed strong hydrogen bond interactions with the backbone atoms of active-site amino acid residues (Asp29 and Asp30) of HIV-1 protease. A strong hydrogen bonding formation between the hydroxyl moiety of Tp-THF and a carbonyl oxygen atom of Gly48 was newly identified. The present findings indicate that the three compounds warrant further study as possible therapeutic agents for treating individuals harboring wild-type HIV and/or HIVMDR.

IMPORTANCE

Darunavir (DRV) inhibits the replication of most existing multidrug-resistant HIV-1 strains and has a high genetic barrier. However, the emergence of highly DRV-resistant HIV-1 strains (HIVDRV(R) ) has recently been observed in vivo and in vitro. Here, we identified three novel HIV-1 protease inhibitors (PIs) containing a tetrahydropyrano-tetrahydrofuran (Tp-THF) moiety with a C-5 hydroxyl (GRL-015, -085, and -097) which potently suppress the replication of HIVDRV(R) . Moreover, the emergence of HIV-1 strains resistant to the three compounds was considerably delayed compared to the case of DRV. The C-5 hydroxyl formed a strong hydrogen bonding interaction with the carbonyl oxygen atom of Gly48 of protease as examined in the structural analyses. Interestingly, a compound with Tp-THF lacking the hydroxyl moiety substantially decreased activity against HIVDRV(R) . The three novel compounds should be further developed as potential drugs for treating individuals harboring wild-type and multi-PI-resistant HIV variants as well as HIVDRV(R) .

Aureonitol, a Fungi-Derived Tetrahydrofuran, Inhibits Influenza Replication by Targeting Its Surface Glycoprotein Hemagglutinin

The influenza virus causes acute respiratory infections, leading to high morbidity and mortality in groups of patients at higher risk. Antiviral drugs represent the first line of defense against influenza, both for seasonal infections and pandemic outbreaks. Two main classes of drugs against influenza are in clinical use: M2-channel blockers and neuraminidase inhibitors. Nevertheless, because influenza strains that are resistant to these antivirals have been described, the search for novel compounds with different mechanisms of action is necessary. Here, we investigated the anti-influenza activity of a fungi-derived natural product, aureonitol. This compound inhibited influenza A and B virus replication. This compound was more effective against influenza A(H3N2), with an EC₅₀ of 100 nM. Aureonitol cytoxicity was also very low, with a CC₅₀ value of 1426 μM. Aureonitol inhibited influenza hemagglutination and, consequently, significantly impaired virus adsorption. Molecular modeling studies revealed that aureonitol docked in the sialic acid binding site of hemagglutinin, forming hydrogen bonds with highly conserved residues. Altogether, our results indicate that the chemical structure of aureonitol is promising for future anti-influenza drug design.

The L33F darunavir resistance mutation acts as a molecular anchor reducing the flexibility of the HIV-1 protease 30s and 80s loops

HIV-1 protease (PR) is a 99 amino acid protein responsible for proteolytic processing of the viral polyprotein - an essential step in the HIV-1 life cycle. Drug resistance mutations in PR that are selected during antiretroviral therapy lead to reduced efficacy of protease inhibitors (PI) including darunavir (DRV). To identify the structural mechanisms associated with the DRV resistance mutation L33F, we performed X-ray crystallographic studies with a multi-drug resistant HIV-1 protease isolate that contains the L33F mutation (MDR769 L33F). In contrast to other PR L33F DRV complexes, the structure of MDR769 L33F complexed with DRV reported here displays the protease flaps in an open conformation. The L33F mutation increases noncovalent interactions in the hydrophobic pocket of the PR compared to the wild-type (WT) structure. As a result, L33F appears to act as a molecular anchor, reducing the flexibility of the 30s loop (residues 29-35) and the 80s loop (residues 79-84). Molecular anchoring of the 30s and 80s loops leaves an open S1/S1' subsite and distorts the conserved hydrogen-bonding network of DRV. These findings are consistent with previous reports despite structural differences with regards to flap conformation.