Mechanism and Kinetics of Copper Complexes Binding to the Influenza A M2 S31N and S31N/G34E Channels.

Copper(II) is known to bind in the influenza virus His37 cluster in the homotetrameric M2 proton channel and block the proton current needed for uncoating. Copper complexes based on iminodiacetate also block the M2 proton channel and show reduced cytotoxicity and zebrafish-embryo toxicity. In voltage-clamp oocyte studies using the ubiquitous amantadine-insensitive M2 S31N variant, the current block showed fast and slow phases, in contrast to the single phase found for amantadine block of wild-type M2. Here, we evaluate the mechanism of block by copper adamantyl iminodiacitate and copper cyclooctyl iminodiacitate complexes and address whether the complexes can coordinate with one or more of the His37 imidazoles. The current traces were fitted to parametrized master equations. The energetics of binding and the rate constants suggest that the first step is copper complex binding within the channel, and the slow step in the current block is the formation of a Cu-histidine coordination complex. Solution-phase isothermal titration calorimetry and density functional theory (DFT) calculations indicate that imidazole binds to the copper complexes. Structural optimization using DFT reveals that the complexes fit inside the channel and project the Cu(II) toward the His37 cluster, allowing one imidazole to form a coordination complex with Cu(II). Electrophysiology and DFT studies also show that the complexes block the G34E amantadine-resistant mutant despite some crowding in the binding site by the glutamates.

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 (k1) is somewhat decreased for both AMT and RMT in the S31N, but the dissociation rate constant (k2) 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.

Investigation of a recent rise of dual amantadine-resistance mutations in the influenza A M2 sequence.

BACKGROUND: The S31N amantadine-resistance mutation in the influenza A M2 sequence currently occurs more frequently in nature than the S31 wild type. Overcoming the resistance of the S31N mutation is the primary focus of M2 researchers who aim to develop novel antiviral therapies. Recent studies have noted a possible rise in frequency of the V27A/S31N double amantadine-resistance mutation in recent years. The purpose of this study is to investigate this recent rise in frequency of the double mutation and any possible bias of the other mutations toward co-occurrence with S31N or S31 strains. RESULTS: The primary dataset used for this study was comprised of 24,152 influenza A M2 channel sequences which were downloaded from UniProt. There is an increased frequency for the S31N/V27A dual AR mutation in recent years, especially in swine. A test for difference in two proportions indicates that the V27A mutation is co-occurring with S31N more often than expected (p-value<0.001) when considering individual amino acid frequencies. At the same time, the different propensities for the V27A as compared to the V27T dual mutant may reflect differences in viral fitness or protein energetics, and this information could be exploited to focus drug development so as to reduce further drug insensitivity. CONCLUSIONS: The development of the S31N/V27A variant in the Midwestern US swine may be a harbinger of novel human strain development. V27A/S31N is a possible path forward for the evolution of M2 which may convey a new level of drug resistance and should receive attention in drug design.

Drug sensitivity, drug-resistant mutations, and structures of three conductance domains of viral porins.

Recent controversies associated with the structure of the M2 protein from influenza A virus and the binding site of drug molecules amantadine and rimantadine motivated the comparison here of the drug binding to three viral porins including the M2 proteins from influenza A and B as well as the viral protein 'u' from HIV-1. While the M2 protein from influenza B does not normally bind amantadine, chimeras with the M2 protein from influenza A show blockage by amantadine. Similarly, Vpu does not normally bind rimantadine, but the single site mutation A18H converts a non-specific channel to a selective proton channel that is sensitive to rimantadine. The comparison of structures and amino acid sequences shows that the membrane protein sample environment can have a significant influence on the structural result. While a bilayer surface bound amphipathic helix has been characterized for AM2, such a helix may be possible for BM2 although it has evaded structural characterization in detergent micelles. A similar amphipathic helix seems less likely for Vpu. Even though the A18H Vpu mutant forms rimantadine sensitive proton channels, the binding of drug and its influence on the protein structure appears to be very different from that for the M2 proteins. Indeed, drug binding and drug resistance in these viral porins appears to result from a complex set of factors.

Functional reconstitution of influenza A M2(22-62).

