Affinity of Rimantadine Enantiomers against Influenza A/M2 Protein Revisited

Nascent pharmacological advancement in adamantane derivatives

The adamantane moiety has attracted significant attention since its discovery in 1933 due to its remarkable structural, chemical, and medicinal properties. This molecule has a notable impact in the therapeutic field because of its "add-on" lipophilicity to any pharmacophoric moieties. As in the case of molecular hybridization, in which one pharmacophore is attached to another one(s) with a probability of increasing the biological activity, adding an adamantane unit improves the absorption distribution, metabolism and excretion properties of the resultant hybrid molecule. This review summarizes various reports highlighting the biological activities of adamantane-based synthetic compounds and their structure-activity relationship study. The information presented in this review may open up possible dimensions for adamantane-based drug development and discovery in the pharmaceutical industry after proper structural modifications.

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

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

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.

Dynamic Interactions Between Brilliant Green and MscL Investigated by Solid-State NMR Spectroscopy and Molecular Dynamics Simulations

The mechanosensitive ion channel of large conductance (MscL) is a promising template for the development of new antibiotics due to its high conservation and uniqueness to microbes. Brilliant green (BG), a triarylmethane dye, has been identified as a new antibiotic targeted MscL. However, the detailed binding sites to MscL and the dynamic pathway of BG through the MscL channel remain unknown. Here, the dynamic interactions between BG and MscL were investigated using solid-state NMR spectroscopy and molecule dynamics (MD) simulations. Residue site-specific binding sites of BG to the MscL channel were identified by solid-state NMR. In addition, MD simulations revealed that BG conducts through the MscL channel via residues along the inner surface of the pore sequentially, in which the strong hydrophobic interactions between BG and hydrophobic residues F23 and I27 in the hydrophobic gate region of the MscL channel are major restrictions. Particularly, it was demonstrated that BG activates the MscL channel by reducing the hydrophobicity of the F23 in the gate region by water molecules that are bound to BG. Taken together, these simulations and experimental data provide novel insights into the dynamic interactions between BG and MscL, based on which new hydrophobic antibiotics and adjuvants targeting MscL can be developed.

Bayesian regression and model selection for isothermal titration calorimetry with enantiomeric mixtures

Bayesian regression is performed to infer parameters of thermodynamic binding models from isothermal titration calorimetry measurements in which the titrant is an enantiomeric mixture. For some measurements the posterior density is multimodal, indicating that additional data with a different protocol are required to uniquely determine the parameters. Models of increasing complexity-two-component binding, racemic mixture, and enantiomeric mixture-are compared using model selection criteria. To precisely estimate one of these criteria, the Bayes factor, a variation of bridge sampling is developed.

In Vitro Inhibition of Zika Virus Replication with Amantadine and Rimantadine Hydrochlorides

Zika virus (ZIKV) is a mosquito-borne flavivirus in which human infection became relevant during recent outbreaks in Latin America due to its unrecognized association with fetal neurological disorders. Currently, there are no approved effective antivirals or vaccines for the treatment or prevention of ZIKV infections. Amantadine and rimantadine are approved antivirals used against susceptible influenza A virus infections that have been shown to have antiviral activity against other viruses, such as dengue virus (DENV). Here, we report the in vitro effectiveness of both amantadine and rimantadine hydrochlorides against ZIKV replication, resulting in a dose-dependent reduction in viral titers of a ZIKV clinical isolate and two different ZIKV reference strains. Additionally, we demonstrate similar in vitro antiviral activity of these drugs against DENV-1 and yellow fever virus (YFV), although at higher drug concentrations for the latter. ZIKV replication was inhibited at drug concentrations well below cytotoxic levels of both compounds, as denoted by the high selectivity indexes obtained with the tested strains. Further work is absolutely needed to determine the potential clinical use of these antivirals against ZIKV infections, but our results suggest the existence of a highly conserved mechanism across flavivirus, susceptible to be blocked by modified more specific adamantane compounds.

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.

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.

Put a cork in it: Plugging the M2 viral ion channel to sink influenza

The ongoing threat of seasonal and pandemic influenza to human health requires antivirals that can effectively supplement existing vaccination strategies. The M2 protein of influenza A virus (IAV) is a proton-gated, proton-selective ion channel that is required for virus replication and is an established antiviral target. While licensed adamantane-based M2 antivirals have been historically used, M2 mutations that confer major adamantane resistance are now so prevalent in circulating virus strains that these drugs are no longer recommended. Here we review the current understanding of IAV M2 structure and function, mechanisms of inhibition, the rise of drug resistance mutations, and ongoing efforts to develop new antivirals that target resistant forms of M2.

Metabolic host response and therapeutic approaches to influenza infection

Based on available metabolomic studies, influenza infection affects a variety of cellular metabolic pathways to ensure an optimal environment for its replication and production of viral particles. Following infection, glucose uptake and aerobic glycolysis increase in infected cells continually, which results in higher glucose consumption. The pentose phosphate shunt, as another glucose-consuming pathway, is enhanced by influenza infection to help produce more nucleotides, especially ATP. Regarding lipid species, following infection, levels of triglycerides, phospholipids, and several lipid derivatives undergo perturbations, some of which are associated with inflammatory responses. Also, mitochondrial fatty acid β-oxidation decreases significantly simultaneously with an increase in biosynthesis of fatty acids and membrane lipids. Moreover, essential amino acids are demonstrated to decline in infected tissues due to the production of large amounts of viral and cellular proteins. Immune responses against influenza infection, on the other hand, could significantly affect metabolic pathways. Mainly, interferon (IFN) production following viral infection affects cell function via alteration in amino acid synthesis, membrane composition, and lipid metabolism. Understanding metabolic alterations required for influenza virus replication has revealed novel therapeutic methods based on targeted inhibition of these cellular metabolic pathways.

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.

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.

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

Recently, the binding kinetics of a ligand-target interaction, such as the residence time of a small molecule on its protein target, are seen as increasingly important for drug efficacy. Here, we investigate these concepts to explain binding and proton blockage of rimantadine variants bearing progressively larger alkyl groups to influenza A virus M2 wild type (WT) and M2 S31N protein proton channel. We showed that resistance of M2 S31N to rimantadine analogues compared to M2 WT resulted from their higher k_off rates compared to the k_on rates according to electrophysiology (EP) measurements. This is due to the fact that, in M2 S31N, the loss of the V27 pocket for the adamantyl cage resulted in low residence time inside the M2 pore. Both rimantadine enantiomers have similar channel blockage and binding k_on and k_off against M2 WT. To compare the potency between the rimantadine variants against M2, we applied approaches using different mimicry of M2, i.e., isothermal titration calorimetry and molecular dynamics simulation, EP, and antiviral assays. It was also shown that a small change in an amino acid at site 28 of M2 WT, which does not line the pore, seriously affects M2 WT blockage kinetics.