Interaction of aminoadamantane derivatives with the influenza A virus M2 channel-Docking using a pore blocking model

Interaction of aminoadamantanes with the influenza A virus M2 proton channel was analyzed by docking simulations of a series of synthetic aminoadamantane derivatives, of differing binding affinity, into the crystal structure of the transmembrane (M2TM) pore. The pore blocking model tested in the ‘gas phase’ describes qualitatively the changes on the relative binding affinities of the compounds (although a series of highly hydrophobic ligands which seem to have little capacity for different specific interactions with their receptor). The docking calculations predicted poses in which the adamantane ring is surrounded mainly by the alkyl side chains of Val27 or Ala30 and the ligand’s amino group is generally hydrogen bonded with hydroxyls of Ser31 or carbonyls of Val27 or carbonyls of Ala30, the former (Ser31) being the most stable and most frequently observed. The binding of the ligand is a compromise between hydrogen bonding ability, which is elevated by a primary amino group, and apolar interactions, which are increased by the ability of the lipophilic moiety to adequately fill a hydrophobic pocket within the M2TM pore. A delicate balance of these hydrophobic contributions is required for optimal interaction.     

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.     

Structural basis for proton conduction and inhibition by the influenza M2 protein

The influenza M2 protein forms an acid‐activated and drug‐sensitive proton channel in the virus envelope that is important for the virus lifecycle. The functional properties and high‐resolution structures of this proton channel have been extensively studied to understand the mechanisms of proton conduction and drug inhibition. We review biochemical and electrophysiological studies of M2 and discuss how high‐resolution structures have transformed our understanding of this proton channel. Comparison of structures obtained in different membrane‐mimetic solvents and under different pH using X‐ray crystallography, solution NMR, and solid‐state NMR spectroscopy revealed how the M2 structure depends on the environment and showed that the pharmacologically relevant drug‐binding site lies in the transmembrane (TM) pore. Competing models of proton conduction have been evaluated using biochemical experiments, high‐resolution structural methods, and computational modeling. These results are converging to a model in which a histidine residue in the TM domain mediates proton relay with water, aided by microsecond conformational dynamics of the imidazole ring. These mechanistic insights are guiding the design of new inhibitors that target drug‐resistant M2 variants and may be relevant for other proton channels.     

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.     

Detection of drug‐induced conformational change of a transmembrane protein in lipid bilayers using site‐directed spin labeling

As a target of antiviral drugs, the influenza A M2 protein has been the focus of numerous structural studies and has been extensively explored as a model ion channel. In this study, we capitalize on the expanding body of high‐resolution structural data available for the M2 protein to design and interpret site‐directed spin‐labeling electron paramagnetic resonance spectroscopy experiments on drug‐induced conformational changes of the M2 protein embedded in lipid bilayers. We obtained data in the presence of adamantane drugs for two different M2 constructs (M2TM 22–46 and M2TMC 23–60). M2TM peptides were spin labeled at the N‐terminal end of the transmembrane domain. M2TMC peptides were spin labeled site specifically at cysteine residues substituted for amino acids within the transmembrane domain (L36, I39, I42, and L43) and the C‐terminal amphipathic helix (L46, F47, F48, C50, I51, Y52, R53, F54, F55, and E56). Addition of adamantane drugs brought about significant changes in measured electron paramagnetic resonance spectroscopy environmental parameters consistent with narrowing of the transmembrane channel pore and closer packing of the C‐terminal amphipathic helices.     

Solution NMR structure of the V27A drug resistant mutant of influenza A M2 channel

The M2 protein of influenza A virus forms a proton-selective channel that is required for viral replication. It is the target of the anti-influenza drugs, amantadine and rimantadine. Widespread drug resistant mutants, however, has greatly compromised the effectiveness of these drugs. Here, we report the solution NMR structure of the highly pathogenic, drug resistant mutant V27A. The structure reveals subtle structural differences from wildtype that maybe linked to drug resistance. The V27A mutation significantly decreases hydrophobic packing between the N-terminal ends of the transmembrane helices, which explains the looser, more dynamic tetrameric assembly. The weakened channel assembly can resist drug binding either by destabilizing the rimantadine-binding pocket at Asp44, in the case of the allosteric inhibition model, or by reducing hydrophobic contacts with amantadine in the pore, in the case of the pore-blocking model. Moreover, the V27A structure shows a substantially increased channel opening at the N-terminal end, which may explain the faster proton conduction observed for this mutant. Furthermore, due to the high quality NMR data recorded for the V27A mutant, we were able to determine the structured region connecting the channel domain to the C-terminal amphipathic helices that was not determined in the wildtype structure. The new structural data show that the amphipathic helices are packed much more closely to the channel domain and provide new insights into the proton transfer pathway.     

