How pH opens a H+ channel: the gating mechanism of influenza A M2

Conformational heterogeneity and structural features for function of the prototype viroporin influenza AM2

The 97-residue influenza A matrix 2 (ΑM2) protein, a prototype for viroporins, transports protons through water molecules and His37. We discuss structural biology and molecular biophysics experiments and some functional assays that have transformed over 40 years our understanding of the structure and function of AM2. The structural studies on ΑM2 have been performed with different conditions (pH, temperature, lipid, constructs) and using various protein constructs, e.g., AM2 transmembrane (AM2TM) domain, AM2 conductance domain (AM2CD), ectodomain-containing or ectodomain-truncated, AM2 full length (AM2FL) and aimed to describe the different conformations and structural details that are necessary for the stability and function of AM2. However, the conclusions from these experiments appeared sometimes ambiguous and caused exciting debates. This was not due to inaccurate measurements, but instead because of the different membrane mimetic environment used, e.g., detergent, micelles or phospholipid bilayer, the method (e.g., X-ray crystallography, solid state NMR, solution NMR, native mass spectrometry), the used protein construct (e.g., AM2TM or AM2CD), or the amino acids residues to follow observables (e.g., NMR chemical shifts). We present these results according to the different used biophysical methods, the research groups and often by keeping a chronological order for presenting the progress in the research. We discuss ideas for additional research on structural details of AM2 and how the present findings can be useful to explore new routes of influenza A inhibition. The AM2 research can provide inspiration to study other viroporins as drug targets.

The Influenza A Virus Replication Cycle: A Comprehensive Review

Influenza A virus (IAV) is the primary causative agent of influenza, colloquially called the flu. Each year, it infects up to a billion people, resulting in hundreds of thousands of human deaths, and causes devastating avian outbreaks with worldwide losses worth billions of dollars. Always present is the possibility that a highly pathogenic novel subtype capable of direct human-to-human transmission will spill over into humans, causing a pandemic as devastating if not more so than the 1918 influenza pandemic. While antiviral drugs for influenza do exist, they target very few aspects of IAV replication and risk becoming obsolete due to antiviral resistance. Antivirals targeting other areas of IAV replication are needed to overcome this resistance and combat the yearly epidemics, which exact a serious toll worldwide. This review aims to summarise the key steps in the IAV replication cycle, along with highlighting areas of research that need more focus.

Understanding the conformational changes in the influenza B M2 ion channel at various protonation states

The characterization of influenza (A/B M2) ion channels is very important as they are potential binding sites for the drugs. We report the all-atom molecular dynamics study of the influenza B M2 ion channel in the presence of explicit solvent and lipid bilayers using the high resolution solid-state NMR structures. The importance of the various protonation states of histidine in the activation of the ion channel is discussed. The conformational changes at the closed and the open structures clearly show that the increase in tilt angle is necessary for the activation of the ion channel. Additionally, the free energy surfaces of the eight systems show the importance of the protonation state of the histidine residues in the activation of the influenza B M2 ion channel. The protonation of the histidine residues increases the tilt angle and the intra-helix distance which is evident from the superimposition of the structures corresponding to the maxima and the minima in the free energy landscape. The findings imply differences in the singly protonated and double protonated conformational states of BM2 ion channel and provide insights to help further studies of these ion channels as the drug targets for the influenza virus.

Pore-Bound Water at the Key Residue Histidine 37 in Influenza A M2

Atomic details of structured water molecules are indispensable to understand the thermodynamics of important biological processes including the proton conduction mechanism of the M2 protein. Despite the expectation of structured water molecules based on crystal structures of Influenza A M2, only two water populations have been observed by NMR in reconstituted lipid bilayer samples. These are the bulk‐ and lipid‐associated water populations typically seen in membrane samples. Here, we detect a bound water molecule at a chemical shift of 11 ppm, located near the functional histidine 37 residue in the M2 conductance domain, which comprises residues 18 to 60. Combining 100 kHz magic‐angle spinning NMR, dynamic nuclear polarization and density functional theory calculations, we show that the bound water forms a hydrogen bond to the δ1 nitrogen of histidine 37.

Influenza Viruses: Harnessing the Crucial Role of the M2 Ion-Channel and Neuraminidase toward Inhibitor Design

As a member of the Orthomyxoviridae family of viruses, influenza viruses (IVs) are known causative agents of respiratory infection in vertebrates. They remain a major global threat responsible for the most virulent diseases and global pandemics in humans. The virulence of IVs and the consequential high morbidity and mortality of IV infections are primarily attributed to the high mutation rates in the IVs' genome coupled with the numerous genomic segments, which give rise to antiviral resistant and vaccine evading strains. Current therapeutic options include vaccines and small molecule inhibitors, which therapeutically target various catalytic processes in IVs. However, the periodic emergence of new IV strains necessitates the continuous development of novel anti-influenza therapeutic options. The crux of this review highlights the recent studies on the biology of influenza viruses, focusing on the structure, function, and mechanism of action of the M2 channel and neuraminidase as therapeutic targets. We further provide an update on the development of new M2 channel and neuraminidase inhibitors as an alternative to existing anti-influenza therapy. We conclude by highlighting therapeutic strategies that could be explored further towards the design of novel anti-influenza inhibitors with the ability to inhibit resistant strains.

