Drug Repurposing Identifies Inhibitors of Oseltamivir-Resistant Influenza Viruses
The inhibitor of neuraminidase (NA), oseltamivir, is a widely used anti-influenza drug. However, oseltamivir-resistant H1N1 influenza viruses carrying the H275Y NA mutation spontaneously emerged as a result of natural genetic drift and drug treatment. Because this and other potential mutations may generate a future pandemic influenza strain that is oseltamivir-resistant, alternative therapy options are needed. Here show that a structure-based computational method can be used to identify existing drugs that inhibit resistant viruses, thereby providing a first line of pharmaceutical defense against this very possible scenario. We identified two drugs, nalidixic acid and dorzolamide, that potently inhibit the NA activity of oseltamivir-resistant H1N1 viruses with the H275Y NA mutation at very low concentrations, but have no effect on wild-type H1N1 NA even at a much higher concentration, suggesting that the oseltamivir-resistance mutation itself caused susceptibility to these drugs.
<p>During the early stages of any influenza pandemic, specific vaccines are unavailable; therefore, effective drugs are the only means of preventing the spread of infection.[1] The 2009 H1N1 pandemic, for example, reached its peak well before a vaccine became available in the United States.[2] The potent neuraminidases (NA) inhibitor oseltamivir (Tamiflu) is the most widely used anti-influenza drug. However, oseltamivir-resistant H1N1 influenza viruses carrying the H275Y NA mutation spontaneously emerged worldwide as a result of natural genetic drift and drug treatment.[3] NA catalyzes the hydrolysis of sialic acid residues between newly formed virions and host glycoproteins to allow virion release. Oseltamivir was designed to fit perfectly into the NA enzyme active site of influenza virus. This site is altered in oseltamivir-resistant viruses, inhibiting binding of the drug. Most oseltamivir-resistant strains of seasonal and pandemic H1N1 influenza virus to date have carried the H275Y NA mutation.[4] This mutant has a bulky hydrophobic side chain at Y275 that forces the E277 side chain to rotate, disrupting the intramolecular E277-R225 electrostatic interactions normally formed upon binding to oseltamivir. The E277 side chain is also forced toward the center of the NA sialic acid binding site by the side chain of Y275. These changes reduce the binding affinity of the site for the large aliphatic side chain of oseltamivir.[5] Had oseltamivir-resistant viruses dominated during the 2009 H1N1 pandemic mortality would have been much greater.</p><p>We reasoned that available drugs that are structurally similar to oseltamivir may bind to the mutant NA active site and that sensitivity to these drugs may be enhanced by the structural changes in the mutant NA binding site. To search for such drugs, we used hierarchical computational screening.[6] We first conducted combined 2D/3D shape screening[7] of the library of current FDA-approved drugs.[8] Using a 0 (no overlap) to 1 (identical shape) scale, we selected approximately 300 compounds whose structural shapes have similarity > 0.7 to that of oseltamivir. We then docked each of these compounds to the active site of the H275Y-mutant N1 NA, using the conformation of oseltamivir bound to the wild-type N1 NA at this site[5a, 5c] as the reference. The docking results were visually inspected to ensure that the key ("hotspot"[9]) interactions between the bound drugs and H257Y-mutant N1 NA in the docked complexes had a maximum overlap with the key interactions between the bound oseltamivir and wild-type N1 NA. The 20 compounds with the highest docking scores were selected, and their structures in complex with H275Y N1 NA were refined by molecular dynamics (MD) simulation studies[10] using 5-ns unrestrained simulation with explicit solvation and calculation of molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy.[11] The 10 compounds with the highest calculated binding free energy were selected for further testing.</p><p>The 10 compounds were tested against two oseltamivir-sensitive (A/Brisbane/59/2007 and A/Netherlands/049/2008) and two oseltamivir-resistant (A/New Jersey/15/2007 and A/Netherlands/028/2008) H1N1 influenza viruses. The NA sequences of all four strains were confirmed in our laboratory and showed >85% similarity to those of the 1918 pandemic H1N1, 2009 pandemic H1N1, and highly pathogenic H5N1 viruses. The NA activity of the viruses and the inhibitory constants of the selected compounds were measured by using an established fluorescence–based NA enzyme inhibition assay.[12]</p><p>Although several drugs inhibited the NA activity of the resistant and nonresistant viruses to some extent (see Table S1 in the Supporting Information), nalidixic acid and dorzolamide exerted a selective inhibitory effect against the two oseltamivir-resistant strains. In additional confirmatory testing (four independent measurements at each concentration), both compounds stably inhibited the NA activity of both oseltamivir-resistant strains in a dose-dependent manner and at a sub-micromolar 50% inhibitory concentration (IC50) but did not inhibit the NA activity of the oseltamivir-susceptible strains (Figure 1). Of the two compounds, nalidixic acid is of particular interest, as it was the first synthetic quinolone antibiotic and is widely used in treating bacterial infections,[13] while dorzolamide is most used in eye drop for glaucoma treatment.</p><p>We used MD simulation studies to examine the binding of the two compounds to the H275Y NA mutant in detail. In both complexes, a salt bridge (observed in neither starting structure) formed between E277 and R225 of NA after approximately 500 ps; this salt bridge was previously observed in wild-type NA (but not H275Y mutant NA) bound to oseltamivir[5b ] (Figure 2). The bound conformations of nalidixic acid and dorzolamide closely resembled that of oseltamivir with the exception of the oseltamivir hydrophobic side chain, which sterically hindered its binding to H275Y N1 NA (Figure 2 and 3).</p><p>Oseltamivir interacts with N1 NA at several key hotspots; it forms salt bridges with three arginines (R119, R293 and R372) and E120 and forms a hydrogen bond with R153 (Figure 3). The acidic group of nalidixic acid and the sulfonamide group of dorzolamide bind to the H275Y N1 NA in a manner similar to that of the oseltamivir acidic group, interacting with three arginine residues (R119, R293 and R372) (Figure 2). In addition, the secondary amine group of dorzolamide forms salt bridges with D152 and E120 (which interact with the primary amine group of oseltamivir), and the dorzolamide sulfone group forms a hydrogen bond with the side chain of R153, as does the amide group of oseltamivir (Figure 2b). Although some of the interactions observed between oseltamivir and wild-type NA are lacking in the binding of nalidixic acid and dorzolamide to the H275Y N1 NA mutant active site, restoration of the E277-R225 salt bridge appears to generate a similar binding affinity.</p><p>To further investigate why the two drugs do not bind to wild-type N1 NA, we used MD simulations and MM-PBSA binding calculations. Nalidixic acid remained stably within the sialic acid binding site of wild-type NA, but its calculated binding free energy, which was 5 kcal/mol less than that of binding to H275Y N1 NA, demonstrated ineffective binding. This finding may be explained by the absence of specific electrostatic interactions, which are important in determining the binding orientation of wild-type N1 NA despite their small contribution to binding energetics.[14] In the case of dorzolamide, its hydrophilic nature caused wide fluctuation of its position in wild-type N1 NA binding sites during MD simulations, indicating the absence of binding to wild-type N1 NA.</p><p>Influenza-specific NA inhibitors are a unique class of drugs[15] designed to mimic sialic acid, the influenza NA substrate. Our 2D/3D shape screening results suggest that the chemical space and the spatial shape of this sialic acid mimetic overlaps with those of the other two drugs. Our structural studies suggest that off-targets of current NA inhibitors and other scaffolds for potent NA inhibitors are rarely reported because wild-type N1 NA is markedly sensitive to substrate shape and specific functional groups. Indeed, the fact that nalidixic acid and dorzolamide inhibit only oseltamivir-resistant N1 NA demonstrates that the slight alteration of the substrate recognition profile of N1 NA caused by the H275Y mutation enhances its vulnerability to some compounds that do not bind efficiently to wild-type N1 NA. On the other hand, because the N1 NA active site is highly conserved in both the wild-type and oseltamivir-resistant strains, the structure of an oseltamivir-resistant mutant NA can be modeled on the basis of the wild-type structure and the model can then be used in structure-based ligand screening to identify small-molecule drugs (e.g., dorzolamide and nalidixic acid) that target the mutant strain.</p><p>This proof of principle study demonstrates that drug resistance caused by a specific influenza N1 NA mutation (H275Y in this case) can be circumvented by using drug repositioning and structure-based computational methods to identify alternative agents. Although it is difficult to predict the site(s) of future NA mutations that may cause oseltamivir resistanc,[16] the conformational changes induced by these mutations can be predicted after the NA is sequenced.[5b] Therefore, our protocol can be used to circumvent future NA mutations that cause oseltamivir resistance and may be useful in rapidly controlling the spread of oseltamivir-resistant pandemic influenza viruses.[17] This approach also allows the identification of existing drugs that bind and inhibit various viruses (or microorganisms) that have acquired drug-resistance mutations.</p>