Cu(II) inhibition of the proton translocation machinery of the influenza A virus M2 protein

The homotetrameric M2 integral membrane protein of influenza virus forms a proton-selective ion channel. An essential histidine residue (His-37) in the M2 transmembrane domain is believed to play an important role in the conduction mechanism of this channel. Also, this residue is believed to form hydrogen-bonded interactions with the ammonium group of the anti-viral compound, amantadine. A molecular model of this channel suggests that the imidazole side chains of His-37 from symmetry-related monomers of the homotetrameric pore converge to form a coordination site for transition metals. Thus, membrane currents of oocytes of Xenopus laevis expressing the M2 protein were recorded when the solution bathing the oocytes contained various transition metals. Membrane currents were strongly and reversibly inhibited by Cu2+ with biphasic reaction kinetics. The biphasic inhibition curves may be explained by a two-site model involving a fast-binding peripheral site with low specificity for divalent metal ions, as well as a high affinity site (Kdiss approximately 2 microM) that lies deep within the pore and shows rather slow-binding kinetics (kon = 18.6 +/- 0.9 M-1 s-1). The pH dependence of the interaction with the high affinity Cu2+-binding site parallels the pH dependence of inhibition by amantadine, which has previously been ascribed to protonation of His-37. The voltage dependence of the inhibition at the high affinity site indicates that the binding site lies within the transmembrane region of the pore. Furthermore, the inhibition by Cu2+ could be prevented by prior application of the reversible blocker of M2 channel activity, BL-1743, providing further support for the location of the site within the pore region of M2. Finally, substitutions of His-37 by alanine or glycine eliminated the high affinity site and resulted in membrane currents that were only partially inhibited at millimolar concentrations of Cu2+. Binding of Cu2+ to the high affinity site resulted in an approximately equal inhibition of both inward and outward currents. The wild-type protein showed very high specificity for Cu2+ and was only partially inhibited by 1 mM Ni2+, Pt2+, and Zn2+. These data are discussed in terms of the functional role of His-37 in the mechanism of proton translocation through the channel.     

Energetic analysis of the two controversial drug binding sites of the M2 proton channel in influenza A virus

Understanding the mechanism of the M2 proton channel of influenza A is crucially important to both basic research and drug discovery. Recently, the structure was determined independently by high-resolution NMR and X-ray crystallography. However, the two studies lead to completely different drug-binding mechanisms: the X-ray structure shows the drug blocking the pore from inside; whereas the NMR structure shows the drug inhibiting the channel from outside by an allosteric mechanism. Which one of the two is correct? To address this problem, we conducted an in-depth computational analysis. The conclusions drawn from various aspects, such as energetics, the channel-gating dynamic process, the pK_a shift and its impact on the channel, and the consistency with the previous functional studies, among others, are all in favour to the allosteric mechanism revealed by the NMR structure. The findings reported here may stimulate and encourage new strategies for developing effective drugs against influenza A, particularly in dealing with the drug-resistant problems.     

Analysis by metadynamics simulation of binding pathway of influenza virus M2 channel blockers

M2 protein of influenza A virus is a proton channel spanning the viral envelope. Activity of this proton channel is required for uncoating of viral particles and equilibrating the pH across the trans Golgi apparatus, which prevents conformational change in hemagglutinin. Amantadine, an anti‐influenza A virus drug, inhibits M2 proton channel activity by binding to the channel pore; however, most currently circulating influenza A viruses are amantadine‐resistant. The most prevalent resistant mutation is a substitution from Ser31 to Asn31 in M2. Further atomistic analysis of ligand‐M2 complexes is needed to provide new approaches for the design of novel M2 channel blockers. Here, the free energy profiles of the binding kinetics of M2 channel blockers were examined by well‐tempered metadynamics simulations and it was found that amantadine first binds to Asp24 of S31 M2 and forms a metastable conformation. In contrast, the free energy profiles of adamantyl bromothiophene dual inhibitor with either S31 M2 or N31 M2 are broad funnel‐shaped curves, suggesting that adamantyl bromothiophene does not form metastable complexes with M2. The trajectory of well‐tempered metadynamics simulations revealed that steric hindrance between adamantyl bromothiophene and S31 M2 interrupts formation of a metastable conformation at Asp24 and that a halogen bond between the bromine atom and N31 is responsible for pulling down the ligand to the channel pore of N31 M2 in the absence of a metastable state. Binding pathways of M2 channel blockers to M2 are here proposed on the basis of these findings; they may provide new approaches to designing further M2 channel blockers.     

