Flu channel drug resistance: a tale of two sites

The M2 proteins of influenza A and B virus, AM2 and BM2, respectively, are transmembrane proteins that oligomerize in the viral membrane to form proton-selective channels. Proton conductance of the M2 proteins is required for viral replication; it is believed to equilibrate pH across the viral membrane during cell entry and across the trans-Golgi membrane of infected cells during viral maturation. In addition to the role of M2 in proton conductance, recent mutagenesis and structural studies suggest that the cytoplasmic domains of the M2 proteins also play a role in recruiting the matrix proteins to the cell surface during virus budding. As viral ion channels of minimalist architecture, the membrane-embedded channel domain of M2 has been a model system for investigating the mechanism of proton conduction. Moreover, as a proven drug target for the treatment of influenza A infection, M2 has been the subject of intense research for developing new anti-flu therapeutics. AM2 is the target of two anti-influenza A drugs, amantadine and rimantadine, both belonging to the adamantane class of compounds. However, resistance of influenza A to adamantane is now widespread due to mutations in the channel domain of AM2. This review summarizes the structure and function of both AM2 and BM2 channels, the mechanism of drug inhibition and drug resistance of AM2, as well as the development of new M2 inhibitors as potential anti-flu drugs.     

The oligomeric state of the active BM2 ion channel protein of influenza B virus

Influenza A virus and influenza B virus particles both contain small integral membrane proteins (A/M2 and BM2, respectively) that function as a pH-sensitive proton channel and are essential for virus replication. The mechanism of action of the M2 channels is a subject of scientific interest particularly as A/M2 channel was shown to be a target for the action of the antiviral drug amantadine. Unfortunately, an inhibitor of the BM2 channel activity is not known. Thus, knowledge of the structural and functional properties of the BM2 channel is essential for the development of potent antiviral drugs. The characterization of the oligomeric state of the BM2 channel is an essential first step in the understanding of channel function. Here we describe determination of the stoichiometry of the BM2 proton channel by utilizing three different approaches. 1) We demonstrated that BM2 monomers can be chemically cross-linked to yield species consistent with dimers, trimers, and tetramers. 2) We studied electrophysiological and biochemical properties of mixed oligomers consisting of wild-type and mutated BM2 subunits and related these data to predicted binomial distribution models. 3) We used fluorescence resonance energy transfer (FRET) in combination with biochemical measurements to estimate the relationships between BM2 channel subunits expressed in the plasma membrane. Our experimental data are consistent with a tetrameric structure of the BM2 channel. Finally, we demonstrated that BM2 transmembrane domain is responsible for the channel oligomerization.     

Designing Inhibitors of M2 Proton Channel against H1N1 Swine Influenza Virus

Background

M2 proton channel of H1N1 influenza A virus is the target protein of anti-flu drugs amantadine and rimantadine. However, the two once powerful adamantane-based drugs lost their 90% bioactivity because of mutations of virus in recent twenty years. The NMR structure of the M2 channel protein determined by Schnell and Chou (Nature, 2008, 451, 591–595) may help people to solve the drug-resistant problem and develop more powerful new drugs against H1N1 influenza virus.

Methodology

Docking calculation is performed to build the complex structure between receptor M2 proton channel and ligands, including existing drugs amantadine and rimantadine, and two newly designed inhibitors. The computer-aided drug design methods are used to calculate the binding free energies, with the computational biology techniques to analyze the interactions between M2 proton channel and adamantine-based inhibitors.

Conclusions

1) The NMR structure of M2 proton channel provides a reliable structural basis for rational drug design against influenza virus. 2) The channel gating mechanism and the inhibiting mechanism of M2 proton channel, revealed by the NMR structure of M2 proton channel, provides the new ideas for channel inhibitor design. 3) The newly designed adamantane-based inhibitors based on the modeled structure of H1N1-M2 proton channel have two pharmacophore groups, which act like a “barrel hoop”, holding two adjacent helices of the H1N1-M2 tetramer through the two pharmacophore groups outside the channel. 4) The inhibitors with such binding mechanism may overcome the drug resistance problem of influenza A virus to the adamantane-based drugs.     

An in-depth analysis of the biological functional studies based on the NMR M2 channel structure of influenza A virus

