Computational and experimental analysis of drug binding to the Influenza M2 channel

Accurate calculation of affinity changes to the close state of influenza A M2 transmembrane domain in response to subtle structural changes of adamantyl amines using free energy perturbation methods in different lipid bilayers

Experimental binding free energies of 27 adamantyl amines against the influenza M2(22-46) WT tetramer, in its closed form at pH 8, were measured by ITC in DPC micelles. The measured Kd's range is ~44 while the antiviral potencies (IC50) range is ~750 with a good correlation between binding free energies computed with Kd and IC50 values (r = 0.76). We explored with MD simulations (ff19sb, CHARMM36m) the binding profile of complexes with strong, moderate and weak binders embedded in DMPC, DPPC, POPC or a viral mimetic membrane and using different experimental starting structures of M2. To predict accurately differences in binding free energy in response to subtle changes in the structure of the ligands, we performed 18 alchemical perturbative single topology FEP/MD NPT simulations (OPLS2005) using the BAR estimator (Desmond software) and 20 dual topology calculations TI/MD NVT simulations (ff19sb) using the MBAR estimator (Amber software) for adamantyl amines in complex with M2(22-46) WT in DMPC, DPPC, POPC. We observed that both methods with all lipids show a very good correlation between the experimental and calculated relative binding free energies (r = 0.77-0.87, mue = 0.36-0.92 kcal mol-1) with the highest performance achieved with TI/MBAR and lowest performance with FEP/BAR in DMPC bilayers. When antiviral potencies are used instead of the Kd values for computing the experimental binding free energies we obtained also good performance with both FEP/BAR (r = 0.83, mue = 0.75 kcal mol-1) and TI/MBAR (r = 0.69, mue = 0.77 kcal mol-1).

Zika M—A Potential Viroporin: Mutational Study and Drug Repurposing

Genus Flavivirus contains several important human pathogens. Among these, the Zika virus is an emerging etiological agent that merits concern. One of its structural proteins, prM, plays an essential role in viral maturation and assembly, making it an attractive drug and vaccine development target. Herein, we have characterized ZikV-M as a potential viroporin candidate using three different bacteria-based assays. These assays were subsequently employed to screen a library of repurposed drugs from which ten compounds were identified as ZikV-M blockers. Mutational analyses of conserved amino acids in the transmembrane domain of other flaviviruses, including West Nile and Dengue virus, were performed to study their role in ion channel activity. In conclusion, our data show that ZikV-M is a potential ion channel that can be used as a drug target for high throughput screening and drug repurposing.

Effect of cholesterol on the membrane partitioning dynamics of hepatitis A virus-2B peptide

Understanding viral peptide detection and partitioning and the subsequent host membrane composition-based response is essential for gaining insights into the viral mechanism. Here, we probe the crucial role of the presence of membrane lipid packing defects, depending on the membrane composition, in allowing the viral peptide belonging to C-terminal Hepatitis A Virus-2B (HAV-2B) to detect, attach and subsequently partition into host cell membrane mimics. Using molecular dynamics simulations, we conclusively show that the hydrophobic residues in the viral peptide detect transiently present lipid packing defects, insert themselves into such defects, form anchor points and facilitate the partitioning of the peptide, thereby inducing membrane disruption. We also show that the presence of cholesterol significantly alters such lipid packing defects, both in size and in number, thus mitigating the partitioning of the membrane active viral peptide into cholesterol-rich membranes. Our results are in excellent agreement with previously published experimental data and further explain the role of lipid defects in understanding such data. These results show differential ways in which the presence and absence of cholesterol can alter the permeability of the host membranes to the membrane active peptide component of HAV-2B virus, via lipid packing defects, and can possibly be a part of the general membrane detection mechanism for viroporins.

Identification of SARS-CoV-2 E Channel Blockers from a Repurposed Drug Library

SARS-CoV-2, the etiological agent of the COVID-19 pandemic, is a member of the Coronaviridae family. It is an enveloped virus with ion channels in its membrane, the most characterized of which is the E protein. Therefore, in an attempt to identify blockers of the E channel, we screened a library of 2839 approved-for-human-use drugs. Our approach yielded eight compounds that exhibited appreciable activity in three bacteria-based channel assays. Considering the fact that the E channel is the most conserved of all SARS-CoV-2 proteins, any inhibitor of its activity may provide an option to curb the viral spread. In addition, inhibitors can also enhance our ability to understand the exact role played by the E protein during the infectivity cycle. Finally, detailed electrophysiological analyses, alongside in vitro and in vivo studies will be needed to establish the exact potential of each of the blockers identified in our study.