The genesis and spread of reassortment human influenza A/H3N2 viruses conferring adamantane resistance

A dramatic rise in the frequency of resistance to adamantane drugs by influenza A (H3N2) viruses has occurred in recent years -- from approximately 2% to approximately 90% in multiple countries worldwide-and associated with a single S31N amino acid replacement in the viral matrix M2 protein. To explore the emergence and spread of these adamantane resistant viruses we performed a phylogenetic analysis of recently sampled complete A/H3N2 genome sequences. Strikingly, all adamantane resistant viruses belonged to a single lineage (the "N-lineage") characterized by 17 amino acid replacements across the viral genome. Further, our analysis revealed that the genesis of the N-lineage was due to a 4+4 segment reassortment event involving 2 distinct lineages of influenza A/H3N2 virus. A subsequent study of hemagglutinin HA1 sequences suggested that the N-lineage was circulating widely in Asia during 2005, and then dominated the Northern hemisphere 2005-2006 season in Japan and the USA. Given the infrequent use of adamantane drugs in many countries, as well as the decades of use in the US associated with little drug resistance, we propose that the globally increasing frequency of adamantane resistance is more likely attributable to its interaction with fitness-enhancing mutations at other genomic sites rather than to direct drug selection pressure. This implies that adamantanes may not be useful for treatment and prophylaxis against influenza viruses in the long term. More generally, these findings illustrate that drug selection pressure is not the sole factor determining the evolution and maintenance of drug resistance in human pathogens.     

Influenza virus M2 protein has ion channel activity.

The influenza virus M2 protein was expressed in Xenopus laevis oocytes and shown to have an associated ion channel activity selective for monovalent ions. The anti-influenza virus drug amantadine hydrochloride significantly attenuated the inward current induced by hyperpolarization of oocyte membranes. Mutations in the M2 membrane-spanning domain that confer viral resistance to amantadine produced currents that were resistant to the drug. Analysis of the currents of these altered M2 proteins suggests that the channel pore is formed by the transmembrane domain of the M2 protein. The wild-type M2 channel was found to be regulated by pH. The wild-type M2 ion channel activity is proposed to have a pivotal role in the biology of influenza virus infection.     

Characterization of inhibition of M2 ion channel activity by BL-1743, an inhibitor of influenza A virus.

The influenza A virus M2 integral membrane protein has ion channel activity that can be inhibited by the antiviral drug amantadine. Recently, a spirene-containing compound, BL-1743 (2-[3-azaspiro (5,5)undecanol]-2-imidazoline), that inhibits influenza virus growth was identified (S. Kurtz, G. Lao, K. M. Hahnenberger, C. Brooks, O. Gecha, K. Ingalls, K.-I. Numata, and M. Krystal, Antimicrob. Agents Chemother. 39:2204-2209, 1995). We have examined the ability of BL-1743 to inhibit the M2 ion channel when expressed in oocytes of Xenopus laevis. BL-1743 inhibition is complete as far as can be measured by electrophysiological methods and is reversible, with a reverse reaction rate constant of 4.0 x 10(-3) s(-1). In contrast, amantadine inhibition is irreversible within the time frame of the experiment. However, BL-1743 inhibition and amantadine inhibition have similar properties. The majority of isolated influenza viruses resistant to BL-1743 are also amantadine resistant. In addition, all known amino acid changes which result in amantadine resistance also confer BL-1743 resistance. However, one BL-1743-resistant virus isolated, designated M2-I35T, contained the change Ile-35-->Thr. This virus is >70-fold more resistant to BL-1743 and only 10-fold more resistant to amantadine than the wild-type virus. When the ion channel activity of M2-I35T was examined in oocytes, it was found that M2-I35T is BL-1743 resistant but is reversibly inhibited by amantadine. These findings suggest that these two drugs interact differently with the M2 protein transmembrane pore region.     

Ion channel activity of influenza A virus M2 protein: characterization of the amantadine block.