The long-sought three-dimensional structure of the M2 proton channel of influenza A virus was successfully determined recently by the high-resolution NMR [J.R. Schnell, J.J. Chou, Structure and mechanism of the M2 proton channel of influenza A virus, Nature 451 (2008) 591-595]. Such a milestone work has provided a solid structural basis for studying drug-resistance problems. However, the action mechanism revealed from the NMR structure is completely different from the traditional view and hence prone to be misinterpreted as “conflicting” with some previous biological functional studies. To clarify this kind of confusion, an in-depth analysis was performed for these functional studies, particularly for the mutations D44N, D44A and N44D on position 44, and the mutations on positions 27-38. The analyzed results have provided not only compelling evidences to further validate the NMR structure but also very useful clues for dealing with the drug-resistance problems and developing new effective drugs against H5N1 avian influenza virus, an impending threat to human beings.     

Kir6.2 Channel Gating by Intracellular Protons: Subunit Stoichiometry for Ligand Binding and Channel Gating

The adenosine triphosphate-sensitive K⁺ (K_ATP) channels are gated by several metabolites, whereas the gating mechanism remains unclear. Kir6.2, a pore-forming subunit of the K_ATP channels, has all machineries for ligand binding and channel gating. In Kir6.2, His175 is the protonation site and Thr71 and Cys166 are involved in channel gating. Here, we show how individual subunits act in proton binding and channel gating by selectively disrupting functional subunits using these residues. All homomeric dimers and tetramers showed pH sensitivity similar to the monomeric channels. Concatenated construction of wild type with disrupted subunits revealed that none of these residues had a dominant-negative effect on the proton-dependent channel gating. Subunit action in proton binding was almost identical to that for channel gating involving Cys166, suggesting a one-to-one coupling from the C terminus to the M2 helix. This was significantly different from the effect of T71Y heteromultimers, suggesting distinct contributions of M1 and M2 helices to channel gating. Subunits underwent concerted rather than independent action. Two wild-type subunits appeared to act as a functional dimer in both cis and trans configurations. The understanding of K_ATP channel gating by intracellular pH has a profound impact on cellular responses to metabolic stress as a significant drop in intracellular pH is more frequently seen under a number of physiological and pathophysiological conditions than a sole decrease in intracellular ATP levels. Runping Wang, Junda Su contributed equally to this work.     

Design of an Efficient Inhibitor for the Influenza A Virus M2 Ion Channel

Influenza A virus is capable of rapidly infecting large human populations, warranting the development of novel drugs to efficiently inhibit virus replication. A transmembrane ion channel formed by the M2 protein plays an important role in influenza virus replication. A reasonable approach to designing an effective antivirus drug is constructing a molecule that binds in the M2 transmembrane proton channel, blocks H+ proton diffusion through the channel, and thus the influenza A virus cycle. The known anti-influenza drugs amantadine and rimantadine have a weak effect on influenza A virus replication. A new class of positively charged molecules, diazabicyclooctane derivatives with a constant charge of +2, was proposed to block proton diffusion through the M2 ion channel. Molecular dynamics simulations were performed to study the temperature fluctuations in the M2 structure, and ionization states of histidine residues were established at physiological pH values. Two types of diazabicyclooctane derivatives were analyzed for binding with the M2 ion channel. An optimal structure was determined for a blocker to most efficiently bind with the M2 ion channel and block proton diffusion. The new molecule is advantageous over amantadine and rimantadine in having a positive charge of +2, which creates a positive electrostatic potential barrier to proton transport through the M2 ion channel in addition to a steric barrier.     

Conformational plasticity of the influenza A M2 transmembrane helix in lipid bilayers under varying pH, drug binding, and membrane thickness

Membrane proteins change their conformations to respond to environmental cues, thus conformational plasticity is important for function. The influenza A M2 protein forms an acid-activated proton channel important for the virus lifecycle. Here we have used solid-state NMR spectroscopy to examine the conformational plasticity of membrane-bound transmembrane domain of M2 (M2TM). (13)C and (15)N chemical shifts indicate coupled conformational changes of several pore-facing residues due to changes in bilayer thickness, drug binding, and pH. The structural changes are attributed to the formation of a well-defined helical kink at G34 in the drug-bound state and in thick lipid bilayers, nonideal backbone conformation of the secondary-gate residue V27 in the presence of drug, and nonideal conformation of the proton-sensing residue H37 at high pH. The chemical shifts constrained the (?, ψ) torsion angles for three "basis" states, the equilibrium among which explains the multiple resonances per site in the NMR spectra under different combinations of bilayer thickness, drug binding, and pH conditions. Thus, conformational plasticity is important for the proton conduction and inhibition of M2TM. The study illustrates the utility of NMR chemical shifts for probing the structural plasticity and folding of membrane proteins.     

Toxicity of Influenza A Virus Matrix Protein 2 for Mammalian Cells is Associated with its Intrinsic Proton-Channeling Activity