A pH-Dependent Conformational Switch Controls N. meningitidis ClpP Protease Function

ClpPs are a conserved family of serine proteases that collaborate with ATP-dependent translocases to degrade protein substrates. Drugs targeting these enzymes have attracted interest for the treatment of cancer and bacterial infections due to their critical role in mitochondrial and bacterial proteostasis, respectively. As such, there is significant interest in understanding structure-function relationships in this protein family. ClpPs are known to crystallize in extended, compact, and compressed forms; however, it is unclear what conditions favor the formation of each form and whether they are populated by wild-type enzymes in solution. Here, we use cryo-EM and solution NMR spectroscopy to demonstrate that a pH-dependent conformational switch controls an equilibrium between the active extended and inactive compressed forms of ClpP from the Gram-negative pathogen Neisseria meningitidis. Our findings provide insight into how ClpPs exploit their rugged energy landscapes to enable key conformational changes that regulate their function.

Influenza A Virus M2 Protein: Roles from Ingress to Egress

Influenza A virus (IAV) matrix protein 2 (M2) is among the smallest bona fide, hence extensively studied, ion channel proteins. The M2 ion channel activity is not only essential for virus replication, but also involved in modulation of cellular homeostasis in a variety of ways. It is also the target for ion channel inhibitors, i.e., anti-influenza drugs. Thus far, several studies have been conducted to elucidate its biophysical characteristics, structure-function relationships of the ion channel, and the M2-host interactome. In this review, we discuss M2 protein synthesis and assembly into an ion channel, its roles in IAV replication, and the pathophysiological impact on the host cell.

Activation pH and Gating Dynamics of Influenza A M2 Proton Channel Revealed by Single-Molecule Spectroscopy

Because of its importance in viral replication, the M2 proton channel of the influenza A virus has been the focus of many studies. Although we now know a great deal about the structural architecture underlying its proton conduction function, we know little about its conformational dynamics, especially those controlling the rate of this action. Herein, we employ a single-molecule fluorescence method to assess the dynamics of the inter-helical channel motion of both full-length M2 and the transmembrane domain of M2. The rate of this motion depends not only on the identity of the channel and membrane composition but also on the pH in a sigmoidal manner. For the full-length M2 channel, the rate is increased from approximately 190 μs-1 at high pH to approximately 80 μs-1 at low pH, with a transition midpoint at pH 6.1. Because the latter value is within the range reported for the conducting pKa value of the His37 tetrad, we believe that this inter-helical motion accompanies proton conduction.

Infrared and fluorescence assessment of the hydration status of the tryptophan gate in the influenza A M2 proton channel

The M2 proton channel of the influenza A virus has been the subject of extensive studies because of its critical role in viral replication. As such, we now know a great deal about its mechanism of action, especially how it selects and conducts protons in an asymmetric fashion. The conductance of this channel is tuned to conduct protons at a relatively low biologically useful rate, which allows acidification of the viral interior of a virus entrapped within an endosome, but not so great as to cause toxicity to the infected host cell prior to packaging of the virus. The dynamic, structural and chemical features that give rise to this tuning are not fully understood. Herein, we use a tryptophan (Trp) analog, 5-cyanotryptophan, and various methods, including linear and nonlinear infrared spectroscopies, static and time-resolved fluorescence techniques, and molecular dynamics simulations, to site-specifically interrogate the structure and hydration dynamics of the Trp41 gate in the transmembrane domain of the M2 proton channel. Our results suggest that the Trp41 sidechain adopts the t90 rotamer, the χ2 dihedral angle of which undergoes an increase of approximately 35° upon changing the pH from 7.4 to 5.0. Furthermore, we find that Trp41 is situated in an environment lacking bulk-like water, and somewhat surprisingly, the water density and dynamics do not show a measurable difference between the high (7.4) and low (5.0) pH states. Since previous studies have shown that upon channel opening water flows into the cavity above the histidine tetrad (His37), the present finding thus provides evidence indicating that the lack of sufficient water molecules near Trp41 needed to establish a continuous hydrogen bonding network poses an additional energetic bottleneck for proton conduction.

Biomimetic heterogeneous multiple ion channels: a honeycomb structure composite film generated by breath figures

We design a novel type of artificial multiple nanochannel system with remarkable ion rectification behavior via a facile breath figure (BF) method. Notably, even though the charge polarity in the channel wall reverses under different pH values, this nanofluidic device displays the same ionic rectification direction. Compared with traditional nanochannels, this composite multiple ion channel device can be more easily obtained and has directional ionic rectification advantages, which can be applied in many fields.

