Unraveling the Binding, Proton Blockage, and Inhibition of Influenza M2 WT and S31N by Rimantadine Variants

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).

Antiviral Compounds to Address Influenza Pandemics: An Update from 2016-2022

In recent decades, the world has gained experience of the dangerous effects of pandemic events caused by emerging respiratory viruses. In particular, annual epidemics of influenza are responsible for severe illness and deaths. Even if conventional influenza vaccines represent the most effective tool for preventing virus infections, they are not completely effective in patients with severe chronic disease and immunocompromised and new small molecules have emerged to prevent and control the influenza viruses. Thus, the attention of chemists is continuously focused on the synthesis of new antiviral drugs able to interact with the different molecular targets involved in the virus replication cycle. To date, different classes of influenza viruses inhibitors able to target neuraminidase enzyme, hemagglutinin protein, Matrix-2 (M2) protein ion channel, nucleoprotein or RNAdependent RNA polymerase have been synthesized using several synthetic strategies comprising the chemical modification of currently used drugs. The best results, in terms of inhibitory activity, are in the nanomolar range and have been obtained from the chemical modification of clinically used drugs such as Peramivir, Zanamivir, Oseltamir, Rimantadine, as well as sialylated molecules, and hydroxypyridinone derivatives. The aim of this review is to report, covering the period 2016-2022, the most recent routes related to the synthesis of effective influenza virus inhibitors.

A Study of the Activity of Adamantyl Amines against Mutant Influenza A M2 Channels Identified a Polycyclic Cage Amine Triple Blocker, Explored by Molecular Dynamics Simulations and Solid-State NMR

We compared the anti-influenza potencies of 57 adamantyl amines and analogs against influenza A virus with serine-31 M2 proton channel, usually termed as WT M2 channel, which is amantadine sensitive. We also tested a subset of these compounds against viruses with the amantadine-resistant L26F, V27A, A30T, G34E M2 mutant channels. Four compounds inhibited WT M2 virus in vitro with mid-nanomolar potency, with 27 compounds showing sub-micromolar to low micromolar potency. Several compounds inhibited L26F M2 virus in vitro with sub-micromolar to low micromolar potency, but only three compounds blocked L26F M2-mediated proton current as determined by electrophysiology (EP). One compound was found to be a triple blocker of WT, L26F, V27A M2 channels by EP assays, but did not inhibit V27A M2 virus in vitro, and one compound inhibited WT, L26F, V27A M2 in vitro without blocking V27A M2 channel. One compound blocked only L26F M2 channel by EP, but did not inhibit virus replication. The triple blocker compound is as long as rimantadine, but could bind and block V27A M2 channel due to its larger girth as revealed by molecular dynamics simulations, while MAS NMR informed on the interaction of the compound with M2(18-60) WT or L26F or V27A.

Kolocouris A. Synthesis and Testing of Analogs of the Tuberculosis Drug Candidate SQ109 against Bacteria and Protozoa: Identification of Lead Compounds against Mycobacterium abscessus and Malaria Parasites

SQ109 is a tuberculosis drug candidate that has high potency against Mycobacterium tuberculosis and is thought to function at least in part by blocking cell wall biosynthesis by inhibiting the MmpL3 transporter. It also has activity against bacteria and protozoan parasites that lack MmpL3, where it can act as an uncoupler, targeting lipid membranes and Ca2+ homeostasis. Here, we synthesized 18 analogs of SQ109 and tested them against M. smegmatis, M. tuberculosis, M. abscessus, Bacillus subtilis, and Escherichia coli, as well as against the protozoan parasites Trypanosoma brucei, T. cruzi, Leishmania donovani, L. mexicana, and Plasmodium falciparum. Activity against the mycobacteria was generally less than with SQ109 and was reduced by increasing the size of the alkyl adduct, but two analogs were ∼4-8-fold more active than SQ109 against M. abscessus, including a highly drug-resistant strain harboring an A309P mutation in MmpL3. There was also better activity than found with SQ109 with other bacteria and protozoa. Of particular interest, we found that the adamantyl C-2 ethyl, butyl, phenyl, and benzyl analogs had 4-10× increased activity against P. falciparum asexual blood stages, together with low toxicity to a human HepG2 cell line, making them of interest as new antimalarial drug leads. We also used surface plasmon resonance to investigate the binding of inhibitors to MmpL3 and differential scanning calorimetry to investigate binding to lipid membranes. There was no correlation between MmpL3 binding and M. tuberculosis or M. smegmatis cell activity, suggesting that MmpL3 is not a major target in mycobacteria. However, some of the more active species decreased lipid phase transition temperatures, indicating increased accumulation in membranes, which is expected to lead to enhanced uncoupler activity.