Molecules of influenza matrix protein 2 (M2) are organized in tetramers that constitute a well-conserved virion component and also form proton channels in the plasma membrane of infected cells. In this report we demonstrate that influenza M2 protein is cytopathic in vitro for mammalian cells. An M2 point-mutant (M2pm) protein was constructed that contained amino acid changes designed to block the proton channel via introduction of large hydrophobic residues. This mutant was significantly less toxic upon transient transfection in vitro than the wild-type M2 (M2wt). To assess the possible correlation between M2 cytotoxicity and its proton channel activity, we monitored changes in mitochondria membrane potential induced by M2wt and M2pm. M2wt rapidly decreased mitochondria membrane potential reflecting the transmembrane proton gradient, while M2pm was markedly less efficient. Thus, M2 is cytotoxic for mammalian cells, likely via its proton channel activity and may therefore contribute to influenza pathogenesis through this previously unknown mechanism.     

Computational Investigation of Drug-Resistant Mutant of M2 Proton Channel (S31N) Against Rimantadine

M2 proton channel is the target for treating the patients who ere suffering from influenza A infection, which facilitates the spread of virions. Amantadine and rimantadine are adamantadine-based drugs, which target M2 proton channel and inhibit the viral replication. Preferably, rimantadine drug is used more than amantadine because of its fewer side effects. However, S31N mutation in the M2 proton channel was highly resistant to the rimantadine drug. Therefore, in the present study, we focused to understand the drug-resistance mechanism of S31N mutation with the aid of molecular docking and dynamics approach. The docking analysis undoubtedly indicates that affinity for rimantadine with mutant-type M2 proton channel is significantly lesser than the native-type M2 proton channel. In addition, RMSD, RMSF, and principal component analysis suggested that the mutation shows increased flexibility. Furthermore, the intermolecular hydrogen bonds analysis showed that there is a complete loss of hydrogen bonds in the mutant complex. On the whole, we conclude that the intermolecular contact was maintained by D-44, a key residue for stable binding of rimantadine. These findings are certainly helpful for better understanding of drug-resistance mechanism and also helpful for designing new drugs for treating influenza infection against drug-resistance target.     

The mode of interaction of the relaxin-like factor (RLF) with the leucine-rich repeat G protein-activated receptor 8

The relaxin-like factor (RLF, also named INSL3) is a critical component in the chain of events that lead to the normal positioning of the gonads in the male fetus. RLF and relaxin share features of the secondary structure to the extent that relaxin cross-reacts with the LGR8, the RLF receptor. Although both hormones interact with their receptors essentially via the B chain, the sharply defined binding cassette of relaxin is not present in RLF. Structure and function analysis of RLF derivatives with single amino acid replacements revealed that the most important binding residues are tryptophan B27, followed by arginine B16 and valine B19. Single alanine replacements for each individual position resulted in a relative receptor affinity of 4.0% (B16), 6.1% (B19), and 0.5% (B27). Tryptophan B27 is located on an extended structure, and arginine B16 and valine B19 are positioned on the exposed surface of the B chain helix. The 3 residues could be brought together to form a contiguous binding area if the C-terminal end of the B chain were free to fold back against the central portion of the B chain helix. Such a movement depends critically on the flexibility of the C-terminal end, which is controlled by positions B23-25. In as much as these major binding residues seem hardly sufficient to explain the strong binding of RLF to LGR8 we searched for and found an extended region where little contributions by individual residues added up to a strong receptor affinity. This mode of interaction could drive the binding energy sufficiently high to account for the picomolar binding constant of RLF and its receptor.     

Mechanistic Insight into the H<Subscript>2</Subscript>O/D<Subscript>2</Subscript>O Isotope Effect in the Proton Transport of the Influenza Virus M2 Protein

The M2 proton channel is essential for the replication of the flu virus and is a known drug target. The functional mechanism of channel activation and conductance is key to both the basic biology of viral replication and the design of drugs that can withstand mutations. A quantitative model was previously developed for calculating the rate of proton transport through the M2 channel. The permeant proton was assumed to diffuse to the pore, obligatorily bind to the His37 tetrad, and then dissociate and be released to either side of the tetrad. Here the model is used to calculate the effect of a change in solvent from H₂O to D₂O on the rate of proton transport. The solvent substitution affects two parameters in the model: the proton diffusion constant and the pK _a for proton binding to the His37 tetrad. When the known effects on these two parameters are included, the deuterium isotope effect calculated from the model is in quantitatively agreement with experimental results. This strict test of the theoretical model provides strong support for the hypothesis that the permeant proton obligatorily binds to and then unbinds from the His37 tetrad. This putatively essential role of the His37 tetrad in the functional mechanism of the M2 channel makes it a promising target for designing mutation-tolerant drugs.     