Modelling of Thyroid Peroxidase Reveals Insights into Its Enzyme Function and Autoantigenicity

Thyroid peroxidase (TPO) catalyses the biosynthesis of thyroid hormones and is a major autoantigen in Hashimoto's disease--the most common organ-specific autoimmune disease. Epitope mapping studies have shown that the autoimmune response to TPO is directed mainly at two surface regions on the molecule: immunodominant regions A and B (IDR-A, and IDR-B). TPO has been a major target for structural studies for over 20 years; however, to date, the structure of TPO remains to be determined. We have used a molecular modelling approach to investigate plausible modes of TPO structure and dimer organisation. Sequence features of the C-terminus are consistent with a coiled-coil dimerization motif that most likely anchors the TPO dimer in the apical membrane of thyroid follicular cells. Two contrasting models of TPO were produced, differing in the orientation and exposure of their active sites relative to the membrane. Both models are equally plausible based upon the known enzymatic function of TPO. The "trans" model places IDR-B on the membrane-facing side of the myeloperoxidase (MPO)-like domain, potentially hindering access of autoantibodies, necessitating considerable conformational change, and perhaps even dissociation of the dimer into monomers. IDR-A spans MPO- and CCP-like domains and is relatively fragmented compared to IDR-B, therefore most likely requiring domain rearrangements in order to coalesce into one compact epitope. Less epitope fragmentation and higher solvent accessibility of the "cis" model favours it slightly over the "trans" model. Here, IDR-B clusters towards the surface of the MPO-like domain facing the thyroid follicular lumen preventing steric hindrance of autoantibodies. However, conformational rearrangements may still be necessary to allow full engagement with autoantibodies, with IDR-B on both models being close to the dimer interface. Taken together, the modelling highlights the need to consider the oligomeric state of TPO, its conformational properties, and its proximity to the membrane, when interpreting epitope-mapping data.

Drug inhibition and proton conduction mechanisms of the influenza a M2 proton channel

The influenza A virus matrix protein 2 (M2 protein) is a pH-regulated proton channel embedded in the viral membrane. Inhibition of the M2 proton channel has been used to treat influenza infections for decades due to the crucial role of this protein in viral infection and replication. However, the widely-used M2 inhibitors, amantadine and rimantadine, have gradually lost their efficiencies because of naturally-occurring drug resistant mutations. Therefore, investigation of the structure and function of the M2 proton channel will not only increase our understanding of this important biological system, but also lead to the design of novel and effective anti-influenza drugs. Despite the simplicity of the M2 molecular structure, the M2 channel is highly flexible and there have been controversies and arguments regarding the channel inhibition mechanism and the proton conduction mechanism. In this book chapter, we will first carefully review the experimental and computational studies of the two possible drug binding sites on the M2 protein and explain the mechanisms regarding how inhibitors prevent proton conduction. Then, we will summarize our recent molecular dynamics simulations of the drug-resistant mutant channels and propose mechanisms for drug resistance. Finally, we will discuss two existing proton conduction mechanisms and talk about the remaining questions regarding the proton-relay process through the channel. The studies reviewed here demonstrate how molecular modeling and simulations have complemented experimental work and helped us understand the M2 channel structure and function.

Fast Atomic Charge Calculation for Implementation into a Polarizable Force Field and Application to an Ion Channel Protein

Polarization of atoms plays a substantial role in molecular interactions. Class I and II force fields mostly calculate with fixed atomic charges which can cause inadequate descriptions for highly charged molecules, for example, ion channels or metalloproteins. Changes in charge distributions can be included into molecular mechanics calculations by various methods. Here, we present a very fast computational quantum mechanical method, the Bond Polarization Theory (BPT). Atomic charges are obtained via a charge calculation method that depend on the 3D structure of the system in a similar way as atomic charges ofab initiocalculations. Different methods of population analysis and charge calculation methods and their dependence on the basis set were investigated. A refined parameterization yielded excellent correlation ofR=0.9967. The method was implemented in the force field COSMOS-NMR and applied to the histidine-tryptophan-complex of the transmembrane domain of the M2 protein channel of influenza A virus. Our calculations show that moderate changes of side chain torsion angleχ1and small variations ofχ2of Trp-41 are necessary to switch from the inactivated into the activated state; and a rough two-side jump model of His-37 is supported for proton gating in accordance with a flipping mechanism.

Activation and proton transport mechanism in influenza A M2 channel