Increased Dissociation of Adamantanamines in Influenza A M2 S31N with Partial Block by Rimantadine

The ubiquitous mutation from serine (WT) to asparagine at residue 31 (S31N) in the influenza A M2 channel renders it insensitive to amantadine (AMT) and rimantadine (RMT) block, but it is unknown whether the inhibition results from weak binding or incomplete block. Two-electrode voltage clamp (TEVC) of transfected Xenopus oocytes revealed that the M2 S31N channel is essentially fully blocked by AMT at 10 mM, demonstrating that, albeit weak, AMT binding in a channel results in complete block of its proton current. In contrast, RMT achieves only a modest degree of block in the M2 S31N channel at 1 mM, with very little increase in block at 10 mM, indicating that the RMT binding site in the channel saturates with only modest block. From exponential curve fits to families of proton current wash-in and wash-out traces, the association rate constant (k₁) is somewhat decreased for both AMT and RMT in the S31N, but the dissociation rate constant (k₂) is dramatically increased compared with WT. The potentials of mean force (PMF) from adaptive biasing force (ABF) molecular dynamics simulations predict that rate constants should be exquisitely sensitive to the charge state of the His37 selectivity filter of M2. With one exception out of eight cases, predictions from the simulations with one and three charged side chains bracket the experimental rate constants, as expected for the acidic bath used in the TEVC assay. From simulations, the weak binding can be accounted for by changes in the potentials of mean force, but the partial block by RMT remains unexplained.

Virtual Screening Identifies Chebulagic Acid as an Inhibitor of the M2(S31N) Viral Ion Channel and Influenza A Virus

The increasing prevalence of drug-resistant influenza viruses emphasizes the need for new antiviral countermeasures. The M2 protein of influenza A is a proton-gated, proton-selective ion channel, which is essential for influenza replication and an established antiviral target. However, all currently circulating influenza A virus strains are now resistant to licensed M2-targeting adamantane drugs, primarily due to the widespread prevalence of an M2 variant encoding a serine to asparagine 31 mutation (S31N). To identify new chemical leads that may target M2(S31N), we performed a virtual screen of molecules from two natural product libraries and identified chebulagic acid as a candidate M2(S31N) inhibitor and influenza antiviral. Chebulagic acid selectively restores growth of M2(S31N)-expressing yeast. Molecular modeling also suggests that chebulagic acid hydrolysis fragments preferentially interact with the highly-conserved histidine residue within the pore of M2(S31N) but not adamantane-sensitive M2(S31). In contrast, chebulagic acid inhibits in vitro influenza A replication regardless of M2 sequence, suggesting that it also acts on other influenza targets. Taken together, results implicate chebulagic acid and/or its hydrolysis fragments as new chemical leads for M2(S31N) and influenza-directed antiviral development.

Investigation of the Drug Resistance Mechanism of M2-S31N Channel Blockers through Biomolecular Simulations and Viral Passage Experiments

Recent efforts in drug development against influenza A virus (IAV) M2 proton channel S31N mutant resulted in conjugates of amantadine linked with aryl head heterocycles. To understand the mechanism of drug resistance, we chose a representative M2-S31N inhibitor, compound 3, as a chemical probe to identify resistant mutants. To increase the possibility of identifying novel resistant mutants, serial viral passage experiments were performed with multiple strains of H1N1 and H3N2 viruses in different cell lines. This approach not only identified M2 mutations around the drug-binding site, including the pore-lining residues (V27A, V27F, N31S, and G34E) and an interhelical residue (I32N), but also a new allosteric mutation (R45H), in addition to L46P previously identified, located at the C-terminus of M2 that is more than 10 Å away from the drug-binding site. The effects of each mutation were next investigated using electrophysiology, recombinant viruses, and molecular dynamics (MD) simulations. The reduced sensitivity in channel blockage correlated with increased drug resistance in antiviral assays using recombinant viruses. The MD simulations show that the V27A, V27F, G34E, and R45H mutations increase the diameter and hydration state of the pore in complex with compound 3. The Molecular Mechanics Generalized Born (MM-GBSA) calculations result in more positive binding free energies for the complexes of resistant M2 (V27A, V27F, G34E, R45H) with compound 3 compared to the stable complexes (S31N and I32N). Overall, this is the first systematic study of the drug resistance mechanism of M2-S31N channel blockers using multiple viruses in different cell lines.