Site-directed glycosylation tagging of functional Kir2.1 reveals that the putative pore-forming segment is extracellular

Inwardly rectifying K+ channels or Kirs are a large gene family and have been predicted to have two transmembrane segments, M1 and M2, intracellular N and C termini, and two extracellular loops, E1 and E2, separated by an intramembranous pore-forming segment, H5. H5 contains a stretch of eight residues that are similar in voltage-dependent K+ channels, Kvs, and this stretch is called the signature sequence of K+ channels. Because mutations in this sequence altered selectivity in Kvs, it has been designated as the selectivity filter. Previously, we used N-glycosylation substitution mutants to map the extracellular topology of a weak inwardly rectifying K+ channel, Kir1.1 or ROMK1, and found that the entire H5 segment was extracellular. We now report utilization of introduced N-glycosylation sites, NX(S/T), at positions Ser(128) in E1, and Gln(140), Ileu(143), and Phe(147) in the H5 sequence of a strong inwardly rectifying K+ channel, Kir2.1. Furthermore, we show that biotinylated channel proteins with N-linked oligosaccharides attached at positions 140 and 143 in the signature sequence are located at the cell surface. Mutant channels were functional as detected by whole-cell and single-channel recordings. Unlike Kir1.1, position Lys(117) was not occupied. We conclude that, for yet another K+ channel, the invariant G(Y/F)G sequence is extracellular rather than intramembranous.     

Independent Contribution of Extracellular Proton Binding Sites to ASIC1a Activation

Acid-sensing ion channels (ASICs) are a group of trimeric cation permeable channels gated by extracellular protons that are mainly expressed in the nervous system. Despite the structural information available for ASIC1, there is limited understanding of the molecular mechanism that allows these channels to sense and respond to drops in extracellular pH. In this report, we employed the substituted cysteine accessibility method and site-directed mutagenesis to examine the mechanism of activation of ASIC1a by extracellular protons. We found that the modification of E238C and D345C channels by MTSET reduced proton apparent affinity for activation. Furthermore, the introduction of positively charged residues at position 345 rendered shifted biphasic proton activation curves. Likewise, channels bearing mutations at positions 79 and 416 in the palm domain of the channel showed reduced proton apparent affinity and biphasic proton activation curves. Of significance, the effect of the mutations at positions 79 and 345 on channel activation was additive. E79K-D345K required a change to a pH lower than 2 for maximal activation. In summary, this study provides direct evidence for the presence of two distinct proton coordination sites in the extracellular region of ASIC1a, which jointly facilitate pore opening in response to extracellular acidification.     

Multiscale simulation of transmembrane proteins

Multiscale simulation is employed to examine changes in atomistic-level protein structure due to long wavelength membrane undulations and plane stress fields. An ensemble of atomistic-level simulations of a model of a transmembrane influenza A virus M2 proton channel in a dimyristoylphosphatidylcholine (DMPC) bilayer is coupled to a corresponding mesoscopic model of a DMPC bilayer in an explicit mesoscopic solvent. Structural variations in the key proton gating His37 residues of the M2 channel are examined. Small, but distinct variations in the structure of the His37 residues are observed in both the open and closed states of the channel as a result of the coupling to mesoscopic-level membrane motions.     

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.     

Roles of the histidine and tryptophan side chains in the M2 proton channel from influenza A virus

The M2 protein form influenza A virus forms a tetrameric ion channel, which enables proton passage across biological membranes when the N-terminal side is acidified. Among the amino acid residues in the transmembrane domain of the M2 protein, His37 and Trp41 are essential for the pH-regulated proton conductance. Current knowledge about the structures and interactions of His37 and Trp41 suggests a model for the M2 ion channel, in which the channel is closed by a network of His37 hydrogen bonds at neutral pH and is opened by a His37-Trp41 cation-pi interaction at acidic pH.     

Structural Influences: Cholesterol,Drug, and Proton Binding to Full-Length Influenza A M2 Protein