Molecular dynamics trajectories 2 μs in length have been generated for the pH-activated, tetrameric M2 proton channel of the influenza A virus in all protonation states of the pH sensor located at the His(37) tetrad. All simulated structures are in very good agreement with high-resolution structures. Changes in the channel caused by progressive protonation of His(37) provide insight into the mechanism of proton transport. The channel is closed at both His(37) and Trp(41) sites in the singly and doubly protonated states, but it opens at Trp(41) upon further protonation. Anions access the charged His(37) and by doing so stabilize the protonated states of the channel. The narrow opening at the His(37) site, further blocked by anions, is inconsistent with the water-wire mechanism of proton transport. Instead, conformational interconversions of His(37) correlated with hydrogen bonding to water molecules indicate that these residues shuttle protons in high-protonation states. Hydrogen bonds between charged and uncharged histidines are rare. The valve at Val(27) remains on average quite narrow in all protonation states but fluctuates sufficiently to support water and proton transport. A proton transport mechanism in which the channel, depending on pH, opens at either the histidine or valine gate is only partially supported by the simulations.

Philosophy of voltage-gated proton channels

In this review, voltage-gated proton channels are considered from a mainly teleological perspective. Why do proton channels exist? What good are they? Why did they go to such lengths to develop several unique hallmark properties such as extreme selectivity and ΔpH-dependent gating? Why is their current so minuscule? How do they manage to be so selective? What is the basis for our belief that they conduct H(+) and not OH(-)? Why do they exist in many species as dimers when the monomeric form seems to work quite well? It is hoped that pondering these questions will provide an introduction to these channels and a way to logically organize their peculiar properties as well as to understand how they are able to carry out some of their better-established biological functions.

Ab initio calculations and validation of the pH-dependent structures of the His37-Trp41 quartet, the heart of acid activation and proton conductance in the M2 protein of Influenza A virus

The M2 protein of Influenza A virus forms a homotetrameric proton channel activated by low pH. The His37-Trp41 quartet is the heart of acid activation and proton conductance, but the functional mechanism is still controversial. We carried out ab initio calculations to model the pH-dependent structures of the His37-Trp41 quartet. In our model at neutral pH, the four His37 residues are configured into a pair of dimers; in each dimer, a proton is shared between Nδ1 on one residue and Nε2 on the other, and, under the restraint of the backbone, the two imidazole rings are nearly parallel, in contrast to a perpendicular arrangement for a free imidazole-imidazolium dimer. Within each dimer the +1 charge is highly delocalized, contributing to its stabilization in a low dielectric environment. The Nδ1-H-Nε2 strong hydrogen bonds result in significantly downfield shifted Nδ1 and Nε2 chemical shifts (at 169.7 and 167.6 ppm, respectively), in good agreement with experiments. In our model at acidic pH (where the channel becomes activated), a third proton binds to an imidazole-imidazolium dimer; the imidazole rings rotate away (each by ~55°) from each other, destroying the dimer structure. The two imidazoliums are stabilized by hydrogen bonds with water molecules and a cation-π interaction with Trp41. The Raman spectra calculated for the His37-Trp41 quartet at neutral and acidic pH are in agreement with experiments. Our calculations support an activation and conductance mechanism in which a hydronium ion from the N-terminal side passes a proton to an imidazole-imidazolium dimer; when the Trp41 gate is open, relaying of a proton onto a water molecule from the C-terminal side then allows the imidazole-imidazolium dimer to reform and be ready for the next round of proton conductance.

Modeling the membrane environment has implications for membrane protein structure and function: influenza A M2 protein

The M2 protein, a proton channel, from Influenza A has been structurally characterized by X-ray diffraction and by solution and solid-state NMR spectroscopy in a variety of membrane mimetic environments. These structures show substantial backbone differences even though they all present a left-handed tetrameric helical bundle for the transmembrane domain. Variations in the helix tilt influence drug binding and the chemistry of the histidine tetrad responsible for acid activation, proton selectivity and transport. Some of the major structural differences do not arise from the lack of precision, but instead can be traced to the influences of the membrane mimetic environments. The structure in lipid bilayers displays unique chemistry for the histidine tetrad, which binds two protons cooperatively to form a pair of imidazole-imidazolium dimers. The resulting interhistidine hydrogen bonds contribute to a three orders of magnitude enhancement in tetramer stability. Integration with computation has provided detailed understanding of the functional mechanism for proton selectivity, conductance and gating of this important drug target.

Voltage-gated proton channels