Amantadine-sensitive proton uptake by liposomes is currently the preferred method of demonstrating M2 functionality after reconstitution, to validate structural determination with techniques such as solid-state NMR. With strong driving forces (two decades each of both [K(+)] gradient-induced membrane potential and [H(+)] gradient), M2(22-62) showed a transport rate of 78 H(+)/tetramer-s (pH(o) 6.0, pH(i) 8.0, nominal V(m)=-114 mV), higher than previously measured for similar, shorter, and full-length constructs. Amantadine sensitivity of the conductance domain at pH 6.8 was also comparable to other published reports. Proton flux rate was optimal at protein densities of 0.05-1.0% (peptide wt.% in lipid). Rundown of total proton uptake after addition of valinomycin and CCCP, as detected by delayed addition of valinomycin, indicated M2-induced K(+) flux of 0.1K(+)/tetramer-s, and also demonstrated that the K(+) permeability, relative to H(+), was 2.8 × 10(-6). Transport rate, amantadine and cyclooctylamine sensitivity, acid activation, and H(+) selectivity were all consistent with full functionality of the reconstituted conductance domain. Decreased external pH increased proton uptake with an apparent pK(a) of 6.

Gramicidin channels are internally gated.

Gramicidin channels are archetypal molecular subjects for solid-state NMR studies and investigations of single-channel or cation conductance. Until now, the transitions between on and off conductance states have been thought, based on multichannel studies, to represent monomer <--> dimer reactions. Here we use a single-molecule deposition method (vesicle fusion to a planar bilayer) to show that gramicidin dimer channels do not normally dissociate when conductance terminates. Furthermore, the observation of two 13C peaks in solid-state NMR indicates very stable dichotomous conformations for both the first and second peptide bonds in the monomers, and a two-dimensional chemical exchange spectrum with a 12-s mixing time demonstrates that the Val1 carbonyl conformations exchange slowly, with lifetimes of several seconds. It is proposed that gramicidin channels are gated by small conformational changes in the channel near the permeation pathway. These studies demonstrate how regulation of conformations governing closed <--> open transitions may be achieved and studied at the molecular level.

Computational studies of gramicidin permeation: an entry way sulfonate enhances cation occupancy at entry sites.

The impact on the cation-transport free-energy profile of replacing the C-terminal ethanolamine in the gramicidin A channel with a taurine residue is studied using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl saline solution. The potential of mean force for ion transport is obtained by umbrella sampling. The presence of a negatively charged sulfonate group at the entrance of the gramicidin channel affects the depth and the location of the binding sites, producing a strong attraction for the cations in the bulk. The potential of mean force by the sulfonate acting directly through electrostatics and van der Waals interactions on the test ion is highly modulated by indirect effects (i.e., sulfonate effects on other components of the system that, in turn, affect the ion free-energy profile in the channel). Because the "entry" sites are located symmetrically at both entry and exit of the channel, the deeper free-energy wells should inhibit exit. Given that the channel has increased conductance experimentally, the simulation results suggest that the channel conductance is normally entry limited.

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

We report on using the synthetic aminoadamantane-CH2-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 kon/low koff and high kon/high koff rate constants, compared to inactive 2-5, which have much lower kon and higher koff. Compounds 1-5 block the M2 WT channel by binding in the longer area from V27-H37, in the inward orientation, with high kon and low koff 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.

The boundary lipid around DMPC-spanning influenza A M2 transmembrane domain channels: Its structure and potential for drug accommodation.

