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

From Acid Activation Mechanisms of Proton Conduction to Design of Inhibitors of the M2 Proton Channel of Influenza A Virus

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.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (N_δ or Ν_ε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having N_ε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 μs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

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

The V27A mutation confers adamantane resistance on the influenza A matrix 2 (M2) proton channel and is becoming more prevalent in circulating populations of influenza A virus. We have used X-ray crystallography to determine structures of a spiro-adamantyl amine inhibitor bound to M2(22-46) V27A and also to M2(21-61) V27A in the Inwardclosed conformation. The spiro-adamantyl amine binding site is nearly identical for the two crystal structures. Compared to the M2 "wild type" (WT) with valine at position 27, we observe that the channel pore is wider at its N-terminus as a result of the V27A mutation and that this removes V27 side chain hydrophobic interactions that are important for binding of amantadine and rimantadine. The spiro-adamantyl amine inhibitor blocks proton conductance in the WT and V27A mutant channels by shifting its binding site in the pore depending on which residue is present at position 27. Additionally, in the structure of the M2(21-61) V27A construct, the C-terminus of the channel is tightly packed relative to that of the M2(22-46) construct. We observe that residues Asp44, Arg45, and Phe48 face the center of the channel pore and would be well-positioned to interact with protons exiting the M2 channel after passing through the His37 gate. A 300 ns molecular dynamics simulation of the M2(22-46) V27A-spiro-adamantyl amine complex predicts with accuracy the position of the ligands and waters inside the pore in the X-ray crystal structure of the M2(22-46) V27A complex.

Decoupled side chain and backbone dynamics for proton translocation - M2 of influenza A

The 97 amino acid bitopic membrane protein M2 of influenza A forms a tetrameric bundle in which two of the monomers are covalently linked via a cysteine bridge. In its tetrameric assembly the protein conducts protons across the viral envelope and within intracellular compartments during the infectivity cycle of the virus. A key residue in the translocation of the protons is His-37 which forms a planar tetrad in the configuration of the bundle accepting and translocating the incoming protons from the N terminal side, exterior of the virus, to the C terminal side, inside the virus. With experimentally available data from NMR spectroscopy of the transmembrane domains of the tetrameric M2 bundle classical MD simulations are conducted with the protein bundle in different protonation stages in respect to His-37. A full correlation analysis (FCA) of the data sets with the His-37 tetrad either in a fully four times unprotonated or protonated state, assumed to mimic high and low pH in vivo, respectively, in both cases reveal asymmetric backbone dynamics. His-37 side chain rotation dynamics is increased at full protonation of the tetrad compared to the dynamics in the fully unprotonated state. The data suggest that proton translocation can be achieved by decoupled side chain or backbone dynamics. Graphical abstract Visualization of the tetrameric bundle of the transmembrane domains of M2 of influenza A after 200 ns of MD simulations (upper left). The four histidine residues 37 are either not protonated as in M2⁰ or fully protonated is in M2⁴⁺. The asymmetric dynamics of the backbones are shown after a full correlation analysis (FCA) in blue (lower left). The rotational dynamics of the χ2 dihedral angles of the histidines in M2⁰ (upper right) are less than those in M2⁴⁺ (lower right).

Resistance-Mutation (N31) Effects on Drug Orientation and Channel Hydration in Amantadine-Bound Influenza A M2

The mechanism of amantadine binding to the S31 variant of the M2 protein of Influenza A is well understood, but the reasons behind N31 M2 amantadine insensitivity remain under investigation. Many molecular dynamics studies have evaluated the influence of amantadine position within the channel pore on its ability to inhibit proton conductance in M2, but little is known about the influence of amantadine rotational orientation. Replica-exchange umbrella sampling, steered, and classic molecular dynamics simulations were performed on amantadine in the solid-state NMR structure of S31 M2 and an N31 M2 homologue, both in the homotetramer configuration, to explore the effects of the position and tilt angle of amantadine on inhibition of the M2 channel. Steered simulations show that amantadine rotates with the amine toward the bulk water as it passes into the hydrophobic entryway lined by Val27 side chains. Results from all simulation types performed indicate that amantadine has a strong, specific orientation with the amine turned inward toward the central cavity in the S31 M2 pore but has variable orientation and a strong propensity to remain outward pointing in N31 M2. Free energy profiles from umbrella sampling, measured relative to bulk water, show amantadine binds more strongly to the S31 M2 pore by 8 kcal/mol in comparison to amantadine in the N31 pore, suggesting that it can escape more readily from the N31 channel through the Val27 secondary gate, whereas it is captured by the S31 channel in the same region. Lower water density and distribution near amantadine in S31 M2 reveal that the drug inhibits proton conductance in S31 M2 because of its inward-pointing configuration, whereas in N31 M2, amantadine forms hydrogen bonds with an N31 side chain and does not widely occlude water occupancy in any configuration. Both amantadine's weaker binding to and weaker water occlusion in N31 M2 might contribute to its inefficacy as an inhibitor of the mutant protein.

Why bound amantadine fails to inhibit proton conductance according to simulations of the drug-resistant influenza A M2 (S31N)

The mechanisms responsible for drug resistance in the Asn31 variant of the M2 protein of influenza A are not well understood. Molecular dynamics simulations were performed on wild-type (Ser31) and S31N influenza A M2 in the homotetramer configuration. After evaluation of 13 published M2 structures, a solid-state NMR structure with amantadine bound was selected for simulations, an S31N mutant structure was developed and equilibrated, and the native and mutant structures were used to determine the binding behavior of amantadine and the dynamics of water in the two channels. Amantadine is stable in the plugging region of wild-type M2, with the adamantane in contact with the Val27 side chains, while amantadine in S31N M2 has more variable movement and orientation, and spontaneously moves lower into the central cavity of the channel. Free energy profiles from umbrella sampling support this observation. In this configuration, water surrounds the drug and can easily transport protons past it, so the drug binds without blocking proton transport in the S31N M2 channel.

Hydrogen-bonded water molecules in the M2 channel of the influenza A virus guide the binding preferences of ammonium-based inhibitors

The tetrameric M2 proton channel of influenza A virus is an integral membrane protein responsible for the acidification of the viral interior. Drugs such as amantadine target the transmembrane region of wild type M2 by acting as pore blockers. However, a number of mutations affecting this domain confer drug resistance, prompting the need for alternative inhibitors. The availability of high-resolution structures of drug-bound M2, paired with computational investigations, revealed that inhibitors can bind at different sites, and provided useful insights in understanding the principles governing proton conduction. Here, we investigated by computation the energetic and geometric factors determining the relative stability of pore blockers at individual sites of different M2 strains. We found that local free energy minima along the translocation pathway of positively charged chemical species correspond to experimentally determined binding sites of inhibitors. Then, by examining the structure of water clusters hydrating each site, as well as of those displaced by binding of hydrophobic scaffolds, we predicted the binding preferences of M2 ligands. This information can be used to guide the identification of novel classes of inhibitors.

Stepwise priming by acidic pH and a high K+ concentration is required for efficient uncoating of influenza A virus cores after penetration

Influenza A virus (IAV) uses the low pH in late endocytic vacuoles as a cue for penetration by membrane fusion. Here, we analyzed the prefusion reactions that prepare the core for uncoating after it has been delivered to the cytosol. We found that this priming process occurs in two steps that are mediated by the envelope-embedded M2 ion channel. The first weakens the interactions between the matrix protein, M1, and the viral ribonucleoprotein bundle. It involves a conformational change in a linker sequence and the C-terminal domain of M1 after exposure to a pH below 6.5. The second step is triggered by a pH of <6.0 and by the influx of K(+) ions. It causes additional changes in M1 as well as a loss of stability in the viral ribonucleoprotein bundle. Our results indicate that both the switch from Na(+) to K(+) in maturing endosomes and the decreasing pH are needed to prime IAV cores for efficient uncoating and infection of the host cell.

IMPORTANCE

The entry of IAV involves several steps, including endocytosis and fusion at late endosomes. Entry also includes disassembly of the viral core, which is composed of the viral ribonucleoproteins and the RNA genome. We have found that the uncoating process of IAV is initiated long before the core is delivered into the cytosol. M2, an ion channel in the viral membrane, is activated when the virus passes through early endosomes. Here, we show that protons entering the virus through M2 cause a conformational change in the matrix protein, M1. This weakens interactions between M1 and the viral ribonucleoproteins. A second change was found to occur when the virus enters late endosomes. The preacidified core is then exposed to a high concentration of K(+), which affects the interactions between the ribonucleoproteins. Thus, when cores are finally delivered to the cytosol, they are already partially destabilized and, therefore, uncoating competent and infectious.

Multiscale molecular dynamics simulations of membrane proteins

The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.

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.

Use of Umbrella Sampling to Calculate the Entrance/Exit Pathway for Z-Pro-Prolinal Inhibitor in Prolyl Oligopeptidase

Prolyl oligopeptidase (POP), a member of the prolyl endopeptidase family, is known to play a role in several neurological disorders. Its primary function is to cleave a wide range of small oligopeptides, including neuroactive peptides. We have used force biased molecular dynamics simulation to study the binding mechanism of POP. We examined three possible binding pathways using Steered Molecular Dynamics (SMD) and Umbrella Sampling (US) on a crystal structure of porcine POP with bound Z-pro-prolinal (ZPP). Using SMD, an exit pathway between the first and seventh blade of the β-propeller domain of POP was found to be a nonviable route. US on binding pathways through the β-propeller tunnel and the TYR190-GLN208 flexible loop at the interface between both POP domains allowed us to isolate the flexible loop pathway as the most probable. Further analysis of that pathway suggests a long-range covariation of the interdomain H-bond network, which indicates the possibility of large-scale domain reorientation observed in bacterial homologues and hypothesized to also occur in human POP.

Multiple Proton Confinement in the M2 Channel from the Influenza A Virus

The tetrameric M2 protein bundle of the influenza A virus is the proton channel responsible for the acidification of the viral interior, a key step in the infection cycle. Selective proton transport is achieved by successive protonation of the conserved histidine amino acids at position 37. A recent X-ray structure of the tetrameric transmembrane (TM) domain of the protein (residues 22-46) resolved several water clusters in the channel lumen, which suggest possible proton pathways to the His37 residues. To explore this hypothesis, we have carried out molecular dynamics (MD) simulations of a proton traveling towards the His37 side chains using MD with classical and quantum force fields. Diffusion through the first half of the channel to the "entry" water cluster near His37 may be hampered by significant kinetic barriers due to electrostatic repulsion. However, once in the entry cluster, a proton can move to one of the acceptor His37 in a nearly barrierless fashion, as evidenced both by MD simulations and a scan of the potential energy surface (PES). Water molecules of the entry cluster, although confined in the M2 pore and restricted in their motions, can conduct protons with a rate very similar to that of bulk water.

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

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

Structural basis for the function and inhibition of an influenza virus proton channel

The M2 protein from influenza A virus is a pH-activated proton channel that mediates acidification of the interior of viral particles entrapped in endosomes. M2 is the target of the anti-influenza drugs amantadine and rimantadine; recently, resistance to these drugs in humans, birds and pigs has reached more than 90% (ref. 1). Here we describe the crystal structure of the transmembrane-spanning region of the homotetrameric protein in the presence and absence of the channel-blocking drug amantadine. pH-dependent structural changes occur near a set of conserved His and Trp residues that are involved in proton gating. The drug-binding site is lined by residues that are mutated in amantadine-resistant viruses. Binding of amantadine physically occludes the pore, and might also perturb the pK(a) of the critical His residue. The structure provides a starting point for solving the problem of resistance to M2-channel blockers.