Voltage-gated proton channels, HV1, have vaulted from the realm of the esoteric into the forefront of a central question facing ion channel biophysicists, namely, the mechanism by which voltage-dependent gating occurs. This transformation is the result of several factors. Identification of the gene in 2006 revealed that proton channels are homologues of the voltage-sensing domain of most other voltage-gated ion channels. Unique, or at least eccentric, properties of proton channels include dimeric architecture with dual conduction pathways, perfect proton selectivity, a single-channel conductance approximately 10(3) times smaller than most ion channels, voltage-dependent gating that is strongly modulated by the pH gradient, ΔpH, and potent inhibition by Zn(2+) (in many species) but an absence of other potent inhibitors. The recent identification of HV1 in three unicellular marine plankton species has dramatically expanded the phylogenetic family tree. Interest in proton channels in their own right has increased as important physiological roles have been identified in many cells. Proton channels trigger the bioluminescent flash of dinoflagellates, facilitate calcification by coccolithophores, regulate pH-dependent processes in eggs and sperm during fertilization, secrete acid to control the pH of airway fluids, facilitate histamine secretion by basophils, and play a signaling role in facilitating B-cell receptor mediated responses in B-lymphocytes. The most elaborate and best-established functions occur in phagocytes, where proton channels optimize the activity of NADPH oxidase, an important producer of reactive oxygen species. Proton efflux mediated by HV1 balances the charge translocated across the membrane by electrons through NADPH oxidase, minimizes changes in cytoplasmic and phagosomal pH, limits osmotic swelling of the phagosome, and provides substrate H(+) for the production of H2O2 and HOCl, reactive oxygen species crucial to killing pathogens.

Proton and calcium-gated ionic mesochannels: phosphate-bearing polymer brushes hosted in mesoporous thin films as biomimetic interfacial architectures

Rational construction of interfaces based on multicomponent responsive systems in which molecular transport is mediated by structures of nanoscale dimensions has become a very fertile research area in biomimetic supramolecular chemistry. Herein, we describe the creation of hybrid mesostructured interfaces with reversible gate-like transport properties that can be controlled by chemical inputs, such as protons or calcium ions. This was accomplished by taking advantage of the surface-initiated polymerization of 2-(methacryloyloxy)ethyl phosphate (MEP) monomer units into and onto mesoporous silica thin films. In this way, phosphate-bearing polymer brushes were used as "gatekeepers" located not only on the outer surface of mesoporous thin films but also in the inner environment of the porous scaffold. Pore-confined PMEP brushes respond to the external triggering chemical signals not only by altering their physicochemical properties but also by switching the transport properties of the mesoporous film. The ion-gate response/operation was based on the protonation and/or chelation of phosphate monomer units in which the polymer brush works as an off-on switch in response to the presence of protons or Ca(2+) ions. The hybrid meso-architectured interface and their functional features were studied by a combination of experimental techniques including ellipso-porosimetry, cyclic voltammetry, X-ray reflectivity, grazing incidence small-angle X-ray scattering, X-ray photoelectron spectroscopy, and in situ atomic force microscopy. In this context, we believe that the integration of stimuli-responsive polymer brushes into nanoscopic supramolecular architectures would provide new routes toward multifunctional biomimetic nanosystems displaying transport properties similar to those encountered in biological ligand-gated ion channels.

Free energy calculations on the two drug binding sites in the M2 proton channel

Two alternative binding sites of adamantane-type drugs in the influenza A M2 channel have been suggested, one with the drug binding inside the channel pore and the other with four drug molecule S-binding to the C-terminal surface of the transmembrane domain. Recent computational and experimental studies have suggested that the pore binding site is more energetically favorable but the external surface binding site may also exist. Nonetheless, which drug binding site leads to channel inhibition in vivo and how drug-resistant mutations affect these sites are not completely understood. We applied molecular dynamics simulations and potential of mean force calculations to examine the structures and the free energies associated with these putative drug binding sites in an M2-lipid bilayer system. We found that, at biological pH (~7.4), the pore binding site is more thermodynamically favorable than the surface binding site by ~7 kcal/mol and, hence, would lead to more stable drug binding and channel inhibition. This result is in excellent agreement with several recent studies. More importantly, a novel finding of ours is that binding to the channel pore requires overcoming a much higher energy barrier of ~10 kcal/mol than binding to the C-terminal channel surface, indicating that the latter site is more kinetically favorable. Our study is the first computational work that provides both kinetic and thermodynamic energy information on these drug binding sites. Our results provide a theoretical framework to interpret and reconcile existing and often conflicting results regarding these two binding sites, thus helping to expand our understanding of M2-drug binding, and may help guide the design and screening of novel drugs to combat the virus.

Computational methods for studying serpin conformational change and structural plasticity

Currently, over a hundred high-resolution structures of serpins are available, exhibiting a wide range of conformations. However, our understanding of serpin dynamics and conformational change is still limited, mainly due to challenges of monitoring structural changes and characterizing transient conformations using experimental methods. Insight can be provided, however, by employing theoretical and computational approaches. In this chapter, we present an overview of such methods, focusing on molecular dynamics and simulation. As serpin conformational dynamics span a wide range of timescales, we discuss the relative merits of each method and suggest which method is suited to specific conformational phenomena.

How does influenza virus a escape from amantadine?