We have investigated the perturbation of influenza A M2TM in DMPC bilayers. We have shown that (a) DSC and SAXS detect changes in membrane organization caused by small changes (micromolar) in M2TM or aminoadamantane concentration and aminoadamantane structure, by comparison of amantadine and spiro[pyrrolidine-2,2'-adamantane] (AK13), (b) that WAXS and MD can suggest details of ligand topology. DSC and SAXS show that at a low M2TM micromolar concentration in DPMC bilayers, two lipid domains are observed, which likely correspond to M2TM boundary lipids and bulk-like lipids. At higher M2TM concentrations, one domain only is identified, which constitutes essentially all of the lipid molecules behaving as boundary lipids. According to SAXS, WAXS, and DSC in the absence of M2TM, both aminoadamantane drugs exert a similar perturbing effect on the bilayer at low concentrations. At the same concentrations of the drug when M2TM is present, amantadine and, to a lesser extent, AK13 cause, according to WAXS, a significant disordering of chain-stacking, which also leads to the formation of two lipid domains. This effect is likely due, according to MD simulations, to the preference of the more lipophilic AK13 to locate closer to the lateral surfaces of M2TM when compared to amantadine, which forms stronger ionic interactions with phosphate groups. The preference of AK13 to concentrate inside the lipid bilayer close to the exterior of the hydrophobic M2TM helices may contribute to its higher binding affinity compared to amantadine.

Free-energy profiles for ions in the influenza M2-TMD channel.

M(2) transmembrane domain channel (M(2)-TMD) permeation properties are studied using molecular dynamics simulations of M(2)-TMD (1NYJ) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. This study allows examination of spontaneous cation and anion entry into the selectivity filter. Three titration states of the M(2)-TMD tetramer are modeled for which the four His(37) residues, forming the selectivity filter, are net uncharged, +2 charged, or +3 charged. M(2)-TMD structural properties from our simulations are compared with the properties of other models extracted from NMR and X-ray studies. During 10 ns simulations, chloride ions occasionally occupy the positively-charged selectivity filter region, and from umbrella sampling simulations, Cl(-) has a lower free-energy barrier in the selectivity-filter region than either Na(+) or NH(4) (+), and NH(4) (+) has a lower free-energy barrier than Na(+). For Na(+) and Cl(-), the free-energy barriers are less than 5 kcal/mol, suggesting that the 1NYJ conformation would probably not be exquisitely proton selective. We also point out a rotameric configuration of Trp(41) that could fully occlude the channel.

Novel influenza inhibitors designed to target PB1 interactions with host importin RanBP5.

In search of novel targets for influenza inhibitors, a site on PB1 was selected for its high conservation and probable interaction with a host protein, RanBP5, that is key to nuclear import of PB1, where it complexes with PB2, PA, and NP to transcribe viral RNA. Docking with libraries of drug-like compounds led to a selection of five candidates that bound tightly and with a pose likely to inhibit protein binding. These were purchased and tested in vitro, found to be active, and then one was synthetically expanded to explore the structure-activity relationship. The top candidates had a carboxylic acid converted to an ester and electron-withdrawing substituents added to a phenyl group in the original structure. Resistance was slow to develop, but cytotoxicity was moderately high. Nuclear localization of PB1 and in vitro polymerase activity were both strongly inhibited.

Proton transport through influenza A virus M2 protein reconstituted in vesicles.

Influenza A virus M2 protein is known to form acid-activated, proton-selective, amantadine-sensitive channels. We directly measured proton uptake in vesicles containing reconstituted M2 by monitoring external pH after addition of valinomycin to vesicles with 100-fold-diluted external [K(+)]. External pH typically increased by a few tenths of a pH unit over a few minutes after valinomycin addition, but proton uptake was not significantly altered by acidification. Under neutral conditions, external addition of 1 mM amantadine produced a reduction in flux consistent with randomly ordered channels; however, experimental variation is high with this method and the block was not statistically significant. Amantadine block was reduced at pH 5.4. In accord with Lin and Schroeder's study of reconstituted M2 using a pH-sensitive dye to monitor intravesicular pH, we conclude that bath pH weakly affects or does not significantly affect proton flow in the pH range 5.4-7.0 for the reconstituted system, contrary to results from electrophysiological studies. Theoretical analysis of the relaxation to Donnan equilibrium utilized for such vesicle uptake assays illuminates the appropriate timescale of the initial slope and an important limitation that must be placed on inferences about channel ion selectivity. The rise in pH over 10 s after ionophore addition yielded time-averaged single-channel conductances of 0.35 +/- 0.20 aS and 0.72 +/- 0.42 aS at pH 5.4 and 7.0, respectively, an order of magnitude lower than previously reported in vesicles. Assuming complete membrane incorporation and tetramerization of the reconstituted protein, such a low time-averaged conductance in the face of previously observed single-channel conductance (6 pS at pH 3) implies an open channel probability of 10(-6)-10(-4). Based on leakage of potassium from M2-containing vesicles, compared to protein-free vesicles, we conclude that M2 exhibits approximately 10(7) selectivity for hydrogen over potassium.

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 koff rates compared to the kon 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 kon and koff 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.

Divalent copper complexes as influenza A M2 inhibitors.

New M2 blockers effective against the ubiquitous amantadine-resistant S31N M2 mutation in influenza A are needed. Six copper complexes, 2, 4, 6, 8, 9, and 10, were synthesized and found to block both wild type and S31N M2. Free Cu2+ also blocks M2 S31N but not S31N/H37A. The copper complexes do not block M2 H37A (either S31 or S31N). The complexes were effective against three influenza A strains in cell-culture assays, but less toxic to cells than CuCl2. For example 4, Cu(cyclooctylamineiminodiacetate), which was stable at pH > 4 in the buffers used, had an EC50 against A/Calif/07/2009 H1N1 of 0.7 ± 0.1 µM with a CC50 of 147 µM (therapeutic index, averaged over three strains, 67.8). In contrast, CuCl2 had an EC50 of 3.8 ± 0.9 µM and CC50 of 19 µM. Because M2 H37 is highly conserved, these complexes show promise for further testing as drugs against all strains of influenza A.

Applied field nonequilibrium molecular dynamics simulations of ion exit from a beta-barrel model of the L-type calcium channel.

We present results of applied field nonequilibrium molecular dynamics simulations (AF NEMD) of a minimal beta-barrel model channel intended to represent an L-type calcium channel that suggests a possible relationship between glutamate side chain conformational changes and ion flux in calcium channels. The beta-barrel is used to provide a scaffolding for glutamate side chains and a confinement for electrolyte of dimensions similar to the expected channel structure. It was preloaded with ions to explore relative rates of ion exit for different occupancy configurations. Our simulations with an asymmetrical flexible selectivity filter represented by four glutamate side chains (EEEE), one of which differs in initial dihedrals from the other three, indicate a plausible mechanism for the observed anomalous mole fraction effect seen in calcium channels. Apparent rates of electric field-induced exit from channels preloaded with three Na+ ions are much higher than for channels with one Ca2+ followed by two Na+ ions, consistent with the common notion that Ca2+ block of Na+ current is due to competition between the Ca2+ and Na+ ions for the negatively charged (EEEE) locus. In our model, the Ca2+ ion ligates simultaneously to the four negatively charged glutamate side chains and sterically blocks the permeation pathway. Ca2+-relief of Ca2+-block is suggested by a much higher rate of exit for channels preloaded with three Ca2+ ions than for channels with two Ca2+ ions.

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.

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.

Spectroscopic analysis of low pH and lipid-induced structural changes in type A botulinum neurotoxin relevant to membrane channel formation and translocation.

Botulinum neurotoxin (BoNT) is an extremely toxic protein to animals and humans. In its mode of action, one of its subunits mediates its translocation by integrating itself into the membrane bilayer. We have examined the membrane channel activity of type A BoNT (BoNT/A) and its heavy (H) chain in planar lipid membrane under various pH conditions to understand the possible role of the channel activity in the translocation of the BoNT/A light (L) chain under physiological conditions. Only BoNT/A H chain, and not the BoNT/A, exhibited membrane channel activity for translocation of ions. The H chain-induced increase in conductance did not require a pH gradient across the lipid membrane, although it was enhanced by a pH gradient. To understand the molecular basis of the membrane channel activity and the translocation of the L chain, the secondary structure of BoNT/A and its H and L chains were analyzed using circular dichroism (CD) and Fourier-transform infrared (FT-IR) spectroscopy at different pH values. BoNT/A showed no structural alternation upon acidifying the buffer pH. However, an increase in beta-sheet content of BoNT/A H chain at low pH was noted when examined by FT-IR. The L chain structure significantly changed with decrease in pH, and the change was mostly reversible. In addition, the neurotoxin and its subunit chains induced a partially reversible aggregation of liposomes at low pH, which indicated their integration into the lipid bilayer. Temperature-induced denaturation studies of BoNT/A H chain indicated major structural reorganization upon its interaction with membrane, especially at low pH.