The influenza A virus M2 integral membrane protein has ion channel activity which can be blocked by the antiviral drug amantadine. The M2 protein transmembrane domain is highly conserved in amino acid sequence for all the human, swine, equine, and avian strains of influenza A virus, and thus, known amino acid differences could lead to altered properties of the M2 ion channel. We have expressed in oocytes of Xenopus laevis the M2 protein of human influenza virus A/Udorn/72 and the avian virus A/chicken/Germany/34 (fowl plague virus, Rostock) and derivatives of the Rostock ion channel altered in the presumed pore region. The pH of activation of the M2 ion channels and amantadine block of the M2 ion channels were investigated. The channels were found to be activated by pH in a similar manner but differed in their apparent Kis for amantadine block.     

Exploring models of the influenza A M2 channel: MD simulations in a phospholipid bilayer

The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (, Biophys. J. 77:2400-2410).     

Modeling of the pore domain of the GLUR1 channel: homology with K+ channel and binding of channel blockers

Molecular models of the M2 segments of the GluR1 channel have been elaborated using a molecular mechanics approach. The models are based on the homology between pore-lining segments of AMPA receptor channels and the KcsA K+ channel and on cyclic H bonds at the Q/R site of the AMPA receptor channel. The N-terminal region of an M2 segment of the channel is assumed, like that of the K+ channel, to adopt a helical conformation. Due to a deletion, the C-terminal end of the M2 segment of the AMPA receptor is more stretched than that of the K+ channel. As a result, only a single oxygen ring may be exposed to the AMPA receptor channel pore. Data on the block of AMPA receptor channels by dicationic adamantane derivatives have been used to select the most relevant model. The model with the oxygen of a Gly residue (position +2 from the Q/R site) exposed to the pore best fits the experimental data. This model also fits experimental data for another class of AMPA receptor antagonists, the polyamine amides. According to the model, the side-chains of the C-terminal residues are involved in intra-receptor interactions that stabilize the structure of the channel rather than in interactions with ions in the pore.     

Binding Hot Spots and Amantadine Orientation in the Influenza A Virus M2 Proton Channel

Structures of truncated versions of the influenza A virus M2 proton channel have been determined recently by x-ray crystallography in the open conformation of the channel, and by NMR in the closed state. The structures differ in the position of the bound inhibitors. The x-ray structure shows a single amantadine molecule in the middle of the channel, whereas in the NMR structure four drug molecules bind at the channel's outer surface. To study this controversy we applied computational solvent mapping, a technique developed for the identification of the most druggable binding hot spots of proteins. The method moves molecular probes—small organic molecules containing various functional groups—around the protein surface, finds favorable positions using empirical free energy functions, clusters the conformations, and ranks the clusters on the basis of the average free energy. The results of the mapping show that in both structures the primary hot spot is an internal cavity overlapping the amantadine binding site seen in the x-ray structure. However, both structures also have weaker hot spots at the exterior locations that bind rimantadine in the NMR structure, although these sites are partially due to the favorable interactions with the interfacial region of the lipid bilayer. As confirmed by docking calculations, the open channel binds amantadine at the more favorable internal site, in good agreement with the x-ray structure. In contrast, the NMR structure is based on a peptide/micelle construct that is able to accommodate the small molecular probes used for the mapping, but has a too narrow pore for the rimantadine to access the internal hot spot, and hence the drug can bind only at the exterior sites.     

Proton conduction through the M2 protein of the influenza A virus; a quantitative,mechanistic analysis of experimental data

The M2 proton channel from influenza A virus forms proton-selective ion channels, which are the target of the drug amantadine. Here, existing experimental data are quantitatively examined for insights into mechanisms to account for the pH- and voltage-dependences of M2 proton conduction. The analysis shows that a model involving protonation equilibria of His37, including pH-dependent changes in the relative rates of diffusion on either side of the pore, is quantitatively able to account for recently reported electrophysiological data examining the pH- and voltage-dependences of Rostock and Weybridge strain M2 proton conduction.     

Substrate flow in catalases deduced from the crystal structures of active site variants of HPII from Escherichia coli.