Antiflu drugs such as amantadine (AMT) were reported to be insensitive to influenza A virus gradually after their marketing. Mutation experiments indicate that the trans-membrane domain of M2 protein plays an essential role in AMT resistance, especially the S31N mutation. To investigate the details of structure and mechanism, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out on both the wild-type protein and its S31N mutant. Our MD simulations reveal AMT can occupy different binding positions in the pore of M2 channel, and the binding modes have also been verified and analyzed by QM/MM calculations. More importantly, we find the formation of a water wire modulated by the binding position of AMT to be essential for the function of M2 protein, and, the block of water wire can inhibit channel function in the WT system. Failure of channel blocking would cause AMT drug resistance in the S31N mutant. These results support one of the conflicting views about M2-drug binding sites: AMT binds to the pore of M2 channel. Our findings help clarify the resistant mechanism of AMT to M2 protein and should facilitate the discovery of new drugs for treating influenza A virus.

Ion channels as antivirus targets

Ion channels are membrane proteins that are found in a number of viruses and which are of crucial physiological importance in the viral life cycle. They have one common feature in that their action mode involves a change of electrochemical or proton gradient across the bilayer lipid membrane which modulates viral or cellular activity. We will discuss a group of viral channel proteins that belong to the viroproin family, and which participate in a number of viral functions including promoting the release of viral particles from cells. Blocking these channel-forming proteins may be "lethal", which can be a suitable and potential therapeutic strategy. In this review we discuss seven ion channels of viruses which can lead serious infections in human beings: M2 of influenza A, NB and BM2 of influenza B, CM2 of influenza C, Vpu of HIV-1, p7 of HCV and 2B of picornaviruses.

Proton-regulated rectified ionic transport through solid-state conical nanopores modified with phosphate-bearing polymer brushes

We describe the use of polyprotic polymer brushes to construct robust signal-responsive chemical devices mimicking the transport properties of proton regulated biological channels.

Simulations of the BM2 proton channel transmembrane domain from influenza virus B

BM2 is a small integral membrane protein from influenza B virus which forms proton-permeable channels. Coarse-grained (CG) molecular dynamics simulations have been used to produce a model of the BM2 channel by self-assembly of a tetrameric bundle of BM2 transmembrane helices in a lipid bilayer. The BM2 channel model is conformationally stable on a 5 mus time scale. This CG model was converted to atomistic resolution to refine interhelix and channel-water interactions. Atomistic molecular dynamics simulations indicate that the BM2 channel is closed when no more than two of the four His19 residues are protonated. Protonating a third His19 side chain initiates a conformational change that opens the channel. In summary, our simulations suggest a common mechanism for BM2 and A/M2, whereby changes in helix packing play a functional role in channel gating.

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

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

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.

Synthetic proton-gated ion channels via single solid-state nanochannels modified with responsive polymer brushes

The creation of switchable and tunable nanodevices displaying transport properties similar to those observed in biological pores poses a major challenge in molecular nanotechnology. Here, we describe the construction of a fully "abiotic" nanodevice whose transport properties can be accurately controlled by manipulating the proton concentration in the surrounding environment. The ionic current switching characteristics displayed by the nanochannels resemble the typical behavior observed in many biological channels that fulfill key pH-dependent transport functions in living organisms, that is, the nanochannel can be switched from an "off" state to an "on" state in response to a pH drop. The construction of such a chemical nanoarchitecture required the integration of stable and ductile macromolecular building blocks constituted of pH-responsive poly(4-vinyl pyridine) brushes into solid state nanopores that could act as gate-keepers managing and constraining the flow of ionic species through the confined environment. In this context, we envision that the integration of environmental stimuli-responsive brushes into solid-state nanochannels would provide a plethora of new chemical alternatives for molecularly design robust signal-responsive "abiotic" devices mimicking the function of proton-gated ion channels commonly encountered in biological membranes.

Why amantadine loses its function in influenza m2 mutants: MD simulations

Molecular dynamics simulations of the drug-resistant M2 mutants, A30T, S31N, and L26I, were carried out to investigate the inhibition of M2 activity using amantadine (AMT). The closed and open channel conformations were examined via non- and triply protonated H37. For the nonprotonated state, these mutants exhibited zero water density in the conducting region, and AMT was still bound to the channel pore. Thus, water transport is totally suppressed, similar to the wild-type channel. In contrast, the triply protonated states of the mutants exhibited a different water density and AMT position. A30T and L26I both have a greater water density compared to the wild-type M2, while for the A30T system, AMT is no longer inside the pore. Hydrogen bonding between AMT and H37 crucial for the bioactivity is entirely lost in the open conformation. The elimination of this important interaction of these mutations is responsible for the lost of AMT's function in influenza A M2. This is different for the S31N mutant in which AMT was observed to locate at the pore opening region and bond with V27 instead of S31.

pH-driven helix rotations in the influenza M2 H+ channel: a potential gating mechanism