From Acid Activation Mechanisms of Proton Conduction to Design of Inhibitors of the M2 Proton Channel of Influenza A Virus

With an estimated 1 billion people affected across the globe, influenza is one of the most serious health concerns worldwide. Therapeutic treatments have encompassed a number of key functional viral proteins, mainly focused on the M2 proton channel and neuraminidase. This review highlights the efforts spent in targeting the M2 proton channel, which mediates the proton transport toward the interior of the viral particle as a preliminary step leading to the release of the fusion peptide in hemagglutinin and the fusion of the viral and endosomal membranes. Besides the structural and mechanistic aspects of the M2 proton channel, attention is paid to the challenges posed by the development of efficient small molecule inhibitors and the evolution toward novel ligands and scaffolds motivated by the emergence of resistant strains.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (N_δ or Ν_ε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having N_ε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 μs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

Rimantadine Binds to and Inhibits the Influenza A M2 Proton Channel without Enantiomeric Specificity

The influenza A M2 wild-type (WT) proton channel is the target of the anti-influenza drug rimantadine. Rimantadine has two enantiomers, though most investigations into drug binding and inhibition have used a racemic mixture. Solid-state NMR experiments using the full length-M2 WT have shown significant spectral differences that were interpreted to indicate tighter binding for (R)- vs (S)-rimantadine. However, it was unclear if this correlates with a functional difference in drug binding and inhibition. Using X-ray crystallography, we have determined that both (R)- and (S)-rimantadine bind to the M2 WT pore with slight differences in the hydration of each enantiomer. However, this does not result in a difference in potency or binding kinetics, as shown by similar values for kon, koff, and Kd in electrophysiological assays and for EC50 values in cellular assays. We concluded that the slight differences in hydration for the (R)- and (S)-rimantadine enantiomers are not relevant to drug binding or channel inhibition. To further explore the effect of the hydration of the M2 pore on binding affinity, the water structure was evaluated by grand canonical ensemble molecular dynamics simulations as a function of the chemical potential of the water. Initially, the two layers of ordered water molecules between the bound drug and the channel's gating His37 residues mask the drug's chirality. As the chemical potential becomes more unfavorable, the drug translocates down to the lower water layer, and the interaction becomes more sensitive to chirality. These studies suggest the feasibility of displacing the upper water layer and specifically recognizing the lower water layers in novel drugs.

Influenza A M2 Inhibitor Binding Understood through Mechanisms of Excess Proton Stabilization and Channel Dynamics

Prevalent resistance to inhibitors that target the influenza A M2 proton channel has necessitated a continued drug design effort, supported by a sustained study of the mechanism of channel function and inhibition. Recent high-resolution X-ray crystal structures present the first opportunity to see how the adamantyl amine class of inhibitors bind to M2 and disrupt and interact with the channel’s water network, providing insight into the critical properties that enable their effective inhibition in wild-type M2. In this work, we examine the hypothesis that these drugs act primarily as mechanism-based inhibitors by comparing hydrated excess proton stabilization during proton transport in M2 with the interactions revealed in the crystal structures, using the Multiscale Reactive Molecular Dynamics (MS-RMD) methodology. MS-RMD, unlike classical molecular dynamics, models the hydrated proton (hydronium-like cation) as a dynamic excess charge defect and allows bonds to break and form, capturing the intricate interactions between the hydrated excess proton, protein atoms, and water. Through this, we show that the ammonium group of the inhibitors is effectively positioned to take advantage of the channel’s natural ability to stabilize an excess protonic charge and act as a hydronium mimic. Additionally, we show that the channel is especially stable in the drug binding region, highlighting the importance of this property for binding the adamantane group. Finally, we characterize an additional hinge point near Val27, which dynamically responds to charge and inhibitor binding. Altogether, this work further illuminates a dynamic understanding of the mechanism of drug inhibition in M2, grounded in the fundamental properties that enable the channel to transport and stabilize excess protons, with critical implications for future drug design efforts.

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.

Chemical Probes for Blocking of Influenza A M2 Wild-type and S31N Channels

We report on using the synthetic aminoadamantane-CH₂-aryl derivatives 1-6 as sensitive probes for blocking M2 S31N and influenza A virus (IAV) M2 wild-type (WT) channels as well as virus replication in cell culture. The binding kinetics measured using electrophysiology (EP) for M2 S31N channel are very dependent on the length between the adamantane moiety and the first ring of the aryl headgroup realized in 2 and 3 and the girth and length of the adamantane adduct realized in 4 and 5. Study of 1-6 shows that, according to molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations, all bind in the M2 S31N channel with the adamantyl group positioned between V27 and G34 and the aryl group projecting out of the channel with the phenyl (or isoxazole in 6) embedded in the V27 cluster. In this outward binding configuration, an elongation of the ligand by only one methylene in rimantadine 2 or using diamantane or triamantane instead of adamantane in 4 and 5, respectively, causes incomplete entry and facilitates exit, abolishing effective block compared to the amantadine derivatives 1 and 6. In the active M2 S31N blockers 1 and 6, the phenyl and isoxazolyl head groups achieve a deeper binding position and high k_on/low k_off and high k_on/high k_off rate constants, compared to inactive 2-5, which have much lower k_on and higher k_off. Compounds 1-5 block the M2 WT channel by binding in the longer area from V27-H37, in the inward orientation, with high k_on and low k_off rate constants. Infection of cell cultures by influenza virus containing M2 WT or M2 S31N is inhibited by 1-5 or 1-4 and 6, respectively. While 1 and 6 block infection through the M2 block mechanism in the S31N variant, 2-4 may block M2 S31N virus replication in cell culture through the lysosomotropic effect, just as chloroquine is thought to inhibit SARS-CoV-2 infection.

Site-directed M2 proton channel inhibitors enable synergistic combination therapy for rimantadine-resistant pandemic influenza

Pandemic influenza A virus (IAV) remains a significant threat to global health. Preparedness relies primarily upon a single class of neuraminidase (NA) targeted antivirals, against which resistance is steadily growing. The M2 proton channel is an alternative clinically proven antiviral target, yet a near-ubiquitous S31N polymorphism in M2 evokes resistance to licensed adamantane drugs. Hence, inhibitors capable of targeting N31 containing M2 (M2-N31) are highly desirable. Rational in silico design and in vitro screens delineated compounds favouring either lumenal or peripheral M2 binding, yielding effective M2-N31 inhibitors in both cases. Hits included adamantanes as well as novel compounds, with some showing low micromolar potency versus pandemic “swine” H1N1 influenza (Eng195) in culture. Interestingly, a published adamantane-based M2-N31 inhibitor rapidly selected a resistant V27A polymorphism (M2-A27/N31), whereas this was not the case for non-adamantane compounds. Nevertheless, combinations of adamantanes and novel compounds achieved synergistic antiviral effects, and the latter synergised with the neuraminidase inhibitor (NAi), Zanamivir. Thus, site-directed drug combinations show potential to rejuvenate M2 as an antiviral target whilst reducing the risk of drug resistance.

"Swine flu" illustrated that the spread of influenza pandemics in the modern era is rapid, making antiviral drugs the best way of limiting disease. One proven influenza drug target is the M2 proton channel, which plays an essential role during virus entry. However, resistance against licensed drugs targeting this protein is now ubiquitous, largely due to an S31N change in the M2 sequence. Understandably, considerable effort has focused on developing M2-N31 inhibitors, yet this has been hampered by controversy surrounding two potential drug binding sites. Here, we show that both sites can in fact be targeted by new M2-N31 inhibitors, generating synergistic antiviral effects. Developing such drug combinations should improve patient outcomes and minimise the emergence of future drug resistance.

N-Isopropylsulfinylimines vs. N-tert-butylsulfinylimines in the stereoselective synthesis of sterically hindered amines: an improved synthesis of enantiopure (R)- and (S)-rimantadine and the trifluoromethylated analogues

An improved fully stereoselective synthesis of both enantiomers of rimantadine and its trifluoromethylated analogues has been developed, using N-isopropylsulfinylimines as a starting chiral material, proving the superiority of the isopropyl group as a chiral inducer over the tert-butyl group in the case of hindered N-sulfinylimines.