The active site of heme catalases is buried deep inside a structurally highly conserved homotetramer. Channels leading to the active site have been identified as potential routes for substrate flow and product release, although evidence in support of this model is limited. To investigate further the role of protein structure and molecular channels in catalysis, the crystal structures of four active site variants of catalase HPII from Escherichia coli (His128Ala, His128Asn, Asn201Ala, and Asn201His) have been determined at approximately 2.0-A resolution. The solvent organization shows major rearrangements with respect to native HPII, not only in the vicinity of the replaced residues but also in the main molecular channel leading to the heme distal pocket. In the two inactive His128 variants, continuous chains of hydrogen bonded water molecules extend from the molecular surface to the heme distal pocket filling the main channel. The differences in continuity of solvent molecules between the native and variant structures illustrate how sensitive the solvent matrix is to subtle changes in structure. It is hypothesized that the slightly larger H(2)O(2) passing through the channel of the native enzyme will promote the formation of a continuous chain of solvent and peroxide. The structure of the His128Asn variant complexed with hydrogen peroxide has also been determined at 2.3-A resolution, revealing the existence of hydrogen peroxide binding sites both in the heme distal pocket and in the main channel. Unexpectedly, the largest changes in protein structure resulting from peroxide binding are clustered on the heme proximal side and mainly involve residues in only two subunits, leading to a departure from the 222-point group symmetry of the native enzyme. An active role for channels in the selective flow of substrates through the catalase molecule is proposed as an integral feature of the catalytic mechanism. The Asn201His variant of HPII was found to contain unoxidized heme b in combination with the proximal side His-Tyr bond suggesting that the mechanistic pathways of the two reactions can be uncoupled.     

Voltage-dependent and -independent block of α-amino-3-hydroxy-5-methylisoxazole-4-propionate receptor channels

Polyamine-containing toxins and synthetic dicationic derivatives of adamantane and phenylcyclohexyl selectively antagonize Ca(2+)-permeable α-amino-3-hydroxy-5-methylisoxazole-4-propionate (AMPA) receptor channels. These compounds demonstrate voltage-dependent open-channel block and are trapped by closed channels. In this study, we describe an alternative mechanism of non-competitive AMPA receptor inhibition caused by 9-aminoacridine and some of its derivatives. These compounds exhibit similar potency against Ca(2+)-permeable and Ca(2+)-impermeable AMPA receptors. The inhibition is largely voltage-independent, binding and unbinding do not require presence of agonist. We conclude that 9-aminoacridine binds to a shallow site in the AMPA receptor, which is located above the activation gate. A comparison of three-dimensional structures of the antagonists suggests that the 'V-like' shape of the hydrophobic headgroup favors voltage-dependent binding to the deep site in the channel pore, whereas the compounds possessing flat aromatic headgroups preferably bind to the shallow site. The characterization of the novel mechanism of AMPA receptor channel antagonism opens a way to develop a new family of pharmacological agents, which can be of scientific and practical importance.     

Role of the lateral channel in catalase HPII of Escherichia coli

The heme-containing catalase HPII of Escherichia coli consists of a homotetramer in which each subunit contains a core region with the highly conserved catalase tertiary structure, to which are appended N- and C-terminal extensions making it the largest known catalase. HPII does not bind NADPH, a cofactor often found in catalases. In HPII, residues 585-590 of the C-terminal extension protrude into the pocket corresponding to the NADPH binding site in the bovine liver catalase. Despite this difference, residues that define the NADPH pocket in the bovine enzyme appear to be well preserved in HPII. Only two residues that interact ionically with NADPH in the bovine enzyme (Asp212 and His304) differ in HPII (Glu270 and Glu362), but their mutation to the bovine sequence did not promote nucleotide binding. The active-site heme groups are deeply buried inside the molecular structure requiring the movement of substrate and products through long channels. One potential channel is about 30 A in length, approaches the heme active site laterally, and is structurally related to the branched channel associated with the NADPH binding pocket in catalases that bind the dinucleotide. In HPII, the upper branch of this channel is interrupted by the presence of Arg260 ionically bound to Glu270. When Arg260 is replaced by alanine, there is a threefold increase in the catalytic activity of the enzyme. Inhibitors of HPII, including azide, cyanide, various sulfhydryl reagents, and alkylhydroxylamine derivatives, are effective at lower concentration on the Ala260 mutant enzyme compared to the wild-type enzyme. The crystal structure of the Ala260 mutant variant of HPII, determined at 2.3 A resolution, revealed a number of local structural changes resulting in the opening of a second branch in the lateral channel, which appears to be used by inhibitors for access to the active site, either as an inlet channel for substrate or an exhaust channel for reaction products.