The pH activated M2 H(+) channel from influenza A has been a subject of numerous studies due to following: (1) It serves as a target for the aminoadamantane drugs that block its channel activity. (2) M2's small size makes it amenable to biophysical scrutiny. (3) A single histidine residue is thought to control the pH gating of the channel. Recent FTIR analysis proposed that the helices of the channel rotate about their directors during pH activation. Herein, we report on molecular dynamics simulations of the X-ray structure of the protein with three charged histidine residues, representing the open form of the protein and two rotated forms with neutral histidines, representing its closed form. We compare the channel stability, convergence, interaction with water and hydration of the histidine residues that have been implicated in channel gating. Taken together, we show that both forms of the protein are stable during the course of the MD simulation and that indeed a rotation of the helices leads to channel closure. Finally, we propose a mechanism for channel gating that involves protonation of the histidine residues that necessities their increased solvation.

Molecular dynamics calculations suggest a conduction mechanism for the M2 proton channel from influenza A virus

The M2 protein of the influenza A virus is activated by low endosomal pH and performs the essential function of proton transfer into the viral interior. The resulting decrease in pH within the virion is essential for the uncoating and further replication of the viral genetic material. The x-ray crystal [Stouffer AL, et al. (2008) Nature 451:596-599] and solution NMR [Schnell JR, Chou JJ (2008) Nature 451:591-595] structures of the transmembrane region of the M2 homo-tetrameric bundle both revealed pores with narrow constrictions at one end, leaving a question as to how protons enter the channel. His-37, which is essential for proton-gating and selective conduction of protons, lies in the pore of the crystallographic and NMR structures. Here, we explore the different protonation states of the His-37 residues of the M2 bundle in a bilayer using molecular dynamics (MD) simulations. When the His-37 residues are neutral, the protein prefers an Open(out)-Closed(in) conformation in which the channel is open to the environment on the outside of the virus but closed to the interior environment of the virus. Diffusion of protons into the channel from the outside of the virus and protonation of His-37 residues in the tetramer stabilizes an oppositely gated Closed(out)-Open(in) conformation. Thus, protons might be conducted through a transporter-like mechanism, in which the protein alternates between Open(out)-Closed(in) and Closed(out)-Open(in) conformations, and His-37 is protonated/deprotonated during each turnover. The transporter-like mechanism is consistent with the known properties of the M2 bundle, including its relatively low rate of proton flux and its strong rectifying behavior.

Voltage-gated proton channels: what's next?

This review is an attempt to identify and place in context some of the many questions about voltage-gated proton channels that remain unsolved. As the gene was identified only 2 years ago, the situation is very different than in fields where the gene has been known for decades. For the proton channel, most of the obvious and less obvious structure-function questions are still wide open. Remarkably, the proton channel protein strongly resembles the voltage-sensing domain of many voltage-gated ion channels, and thus offers a novel approach to study gating mechanisms. Another surprise is that the proton channel appears to function as a dimer, with two separate conduction pathways. A number of significant biological questions remain in dispute, unanswered, or in some cases, not yet asked. This latter deficit is ascribable to the intrinsic difficulty in evaluating the importance of one component in a complex system, and in addition, to the lack, until recently, of a means of performing an unambiguous lesion experiment, that is, of selectively eliminating the molecule in question. We still lack a potent, selective pharmacological inhibitor, but the identification of the gene has allowed the development of powerful new tools including proton channel antibodies, siRNA and knockout mice.

Voltage-gated proton channels

The history of research on voltage-gated proton channels is recounted, from their proposed existence in dinoflagellates by Hastings in 1972 and their demonstration in snail neurons by Thomas and Meech in 1982 to the discovery in 2006 (after a decade of controversy) of genes that unequivocally code for proton channels. Voltage-gated proton channels are perfectly selective for protons, conduct deuterons half as well, and the conductance is strongly temperature dependent. These properties are consistent with a conduction mechanism involving hydrogen-bonded-chain transfer, in which the selectivity filter is a titratable amino acid residue. Channel opening is regulated stringently by pH such that only outward current is normally activated. Main functions of proton channels include acid extrusion from cells and charge compensation for the electrogenic activity of the phagocyte NADPH oxidase. Genetic approaches hold the promise of rapid progress in the near future.

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of alpha-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 x 5 micros) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2x rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of approximately 16 degrees. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.

How amantadine and rimantadine inhibit proton transport in the M2 protein channel

To understand how antiviral drugs inhibit the replication of influenza A virus via the M2 ion channel, molecular dynamics simulations have been applied to the six possible protonation states of the M2 ion channel in free form and its complexes with two commercial drugs in a fully hydrated lipid bilayer. Among the six different states of free M2 tetramer, water density was present in the pore of the systems with mono-protonated, di-protonated at adjacent position, tri-protonated and tetra-protonated systems. In the presence of inhibitor, water density in the channel was considerably better reduced by rimantadine than amantadine, agreed well with the experimental IC(50) values. With the preferential position and orientation of the two drugs in all states, two mechanisms of action, where the drug binds to the opening pore and the histidine gate, were clearly explained, i.e., (i) inhibitor was detected to localize slightly closer to the histidine gate and can facilitate the orientation of His37 imidazole rings to lie in the close conformation and (ii) inhibitor acts as a blocker, binding at almost above the opening pore and interacts slightly with the three pore-lining residues, Leu26, Ala30 and Ser31. Here, the inhibitors were found to bind very weakly to the channel due to their allosteric hindrance while theirs side chains were strongly solvated.