The structure and functions of the M2 protein from Influenza A are sensitive to pH, cholesterol, and the antiinfluenza drug Amantadine. This is a tetrameric membrane protein of 97 amino-acid residues that has multiple functions, among them as a proton-selective channel and facilitator of viral budding, replacing the need for the ESCRT proteins that other viruses utilize. Here, various amino-acid-specific-labeled samples of the full-length protein were prepared and mixed, so that only interresidue ¹³C-¹³C cross peaks between two differently labeled proteins representing interhelical interactions are observed. This channel is activated at slightly acidic pH values in the endosome when the His³⁷ residues in the middle of the transmembrane domain take on a +2 or +3 charged state. Changes observed here in interhelical distances in the N-terminus can be accounted for by modest structural changes, and no significant changes in structure were detected in the C-terminal portion of the channel upon activation of the channel. Amantadine, which blocks proton conductance by binding in the aqueous pore near the N-terminus, however, significantly modifies the tetrameric structure on the opposite side of the membrane. The interactions between the juxtamembrane amphipathic helix of one monomer and its neighboring monomer observed in the absence of drug are disrupted in its presence. However, the addition of cholesterol prevents this structural disruption. In fact, strong interactions are observed between cholesterol and residues in the amphipathic helix, accounting for cholesterol binding adjacent to a native palmitoylation site and near to an interhelix crevice that is typical of cholesterol binding sites. The resultant stabilization of the amphipathic helix deep in the bilayer interface facilitates the bilayer curvature that is essential for viral budding.     

Crystal structure of the drug‐resistant S31N influenza M2 proton channel

The M2 protein is a small proton channel found in the influenza A virus that is necessary for viral replication. The M2 channel is the target of a class of drugs called the adamantanes, which block the channel pore and prevent the virus from replicating. In recent decades mutations have arisen in M2 that prevent the adamantanes from binding to the channel pore, with the most prevalent of these mutations being S31N. Here we report the first crystal structure of the S31N mutant crystallized using lipidic cubic phase crystallization techniques and solved to 1.59 Å resolution. The Asn31 residues point directly into the center of the channel pore and form a hydrogen‐bonded network that disrupts the drug‐binding site. Ordered waters in the channel pore form a continuous hydrogen bonding network from Gly34 to His37.     

Divergent CFTR orthologs respond differently to the channel inhibitors CFTRinh-172, glibenclamide, and GlyH-101

Comparison of diverse orthologs is a powerful tool to study the structure and function of channel proteins. We investigated the response of human, killifish, pig, and shark cystic fibrosis transmembrane conductance regulator (CFTR) to specific inhibitors of the channel: CFTR(inh)-172, glibenclamide, and GlyH-101. In three systems, including organ perfusion of the shark rectal gland, primary cultures of shark rectal gland tubules, and expression studies of each ortholog in cRNA microinjected Xenopus laevis oocytes, we observed fundamental differences in the sensitivity to inhibition by these channel blockers. In organ perfusion studies, shark CFTR was insensitive to inhibition by CFTR(inh)-172. This insensitivity was also seen in short-circuit current experiments with cultured rectal gland tubular epithelial cells (maximum inhibition 4 ± 1.3%). In oocyte expression studies, shark CFTR was again insensitive to CFTR(inh)-172 (maximum inhibition 10.3 ± 2.5% at 25 μM), pig CFTR was insensitive to glibenclamide (maximum inhibition 18.4 ± 4.4% at 250 μM), and all orthologs were sensitive to GlyH-101. The amino acid residues considered responsible by previous site-directed mutagenesis for binding of the three inhibitors are conserved in the four CFTR isoforms studied. These experiments demonstrate a profound difference in the sensitivity of different orthologs of CFTR proteins to inhibition by CFTR blockers that cannot be explained by mutagenesis of single amino acids. We believe that the potency of the inhibitors CFTR(inh)-172, glibenclamide, and GlyH-101 on the CFTR chloride channel protein is likely dictated by the local environment and the three-dimensional structure of additional residues that form the vestibules, the chloride pore, and regulatory regions of the channel.     

Azide reduces the hydrophobic barrier of the bacteriorhodopsin proton channel

The sensitivity of a nitroxide spin label to the polarity of its environment has been used to estimate the hydrophobic barrier of the proton channel of the transmembrane proton pump bacteriorhodopsin. By means of site-specific mutagenesis, single cysteine residues were introduced at 10 positions located at the protein surface, in the protein interior, and along the proton pathway. After reaction with a methanethiosulfonate spin label, the principle values of the hyperfine tensor A and the g-tensor were determined from electron paramagnetic resonance spectra measured at 170 K. The shape of the hydrophobic barrier of the proton channel is characterized in terms of a polarity index, DeltaA, determined from the variation of the hyperfine coupling constant Azz. The maximum of the hydrophobic barrier is found to be close to the retinal chromophore in the proton uptake pathway. The effect of the asymmetric distribution of charged and polar residues in the proton release and uptake pathways is clearly reflected in the behavior of the hydrophobic barrier. The presence of azide reduces the barrier height of both the cytoplasmic and extracellular channels. This finding supports the view of azide and other weakly acidic anions as catalysts for the formation of hydrogen-bonded networks in proton pathways of proteins.