Elnaz Aledavood
Beatrice Selmi
Carolina Estarellas
Matteo Masetti
F. Javier Luque- 1Departament de Nutrició, Ciències de l'Alimentació i Gastronomia, Institut de Biomedicina and Institut de Química Teòrica i Computacional, University of Barcelona, Barcelona, Spain
- 2Department of Pharmacy and Biotechnology, Alma Mater Studiorum – Università di Bologna, Bologna, Italy
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.
Introduction
Influenza A virus is an important pathogen that still causes the death of 290,000–650,000 people in seasonal outbreaks worldwide (Iuliano et al., 2018). It pertains to the Orthomyxoviridae family, which is characterized by the presence of a negative-sense, single-stranded, enveloped ribonucleic acid with a segmented genome (Cheung and Poon, 2007). There are four genera in the Orthomyxoviridae family: Influenza A, Influenza B, Influenza C, and Thogotovirus. Influenza A viruses are found in many different animals, and Influenza B viruses circulate widely only among humans. Influenza A viruses are divided into subtypes based on two proteins on the surface of the virus: hemagglutinin (HA), which is responsible for the receptor binding and membrane fusion, and neuraminidase (NA), which assists the release of the viral progeny (Cheung and Poon, 2007). Up to 18 variants of HA (from type H1 to H18) have been described till now, whereas 11 subtypes have been disclosed for NA. All known subtypes of influenza A viruses have been found among birds, except subtypes H17N10 and H18N11, which have only been found in bats. Only some subtypes (H1N1, H2N2, H3N2, H5N1, H7N7, and H9N2) have been isolated from humans, suggesting that there are restrictions to host viruses in humans.
While infections in humans directly from animals are unusual, sporadic infections and outbreaks caused by influenza A viruses have occurred in the past. The most severe pandemic was the “Spanish Flu” of 1918–1919, which caused 20–50 million deaths by influenza A H1N1 virus strain. The subsequent health challenges were the “Asian flu” in 1957 and the “Hong Kong flu” in 1968, which were caused respectively by H2N2 and H3N2 strains, and more recently the 2009 swine influenza pandemics, which pertains to the H1N1 strain (Cox and Subbarao, 2000; Hayward et al., 2014; Watson et al., 2015; Saunders-Hastings and Krewski, 2016). The latent risk of facing future pandemics underlined the Global Influenza Strategy initiative launched by the World Health Organization in 2019, which has the goal of preventing seasonal influenza, controlling the spread of influenza from animals to humans, and preparing for the next influenza pandemic (World Health Organization, 2019). This threatening scenario underscores the urgency in developing better diagnostic tools and promoting the discovery of effective influenza treatments, including vaccines and antiviral drugs.
This review aims to discuss experimental and computational efforts invested in elucidating the structural basis of the proton conducting activity of the M2 protein channel and the discovery of antiviral compounds targeting the M2 transmembrane domain. To this end, the manuscript first describes the role of the M2 proton channel in the replicative life cycle of the influenza A virus. Then, attention is focused on the structure of the proton channel, which is directly implicated in the proton transfer across the viral membrane as a preliminary event leading to the fusion of viral and host membranes and the subsequent release of the viral genetic material. Emphasis is made on the molecular mechanisms that underlie the proton conduction through the lumen of the M2 proton channel. The discussion also highlights the role played by key residues located in the transmembrane region, and especially the efforts spent in the determination of the pKa of the imidazole rings present in the His37 tetrad. Finally, the manuscript examines the mechanism of action of selected inhibitors, and the evolution followed in the design of novel compounds introduced to increase the inhibition potency and to mitigate the emergence of resistance to current treatments.
The Influenza A Virus Replicative Cycle
Influenza A virus replicates in the epithelial cells of the upper respiratory tract of humans, pigs, and horses. Regarding in vitro models, the virus can also affect different cells that contain targets decorated with sialic acid moieties (Julkunen et al., 2001; Sidorenko and Reichl, 2004), which act as the entrance anchor point into the cells.
Enclosed by a lipid-protein envelope, the genetic material of the influenza A virus contains 8 single-stranded RNA segments that encode the genetic information for viral proteins (Dou et al., 2018; Krammer et al., 2018; Jung and Lee, 2020): the membrane-bound hemagglutinin (HA) and neuraminidase (NA), the transmembrane M2 proton channel, the matrix protein M1, the non-structural proteins (NS1, NS2/NEP), and the heterotrimeric RNA-dependent RNA polymerase (RdRp), formed by three domains (PA, PB1, and PB2), and nucleoprotein (NP) components.
The replicative life cycle of the virus is schematically shown in Figure 1. The first step of the cycle is the entry of the virus inside the host cell (step 1). Viral entry into the cell occurs via HA, which recognizes and binds the sialic acid of the host cell receptors. The virions penetrate inside the host cell through clathrin-dependent receptor-mediated endocytosis (step 2). The acidification of the endosome is facilitated by the M2 proton channel, which transfers protons to the interior of the viral particles. This event promotes a large-scale conformational remodeling in HA that facilitates the release of the fusion peptide, and the concomitant structural rearrangement of HA triggers the fusion of the viral and host membranes, leading to conformational changes in the M1 protein that facilitate the release of the viral ribonucleoproteins (vRNPs) to the cellular cytoplasm (step 3) (Honda et al., 2002). These vRNPs comprise complexes formed by the individual viral RNA segments and proteins, including the RNA-dependent RNA polymerase (RdRp), which is associated with the 3′- and 5′-ends of viral RNA segments, and nucleoproteins (NPs) that are coating the remainder of the viral RNA. The vRNP complexes are then transported into the nucleus to initiate transcription and replication by RdRp, which is formed by the polymerase basic protein 1 (PB1), polymerase basic protein 2 (PB2), and polymerase acidic protein (PA). The first step is known as cap-snatching, a process where RdRp (PB2 and PA subunits) interacts with the C-terminal domain of the host RNA polymerase II and removes the first 10–13 nucleotides of a nascent transcript, which is then used as a primer to initiate transcription from the viral RNA template by PB1 (De Vlugt et al., 2018). This is followed by a viral genome replication phase that involves the synthesis of full-length viral RNA and cRNA strands (steps 4 and 5) (Pielak and Chou, 2011). With a late stage of packaging, vRNPs are formed (step 6). The synthesis of M2 protein, HA and NA is carried out by ribosomes bound to the membranes of the endoplasmic reticulum, which are then glycosylated and transported to the Golgi apparatus (step 7) (Sidorenko and Reichl, 2004). Finally, a complex between the membrane-embedded proteins and vRNPs is formed and delivered to the cell membrane, leading to the assembly of new virions that will be released to the extracellular side of the host cell membrane (step 8).
FIGURE 1. Schematic pattern of the different phases of the replicative cycle in influenza A virus.
The preceding discussion suffices to reveal the complexity encoded in the multi-step processes - cell entry, replication, intracellular trafficking, virion assembly, and release - involved in the generation of new virions. In this scenario, the search for effective antiviral treatments has been inspired by the identification of druggable viral targets with a relevant role in the life cycle of virus infection. Early drug discovery projects targeting the M2 proton channel primarily relied on adamantane-related compounds (Moorthy et al., 2014; Wang et al., 2015). However, the emergence of resistant strains motivated the search of novel therapeutic approaches, including the development of 1) NA inhibitors such as zanamivir and oseltamivir, although drug-resistant variants have also been reported, 2) HA inhibitors, as exemplified by the discovery of arbidol and JNJ7918, and, more recently, 3) RdRp inhibitors, such as pimodivir and baloxavir (Loregian et al., 2014; Wu et al., 2017; Han et al., 2018; Mifsud et al., 2019; Ginex and Luque, 2021). The discovery of small-molecule antiviral compounds against these targets has benefited from the progress made in disclosing the fine structural details of these proteins, the molecular mechanisms that underlie their biological function, and the characterization of the mode of action of inhibitory compounds.
Structure of the M2 Proton Channel
The M2 proton channel is a homotetrameric integral transmembrane (TM) protein located in the viral envelope of the influenza A virus. It is a selective, pH-dependent channel that regulates the acidification of the interior of the virion, leading to dissociation of the viral RNA from its bound matrix proteins and release of the viral genetic material for replication (Pinto et al., 1992; Pielak and Chou, 2011; Manzoor et al., 2017).
From a structural point of view, the M2 channel is a homotetramer, each monomer consisting of 97 amino acids and can be divided into four main regions: 1) the highly conserved, unstructured N-terminal domain (residues 1–21) located in the viral exterior, which assists the incorporation of the M2 channel into the virion, 2) the TM domain (residues 22–46) (Figure 2), which assembles into a tetrameric domain and is involved in the proton flux, 3) an amphiphilic membrane-anchored helical domain (residues 47–67), which is located in the interface with the membrane and stabilizes the channel by inducing membrane curvature and mediates membrane scission, and 4) the C-terminal domain (residues 68–97), which is located in the viral interior and interacts with the matrix protein M1 and contributes to virus packaging and budding (McCown and Pekosz, 2006; Rossman et al., 2010; Wang and Hong, 2015a).
FIGURE 2. (Left) Representation of selected residues that shape the inner pore of the M2 proton channel (shown as sticks). (Right) Representation of the four transmembrane helices (residues 22–46; shown as cartoon in magenta, blue, pink and violet), where the amphipathic helices (residues 47–67) that are located near the inner side of the membrane are shown as grey cartoon. Polar heads of membrane lipids in both outer and inner sides of the viral membrane are shown as green spheres (adapted from a Molecular Dynamics simulations of the membrane-embedded M2 proton channel).
The TM region contains the binding site for M2 inhibitors in the interior of the pore, filling a pocket located between the tetrads formed by residues Val27, Ser31, Gly34, and His37 (Schnell and Chou, 2008; Acharya et al., 2010; Cady et al., 2010; Sharma et al., 2010). This pore is implicated in the proton conduction, and the proton current is regulated by the tetrads formed by His37 and Trp41. His37 is a pH sensor and can conduct the proton inside the channel by changing the protonation state of the imidazole side chain, whereas Trp41 acts as a gate of the channel. Accordingly, binding of compounds that fill the binding site in the interior of the M2 channel would sterically impede the passage of protons (Figure 2).
The pH-sensing gating mechanism regulates the proton flux through the channel, which is modulated by the transition from an open state at an acidic pH, where the imidazole rings of the His37 tetrad are protonated, to a closed state at alkaline pH environment, which promotes the deprotonated state of His37 residues. In the open state, the M2 channel is also permeable to some cations, such as Na+ and K+, through an antiporter-like mechanism, carrying out the metal cation efflux in conjunction to proton influx, but the permeability to protons has been estimated to be 106–107 higher than the conduction of alkaline cations (Chizhmakov et al., 1996; Shimbo et al., 1996; Mould et al., 2000a; Manzoor et al., 2017). There is an asymmetric conductance as the proton flux is higher when the pH outside the virus is low, consistently with an acidic endosome. This behavior has been attributed to the presence of the Trp41 gate, which reduces the conduction in the reverse direction (Ma et al., 2013). Thus, when the His37 tetrad is protonated by protons coming from the endosome, the indole moiety of Trp41 undergoes a conformational change and the channel enters an open state that would allow the passage of protons to the interior of the virus.
The functional relevance of His37 and Trp41 is revealed by the drastic effects triggered by mutations in these positions, which cause a reduction in the pH sensitivity of the channel, and a disruption in the proton conduction (Wang et al., 1995; Tang et al., 2002; Venkataraman et al., 2005). Thus, mutations in His37 have a drastic effect on the functional properties of the M2 channel, since the mutated proteins not only increase the proton conduction of the channel, but also suppress the selectivity for proton conduction. As an example, the specific activity determined for the voltage-independent component of the conductance at pH 6.2 for the wild type channel was 0.16 μA/ng, but increased to 1.36 and 30 μA/ng for the H37G and H37E mutated channels, respectively (Holsinger et al., 1994). Furthermore, the mutated ion channels H37G, H37S, and H37T were found to be deficient in ion selectivity (Venkataraman et al., 2005). On the other hand, the expression of the M2 protein in oocytes leads to an inward H+ current upon bathing in acidic media (pH 5.9), which is lost upon return to alkaline conditions (pH 8.5). However, whereas no currents were observed for mutated channels containing Phe41, Cys41, or Ala41 prior to acidification, outward H+ currents were detected by changing acidic conditions to alkaline solution.
Mechanism of Proton Conduction
The specific role played by His37 and Trp41 in assisting the proton flux through the pore of the M2 channel has been the subject of intense efforts, mainly focused on the TM domain, which reproduces most of the biophysical, site-directed mutagenesis and electrophysiological features of the full-length protein. These studies, further assisted by X-ray crystallography of structures solved in membrane-mimetic solvents at different pH and temperature, solution NMR, solid-state NMR (ss-NMR) spectroscopy (see below), in conjunction with computational simulations, have crystallized in several models of proton conduction.
An early mechanistic hypothesis was the water wire model, also known as the shutter mechanism (Figure 3), which assumes the presence of a discontinuous water column when the His37 tetrad is in the deprotonated state, where the channel adopts a closed conformation that would be populated at alkaline pH (Sansom et al., 1997). When the pH decreases, protonation of the His37 tetrad leads to electrostatic repulsion between the charged imidazole rings. This would facilitate the conformational transition to an open state, which enables the diffusion of excess proton through the pore by means of a Grotthuss-type mechanism (Agmon, 1995). By merely working as a gatekeeper, in this model, the His37 tetrad would play a passive role as long as it is not directly involved in the proton transfer.
FIGURE 3. Schematic representation of the shutter mechanism. In the closed state the channel is in the neutral state (imidazole ring in yellow), and there is not an electrostatic repulsion between the imidazole rings of His37. However, when His37 are fully protonated (coloured in blue), the electrostatic repulsion between imidazole rings triggers a conformational change that allows water and protons to pass through the gate.
In contrast to the shutter model, the currently accepted proton conduction mechanism assumes a proton relay model (Figure 4). According to this model, protons diffuse along a water wire until reaching the tetrad formed by His37, generating an imidazolium intermediate, which in turn releases a proton to the inner side (Pinto et al., 1997; Okada et al., 2001). This mechanism was first proposed by Pinto et al. on the basis of conductance measurements in conjunction with Cys scanning mutagenesis (Pinto et al., 1997). However, a number of questions remained to be elucidated, such as the number of His residues required to be protonated for proton conduction, the stabilization of the protonated imidazolium via interactions with other residues in the channel or with water molecules filling the inner pore, and whether translocation of the proton into the viral interior would be mediated by prototropic tautomerism of the imidazole ring or through conformational flip of the His37 side chain.
FIGURE 4. Schematic representation of the proton transfer flow associated to conformational rearrangement of the His37 imidazole ring.
Deeper knowledge into the molecular events implicated in proton conduction has been gained from ss-NMR studies. Hu et al. refined the proton relay model suggesting that proton conduction is concomitant with His37 conformational changes (Figure 4) (Hu et al., 2010). At high external pH the His37 imidazole are tightly packed through CH-π interactions, preventing the formation of hydrogen-bonded water chains in the pore. Acidification protonates the imidazole ring at position Nδ1, and the electrostatic repulsion between protonated histidines would widen the pore (Khurana et al., 2009), enabling the reorientation of the imidazolium ring at the microsecond time scale. Proton conduction would involve imidazolium deprotonation facilitated by Cε1-Hε1 hydrogen bonding and continuous ring flips, which would contribute to the proper alignment of the charged imidazolium with C-terminal water molecules and rearrangement of the neutral imidazole to the N-terminal side of the pore to start another protonation cycle.
This proton transfer model dynamically couples the interaction of His37 with water molecules, which are required for delivering protons to the imidazole ring, with the conformational flips of His37 side chains. The barrier of the conformational rearrangement of the His37 residues was estimated to be close to 14 kcal mol−1, which would justify the temperature dependence of proton transport, though possibly this value may contain a contribution arising from the barrier for proton transfer (Lin and Schroeder, 2001). This model is consistent with the deuterium isotope effects (Mould et al., 2000b), and support comes from 15N NMR data, which yielded an estimated rate of 105 s−1 for the protonation and deprotonation of the imidazole nitrogens (Hu et al., 2012). This rate was consistent with the ring reorientation rate estimated from motionally averaged dipolar couplings, supporting the proposal that ring reorientation is synchronized with, and facilitates, proton transfer. Notably, only a single histidine residue is required to be actively involved in proton transfer according to this mechanism.
Other ssNMR studies suggested that proton conduction may also involve the formation and breaking of hydrogen bonds between adjacent pairs of histidines, in conjunction with the assistance of tryptophan residues, to guide the proton through the channel (Figure 5) (Sharma et al., 2010). Thus, two His37 residues are protonated, each proton being shared with an adjacent His residue through a strong hydrogen bond. Accordingly, the His37 tetrad can be viewed as a pair of imidazole-imidazolium dimer (also denoted dimer-of-dimers). This would define the so-called histidine-locked state, which is further stabilized by a cation-π interaction of the protonated histidine with a tryptophan residue. When the pH decreases, a proton is transferred from the hydronium cation to the interresidue hydrogen bond formed between Nδ1 and Nε2 of the hydrogen-bonded histidines, leading to the activated state that contains a third protonated His37. The two imidazolium rings can then rotate, so that a protonated histidine may form a hydrogen bond with water in the N-terminal pore, while a conformational change in the other histidine enables the formation of a cation-π interaction with an indole of Trp41, which would block water access from the C-terminal pore. The conducting state is obtained when this indole moves aside to expose the Nε2 proton to a water molecule in the C-terminal pore. Finally, the histidine-locked state is recovered after the release of the Nε2 proton to a water molecule in the C-terminus. Recent studies have also suggested that the hydrogen bonding between imidazole rings could be formed even at a pH of 7.8 in the neutral charge state (Movellan et al., 2020), and that a water molecule is hydrogen-bonded to the deprotonated nitrogen of the histidine imidazole (Movellan et al., 2021).
FIGURE 5. Schematic representation of the proton transfer flow associated to formation of hydrogen-bonded interactions between His37 residues.
According to this mechanism, the His37 tetrad would only sense acidification in the N-terminal side, reflecting previous observations obtained upon application of a voltage to drive protons outward in M2-transformed cells (Chizhmakov et al., 2003). The selectivity of proton selective transport is also justified as long as the proton must bind to and unbind from the histidine imidazole ring, which would explain the saturation attained at a moderate pH close to the histidine-tetrad pKa on the N-terminal side (Yi et al., 2009). In fact, the proton flux determined for the H37A mutant exhibits no pH dependence (Venkataraman et al., 2005). Finally, the low conductance would reflect the conformational transition of the Trp41 gate, which opens occasionally to form the conducting state that enables the proton to be released to the C-terminal pore.
More recently, measurement of 15N–15N J-couplings of 15N His37-labeled full-length M2 proton channel have been used to examine the occurrence of various imidazolium-imidazole hydrogen-bonding arrangements in the channel at low pH (Miao et al., 2015; Fu et al., 2020). According to these results, the proton conduction takes place through a dynamical process where low-barrier hydrogen bonds between pairs of His37 residues are broken and reformed.
From a computational point of view, the translocation of the proton has been examined by Carnevale et al. by combining classical Molecular Dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) methods. MD simulations were performed using methylammonium cation to examine the response of water molecules lying in the pore to an incoming positively charged moiety for the +2 state of the M2 channel (Carnevale et al., 2010). The results disclosed a broad minimum in the free energy landscape in the region below Val27 side chains, and the existence of a second local minimum in the region near Ala34, where the positive charge is stabilized by interactions with water molecules located above the His37 tetrad, and with the local dipoles of the backbone carbonyl groups. QM/MM computations aimed to determine the energetics of proton transport through the set of water molecules located above His37. In this case, the results supported a nearly free diffusion of the excess proton, as the transport of the proton from the water molecules to a neutral His37 occurs through several local minima separated by small barriers. Therefore, these findings suggest that the confining environment in the pore has little impact on the kinetics of proton diffusion, which should resemble proton transport in bulk water.
As a final remark, it is worth noting the computational study of proton conduction reported by Liang et al. (2016), which resorted to mutiscale simulations to follow the proton translocation through the channel considering distinct protonation states (fully deprotonated, +1 and +2) of the channel in models containing the TM region without (M2TM) and with (M2CD) the amphipatic helix, which together define the conductance domain. To this end, the free energy profile for proton migration was determined by combining multiscale reactive molecular dynamics (MS-RMD) and QM/MM MD simulations, which enabled to follow the proton diffusion from the N-terminal and C-terminals ends of the pore as well as (de)protonation of His37 tetrad. The free energy profile determined for the fully neutral and +1 states reveal a high barrier (>15 kcal mol−1) for His37 deprotonation to the viral interior, leading to a predicted conductance well below the experimental value, suggesting that the ion channel is inactive at high pH conditions. In the +2 state, however, although the electrostatic repulsion between the incoming proton increases the barrier from viral exterior to His37, deprotonation of His37 is facilitated by the increased electrostatic repulsion between the leaving proton and the His37 tetrad, this process being further assisted by larger exposure to water molecules and increased flexibility of Trp41residues. Thus, the His37 deprotonation barrier is decreased to ∼12 kcal mol−1, and the predicted conductance (1.0 and 7.7 for the M2CD and M2TM models, respectively) is in agreement with the measured values (ranging from 0.4 to 4.0; Mould et al., 2000b; Leiding et al., 2010; Sharma et al., 2010).
The whole of these studies illustrate the complexity of gaining a detailed knowledge of the mechanism of proton conduction in spite of the apparent simplicity encoded in the four-helix bundle of the M2 proton channel. The proton conductance determined from electrophysiological assays encompass a variety of molecular events that take place at different spatiotemporal scales, including the formation and breaking of hydrogen-bond complexes implicated in the proton transfer between water molecules and to/from the His37 imidazole-imidazolium ring, the coupling with conformational changes in the side chains of specific residues, and the structural rearrangements observed in the TM helices depending on the protonation state of the His37 tetrad.
The pKa of Histidines
The mechanistic models described above emphasize the relevance of the microscopic pKa of the His37 tetrad (Table 1) in assisting the proton flow through the inner pore of the M2 channel, which in turn would be associated with the pH-dependent equilibrium between open and closed states of the channel.
TABLE 1. Estimated pKa values determined for the His37 tetrad in different constructs of the M2 channel.
Using ultraviolet Resonance Raman spectroscopy and the transmembrane domain of the M2 proton channel (residues 22–46), Okada et al. reported that the proton channel is activated by low pH with a transition midpoint at pH 5.7 (Okada et al., 2001). Using the same construct, studies by Hu et al. measured pKa values of <5.0, 6.3, 8.2, and 8.2, leading to a +2 protonation state at neutral pH (Hu et al., 2006). Recent studies by Hu et al. supported pKa values of 4.0, 4.2, 6.8, and 7.6, which would support a +3 state at neutral pH, revealing the difficulty of elucidating the protonation preferences of the His37 tetrad and the influence of the experimental conditions used in these assays (Hu et al., 2012). On the other hand, NMR studies performed for a larger construct consisting of residues 21–97, which include the transmembrane helix but also the full cytoplasmic domain, led to pKa values of 7.1 and 5.4 for the second and third protonation (Liao et al., 2015). Likewise, a construct formed by residues 18–60, which encompass both the TM helix and the amphipathic helix that interacts with the interfacial region of the lipid bilayer, yielded estimated pKa values of 7.6 and 4.5 for the second and third protonations (Colvin et al., 2014), whereas pKas of 6.3 and 5.5 for the second and third protonations were determined by Miao et al. for the full-length protein (Miao et al., 2015). Even though the first two pKas are lower than the values estimated for shorter constructs, the authors indicate that at pH 6.2, where the channel should become activated, the M2 channel should be primarily in the +2 state. Moreover, the authors point out that there is little opportunity for a singly charged His37 tetrad.
The differences found in the pKas reported in Table 1 can be ascribed to the usage of different constructs, membrane-mimetic lipid environments, and conformational heterogeneity of the tetrameric helical bundle. In spite of these differences, there is consensus that the first two pKas appear to be above the endosomal pH that activates the M2 channel (pH ∼ 6) (Mould et al., 2000b), suggesting that the physiologically active state for the early activation of the M2 channel corresponds to the protonation state of +2 for the His tetrad.
Computational studies have also been performed with the aim to estimate the pKa of the His37 tetrad, while providing atomistic details of the pH-dependent protonation process. In this regard, Dong et al. performed QM/MM calculations to characterize the structural features of the +3 state and its subsequent deprotonation to the +2 species (Dong et al., 2013; Dong et al., 2014). Calculations were performed considering both a 4-fold symmetric arrangement of the His37 residues (histidine-box) and a 2-fold symmetric configuration (i.e., the dimer-of-dimers arrangement). The results showed that in the triply protonated state the two alternative models (histidine-box and dimer-of-dimers arrangements) converge to a single deprotonation mechanism, and exhibit similar free energy profiles with a barrier height of ∼6.5 kcal/mol to release a proton, supporting deprotonation as a mechanism for proton conduction.
Chen et al. have used constant pH replica-exchange MD simulations to determine the pKa values using a simulation model consisting of the TM helix (Chen et al., 2016). In this study, the simulated system consisted of the TM helix embedded in an explicit DMPC lipid bilayer and surrounded by an explicit solvent model (CHARM22 force field). However, the forces on titration coordinates were determined using the Generalized Born model in conjunction with an implicit membrane. Furthermore, a high-dielectric cylinder that encompass the channel was used to account for the continuous water wire in the pore. Finally, 12 replicas were used in simulations, covering a pH range of 3.5–9.0. The pKa values were estimated to be 8.3, 7.1, 6.2, and 5.7, which are in general close to the first three pKa’s obtained in previous ssNMR studies (see Table 1). The largest deviation is limited to the pKa of the last protonation, likely reflecting limitations arising from the Generalized Born model employed in the hybrid-solvent CpHMD, and to the conformational sensitivity of helices related to differences in the lipid environment used in simulations and experiments.
More recently, Torabifard et al. have used CpHMD simulations performed using explicit solvent and the multisite λ-dynamics to estimate the pKa values of the His37 tetrad (Torabifard et al., 2020). They used two models consisting of the TM helix alone (M2TM) and another one that combines the TM region with the amphipatic helix at the interface of the cytoplasmic domain (M2CD). These models were inserted in bilayers composed of a 4:1 ratio of DOPC:DOPE and solvated by water with an ionic force of 0.15 M (CHARMM36 force field). CpHMD simulations were performed using the multisite λ-dynamics approach (Knight and Brooks, 2011). This computational strategy yielded pKa values of 12.5, 10.5, 7.0, and 5.4 for the M2CD model, which are within the range of uncertainty for measured pKa values, especially regarding the first and second pKa values (larger deviations from the experimental data, however, were observed for the M2TM model).
Although the deviations observed between predicted and experimental pKa values may be likely ascribed to inaccuracies in the biomolecular force field, the approximated nature of the physical models used in calculations, incomplete convergence of the CpHMD simulations, and differences in the nature of lipid environments used in simulations and experiments, the results are encouraging and should provide a basis for gaining insight into the molecular factors that underlie the proton conduction in the M2 channel.
Other Factors That Influence the Proton Conduction of the M2 Channel
Besides the protonation state of the imidazole rings of the His37 tetrad, the efficiency of the M2 channel to perform the proton conduction may be influenced by a number of factors related to both the global structure of the ion channel and to the role of key residues other than His37 and Trp41. In this regard, Asp44 (see Figure 2) may contribute to modulate the proton conduction mechanism, as this residue forms direct or water-mediated hydrogen bonds with the indole moiety of Trp41 in the closed state (Acharya et al., 2010). In fact, Asp44 mutants tend to increase the population of the open state, as there is a loss of a stabilizing interaction of the closed state (Ma et al., 2013).
The influence exerted by the N-terminal ectodomain has been examined by Hong and coworkers (Kwon and Hong, 2016). The ectodomain is highly dynamic, although the motional flexibility is reduced for residues closer to the TM domain, possibly reflecting a tethering effect of the TM helix and the influence of lipid headgroups on the membrane surface. Furthermore, the electrostatic repulsion experienced by acidic residues (i.e., the ectodomain contains four acidic residues and only two cationic residues) appear to promote the adoption of TM conformations that would favor the binding of drugs in the inner pore (see below) even in the absence of the drug. This may be a factor that justifies the lower inhibitory concentration of full-length M2 compared to that of the ectodomain-truncated M2.
In this context, the potential role of Asp21 (see Figure 2) has also been recently highlighted by Jeong and Dyer (2017). Thus, by using a laser-induced pH jump coupled with time-resolved Trp fluorescence spectroscopy, protonation of His37 was estimated to occur in an unusually fast process, as reflected in a protonation rate of 1.6 ± 0.4 × 1010 M−1 s−1, suggesting that Asp21 at the end of the ectodomain and Asp24 at the beginning of the TM helix may act synergistically as proton-collecting antenna residues. The electrostatic field created by these residues would thus create a proton-capturing funnel at the entrance of the ion channel, enhancing proton harvesting from the surrounding aqueous phase. Protons would then be captured by hydrogen-bonded wires of water molecules within the pore, and the hydronium cation could be stabilized by interactions with pore-lining carbonyl groups as well as through bridging water molecules (Thomaston et al., 2015). Furthermore, they concluded that protonation of the His37 tetrad promotes opening of the C-terminal region, enhancing the solvent-exposure of Trp41, with a rate of (4 ± 2) × 103 s−1. The temporal decoupling between His37 protonation and this conformational change suggests that probably each proton transport cycle does not require a further conformational change after M2 activation.
Additional studies have examined the effect of the cytoplasmic tail on the conformational properties of the TM domain (Chizhmakov et al., 2003). The results obtained for the cytoplasmic-containing M2 channel reveal that even at neutral pH cationic histidines are present in the interior of the pore, which is in contrast to the results obtained for the TM peptide alone. This effect might be attributed to the acidic character of the cytoplasmic domain, which could facilitate opening of the TM pore at the His37 constriction. Furthermore, the presence of the cytoplasmic domain favors the adoption of a more helical conformation in the His37 backbone, suggesting an ordering effect on the four-helix bundle. Overall, these results provide a basis for the higher proton conductance of full-length M2 channel compared to the TM peptide.
Finally, the stability of the protonated states of the His37 tetrad may be affected by the presence of chloride anions in the interior of the pore. Inspection of the X-ray crystallographic data available for the M2 proton channel in the Protein Data Bank (Berman et al., 2000; Burley et al., 2021) reveals the occurrence of a subset of structures where chloride anions are found in the interior of the pore (this subset is collected in Table 2) (Thomaston and DeGrado, 2016; Thomaston et al., 2018; Thomaston et al., 2020; Thomaston et al., 2021). In four out of the five cases, the chloride anion occupies a well-defined position close to the plane formed by the Trp41 residues. Indeed, the chloride anion is stabilized by hydrogen-bond interaction with the indole NH groups, as well as by electrostatic interactions with the positive charge of Arg45 residues, though this latter effect is counterbalanced by the repulsion with the negative charge of Asp44 (Figure 6A). More strikingly, a chloride anion was observed along the four-fold axis in the plane defined by Gly34 residues in the X-ray structure of the S31N mutated variant (PDB entry 5C02).
TABLE 2. X-ray structures containing chloride anions at the C-terminus of the M2 proton channel.
FIGURE 6. (A) Representation of chloride anions (green spheres) in selected X-ray crystallographic structures of the transmembrane region of the M2 proton channel (adapted from PDB entries 5C02, 6BMZ, 6NV1, 6US8 and 6US9). (B) Representation of the electrostatic potential (kcal mol−1) in a plane passing through the interior of the M2 channel in (left) absence and (right) presence of chloride anions. Isocontours range from +12 (blue) to −12 kcal/mol (red). His37 and Trp41 residues are shown as sticks. Adapted with permission from (Llabrés et al., 2016) (Copyright 2016 American Chemical Society). For the sake of clarify one of the helices is not shown.
MD simulations have also revealed the presence of chloride anions in the interior of the channel pore, though the presence of chloride anions depends on the protonation state of the M2 proton channel. Thus, Mustafa et al. reported the occurrence of chloride anions at the level of the Trp41 tetrad for the +2 and +3 states in short simulations performed for the 22–46 construct (Mustafa et al., 2009). Wei and Pohorille also indicated that chloride anions penetrate the pore in all charged states of the M2 channel, but not for the fully unprotonated form (Wei and Pohorille, 2013). The position of the chloride anions was also affected by the charge state, as they were located in the space defined by the His37 and Trp41 planes for the diprotonated form of the channel. Gleed and Busath (2015), Gkeka et al. (2013) and Llabrés et al. (2016) also reported the presence of chloride anions close to the protonated His37 residues.
The presence of chloride anions in the pore may contribute to modulate the structural stability of the helical bundle, as suggested by Wei and Pohorille (2013). The MD simulations performed in presence of chloride ions revealed that they reside near Trp41, but their presence in this location of the pore depends on the charged state of the channel. Thus, anions do penetrate the pore, but to different degree, in all charged states, and they were absent only in the unprotonated form of M2 channel. Furthermore, simulations performed in the presence of phosphate anions confirmed their presence in the pore, suggesting that the size and specific chemical nature of counterions are not essential for the structural stabilization of the channel. Remarkably, additional simulations carried out in the presence of a uniform electrostatic field instead of explicit ions revealed a destabilization of the helical bundle in the +3 and +4 states, reflecting the electrostatic stabilization afforded by counterions that would balance the electrostatic repulsion between protonated His37 residues.
Besides the structural stabilization, the presence of anions in the pore has a strong influence on the nature of the electrostatic potential in the interior of the channel, which is highly positive in the absence of anions, but only slightly positive when chloride anions are present near the His37 sites, according to the MD simulations reported for the channel in the +2 state (Figure 6B) (Llabrés et al., 2016). Keeping in mind the functional relevance of electric fields in enzyme catalysis as well as in (un)binding of ligands to their macromolecular targets (Tan et al., 1993; Dillon et al., 2006; Fried and Boxer, 2017; Vaissier Welborn and Head-Gordon, 2019), the presence of negatively charged ions in the interior of the pore might be valuable to facilitate the diffusion of the proton from the bulk solvent along the N-terminal side of the luminal pore, and hence the transition from the +2 charged species to the triply protonated state. Furthermore, counterions may also facilitate the binding of Amt and related inhibitors by screening the electrostatic repulsion of the protonated amine of inhibitors with the charged His37 residues in the +2 state.
Finally, it is unclear whether the presence of counterions may have a functional role in assisting the translocation of the proton from the protonated His37 residues to the interior of the virion. At this point, it is worth noting that previous QM/MM calculations have suggested that the presence of chloride anions may increase the barrier for deprotonation (Dong et al., 2014), which would reduce the proton conduction. This effect, however, can be largely dependent upon the precise location of the chloride anions relative to the His37-Trp41 pair. Moreover, the functional impact might also be alleviated if the residence time of the chloride anion is affected by the structural fluctuations of the C-terminus in the triply charged state of the channel due to the larger electrostatic repulsion between the imidazolium rings. In particular, one can speculate that a fast exchange of chloride anions might transiently induce an electric field that would ease the transfer of the proton from the protonated His37 tetrad to the C-terminal side, thus facilitating the restoration of the +2 state and hence enabling the entry of new protons to the interior of the virion. In our view, this is an issue that deserves more attention.
Inhibition of the M2 Proton Channel
Beyond the definition of the molecular events implicated in proton conduction, a major challenge has been the identification of small molecules that can block the proton flow through the inner pore, and hence be used for therapeutic treatments against flu. Amantadine (Amt) and its ethyl analog, rimantadine (Rmt), were patented in 1961 and 1963, respectively (Du Pont Patent, 1961; Prichard, 1967), but they are no longer recommended for the treatment of flu infection. This obeys to several reasons, such as the limited effectiveness against influenza B virus, unwanted side effects, and the emergence of adamantane-resistant influenza A viral strains, primarily single mutated variants V27A, L26F and S31N. These mutations are found in a specific area facing the interior of the pore, suggesting therefore the location of the binding site in the M2 channel (Grambas et al., 1992; Hoyden, 2006; Deyde et al., 2007). For instance, in contrast to the large and widespread chemical shifts observed upon binding of Rmt to the M2 channel, only minor changes in chemical shifts were detected upon addition of a 4-fold molar excess of Rmt to the S31N mutated channel (Cady et al., 2009; Andreas et al., 2010; Andreas et al., 2012). For our purposes here, the discussion of the molecular determinants implicated in drug binding will be limited to the subset of structures available in the Protein Data Bank that contains drug-like inhibitors (see Tables 3, 4).
TABLE 3. X-ray and ss-NMR structures of the complexes formed by the wild type M2 proton channel and mutated variants with Amt and Rmt.
TABLE 4. X-ray and solution NMR structures of the complexes formed by the wild type M2 proton channel and mutated variants with adamantane-related compounds.
As noted above, the binding site for M2 inhibitors is located in the interior of the pore, filling a pocket shaped by residues Val27, Ser31, Gly34, and His37. Early neutron difraction (Duff et al., 1994) and computational (Samson and Kerr, 1993) studies already identified this pocket as the binding site that mediates the blockade of proton conduction by Amt. Detailed information about the binding mode of Amt and Rmt has subsequently been gained from structural studies based on X-ray diffraction, solution NMR and ss-NMR (Table 3; Schnell and Chou, 2008; Stouffer et al., 2008; Cady et al., 2009; Cady et al., 2010; Thomaston et al., 2018). The inspection of these structures reveals the existence of two binding sites (Figure 7). In most cases the drug is located in the interior of the pore filling the space located between Val27 and Gly34 (pore-binding model), but in two structures (PDB ID 2RLF and 6US8) the drug is also found outside the channel located in a peripheral site close to Leu43 and Asp44 in the C-terminus of the TM helical region (interface-binding model). The high-affinity site lies in the pore, as revealed from the analysis of REDOR dipolar dephasing between 13C-labeled M2TM and perdeuterated Amt in ss-NMR studies performed in a lipid bilayer (Cady et al., 2010). From a computational point of view, free energy calculations also pointed out that binding of Rmt to the pore-binding site is ∼7 kcal/mol more favorable relative to the interface-binding site, leading to more stable drug binding and channel inhibition (Li et al., 2008). The presence of the peripheral binding site can be attributed to the partitioning of the drug in the membrane-like lipidic environment (Duff et al., 1994; Gu et al., 2011). The studies reported by Cady et al. (2010) showed that the highest affinity binding site of amantadine is the N-terminal pore lumen, which is consistent with the known stoichiometry of binding (1:1 drug:channel) and the location of resistant mutations, such as V27A and S31N. However, increasing the drug:channel concentration up to 4:1, the drug contributes ∼7% of the amphiphiles composing the DMPC bilayer, and binding to the low-affinity, peripheral site is observed in NMR studies. Thus, when free amantadine is a major component of the membrane, Amt contacts the C-terminus of the protein, though the affinity for the peripheral site has been estimated to be ∼40-fold lower (see (Wang et al., 2015) for a detailed discussion of biophysical, computational and functional assays that addressed the binding mode of Amt and Rmt).
FIGURE 7. Representation of the binding mode of Rmt (shown as sticks with C-atoms in pink and yellow, respectively) to the TM helical region of the M2 proton channel (residues Val27 and His37 are shown as sticks; obtained from superposition of PDB structures 2RLF and 6BKL).
Inspection of these structures shows that in the pore the drug is oriented along the axis of the channel. Drug binding causes a dehydration of the channel, thus preventing the formation of a water wire from the N-terminus, although the amino group of Amt and Rmt is bound to water molecules that fill the space located between the drug and the His37 tetrad. In the specific case of Rmt, the two enantiomers bind to the pore with only slight differences in the hydration pattern, and exhibit similar values for association and dissociation kinetic rate constants and binding affinities, which was in agreement with the small differences in the relative binding free energy (around 0.3 kcal mol−1) determined from free energy calculations (Thomaston et al., 2020). Furthermore, Rmt was found to be more potent than Amt, as noted in a ∼14-fold ratio of their binding affinities to the M2 channel. On the other hand, it has been shown that binding of Amt reduces the pKa of His37 and alters the propensity to form hydrogen-bond interactions in the His37 tetrad (Hu et al., 2007; Cady et al., 2011).
From a computational standpoint, fine details of the binding mode of adamantane blockers have been extensively investigated. Yi et al. paid attention to the changes in conformational states of the M2 channel upon binding of Amt through comparison of the ensembles collected for both apo and bound forms (Yi et al., 2008; Yi et al., 2009). In particular, they observed a reduction in the conformational heterogeneity upon Amt binding, as reflected in a narrower distribution of kink angles around Gly34 in the transmembrane helices compared to the apo species, which agrees with the broadened resonances observed from ssNMR studies of this latter form. In contrast, large-kink angles were not observed in simulations. The occurrence of kinks also allows the access of water molecules that may stabilize the kinked helices via hydrogen bonds with the backbone carbonyl and amide groups around the kink, which in turn may also influence the proton conduction.
On the other hand, efforts have been conducted to predict the differences in binding affinity for adamantane inhibitors using free energy calculations, including implicit solvent/implicit membrane molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach (Homeyer et al., 2016) and free energy perturbation techniques (Gkeka et al., 2013; Ioannidis et al., 2016). Gratifyingly, alchemical free energy calculations performed for a set of 10 adamantane ligands reproduced satisfactorily the binding potency determined experimentally using isothermal titration calorimetry (ITC) against the M2 channel at high pH (pH = 8), especially keeping in mind that they cover a binding affinity range of only ∼2 kcal mol−1. Nevertheless, it is worth noting that the agreement between computational and experimental data was affected by the nature of the lipids used to model the membrane. In particular, a higher correlation between experimental and computed relative binding free energies was obtained when calculations were performed for the tetramer in a 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer, possibly due to a larger mimetic resemblance with the dodecylphosphocholine micelles used for ITC measurements.
Furthermore, the entire recognition process of Amt to the M2 channel has been explicitly explored using multiple-walkers well-tempered metadynamics calculations (Llabrés et al., 2016). The results suggested that Amt retains some degree of conformational motion in the pore allowing the adoption of two main orientations where the amino group is oriented to the N-terminus and C-terminus (denoted up and down, respectively). Binding of Amt would follow a sequential mechanism that would involve trapping of the protonated Amt by the negative electrostatic field created by the tetrad of Asp24, followed by the adoption of a transiently populated intermediate with Amt in the up state, and finally the reorientation in the interior of the pore to the down state, which is the thermodynamically favored orientation (Figure 8). Release of Amt from the down state to the channel mouth is destabilized by ∼12 kcal mol−1. The agreement with the experimental binding affinity (∼8–9 kcal mol−1) is reasonable, since the theoretical estimate omits the free energy contribution due to the release of Amt to the aqueous solvent, and the fact that the apparent Ki might not properly correspond to an equilibrium measurement of the inhibitory constant (Wang et al., 1993; Ma et al., 2009; Rosenberg and Casarotto, 2010). As a particular remark, this study also pointed out that the down → up transition depends on the presence of chloride anions in the C-terminus of the channel, since the stabilization of the down state can be enhanced by ∼6 kcal mol−1 due to the modulation of the electrostatic field in the interior of the pore. Finally, the barrier from Amt dissociation was estimated to be ∼19 kcal mol−1, which compares a value of 22 kcal mol−1 determined from electrophysiological assays (Wang et al., 1993).
FIGURE 8. Representation of the stepwise mechanism for Amt (carbon/nitrogen atoms shown as orange/blue-colored spheres.) binding to the M2 proton channel. (i) Amt is trapped by electrostatic interactions with Asp24. (ii) Thermal fluctuations of the helices facilitate crossing through the Val27 filter, and Amt fills the channel lumen in the up binding mode (iii) Finally, rotation of Amt leads to the down state, which is further stabilized by chloride anions. Chloride anions are shown as green spheres. For the sake of clarity, only three helices of the tetrameric channel are shown as green cartoon embedded in the lipid membrane. Reprinted with permission from (Llabrés et al., 2016) (Copyright 2016 American Chemical Society).
Mutation of Val27 to Ala causes a drastic alteration in the free energy surface for Amt binding, as noted in a stabilization of the up state less than 1 kcal mol−1, which agrees with the experimentally observed lack of inhibitory potency for the Amt-resistant V27A channel. This effect can be attributed to the increased accesible volume in the inner pore due to the V27A mutation, which encouraged the search of compounds with an expanded hydrophobic cage (Duque et al., 2011; Wang et al., 2011; Rey-Carrizo et al., 2014; Rey-Carrizo et al., 2015; Thomaston et al., 2018). In the case of the spiroadamantane 1 ((1r,1′S,3′S,5′S,7′S)-spiro [cyclohexane-1,2′-tricyclo [3.3.1.1∼3,7∼]decan]-4-amine) shown in Table 4 (Thomaston et al., 2018; Thomaston et al., 2020), comparison of the X-ray structures solved for the wild type and the V27A mutant shows that the compound exhibits a similar arrangement in the ion channel, which is found in the closed conformation (Figure 9) (Thomaston et al., 2020). The ligand exhibits only a slight shift depending on the nature of the residue at position 27. In the V27A complex, the amino group of 1 occupies approximately the same position as the ammonium group of Amt in complex with the wild type channel (PDB entry 6BKK) (Thomaston et al., 2018). This is facilitated by the larger free volume enabled by the presence of Ala27. Nevertheless, in the complex with the wild type channel, the adamantyl group of 1 and Amt overlap, but the ammonium group binds deeper in the pore at the expense of displacing few water molecules. MD simulations performed for the V27A complex with 1 reflected the shift of the ligand toward the N-terminus.
FIGURE 9. Representation of the binding mode of the spiroadamantane 1 bound to (left) the wild type M2 proton channel (construct 22–46; PDB ID 6BMZ) and (right) the V27A mutant (construct 21–61; PDB ID 6NV1 and 6OUG). The ligand is shown as sticks with C-atoms in pink, and the TM helical region of the M2 proton channel as cartoon (residues Val27, Ala27 and His37 are shown as sticks). For the sake of clarify one of the helices is not shown.
A combination of computational techniques, biophysical studies, functional assays and classical medicinal chemistry approaches has been used also to design potent inhibitors of V27A. This mutation imposes a reduction in the steric constraint found in the wild type channel, suggesting the compounds with a size-expanded hydrophobic cage would be better suited to fit the larger volume of the pocket in the V27A channel (Wang et al., 2011). An example is the spiroadamantane 4, which has an IC50 value of 18.7 μM against the wild type channel and of 0.3 μM against the V27A species. MD simulations showed that the ammonium group occupies the lower aqueous site forming solvent-mediated hydrogen bonds with the His37 residues in both wild type and V27A, but in this latter case the compound was slight shifted toward the N-terminus due to the extra space afforded by Ala27.
Other efforts to develop size-expanded analogs have led to the pyrrolidine derivative 5 reported by Vázquez and coworkers (Figure 10) (Duque et al., 2011; Rey-Carrizo et al., 2014). This compound inhibits the wild type M2 channel with an IC50 of 18 μM, while being ca. 26-fold more potent against the V27A mutant. MD simulations revealed an orientation consistent with the down state of Amt, with the amine nitrogen pointing toward the His37 tetrad with an average tilt angle (i.e., the deviation of the amine nitrogen from the pore axis) of 16°. Nevertheless, in the V27A complex 5 adopted the down orientation, but was also found in a up-like arrangement, presumably facilitated not only by the larger volume in the pore, but also by the widening of the helices at the location of Ala27, as the cross-diagonal distance between Cα atoms at position 27 was enlarged by 1.1 Å compared to the wild type channel. The effect of including other polycyclic scaffolds as the hydrophobic cage have led to analogs 6 and 7 (Figure 10), which exhibit a micromolar potency against the wild type channel (IC50 ∼ 2 μM), but a lower inhibition in the V27A mutant (IC50 ∼ 17.2 and 184.6 μM for 7 and 6, respectively) (Rey-Carrizo et al., 2015), suggesting a limiting effect in the size of the hydrophobic cage even for the V27A channel. At this point, MD simulations suggested that even the center-of-mass of 6 and 7 was located between the planes formed by S31 and A27 tetrads, the protonated amine was pointing to the N-terminus in five out of six MD simulations.
FIGURE 10. Representation of selected size-expanded analogs of Amt.
The resistance of the S31N channel, which is the most prevalent in currently circulating influenza viruses, to inhibition by Amt is reflected in the more variable motion and orientation of the ligand in the mutated channel. Thus, whereas Amt occupies a stable position in the wild type channel, it showed larger positional fluctuations in the pore of the mutated channel. This trait is reflected in a decreased energetic stability of the ligand in the proximity to the Asn31 side chains, as reflected in potential of mean force calculations. Finally, this destabilization was accompanied by an increase in the density of water molecules around and above Amt in the S31N channel compared to the wild type species, enabling water molecules to fill laterally the available space in the channel pore (Gleed and Busath, 2015). Furthermore, studies conducted for the binding of Rmt to the wild type M2 channel and the S31N variant have related the resistance with the higher dissociation rate constant determined for the mutated channel (Drakopoulos et al., 2018). In turn, this effect was attributed to the reduction in van der Waals interactions due to the shift of the ligand toward the C-terminus due to due to enhanced repulsive forces of the Asn side chains with the adamantyl ring in the mutated channel, and the concomitant weakening of the stabilizing contacts with Val27.
While these features reveal the challenges for developing inhibitors against the M2 S31N variant, significant progresses have been made in the last years. An example is compound 2 ((3S,5S,7S)-N-{[5-(thiophen-2-yl)-1,2-oxazol-3-yl]methyl}tricyclo [3.3.1.1∼3,7∼]decan-1-aminium), also denoted M2WJ332, Table 4). This compound inhibits the proton conduction in electrophysiological assays of the S31N channel with an IC50 of 16 μM, being however much less effective against the wild type ion channel (Rey-Carrizo et al., 2014). Remarkably, the NMR structure reveals that the compound binds the S31N channel in a different orientation compared to Amt and Rmt (Figure 10), since the adamantly ring is located facing the His37 tetrad, whereas the thienyl group is located at the level of the Val27 tetrad. This arrangement is facilitated by the slight expansion (∼1.0 Å) observed at the N-terminus of the TM helices relative to the wild type M2 channel. As noted in the analysis of MD simulations performed for the complex with 2, this binding mode is reinforced by hydrogen-bonds formed by the protonated amine and the isoxazole ring with three Asn31 residues (Wang et al., 2013).
Starting from compound 2, a successful rational design led to compound 3 ((3s,5s,7s)-N-[(5-bromothiophen-2-yl)methyl]tricyclo [3.3.1.1∼3,7∼]decan-1-aminium; Table 4) (Drakopoulos et al., 2018), which exhibits a similar potency in inhibiting both the wild type and the S31N channel (around 77% inhibition in two-electrode voltage clamp assays, and EC50 values of 4.6 and 1.8 μM, respectively). The chemical scaffold of compound 3 combines an adamantane moiety with a bromothiophene unit. However, the most remarkable structural feature is that 3 adopts completely flipped arrangements in the wild type and mutated channels (Figure 11A). Thus, the bromothiophene unit of the ligand faces the His37 tetrad in the wild type channel, but it is pointing toward the N-terminus in the S31N variant. MD simulations confirmed the structural stability of the flipped orientations in the wild type and mutated channels, enabling an interpretation of the structure-activity relationships (Wu et al., 2014). Inspection of the NMR structures (2MUW and 2MUV) does not support the involvement of halogen bonding, as this noncovalent interaction requires specific geometric features between the halogen atom and an electron-rich atom that are not fulfilled in these structures (see Kolár and Hobza, 2016 for details of halogen bonding). Furthermore, MD simulations revealed the formation of weak interactions between the bromine atom and water molecules, supplemented with transient, nonspecific contacts with the imidazole ring of His37 and the carbonyl group of Gly34 in the wild type channel, and with the Val27 side chains in the S31N channel (Wu et al., 2014).
FIGURE 11. (A) Representation of the binding mode of (left) compound 2 bound to the S31N channel (construct 19–49; PDB ID 2LY0), and compound 3 bound to (middle) the wild type M2 proton channel (construct 19–49; PDB ID 2MUW) and (right) the S31N mutant (construct 19–49; PDB ID 2MUV). The ligand is shown as sticks with C-atoms in pink, and the TM helical region of the M2 proton channel as cartoon (residues Val27, Ala27 and His37 are shown as sticks). For the sake of clarify one of the helices is not shown. (B) Superposition of selected inhibitors of wild type M2 ion channel and its V27A and S31N variants (taken from PBD ID 2KQT, 6US8, 6BMZ, 6NV1, 2LY0, 2MUW, and 2MUV). (C) Representation of the location of the protonated amine nitrogen atom along the channel axis taken from the previous subset of structures. For the sake of simplicity, only three helices of the tetrameric bundle are shown using the backbone skeleton of PDB ID 2KQT as cartoon is shown in plots (B,C).
It is also interesting to highlight the synthesis of aminoadamantane-CH2-aryl derivatives (8 in Figure 10) designed as sensitive probes for blockage of the wild type and S31N ion channels (Tzitzoglaki et al., 2020). These compounds retain the aminoadamantane unit present in 2, but the isoxazole ring has been replaced by a substituted phenyl ring. MD simulations have shown that compound 8 exhibits a stable binding in the S31N mutated channel. The binding mode locates the adamantyl group between Val27 and Gly34, and the phenyl unit fills the space between the side chains of Val27 residues. The results obtained from electrophysiological assays in oocytes evidence that blockage of the proton conduction in the wild type channel is relatively insensitive to chemical changes in the scaffold of 8, such as the insertion of a methylene between the adamantane and amino units or the replacement of adamantane by diamantane or triamantane. However, these changes have a drastic effect on the current blockage of the S31N channel, leading to lower association and higher dissociation rate constants relative to the parent compound 8. These results point out the larger sensitivity of the S31N ion channel to the chemical features of inhibitors upon binding to the pore.
Superposition of the compounds shown in Tables 3, 4 is useful to exemplify the ability of the inhibitors to exploit different sites and binding modes along the interior of the pore in order to prevent proton transport in the wild type M2 channel and the V27A and S31N mutated species. (Figures 11B,C). This feature has been highlighted in previous studies that combined the structural analysis of bound channels with metadynamics calculations (Wang et al., 2013; Gianti et al., 2015), which pointed out the existence of distint locations that can be occupied by the protonated amine along the axis of the channel. These studies also remarked the role played by waters molecules inside the pore, as they contribute to stabilize the inhibitor by facilitating hydrogen-bond bridges with carbonyl groups of the transmembrane helices as well as with the His37 tetrad. From a practical perspective, the presence of multiple interaction sites for the protonated amine and the structure of water clusters that stabilize the inhibitor can be useful hallmarks for guiding the design of new drugs. In this regard, recent studies have indicated that the ammonium group of these inhibitors can act as a hydronium mimic upon binding to the channel pore, suggesting that the adamantane-like compounds function as mechanism-based inhibitors as long as they exploit structural features of the channel implicated in proton transport (Watkins et al., 2020).
As a final remark, it is worth noting a very recent study by Kolokouris and coworkers where a systematic analysis of setup conditions and methodological details, covering from docking calculations to MD simulations, is presented for the M2 channel (Kolocouris et al., 2021). This study examines distinctive traits, such as the usage of M2 (22–46) and M2 (22–62) constructs, lipid content in the bilayer (1,2-dimyristoyl-sn-glycero-3-phosphocholine,1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanol-amine), and the effect of cholesterol, paying particular attention to the location of channel blockers and water structure, the presence of chloride anions, structural properties of the helical bundle, and the binding site of cholesterol. While this study deserves practical interest for the choice of parameters in simulations of the M2 proton channel (and related viroporins), we limit ourselves to highlight three findings. First, the results obtained for M2 complexes with Amt at pH 8 revealed that chloride anions approached Arg45 and occupied a position located between Trp41 and Arg45 (at a distance close to 15.5–17.5 Å from the Amt nitrogen atom) presence of chloride anions stabilized in the interior of the pore. The anions approached in the C-terminus. The anions approached Arg45. Second, the simulations with CHARMM36 preserved the structure of the amphipatic helices, providing a correct description for the wedge shape geometry that might be required to modulate the saddle-shape curvatuve of the membrane for the release of virions (Wang and Hong, 2015b; Martyna et al., 2017). Finally, the ability of cholesterol to interact with the tetrameric channel, forming van der Waals interactions with Ile39, Leu40, Leu43, Leu46, Phe47, Ile51 and Phe54, and a hydrogen bond with Ser50, which is in agreement with previous experimental evidence (Ekanayake et al., 2016; Elkins et al., 2017; Elkins et al., 2018), paving the way to future studies about the role of cholesterol in mediating the M2 protein function.
Conclusion
Resistance to current pharmacological treatments is a severe health challenge, as illustrated by the emergence of viral strains resistant to drugs interfering with the M2 proton channel, such as Amt and Rmt, but also to antiviral agents targeting neuraminidase, such as oseltamivir, as noticed in a resistant strain associated to the His275Tyr mutation in this protein (Meijer et al., 2009; Ginex and Luque, 2021). Currently circulating seasonal H1N1 and H3N2 strains are now resistant to adamantanes, and oseltamivir is no longer effective against the pandemic H1N1 subtypes (Hussain et al., 2017).
Finding new active antiviral compounds is an urgent need, which may be explored resorting to the application of a combinatorial therapy that should potentiate the effect of two or three drugs endowed with distinct mechanisms of action, hopefully being less susceptible to inactivating mutations (Dunning et al., 2014). A long-term strategy is the identification and characterization of the molecular mechanisms associated to the activity of proteins relevant for the life cycle of the virus, as this knowledge may disclose novel strategies for the design of antiviral compounds.
This approach is illustrated by the research efforts invested in the M2 proton channel. In spite of the apparent simplicity posed by the membrane-embedded four-helix bundle, the flow of protons reflects a complex synergy between different molecular events, involving the acidity of the His37 residues, the structural coupling between His37 and the gating Trp41 residue, the structure of the hydrating water molecules in the pore, the pH-dependent conformational arrangement of the TM helices, the electrostatic funnel created by acidic residues in the N-terminus of the channel, and the influence of the amphiphilic helix and cytoplasmic tail on the flexibility of the TM domain, in addition to the potential influence exerted by other factors such as the nature of the lipid environment in the viral capside.
Although a full understanding of the interplay between these factors, which encompass a diverse range of spatial and time processes, has not yet been achieved, the progress consolidated in the last decades through the combined use of structural, biophysical, physiological and computational techniques is impressive, and the successful development of compound 3 is encouraging. This scientific background defines the framework to address novel questions about the mechanisms of drug resistance and the guidance to novel antiviral and treatment approaches, such as the potential dependence between drug resistance and the nature of the viral strain or the infected cell type, which have been the subject of recent studies (Musharrafieh et al., 2019; Musharrafieh et al., 2020). Hopefully, the outcome will be valuable to enrich our current anti-influenza therapeutic arsenal in the form of effective antivirals less susceptible to drug resistance.
Author Contributions
All authors contributed to conception, design and writing of the manuscript. EA and BS wrote the draft of the manuscript and prepared Figures and Tables. CE, MM and FJL revised and completed the draft and prepared the final version of the manuscript.
Funding
We thank the Spanish Ministerio de Ciencia e Innovación (PID2020-117646RB-I00 MCIN/AEI/10.13039/501100011033, and Maria de Maetzu MDM-2017-0767 AEI/FEDER) for financial support and the Barcelona Supercomputing Center (BCV-2019-1-0009) and the Consorci de Serveis Universitaris de Catalunya (CSUC; Molecular Recognition project) for computational resources. E.A. thanks AGAUR (Generalitat of Catalunya; 2018FI-B1-00001) for a fellowship.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Acharya, R., Carnevale, V., Fiorin, G., Levine, B. G., Polishchuk, A. L., Balannik, V., et al. (2010). Structure and Mechanism of Proton Transport through the Transmembrane Tetrameric M2 Protein Bundle of the Influenza A Virus. Proc. Natl. Acad. Sci. 107, 15075–15080. doi:10.1073/pnas.1007071107
Agmon, N. (1995). The Grotthuss Mechanism. Chem. Phys. Lett. 244, 456–462. doi:10.1016/0009-2614(95)00905-j
Andreas, L. B., Eddy, M. T., Pielak, R. M., Chou, J., and Griffin, R. G. (2010). Magic Angle Spinning NMR Investigation of Influenza A M218−60: Support for an Allosteric Mechanism of Inhibition. J. Am. Chem. Soc. 132, 10958–10960. doi:10.1021/ja101537p
Andreas, L. B., Eddy, M. T., Chou, J. J., and Griffin, R. G. (2012). Magic-angle-spinning NMR of the Drug Resistant S31N M2 Proton Transporter from Influenza A. J. Am. Chem. Soc. 134, 7215–7218. doi:10.1021/ja3003606
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The Protein Data Bank. Protein Data Bank Nucleic Acids Res. 28, 235–242. doi:10.1093/nar/28.1.235
Burley, S. K., Bhikadiya, C., Bi, C., Bittrich, S., Chen, L., Crichlow, G., et al. (2021). RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nuc. Acids Res. 49, 437–451. doi:10.1093/nar/gkaa1038
Cady, S. D., Mishanina, T. V., and Hong, M. (2009). Structure of Amantadine-Bound M2 Transmembrane Peptide of Influenza A in Lipid Bilayers from Magic-Angle-Spinning Solid-State NMR: The Role of Ser31 in Amantadine Binding. J. Mol. Biol. 385, 1127–1141. doi:10.1016/j.jmb.2008.11.022
Cady, S. D., Schmidt-Rohr, K., Wang, J., Soto, C. S., DeGrado, W. F., and Hong, M. (2010). Structure of the Amantadine Binding Site of Influenza M2 Proton Channels in Lipid Bilayers. Nature 463, 689–692. doi:10.1038/nature08722
Cady, S. D., Wang, J., Wu, Y., DeGrado, W. F., and Hong, M. (2011). Specific Binding of Adamantane Drugs and Direction of Their Polar Amines in the Pore of the Influenza M2 Transmembrane Domain in Lipid Bilayers and Dodecylphosphocholine Micelles Determined by NMR Spectroscopy. J. Am. Chem. Soc. 133, 4274–4284. doi:10.1021/ja102581n
Carnevale, V., Fiorin, G., Levine, B. G., DeGrado, W. F., and Klein, M. L. (2010). Multiple Proton Confinement in the M2 Channel from the Influenza A Virus. J. Phys. Chem. C 114, 20856–20863. doi:10.1021/jp107431g
Chen, W., Huang, Y., and Shen, J. (2016). Conformational Activation of a Transmembrane Proton Channel from Constant pH Molecular Dynamics. J. Phys. Chem. Lett. 7, 3961–3966. doi:10.1021/acs.jpclett.6b01853
Cheung, T. K. W., and Poon, L. L. M. (2007). Biology of Influenza a Virus. Ann. N. Y. Acad. Sci. 1102, 1–25. doi:10.1196/annals.1408.001
Chizhmakov, I. V., Geraghty, F. M., Ogden, D. C., Hayhurst, A., Antoniou, M., and Hay, A. J. (1996). Selective Proton Permeability and pH Regulation of the Influenza Virus M2 Channel Expressed in Mouse Erythroleukaemia Cells. J. Physiol. 494, 329–336. doi:10.1113/jphysiol.1996.sp021495
Chizhmakov, I. V., Ogden, D. C., Geraghty, F. M., Hayhurst, A., Skinner, A., Betakova, T., et al. (2003). Differences in Conductance of M2 Proton Channels of Two Influenza Viruses at Low and High pH. J. Physiol. 546, 427–438. doi:10.1113/jphysiol.2002.028910
Colvin, M. T., Andreas, L. B., Chou, J. J., and Griffin, R. G. (2014). Proton Association Constants of His 37 in the Influenza-A M218-60 Dimer-Of-Dimers. Biochemistry 53, 5987–5994. doi:10.1021/bi5005393
Cox, N. J., and Subbarao, K. (2000). Global Epidemiology of Influenza: Past and Present. Annu. Rev. Med. 51, 407–421. doi:10.1146/annurev.med.51.1.407
De Vlugt, C., Sikora, D., and Pelchat, M. (2018). Insight into Influenza: A Virus Cap-Snatching. Viruses 10, 641. doi:10.3390/v10110641
Deyde, V. M., Xu, X., Bright, R. A., Shaw, M., Smith, C. B., Zhang, Y., et al. (2007). Surveillance of Resistance to Adamantanes Among Influenza A(H3N2) and A(H1N1) Viruses Isolated Worldwide. J. Infect. Dis. 196, 249–257. doi:10.1086/518936
Dillon, P. F., Root-Bernstein, R. S., and Lieder, C. M. (2006). Molecular Shielding of Electric Field Complex Dissociation. Biophysical J. 90, 1432–1438. doi:10.1529/biophysj.105.071969
Dong, H., Fiorin, G., Degrado, W. F., and Klein, M. L. (2013). Exploring Histidine Conformations in the M2 Channel Lumen of the Influenza A Virus at Neutral pH via Molecular Simulations. J. Phys. Chem. Lett. 4, 3067–3071. doi:10.1021/jz401672h
Dong, H., Fiorin, G., Degrado, W. F., and Klein, M. L. (2014). Proton Release from the Histidine-Tetrad in the M2 Channel of the Influenza A Virus. J. Phys. Chem. B 118, 12644–12651. doi:10.1021/jp5102225
Dou, D., Revol, R., Östbye, H., Wang, H., and Daniels, R. (2018). Influenza A Virus Cell Entry, Replication, Virion Assembly and Movement. Front. Immunol. 9, 1581. doi:10.3389/fimmu.2018.01581
Drakopoulos, A., Tzitzoglaki, C., McGuire, K., Hoffmann, A., Konstantinidi, A., Kolokouris, D., et al. (2018). Unraveling the Binding, Proton Blockage, and Inhibition of Influenza M2 WT and S31N by Rimantadine Variants. ACS Med. Chem. Lett. 9, 198–203. doi:10.1021/acsmedchemlett.7b00458
Du Pont Patent (1961). Pharmaceutical Compositions and Methods Utilizing 1-aminoadamantane and its Derivatives. US3310469A. United States Patent and Trademark Office. . Available at: https://www.uspto.gov.(Accessed September 14, 2021).
Duff, K. C., Gilchrist, P. J., Saxena, A. M., and Bradshaw, J. P. (1994). Neutron Diffraction Reveals the Site of Amantadine Blockade in the Influenza A M2 Ion Channel. Virology 202, 287–293. doi:10.1006/viro.1994.1345
Dunning, J., Baillie, J. K., Cao, B., and Hayden, F. G. (2014). International Severe Acute Respiratory and Emerging Infection Consortium (ISARIC)Antiviral Combinations for Severe Influenza. Lancet Infect. Dis. 14, 1259–1270. doi:10.1016/s1473-3099(14)70821-7
Duque, M. D., Ma, C., Torres, E., Wang, J., Naesens, L., Juárez-Jiménez, J., et al. (2011). Exploring the Size Limit of Templates for Inhibitors of the M2 Ion Channel of Influenza A Virus. J. Med. Chem. 54, 2646–2657. doi:10.1021/jm101334y
Ekanayake, E. V., Fu, R., and Cross, T. A. (2016). Structural Influences: Cholesterol, Drug, and Proton Binding to Full-Length Influenza A M2 Protein. Biophysical J. 110, 1391–1399. doi:10.1016/j.bpj.2015.11.3529
Elkins, M. R., Williams, J. K., Gelenter, M. D., Dai, P., Kwon, B., Sergeyev, I. V., et al. (2017). Cholesterol-binding Site of the Influenza M2 Protein in Lipid Bilayers from Solid-State NMR. Proc. Natl. Acad. Sci. USA 114, 12946–12951. doi:10.1073/pnas.1715127114
Elkins, M. R., Sergeyev, I. V., and Hong, M. (2018). Determining Cholesterol Binding to Membrane Proteins by Cholesterol 13C Labeling in Yeast and Dynamic Nuclear Polarization NMR. J. Am. Chem. Soc. 140, 15437–15449. doi:10.1021/jacs.8b09658
Fried, S. D., and Boxer, S. G. (2017). Electric Fields and Enzyme Catalysis. Annu. Rev. Biochem. 86, 387–415. doi:10.1146/annurev-biochem-061516-044432
Fu, R., Miao, Y., Qin, H., and Cross, T. A. (2020). Observation of the Imidazole-Imidazolium Hydrogen Bonds Responsible for Selective Proton Conductance in the Influenza A M2 Channel. J. Am. Chem. Soc. 142, 2115–2119. doi:10.1021/jacs.9b09985
Gianti, E., Carnevale, V., DeGrado, W. F., Klein, M. L., and Fiorin, G. (2015). Hydrogen-Bonded Water Molecules in the M2 Channel of the Influenza A Virus Guide the Binding Preferences of Ammonium-Based Inhibitors. J. Phys. Chem. B 119, 1173–1183. doi:10.1021/jp506807y
Ginex, T., and Luque, F. J. (2021). Searching for Effective Antiviral Small Molecules against Influenza A Virus: A Patent Review. Expert Opin. Ther. Patents 31, 53–66. doi:10.1080/13543776.2020.1831471
Gkeka, P., Eleftheratos, S., Kolocouris, A., and Cournia, Z. (2013). Free Energy Calculations Reveal the Origin of Binding Preference for Aminoadamantane Blockers of Influenza A/M2TM Pore. J. Chem. Theor. Comput. 9, 1272–1281. doi:10.1021/ct300899n
Gleed, M. L., and Busath, D. D. (2015). Why Bound Amantadine Fails to Inhibit Proton Conductance According to Simulations of the Drug-Resistant Influenza A M2 (S31N). J. Phys. Chem. B 119, 1225–1231. doi:10.1021/jp508545d
Grambas, S., Bennett, M. S., and Hay, A. J. (1992). Influence of Amantadine Resistance Mutations on the pH Regulatory Function of the M2 Protein of Influenza A Viruses. Virology 191, 541–549. doi:10.1016/0042-6822(92)90229-i
Gu, R.-X., Liu, L. A., Wei, D.-Q., Du, J.-G., Liu, L., and Liu, H. (2011). Free Energy Calculations on the Two Drug Binding Sites in the M2 Proton Channel. J. Am. Chem. Soc. 133, 10817–10825. doi:10.1021/ja1114198
Han, J., Perez, J., Schäfer, A., Cheng, H., Peet, N., Rong, L., et al. (2018). Influenza Virus: Small Molecule Therapeutics and Mechanisms of Antiviral Resistance. Curr. Med. Chem. 25, 5115–5127. doi:10.2174/0929867324666170920165926
Hayward, A. C., Fragaszy, E. B., Bermingham, A., Wang, L., Copas, A., Edmunds, W. J., et al. (2014). Comparative Community burden and Severity of Seasonal and Pandemic Influenza: Results of the Flu Watch Cohort Study. Lancet Respir. Med. 2, 445–454. doi:10.1016/s2213-2600(14)70034-7
Holsinger, L. J., Nichani, D., Pinto, L. H., and Lamb, R. A. (1994). Influenza A Virus M2 Ion Channel Protein: A Structure-Function Analysis. J. Virol. 68, 1551–1563. doi:10.1128/jvi.68.3.1551-1563.1994
Homeyer, N., Ioannidis, H., Kolarov, F., Gauglitz, G., Zikos, C., Kolocouris, A., et al. (2016). Interpreting Thermodynamic Profiles of Aminoadamantane Compounds Inhibiting the M2 Proton Channel of Influenza A by Free Energy Calculations. J. Chem. Inf. Model. 56, 110–126. doi:10.1021/acs.jcim.5b00467
Honda, A., Mizumoto, K., and Ishihama, A. (2002). Minimum Molecular Architectures for Transcription and Replication of the Influenza Virus. Proc. Natl. Acad. Sci. 99, 13166–13171. doi:10.1073/pnas.152456799
Hoyden, F. G. (2006). Antiviral Resistance in Influenza Viruses - Implications for Management and Pandemic Response. N. Eng. J. Med. 354, 785–788. doi:10.1056/NEJMp068030
Hu, J., Fu, R., Nishimura, K., Zhang, L., Zhou, H.-X., Busath, D. D., et al. (2006). Histidines, Heart of the Hydrogen Ion Channel from Influenza A Virus: Toward an Understanding of Conductance and Proton Selectivity. Proc. Natl. Acad. Sci. 103, 6865–6870. doi:10.1073/pnas.0601944103
Hu, J., Fu, R., and Cross, T. A. (2007). The Chemical and Dynamical Influence of the Anti-viral Drug Amantadine on the M2 Proton Channel Transmembrane Domain. Biophys. J. 93, 276–283. doi:10.1529/biophysj.106.102103
Hu, F., Luo, W., and Hong, M. (2010). Mechanisms of Proton Conduction and Gating in Influenza M2 Proton Channels from Solid-State NMR. Science 330, 505–508. doi:10.1126/science.1191714
Hu, F., Schmidt-Rohr, K., and Hong, M. (2012). NMR Detection of pH-Dependent Histidine-Water Proton Exchange Reveals the Conduction Mechanism of a Transmembrane Proton Channel. J. Am. Chem. Soc. 134, 3703–3713. doi:10.1021/ja2081185
Hussain, M., Galvin, H., Haw, T. Y., Nutsford, A., and Husain, M. (2017). Drug Resistance in Influenza A Virus: The Epidemiology and Management. Idr 10, 121–134. doi:10.2147/idr.s105473
Ioannidis, H., Drakopoulos, A., Tzitzoglaki, C., Homeyer, N., Kolarov, F., Gkeka, P., et al. (2016). Alchemical Free Energy Calculations and Isothermal Titration Calorimetry Measurements of Aminoadamantanes Bound to the Closed State of Influenza A/M2TM. J. Chem. Inf. Model. 56, 862–876. doi:10.1021/acs.jcim.6b00079
Iuliano, A. D., Roguski, K. M., Chang, H. H., Muscatello, D. J., Palekar, R., Tempia, S., et al. (2018). Estimates of Global Seasonal Influenza-Associated Respiratory Mortality: A Modelling Study. Lancet 391, 1285–1300. doi:10.1016/S0140-6736(17)33293-2
Jeong, B.-S., and Dyer, R. B. (2017). Proton Transport Mechanism of M2 Proton Channel Studied by Laser-Induced pH Jump. J. Am. Chem. Soc. 139, 6621–6628. doi:10.1021/jacs.7b00617
Julkunen, I., Sareneva, T., Pirhonen, J., Ronni, T., Melén, K., and Matikainen, S. (2001). Molecular Pathogenesis of Influenza A Virus Infection and Virus-Induced Regulation of Cytokine Gene Expression. Cytokine Growth Factor. Rev. 12, 171–180. doi:10.1016/s1359-6101(00)00026-5
Jung, H. E., and Lee, H. K. (2020). Host Protective Immune Responses against Influenza A Virus Infection. Viruses 12, 504. doi:10.3390/v12050504
Khurana, E., Peraro, M. D., DeVane, R., Vemparala, S., DeGrado, W. F., and Klein, M. L. (2009). Molecular Dynamics Calculations Suggest a Conduction Mechanism for the M2 Proton Channel from Influenza A Virus. Proc. Natl. Acad. Sci. 106, 1069–1074. doi:10.1073/pnas.0811720106
Knight, J. L., and Brooks, C. L. (2011). Multisite λ Dynamics for Simulated Structure-Activity Relationship Studies. J. Chem. Theor. Comput. 7, 2728–2739. doi:10.1021/ct200444f
Kolár, M. H., and Hobza, P. (2016). Computer Modeling of Halogen Bonds and Other σ-hole Interactions. Chem. Rev. 116, 5155–5187. doi:10.1021/acs.chemrev.5b00560
Kolocouris, D., Kalenderoglou, I. E., and Kolocouris, A. (2021). Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers. J. Chem. Inf. Model. 61, 5550–5568. doi:10.1021/acs.jcim.1c00264
Krammer, F., Smith, G. J. D., Fouchier, R. A. M., Peiris, M., Kedzierska, K., Doherty, P. C., et al. (2018). Influenza. Nat. Rev. Dis. Primers 4, 3. doi:10.1038/s41572-018-0002-y
Kwon, B., and Hong, M. (2016). The Influenza M2 Ectodomain Regulates the Conformational Equilibria of the Transmembrane Proton Channel: Insights from Solid-State Nuclear Magnetic Resonance. Biochemistry 55, 5387–5397. doi:10.1021/acs.biochem.6b00727
Leiding, T., Wang, J., Martinsson, J., DeGrado, W. F., and Arsköld, S. P. (2010). Proton and Cation Transport Activity of the M2 Proton Channel from Influenza A Virus. Proc. Natl. Acad. Sci. 107, 15409–15414. doi:10.1073/pnas.1009997107
Li, C., Yi, M., Hu, J., Zhou, H.-X., and Cross, T. A. (2008). Solid-state NMR and MD Simulations of the Antiviral Drug Amantadine Solubilized in DMPC Bilayers. Biophysical J. 94, 1295–1302. doi:10.1529/biophysj.107.112482
Liang, R., Swanson, J. M. J., Madsen, J. J., Hong, M., DeGrado, W. F., and Voth, G. A. (2016). Acid Activation Mechanism of the Influenza A M2 Proton Channel. Proc. Natl. Acad. Sci. U. S. A. 113, E6955–E6964. doi:10.1073/pnas.1615471113
Liao, S. Y., Yang, Y., Tietze, D., and Hong, M. (2015). The Influenza M2 Cytoplasmic Tail Changes the Proton-Exchange Equilibria and the Backbone Conformation of the Transmembrane Histidine Residue to Facilitate Proton Conduction. J. Am. Chem. Soc. 137, 6067–6077. doi:10.1021/jacs.5b02510
Lin, T.-I., and Schroeder, C. (2001). Definitive Assignment of Proton Selectivity and Attoampere Unitary Current to the M2 Ion Channel Protein of Influenza A Virus. J. Virol. 75, 3647–3656. doi:10.1128/jvi.75.8.3647-3656.2001
Llabrés, S., Juárez-Jiménez, J., Masetti, M., Leiva, R., Vázquez, S., Gazzarrini, S., et al. (2016). Mechanism of the Pseudoirreversible Binding of Amantadine to the M2 Proton Channel. J. Am. Chem. Soc. 138, 15345–15358. doi:10.1021/jacs.6b07096
Loregian, A., Mercorelli, B., Nannetti, G., Compagnin, C., and Palù, G. (2014). Antiviral Strategies against Influenza Virus: Towards New Therapeutic Approaches. Cel. Mol. Life Sci. 71, 3659–3683. doi:10.1007/s00018-014-1615-2
Ma, C., Polishchuk, A. L., Ohigashi, Y., Stouffer, A. L., Schön, A., Magavern, E., et al. (2009). Identification of the Functional Core of the Influenza A Virus A/M2 Proton-Selective Ion Channel. Proc. Natl. Acad. Sci. 106, 12283–12288. doi:10.1073/pnas.0905726106
Ma, C., Fiorin, G., Carnevale, V., Wang, J., Lamb, R. A., Klein, M. L., et al. (2013). Asp44 Stabilizes the Trp41 Gate of the M2 Proton Channel of Influenza A Virus. Structure 21, 2033–2041. doi:10.1016/j.str.2013.08.029
Manzoor, R., Igarashi, M., and Takada, A. (2017). Influenza A Virus M2 Protein: Roles from Ingress to Egress. Int. J. Mol. Sci. 18, 1–16. doi:10.3390/ijms18122649
Martyna, A., Bahsoun, B., Badham, M. D., Srinivasan, S., Howard, M. J., and Rossman, J. S. (2017). Membrane Remodeling by the M2 Amphipathic helix Drives Influenza Virus Membrane Scission. Sci. Rep. 7, 44695. doi:10.1038/srep44695
McCown, M. F., and Pekosz, A. (2006). Distinct Domains of the Influenza A Virus M2 Protein Cytoplasmic Tail Mediate Binding to the M1 Protein and Facilitate Infectious Virus Production. J. Virol. 80, 8178–8189. doi:10.1128/jvi.00627-06
Meijer, A., Lackenby, A., Hungnes, O., Lina, B., van der WerfSchweiger, S. B., Schweiger, B., et al. (2009). Oseltamivir-resistant Influenza Virus A (H1N1), Europe, 2007-08 Season. Emerg. Infect. Dis. 15, 552–560. doi:10.3201/eid1504.181280
Miao, Y., Fu, R., Zhou, H.-X., and Cross, T. A. (2015). Dynamic Short Hydrogen Bonds in Histidine Tetrad of Full-Length M2 Proton Channel Reveal Tetrameric Structural Heterogeneity and Functional Mechanism. Structure 23, 2300–2308. doi:10.1016/j.str.2015.09.011
Mifsud, E. J., Hayden, F. G., and Hurt, A. C. (2019). Antivirals Targeting the Polymerase Complex of Influenza Viruses. Antiviral Res. 169, 104545. doi:10.1016/j.antiviral.2019.104545
Moorthy, N. S., Poongavanam, V., and Pratheepa, V. (2014). Viral M2 Ion Channel Protein: A Promising Target for Anti-influenza Drug Discovery. Mini Rev. Med. Chem. 14, 819–830. doi:10.2174/138955751410141020150822
Mould, J. A., Drury, J. E., Frings, S. M., Kaupp, U. B., Pekosz, A., Lamb, R. A., et al. (2000). Permeation and Activation of the M2 Ion Channel of Influenza A Virus. J. Biol. Chem. 275, 31038–31050. doi:10.1074/jbc.m003663200
Mould, J. A., Li, H.-C., Dudlak, C. S., Lear, J. D., Pekosz, A., Lamb, R. A., et al. (2000). Mechanism for Proton Conduction of the M2 Ion Channel of Influenza A Virus. J. Biol. Chem. 275, 8592–8599. doi:10.1074/jbc.275.12.8592
Movellan, K. T., Wegstroth, M., Overkamp, K., Leonov, A., Becker, S., and Andreas, L. B. (2020). Imidazole-imidazole Hydrogen Bonding in the pH-Sensing Histidine Side Chains of Influenza A M2. J. Am. Chem. Soc. 142, 2704–2708. doi:10.1021/jacs.9b10984
Movellan, K. T., Dervişoğlu, R., Becker, S., and Andreas, L. B. (2021). Pore‐Bound Water at the Key Residue Histidine 37 in Influenza A M2. Angew. Chem. Int. Ed. 60, 24075–24079. doi:10.1002/anie.202103955
Musharrafieh, R., Lagarias, P. I., Ma, C., Tan, G. S., Kolocouris, A., and Wang, J. (2019). The L46P Mutant Confers a Novel Allosteric Mechanism of Resistance toward the Influenza A Virus M2 S31N Proton Channel Blockers. Mol. Pharmacol. 96, 148–157. doi:10.1124/mol.119.116640
Musharrafieh, R., Lagarias, P., Ma, C., Hau, R., Romano, A., Lambrinidis, G., et al. (2020). Investigation of the Drug Resistance Mechanism of M2-S31N Channel Blockers through Biomolecular Simulations and Viral Passage Experiments. ACS Pharmacol. Transl. Sci. 3, 666–675. doi:10.1021/acsptsci.0c00018
Mustafa, M., Henderson, D. J., and Busath, D. D. (2009). Free-Energy Profiles for Ions in the Influenza M2-TMD Channel. Proteins 76, 794–807. doi:10.1002/prot.22376
Okada, A., Miura, T., and Takeuchi, H. (2001). Protonation of Histidine and Histidine−Tryptophan Interaction in the Activation of the M2 Ion Channel from Influenza A Virus. Biochemistry 40, 6053–6060. doi:10.1021/bi0028441
Pielak, R. M., and Chou, J. J. (2011). Influenza M2 Proton Channels. Biochim. Biophys. Acta (Bba) - Biomembr. 1808, 522–529. doi:10.1016/j.bbamem.2010.04.015
Pinto, L. H., Holsinger, L. J., and Lamb, R. A. (1992). Influenza Virus M2 Protein Has Ion Channel Activity. Cell 69, 517–528. doi:10.1016/0092-8674(92)90452-i
Pinto, L. H., Dieckmann, G. R., Gandhi, C. S., Papworth, C. G., Braman, J., Shaughnessy, M. A., et al. (1997). A Functionally Defined Model for the M2 Proton Channel of Influenza A Virus Suggests a Mechanism for its Ion Selectivity. Proc. Natl. Acad. Sci. 94, 11301–11306. doi:10.1073/pnas.94.21.11301
Rey-Carrizo, M., Barniol-Xicota, M., Ma, C., Frigolé-Vivas, M., Torres, E., Naesens, L., et al. (2014). Easily Accessible Polycyclic Amines that Inhibit the Wild-type and Amantadine-Resistant Mutants of the M2 Channel of Influenza A Virus. J. Med. Chem. 57, 5738–5747. doi:10.1021/jm5005804
Rey-Carrizo, M., Gazzarrini, S., Llabrés, S., Frigolé-Vivas, M., Juárez-JiménezFont-Bardia, J. M., Font-Bardia, M., et al. (2015). New Polycyclic Dual Inhibitors of the Wild Type and the V27A Mutant M2 Channel of the Influenza A Virus with Unexpected Binding Mode. Eur. J. Med. Chem. 96, 318–329. doi:10.1016/j.ejmech.2015.04.030
Rosenberg, M. R., and Casarotto, M. G. (2010). Coexistence of Two Adamantane Binding Sites in the Influenza A M2 Ion Channel. Proc. Natl. Acad. Sci. 107, 13866–13871. doi:10.1073/pnas.1002051107
Rossman, J. S., Jing, X., Leser, G. P., and Lamb, R. A. (2010). Influenza Virus M2 Protein Mediates ESCRT-independent Membrane Scission. Cell 142, 902–913. doi:10.1016/j.cell.2010.08.029
Samson, M. S., and Kerr, I. D. (1993). Influenza Virus M2 Protein: A Molecular Modelling Study of the Ion Channel. Protein Eng. 6, 65–74.
Sansom, M. S. P., Kerr, I. D., Smith, G. R., and Son, H. S. (1997). The Influenza A Virus M2 Channel: A Molecular Modeling and Simulation Study. Virology 233, 163–173. doi:10.1006/viro.1997.8578
Saunders-Hastings, P., and Krewski, D. (2016). Reviewing the History of Pandemic Influenza: Understanding Patterns of Emergence and Transmission. Pathogens 5, 66. doi:10.3390/pathogens5040066
Schnell, J. R., and Chou, J. J. (2008). Structure and Mechanism of the M2 Proton Channel of Influenza A Virus. Nature 451, 591–595. doi:10.1038/nature06531
Sharma, M., Yi, M., Dong, H., Qin, H., Peterson, E., Busath, D. D., et al. (2010). Insight into the Mechanism of the Influenza A Proton Channel from a Structure in a Lipid Bilayer. Science 330, 509–512. doi:10.1126/science.1191750
Shimbo, K., Brassard, D. L., Lamb, R. A., and Pinto, L. H. (1996). Ion Selectivity and Activation of the M2 Ion Channel of Influenza Virus. Biophys. J. 70, 1335–1346. doi:10.1016/S0006-3495(96)79690-X
Sidorenko, Y., and Reichl, U. (2004). Structured Model of Influenza Virus Replication in MDCK Cells. Biotechnol. Bioeng. 88, 1–14. doi:10.1002/bit.20096
Stouffer, A. L., Acharya, R., Salom, D., Levine, A. S., Di Costanzo, L., Soto, C. S., et al. (2008). Structural Basis for the Function and Inhibition of an Influenza Virus Proton Channel. Nature 451, 596–599. doi:10.1038/nature06528
Tan, R. C., Truong, T. N., McCammon, J. A., and Sussman, J. L. (1993). Acetylcholinesterase: Electrostatic Steering Increases the Rate of Ligand Binding. Biochemistry 32, 401–403. doi:10.1021/bi00053a003
Tang, Y., Zaitseva, F., Lamb, R. A., and Pinto, L. H. (2002). The Gate of the Influenza Virus M2 Proton Channel Is Formed by a Single Tryptophan Residue. J. Biol. Chem. 277, 39880–39886. doi:10.1074/jbc.m206582200
Thomaston, J. L., and DeGrado, W. F. (2016). Crystal Structure of the Drug-Resistant S31N Influenza M2 Proton Channel. Protein Sci. 25, 1551–1554. doi:10.1002/pro.2937
Thomaston, J. L., Alfonso-Prieto, M., Woldeyes, R. A., Fraser, J. S., Klein, M. L., Fiorin, G., et al. (2015). High-resolution Structures of the M2 Channel from Influenza A Virus Reveal Dynamic Pathways for Proton Stabilization and Transduction. Proc. Natl. Acad. Sci. USA 112, 14260–14265. doi:10.1073/pnas.1518493112
Thomaston, J. L., Polizzi, N. F., Konstantinidi, A., Wang, J., Kolocouris, A., and DeGrado, W. F. (2018). Inhibitors of the M2 Proton Channel Engage and Disrupt Transmembrane Networks of Hydrogen-Bonded Waters. J. Am. Chem. Soc. 140, 15219–15226. doi:10.1021/jacs.8b06741
Thomaston, J. L., Konstantinidi, A., Liu, L., Lambrinidis, G., Tan, J., Caffrey, M., et al. (2020). X-ray crystal Structures of the Influenza M2 Proton Channel Drug-Resistant V27A Mutant Bound to a spiro-adamantyl Amine Inhibitor Reveal the Mechanism of Adamantane Resistance. Biochemistry 59, 627–634. doi:10.1021/acs.biochem.9b00971
Thomaston, J. L., Samways, M. L., Konstantinidi, A., Ma, C., Hu, Y., Bruce Macdonald, H. E., et al. (2021). Rimantadine Binds to and Inhibits the Influenza A M2 Proton Channel without Enantiomeric Specificity. Biochemistry 60, 2471–2482. doi:10.1021/acs.biochem.1c00437
Torabifard, H., Panahi, A., and Brooks, C. L. (2020). M2 Amphipathic Helices Facilitate pH-dependent Conformational Transition in Influenza A Virus. Proc. Natl. Acad. Sci. USA 117, 3583–3591. doi:10.1073/pnas.1913385117
Tzitzoglaki, C., McGuire, K., Lagarias, P., Konstantinidi, A., Hoffmann, A., Fokina, N. A., et al. (2020). Chemical Probes for Blocking of Influenza A M2 Wild-type and S31N Channels. ACS Chem. Biol. 15, 2331–2337. doi:10.1021/acschembio.0c00553
Vaissier Welborn, V., and Head-Gordon, T. (2019). Computational Design of Synthetic Enzymes. Chem. Rev. 119, 6613–6630. doi:10.1021/acs.chemrev.8b00399
Venkataraman, P., Lamb, R. A., and Pinto, L. H. (2005). Chemical rescue of Histidine Selectivity Filter Mutants of the M2 Ion Channel of Influenza A Virus. J. Biol. Chem. 280, 21463–21472. doi:10.1074/jbc.m412406200
Wang, T., and Hong, M. (2015). Investigation of the Curvature Induction and Membrane Localization of the Influenza Virus M2 Protein Using Static and Off-Magic-Angle Spinning Solid-State Nuclear Magnetic Resonance of Oriented Bicelles. Biochemistry 54, 2214–2226. doi:10.1021/acs.biochem.5b00127
Wang, C., Takeuchi, K., Pinto, L. H., and Lamb, R. A. (1993). Ion Channel Activity of Influenza A Virus M2 Protein: Characterization of the Amantadine Block. J. Virol. 67, 5585–5594. doi:10.1128/jvi.67.9.5585-5594.1993
Wang, C., Lamb, R. A., and Pinto, L. H. (1995). Activation of the M2 Ion Channel of Influenza Virus: A Role for the Transmembrane Domain Histidine Residue. Biophys. J. 69, 1363–1371. doi:10.1016/s0006-3495(95)80003-2
Wang, J., Ma, C., Fiorin, G., Carnevale, V., Wang, T., Hu, F., et al. (2011). Molecular Dynamics Simulation Directed Rational Design of Inhibitors Targeting Drug-Resistant Mutants of Influenza A Virus M2. J. Am. Chem. Soc. 133, 12834–12841. doi:10.1021/ja204969m
Wang, J., Wu, Y., Ma, C., Fiorin, G., Wang, J., Pinto, L. H., et al. (2013). Structure and Inhibition of the Drug-Resistant S31N Mutant of the M2 Ion Channel of Influenza A Virus. Proc. Natl. Acad. Sci. USA 110, 1315–1320. doi:10.1073/pnas.1216526110
Wang, J., Li, F., and Ma, C. (2015). Recent Progress in Designing Inhibitors that Target the Drug-Resistant M2 Proton Channels from the Influenza A Viruses. Biopolymers 104, 291–309. doi:10.1002/bip.22623
Watkins, L. C., DeGrado, W. F., and Voth, G. A. (2020). Influenza A M2 Inhibitor Binding Understood through Mechanisms of Excess Proton Stabilization and Channel Dynamics. J. Am. Chem. Soc. 142, 17425–17433. doi:10.1021/jacs.0c06419
Watson, S. J., Langat, P., Reid, S. M., Lam, T. T.-Y., Cotten, M., Kelly, M., et al. (2015). Molecular Epidemiology and Evolution of Influenza Viruses Circulating within European Swine between 2009 and 2013. J. Virol. 89, 9920–9931. doi:10.1128/jvi.00840-15
Wei, C., and Pohorille, A. (2013). Activation and Proton Transport Mechanism in Influenza A M2 Channel. Biophysical J. 105, 2036–2045. doi:10.1016/j.bpj.2013.08.030
World Health Organization (2019). Global Influenza Strategy 2019-2030. Geneva: World Health Organization. Available at: https://apps.who.int/iris/handle/10665/311184 (Accessed September 14, 2021).
Wu, Y., Canturk, B., Jo, H., Ma, C., Gianti, E., Klein, M. L., et al. (2014). Flipping in the Pore: Discovery of Dual Inhibitors that Bind in Different Orientations to the Wild-type versus the Amantadine-Resistant S31N Mutant of the Influenza A Virus M2 Proton Channel. J. Am. Chem. Soc. 136, 17987–17995. doi:10.1021/ja508461m
Wu, X., Wu, X., Sun, Q., Zhang, C., Yang, S., Li, L., et al. (2017). Progress of Small Molecular Inhibitors in the Development of Anti-influenza Virus Agents. Theranostics 7, 826–845. doi:10.7150/thno.17071
Yi, M., Cross, T. A., and Zhou, H.-X. (2008). A Secondary Gate as a Mechanism for Inhibition of the M2 Proton Channel by Amantadine. J. Phys. Chem. B 112, 7977–7979. doi:10.1021/jp800171m
Keywords: M2 proton channel, influenza A virus, proton transport, inhibition, molecular mechanism, drug design
Citation: Aledavood E, Selmi B, Estarellas C, Masetti M and Luque FJ (2022) From Acid Activation Mechanisms of Proton Conduction to Design of Inhibitors of the M2 Proton Channel of Influenza A Virus. Front. Mol. Biosci. 8:796229. doi: 10.3389/fmolb.2021.796229
Received: 16 October 2021; Accepted: 24 December 2021;
Published: 14 January 2022.
Edited by:
Carla Molteni, King’s College London, United KingdomReviewed by:
Mercedes Alfonso-Prieto, Helmholtz Association of German Research Centres (HZ), GermanySuman Chakrabarty, S. N. Bose National Centre for Basic Sciences, India
Giacomo Fiorin, Temple University, United States
Copyright © 2022 Aledavood, Selmi, Estarellas, Masetti and Luque. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Carolina Estarellas, cestarellas@ub.edu; Matteo Masetti, matteo.masetti4@unibo.it; F. Javier Luque, fjluque@ub.edu
†These authors have contributed equally to this work and share first authorship