Potential Viroporin Candidates From Pathogenic Viruses Using Bacteria-Based Bioassays

Viroporins are a family of small hydrophobic proteins found in many enveloped viruses that are capable of ion transport. Building upon the ability to inhibit influenza by blocking its archetypical M2 H⁺ channel, as a family, viroporins may represent a viable target to curb viral infectivity. To this end, using three bacterial assays we analyzed six small hydrophobic proteins from biomedically important viruses as potential viroporin candidates. Our results indicate that Eastern equine encephalitis virus 6k, West Nile virus MgM, Dengue virus 2k, Dengue virus P1, Variola virus gp170, and Variola virus gp151 proteins all exhibit channel activity in the bacterial assays, and as such may be considered viroporin candidates. It is clear that more studies, such as patch clamping, will be needed to characterize the ionic conductivities of these proteins. However, our approach presents a rapid procedure to analyze open reading frames in other viruses, yielding new viroporin candidates for future detailed investigation. Finally, if conductivity is proven vital to their cognate viruses, the bio-assays presented herein afford a simple approach to screen for new channel blockers.

Random Mutagenesis Analysis of the Influenza A M2 Proton Channel Reveals Novel Resistance Mutants

The influenza M2 proton channel is a major drug target, but unfortunately, the acquisition of resistance mutations greatly reduces the functional life span of a drug in influenza treatment. New M2 inhibitors that inhibit mutant M2 channels otherwise resistant to the early adamantine-based drugs have been reported, but it remains unclear whether and how easy resistance could arise to such inhibitors. We have combined a newly developed proton conduction assay with an established method for selection and screening, both Escherichia coli-based, to enable the study of M2 function and inhibition. Combining this platform with two groups of structurally different M2 inhibitors allowed us to isolate drug resistant M2 channels from a mutant library. Two groups of M2 variants emerged from this analysis. A first group appeared almost unaffected by the inhibitor, M_089 (N13I, I35L, and F47L) and M_272 (G16C and D44H), and the single-substitution variants derived from these (I35L, L43P, D44H, and L46P). Functionally, these resemble the known drug resistant M2 channels V27A, S31N, and swine flu. In addition, a second group of tested M2 variants were all still inhibited by drugs but to a lesser extent than wild type M2. Molecular dynamics simulations aided in distinguishing the two groups where drug binding to the wild type and the less resistant M2 group showed a stable positioning of the ligand in the canonical binding pose, as opposed to the drug resistant group in which the ligand rapidly dissociated from the complex during the simulations.

A Robust Proton Flux (pHlux) Assay for Studying the Function and Inhibition of the Influenza A M2 Proton Channel

The M2 protein is an important target for drugs in the fight against the influenza virus. Because of the emergence of resistance against antivirals directed toward the M2 proton channel, the search for new drugs against resistant M2 variants is of high importance. Robust and sensitive assays for testing potential drug compounds on different M2 variants are valuable tools in this search for new inhibitors. In this work, we describe a fluorescence sensor-based assay, which we termed "pHlux", that measures proton conduction through M2 when synthesized from an expression vector in Escherichia coli. The assay was compared to a previously established bacterial potassium ion transport complementation assay, and the results were compared to simulations obtained from analysis of a computational model of M2 and its interaction with inhibitor molecules. The inhibition of M2 was measured for five different inhibitors, including Rimantadine, Amantadine, and spiro type compounds, and the drug resistance of the M2 mutant variants (swine flu, V27A, and S31N) was confirmed. We demonstrate that the pHlux assay is robust and highly sensitive and shows potential for high-throughput screening.

Inhibitors of the M2 Proton Channel Engage and Disrupt Transmembrane Networks of Hydrogen-Bonded Waters