A secondary gate as a mechanism for inhibition of the M2 proton channel by amantadine

The mechanism of inhibition of the influenza A virus M2 proton channel by the antiviral drug amantadine has been under intense investigation. The importance of a mechanistic understanding is heightened by the prevalence of amantadine-resistant mutations. To gain mechanistic insight at the molecular level, we carried out extensive molecular dynamics simulations of the tetrameric M2 proton channel in both apo and amantadine-bound forms in a lipid bilayer. The simulation of the apo form revealed that Val27 from the four M2 subunits can form a secondary gate near the channel entrance and break the water wire in the channel pore. This gate arises from physical occlusion and the elimination of hydrogen-bonding partners for water molecules. In the presence of amantadine, the secondary gate formed by Val27 and the drug molecule lying just below form an extended blockage, which breaks the water wire throughout the simulation. The location and orientation of amantadine inside of the channel pore as found in our simulation are supported by a host of experimental observations. Our study suggests a novel role for Val27 in the inhibition of the M2 proton channel by amantadine.

Solid-state NMR characterization of conformational plasticity within the transmembrane domain of the influenza A M2 proton channel

Membrane protein function within the membrane interstices is achieved by mechanisms that are not typically available to water-soluble proteins. The whole balance of molecular interactions that stabilize three-dimensional structure in the membrane environment is different from that in an aqueous environment. As a result interhelical interactions are often dominated by non-specific van der Waals interactions permitting dynamics and conformational heterogeneity in these interfaces. Here, solid-state NMR data of the transmembrane domain of the M2 protein from influenza A virus are used to exemplify such conformational plasticity in a tetrameric helical bundle. Such data lead to very high resolution structural restraints that can identify both subtle and substantial structural differences associated with various states of the protein. Spectra from samples using two different preparation protocols, samples prepared in the presence and absence of amantadine, and spectra as a function of pH are used to illustrate conformational plasticity.

Proton transport through influenza A virus M2 protein reconstituted in vesicles

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

Proton transport behavior through the influenza A M2 channel: insights from molecular simulation

The structural properties of the influenza A virus M2 transmembrane channel in dimyristoylphosphatidylcholine bilayer for each of the four protonation states of the proton-gating His-37 tetrad and their effects on proton transport for this low-pH activated, highly proton-selective channel are studied by classical molecular dynamics with the multistate empirical valence-bond (MS-EVB) methodology. The excess proton permeation free energy profile and maximum ion conductance calculated from the MS-EVB simulation data combined with the Poisson-Nernst-Planck theory indicates that the triply protonated His-37 state is the most likely open state via a significant side-chain conformational change of the His-37 tetrad. This proposed open state of M2 has a calculated proton permeation free energy barrier of 7 kcal/mol and a maximum conductance of 53 pS compared to the experimental value of 6 pS. By contrast, the maximum conductance for Na(+) is calculated to be four orders of magnitude lower, in reasonable agreement with the experimentally observed proton selectivity. The pH value to activate the channel opening is estimated to be 5.5 from dielectric continuum theory, which is also consistent with experimental results. This study further reveals that the Ala-29 residue region is the primary binding site for the antiflu drug amantadine (AMT), probably because that domain is relatively spacious and hydrophobic. The presence of AMT is calculated to reduce the proton conductance by 99.8% due to a significant dehydration penalty of the excess proton in the vicinity of the channel-bound AMT.

Baculovirus-expressed influenza vaccine. A novel technology for safe and expeditious vaccine production for human use

Effectiveness of the influenza vaccine in persons with high-risk conditions needs to be improved. In this paper, the authors review various vaccination strategies, including repeated doses of the vaccine or the use of higher hemagglutinin (HA) content vaccines that have been shown to result in improved immunogenicity. A recombinant HA vaccine produced in insect cells using the baculovirus vectors system presents the possibility for safe and expeditious vaccine production. The high purity of the antigen enables administration at much higher doses without a significant increase in side effects in human subjects. An overview of the use of this production system for the development of alternative influenza vaccine targets is also provided, such as neuraminidase and possibly M2. However, the role of M2 may be more appropriate as an adjuvant vaccine in combination with standard HA vaccine supplement and needs further evaluation. The conclusion that the insect cell-baculovirus production technology is a modern solution for rapid viral or parasitic antigen production is made and that this technology is particularly suitable for influenza where annual adjustment of the vaccine is required. In addition, a highly purified recombinant protein vaccine results in an improved influenza vaccine response in those with high-risk medical conditions.