Water-mediated interactions play key roles in drug binding. In protein sites with sparse polar functionality, a small-molecule approach is often viewed as insufficient to achieve high affinity and specificity. Here we show that small molecules can enable potent inhibition by targeting key waters. The M2 proton channel of influenza A is the target of the antiviral drugs amantadine and rimantadine. Structural studies of drug binding to the channel using X-ray crystallography have been limited because of the challenging nature of the target, with the one previously solved crystal structure limited to 3.5 Å resolution. Here we describe crystal structures of amantadine bound to M2 in the Inward_closed conformation (2.00 Å), rimantadine bound to M2 in both the Inward_closed (2.00 Å) and Inward_open (2.25 Å) conformations, and a spiro-adamantyl amine inhibitor bound to M2 in the Inward_closed conformation (2.63 Å). These X-ray crystal structures of the M2 proton channel with bound inhibitors reveal that ammonium groups bind to water-lined sites that are hypothesized to stabilize transient hydronium ions formed in the proton-conduction mechanism. Furthermore, the ammonium and adamantyl groups of the adamantyl-amine class of drugs are free to rotate in the channel, minimizing the entropic cost of binding. These drug-bound complexes provide the first high-resolution structures of drugs that interact with and disrupt networks of hydrogen-bonded waters that are widely utilized throughout nature to facilitate proton diffusion within proteins.

Mapping the Resistance Potential of Influenza's H+ Channel against an Antiviral Blocker

The development of drug resistance has long plagued our efforts to curtail viral infections in general and influenza in particular. The problem is particularly challenging since the exact mode of resistance may be difficult to predict, without waiting for untreatable strains to evolve. Herein, a different approach is taken. Using a novel genetic screen, we map the resistance options of influenza's M2 channel against its aminoadamantane antiviral inhibitors. In the process, we could identify clinically known resistant mutations in a completely unbiased manner. Additionally, novel mutations were obtained, which, while known to exist in circulating viruses, were not previously classified as drug resistant. Finally, we demonstrated the approach against an anti-influenza drug that has not seen clinical use, identifying several resistance mutations in the process. In conclusion, we present and employ a method to predict the resistance options of influenza's M2 channel to antiviral agents ahead of clinical use and without medical hazard.

Recent progress in designing inhibitors that target the drug-resistant M2 proton channels from the influenza A viruses

Influenza viruses are the causative agents for seasonal influenza, which results in thousands of deaths and millions of hospitalizations each year. Moreover, sporadic transmission of avian or swan influenza viruses to humans often leads to an influenza pandemic, as there is no preimmunity in the human body to fight against such novel strains. The metastable genome of the influenza viruses, coupled with the reassortment of different strains from a wide range of host origins, leads to the continuous evolution of the influenza virus diversity. Such characteristics of influenza viruses present a grand challenge in devising therapeutic strategies to combat influenza virus infection. This review summarizes recent progress in designing small molecule inhibitors that target the drug-resistant influenza A virus M2 proton channels and highlights the contribution of mechanistic studies of proton conductance to drug discovery. The lessons learned throughout the course of M2 drug discovery might provide insights for designing inhibitors that target other therapeutically important ion channels.

Interpreting Thermodynamic Profiles of Aminoadamantane Compounds Inhibiting the M2 Proton Channel of Influenza A by Free Energy Calculations

The development of novel anti-influenza drugs is of great importance because of the capability of influenza viruses to occasionally cross interspecies barriers and to rapidly mutate. One class of anti-influenza agents, aminoadamantanes, including the drugs amantadine and rimantadine now widely abandoned due to virus resistance, bind to and block the pore of the transmembrane domain of the M2 proton channel (M2TM) of influenza A. Here, we present one of the still rare studies that interprets thermodynamic profiles from isothermal titration calorimetry (ITC) experiments in terms of individual energy contributions to binding, calculated by the computationally inexpensive implicit solvent/implicit membrane molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach, for aminoadamantane compounds binding to M2TM of the avian "Weybridge" strain. For all eight pairs of aminoadamantane compounds considered, the trend of the predicted relative binding free energies and their individual components, effective binding energies and changes in the configurational entropy, agrees with experimental measures (ΔΔG, ΔΔH, TΔΔS) in 88, 88, and 50% of the cases. In addition, information yielded by the MM-PBSA approach about determinants of binding goes beyond that available in component quantities (ΔH, ΔS) from ITC measurements. We demonstrate how one can make use of such information to link thermodynamic profiles from ITC with structural causes on the ligand side and, ultimately, to guide decision making in lead optimization in a prospective manner, which results in an aminoadamantane derivative with improved binding affinity against M2TM(Weybridge).

Resistance-Mutation (N31) Effects on Drug Orientation and Channel Hydration in Amantadine-Bound Influenza A M2

The mechanism of amantadine binding to the S31 variant of the M2 protein of Influenza A is well understood, but the reasons behind N31 M2 amantadine insensitivity remain under investigation. Many molecular dynamics studies have evaluated the influence of amantadine position within the channel pore on its ability to inhibit proton conductance in M2, but little is known about the influence of amantadine rotational orientation. Replica-exchange umbrella sampling, steered, and classic molecular dynamics simulations were performed on amantadine in the solid-state NMR structure of S31 M2 and an N31 M2 homologue, both in the homotetramer configuration, to explore the effects of the position and tilt angle of amantadine on inhibition of the M2 channel. Steered simulations show that amantadine rotates with the amine toward the bulk water as it passes into the hydrophobic entryway lined by Val27 side chains. Results from all simulation types performed indicate that amantadine has a strong, specific orientation with the amine turned inward toward the central cavity in the S31 M2 pore but has variable orientation and a strong propensity to remain outward pointing in N31 M2. Free energy profiles from umbrella sampling, measured relative to bulk water, show amantadine binds more strongly to the S31 M2 pore by 8 kcal/mol in comparison to amantadine in the N31 pore, suggesting that it can escape more readily from the N31 channel through the Val27 secondary gate, whereas it is captured by the S31 channel in the same region. Lower water density and distribution near amantadine in S31 M2 reveal that the drug inhibits proton conductance in S31 M2 because of its inward-pointing configuration, whereas in N31 M2, amantadine forms hydrogen bonds with an N31 side chain and does not widely occlude water occupancy in any configuration. Both amantadine's weaker binding to and weaker water occlusion in N31 M2 might contribute to its inefficacy as an inhibitor of the mutant protein.

Molecular dynamics, monte carlo simulations, and langevin dynamics: a computational review

Macromolecular structures, such as neuraminidases, hemagglutinins, and monoclonal antibodies, are not rigid entities. Rather, they are characterised by their flexibility, which is the result of the interaction and collective motion of their constituent atoms. This conformational diversity has a significant impact on their physicochemical and biological properties. Among these are their structural stability, the transport of ions through the M2 channel, drug resistance, macromolecular docking, binding energy, and rational epitope design. To assess these properties and to calculate the associated thermodynamical observables, the conformational space must be efficiently sampled and the dynamic of the constituent atoms must be simulated. This paper presents algorithms and techniques that address the abovementioned issues. To this end, a computational review of molecular dynamics, Monte Carlo simulations, Langevin dynamics, and free energy calculation is presented. The exposition is made from first principles to promote a better understanding of the potentialities, limitations, applications, and interrelations of these computational methods.

Why bound amantadine fails to inhibit proton conductance according to simulations of the drug-resistant influenza A M2 (S31N)

The mechanisms responsible for drug resistance in the Asn31 variant of the M2 protein of influenza A are not well understood. Molecular dynamics simulations were performed on wild-type (Ser31) and S31N influenza A M2 in the homotetramer configuration. After evaluation of 13 published M2 structures, a solid-state NMR structure with amantadine bound was selected for simulations, an S31N mutant structure was developed and equilibrated, and the native and mutant structures were used to determine the binding behavior of amantadine and the dynamics of water in the two channels. Amantadine is stable in the plugging region of wild-type M2, with the adamantane in contact with the Val27 side chains, while amantadine in S31N M2 has more variable movement and orientation, and spontaneously moves lower into the central cavity of the channel. Free energy profiles from umbrella sampling support this observation. In this configuration, water surrounds the drug and can easily transport protons past it, so the drug binds without blocking proton transport in the S31N M2 channel.

Bacteria-based analysis of HIV-1 Vpu channel activity

HIV-1 Vpu is a small, single-span membrane protein with two attributed functions that increase the virus' pathogenicity: degradation of CD4 and inactivation of BST-2. Vpu has also been shown to posses ion channel activity, yet no correlation has been found between this attribute and Vpu's role in viral release. In order to gain further insight into the channel activity of Vpu we devised two bacteria-based assays that can examine this function in detail. In the first assay Vpu was over-expressed, such that it was deleterious to bacterial growth due to membrane permeabilization. In the second and more sensitive assay, the channel was expressed at low levels in K⁺ transport deficient bacteria. Consequently, Vpu expression enabled the bacteria to grow at otherwise non permissive low K⁺ concentrations. Hence, Vpu had the opposite impact on bacterial growth in the two assays: detrimental in the former and beneficial in the latter. Furthermore, we show that channel blockers also behave reciprocally in the two assays, promoting growth in the first assay and hindering it in the second assay. Taken together, we investigated Vpu's channel activity in a rapid and quantitative approach that is amenable to high-throughput screening, in search of novel blockers.