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

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.

Acid activation mechanism of the influenza A M2 proton channel

The homotetrameric influenza A M2 channel (AM2) is an acid-activated proton channel responsible for the acidification of the influenza virus interior, an important step in the viral lifecycle. Four histidine residues (His37) in the center of the channel act as a pH sensor and proton selectivity filter. Despite intense study, the pH-dependent activation mechanism of the AM2 channel has to date not been completely understood at a molecular level. Herein we have used multiscale computer simulations to characterize (with explicit proton transport free energy profiles and their associated calculated conductances) the activation mechanism of AM2. All proton transfer steps involved in proton diffusion through the channel, including the protonation/deprotonation of His37, are explicitly considered using classical, quantum, and reactive molecular dynamics methods. The asymmetry of the proton transport free energy profile under high-pH conditions qualitatively explains the rectification behavior of AM2 (i.e., why the inward proton flux is allowed when the pH is low in viral exterior and high in viral interior, but outward proton flux is prohibited when the pH gradient is reversed). Also, in agreement with electrophysiological results, our simulations indicate that the C-terminal amphipathic helix does not significantly change the proton conduction mechanism in the AM2 transmembrane domain; the four transmembrane helices flanking the channel lumen alone seem to determine the proton conduction mechanism.

Host Cell Targets in HCV Therapy: Novel Strategy or Proven Practice?

The development of novel antiviral drugs against hepatitis C is a challenging and competitive area of research. Progress of this research has been hampered due to the quasispecies nature of the hepatitis C virus, the absence of cellular infection models and the lack of easily accessible and highly representative animal models. The current combination therapy consisting of interferon-α and ribavirin mainly acts by supporting host cell defence. These therapeutics are the prototypic representatives of indirect antiviral agents as they act on cellular targets. However, the therapy is not a cure, when considered from the long-term perspective, for almost half of the chronically infected patients. This draws attention to the urgent need for more efficient treatments. Novel anti-hepatitis C treatments under study are directed against a number of so-called direct antiviral targets such as polymerases and proteases, which are encoded by the virus. Although such direct antiviral approaches have proven to be successful in several viral indications, there is a risk of resistant viruses developing. In order to avoid resistance, the development of indirect antiviral compounds has to be intensified. These act on host cell targets either by boosting the immune response or by blocking the virus host cell interaction. A particularly interesting approach is the development of inhibitors that interfere with signal transduction, such as protein kinase inhibitors. The purpose of this review is to stress the importance of developing indirect antiviral agents that act on host cell targets. In doing so, a large source of potential targets and mechanisms can be exploited, thus increasing the likelihood of success. Ultimately, combination therapies consisting of drugs against direct and indirect viral targets will most probably provide the solution to fighting and eradicating hepatitis C virus in patients.

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.

Effect of cytosolic pH on inward currents reveals structural characteristics of the proton transport cycle in the influenza A protein M2 in cell-free membrane patches of Xenopus oocytes

Transport activity through the mutant D44A of the M2 proton channel from influenza virus A was measured in excised inside-out macro-patches of Xenopus laevis oocytes at cytosolic pH values of 5.5, 7.5 and 8.2. The current-voltage relationships reveal some peculiarities: 1. "Transinhibition", i.e., instead of an increase of unidirectional outward current with increasing cytosolic H(+) concentration, a decrease of unidirectional inward current was found. 2. Strong inward rectification. 3. Exponential rise of current with negative potentials. In order to interpret these findings in molecular terms, different kinetic models have been tested. The transinhibition basically results from a strong binding of H(+) to a site in the pore, presumably His37. This assumption alone already provides inward rectification and exponential rise of the IV curves. However, it results in poor global fits of the IV curves, i.e., good fits were only obtained for cytosolic pH of 8.2, but not for 7.5. Assuming an additional transport step as e.g. caused by a constriction zone at Val27 resulted in a negligible improvement. In contrast, good global fits for cytosolic pH of 7.5 and 8.2 were immediately obtained with a cyclic model. A "recycling step" implies that the protein undergoes conformational changes (assigned to Trp41 and Val27) during transport which have to be reset before the next proton can be transported. The global fit failed at the low currents at pHcyt = 5.5, as expected from the interference of putative transport of other ions besides H(+). Alternatively, a regulatory effect of acidic cytosolic pH may be assumed which strongly modifies the rate constants of the transport cycle.

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.

Voltage-gated proton channels: molecular biology, physiology, and pathophysiology of the H(V) family

Voltage-gated proton channels (H(V)) are unique, in part because the ion they conduct is unique. H(V) channels are perfectly selective for protons and have a very small unitary conductance, both arguably manifestations of the extremely low H(+) concentration in physiological solutions. They open with membrane depolarization, but their voltage dependence is strongly regulated by the pH gradient across the membrane (ΔpH), with the result that in most species they normally conduct only outward current. The H(V) channel protein is strikingly similar to the voltage-sensing domain (VSD, the first four membrane-spanning segments) of voltage-gated K(+) and Na(+) channels. In higher species, H(V) channels exist as dimers in which each protomer has its own conduction pathway, yet gating is cooperative. H(V) channels are phylogenetically diverse, distributed from humans to unicellular marine life, and perhaps even plants. Correspondingly, H(V) functions vary widely as well, from promoting calcification in coccolithophores and triggering bioluminescent flashes in dinoflagellates to facilitating killing bacteria, airway pH regulation, basophil histamine release, sperm maturation, and B lymphocyte responses in humans. Recent evidence that hH(V)1 may exacerbate breast cancer metastasis and cerebral damage from ischemic stroke highlights the rapidly expanding recognition of the clinical importance of hH(V)1.

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.

Mechanistic insight into the h(2)o/d (2)o isotope effect in the proton transport of the influenza virus m2 protein

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

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.

A theory for the proton transport of the influenza virus M2 protein: extensive test against conductance data

A theory for calculating the proton flux through the influenza virus M2 channel is tested here against an extensive set of conductance data. The flux is determined by the rate constants for binding to the His(37) tetrad from the two sides of the membrane and the corresponding unbinding rate constants. The rate constants are calculated by explicitly treating the structure and dynamics of the protein. Important features revealed by previous studies, such as a gating role for Val(27) at the entrance to the channel pore, and channel activation by viral exterior pH, are incorporated in this theory. This study demonstrates that the conductance function of the M2 proton channel can now be quantitatively rationalized by the structure and dynamics of the protein.

Equivalence of two approaches for modeling ion permeation through a transmembrane channel with an internal binding site

Ion permeation through transmembrane channels has traditionally been modeled using two different approaches. In one approach, the translocation of the permeant ion through the channel pore is modeled as continuous diffusion and the rate of ion transport is obtained from solving the steady-state diffusion equation. In the other approach, the translocation of the permeant ion through the pore is modeled as hopping along a discrete set of internal binding sites and the rate of ion transport is obtained from solving a set of steady-state rate equations. In a recent work [Zhou, J. Phys. Chem. Lett. 1, 1973 (2010)], the rate constants for binding to an internal site were further calculated by modeling binding as diffusion-influenced reactions. That work provided the foundation for bridging the two approaches. Here we show that, by representing a binding site as an energy well, the two approaches indeed give the same result for the rate of ion transport.

A pH-dependent conformational ensemble mediates proton transport through the influenza A/M2 protein

The influenza A/M2 protein exhibits inwardly rectifying, pH-activated proton transport that saturates at low pH. A comparison of high-resolution structures of the transmembrane domain at high and low pH suggests that pH-dependent conformational changes may facilitate proton conduction by alternately changing the accessibility of the N-terminal and C-terminal regions of the channel as a proton transits through the transmembrane domain. Here, we show that M2 functionally reconstituted in liposomes populates at least three different conformational states over a physiologically relevant pH range, with transition midpoints that are consistent with previously reported His37 pK(a) values. We then develop and test two similar, quantitative mechanistic models of proton transport, where protonation shifts the equilibrium between structural states having different proton affinities and solvent accessibilities. The models account well for a collection of experimental data sets over a wide range of pH values and voltages and require only a small number of adjustable parameters to accurately describe the data. While the kinetic models do not require any specific conformation for the protein, they nevertheless are consistent with a large body of structural information based on high-resolution nuclear magnetic resonance and crystallographic structures, optical spectroscopy, and molecular dynamics calculations.

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.

Rate theories for biologists

Some of the rate theories that are most useful for modeling biological processes are reviewed. By delving into some of the details and subtleties in the development of the theories, the review will hopefully help the reader gain a more than superficial perspective. Examples are presented to illustrate how rate theories can be used to generate insight at the microscopic level into biomolecular behaviors. An attempt is made to clear up a number of misconceptions in the literature regarding popular rate theories, including the appearance of Planck's constant in the transition-state theory and the Smoluchowski result as an upper limit for protein-protein and protein-DNA association rate constants. Future work in combining the implementation of rate theories through computer simulations with experimental probes of rate processes, and in modeling effects of intracellular environments so that theories can be used for generating rate constants for systems biology studies is particularly exciting.

Influenza virus morphogenesis and budding

Influenza viruses are enveloped, negative stranded, segmented RNA viruses belonging to Orthomyxoviridae family. Each virion consists of three major sub-viral components, namely (i) a viral envelope decorated with three transmembrane proteins hemagglutinin (HA), neuraminidase (NA) and M2, (ii) an intermediate layer of matrix protein (M1), and (iii) an innermost helical viral ribonucleocapsid [vRNP] core formed by nucleoprotein (NP) and negative strand viral RNA (vRNA). Since complete virus particles are not found inside the cell, the processes of assembly, morphogenesis, budding and release of progeny virus particles at the plasma membrane of the infected cells are critically important for the production of infectious virions and pathogenesis of influenza viruses as well. Morphogenesis and budding require that all virus components must be brought to the budding site which is the apical plasma membrane in polarized epithelial cells whether in vitro cultured cells or in vivo infected animals. HA and NA forming the outer spikes on the viral envelope possess apical sorting signals and use exocytic pathways and lipid rafts for cell surface transport and apical sorting. NP also has apical determinant(s) and is probably transported to the apical budding site similarly via lipid rafts and/or through cortical actin microfilaments. M1 binds the NP and the exposed RNAs of vRNPs, as well as to the cytoplasmic tails (CT) and transmembrane (TM) domains of HA, NA and M2, and is likely brought to the budding site on the piggy-back of vRNP and transmembrane proteins. Budding processes involve bud initiation, bud growth and bud release. The presence of lipid rafts and assembly of viral components at the budding site can cause asymmetry of lipid bilayers and outward membrane bending leading to bud initiation and bud growth. Bud release requires fusion of the apposing viral and cellular membranes and scission of the virus buds from the infected cellular membrane. The processes involved in bud initiation, bud growth and bud scission/release require involvement both viral and host components and can affect bud closing and virus release in both positive and negative ways. Among the viral components, M1, M2 and NA play important roles in bud release and M1, M2 and NA mutations all affect the morphology of buds and released viruses. Disassembly of host cortical actin microfilaments at the pinching-off site appears to facilitate bud fission and release. Bud scission is energy dependent and only a small fraction of virus buds present on the cell surface is released. Discontinuity of M1 layer underneath the lipid bilayer, absence of outer membrane spikes, absence of lipid rafts in the lipid bilayer, as well as possible presence of M2 and disassembly of cortical actin microfilaments at the pinching-off site appear to facilitate bud fission and bud release. We provide our current understanding of these important processes leading to the production of infectious influenza virus particles.

Self-assembly and function of primitive cell membranes

We describe possible pathways for separating amphiphilic molecules from organic material on the early earth to form membrane-bound structures required for the start of cellular life. We review properties of the first membranes and their function as permeability barriers. Finally, we discuss the emergence of protein-mediated ion transport across membranes, which facilitated many other cellular functions.

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.

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.

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.

19F NMR detection of the complex between amantadine and the receptor portion of the influenza A M2 ion channel in DPC micelles

(19)F NMR probes were used to follow interactions between ligands in the aminoadamantane series, amantadine (Am) 1 and 3-F-Am 2, and the 5-F-Trp20 transmembrane fragment of the influenza A M2 proton channel (F-M2TM 3) in dodecylphosphocholine micelles over the pH range 5-8. Above pH 7, when the peptide adopts a tetrameric state that is able to bind channel blocking ligands, (19)F-Trp signals from both the free and bound states of the M2TM tetramer are resolved. This differentiation of bound and unbound states of the M2TM receptor by (19)F NMR may provide a system for SAR studies.

Controlling influenza virus replication by inhibiting its proton channel

Influenza A and viruses encode minimalistic proton-selective ion channels known as A/M2 and BM2, respectively. The A/M2 channel is the target of the antiviral drug, amantadine. The structural and mechanistic aspects of proton conductance of the viral ion channels are described and the review makes a case for the development of more effective antivirals.

Histidines, heart of the hydrogen ion channel from influenza A virus: toward an understanding of conductance and proton selectivity

The heart of the H+ conductance mechanism in the homotetrameric M2 H+ channel from influenza A is a set of four histidine side chains. Here, we show that protonation of the third of these imidazoles coincides with acid activation of this transmembrane channel and that, at physiological pH, the channel is closed by two imidazole-imidazolium dimers, each sharing a low-barrier hydrogen bond. This unique construct succeeds in distributing a pair of charges over four rings and many atoms in a low dielectric environment to minimize charge repulsion. These dimers form with identical pKas of 8.2 +/- 0.2, suggesting cooperative H+ binding and clearly illustrating high H+ affinity for this channel. The protonation behavior of the histidine side chains has been characterized by using solid-state NMR spectroscopy on the M2 transmembrane domain in fully hydrated lipid bilayers where the tetrameric backbone structure is known. Furthermore, electrophysiological measurements of multichannel and single-channel experiments confirm that these protein constructs are functional.

Influenza virus proton channels

The M2 ion channel proteins of influenza A and B viruses are essential to viral replication. The two ion channels share a common motif, HXXXW, that is responsible for proton selectivity and activation. The ion channel for the influenza A virus, but not influenza B virus, is inhibited by the antiviral drug amantadine and amantadine-resistant escape mutants form in treated influenza A patients. The studies reviewed suggest that an antiviral compound directed against the conserved motif would be more useful than amantadine in inhibiting viral replication.

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

The tetrameric M2 protein from influenza A is one of the simplest pH-gated H+ channels known, offering the potential of structurally characterizing its gating mechanism. Since the only ionizable groups in the pore are four histidines, we investigated the stability and dynamics of all six possible protonation states of the protein by using molecular dynamics. We show that while all channel protonation states are surprisingly stable, only systems with two or more charged histidines are appreciably conductive. The structural switch, from a uniprotonated to a biprotonated channel, causes an electrostatic repulsion between the charged histidines that pushes the helices apart. This results in the formation of a continuous water file that conducts protons via a H+ wire. pKa calculations place this transition at a pH of 5.6, in remarkable agreement with the experimental value. Since the conversion from uniprotonation to biprotonation occurs during endosome acidification, this explains how M2 is activated in vivo.

Protein dynamics control proton transfer from bulk solvent to protein interior: a case study with a green fluorescent protein

The kinetics of proton transfer in Green Fluorescent Protein (GFP) have been studied as a model system for characterizing the correlation between dynamics and function of proteins in general. The kinetics in EGFP (a variant of GFP) were monitored by using a laser-induced pH jump method. The pH was jumped from 8 to 5 by nanosecond flash photolysis of the "caged proton," o-nitrobenzaldehyde, and subsequent proton transfer was monitored by following the decrease in fluorescence intensity. The modulation of proton transfer kinetics by external perturbants such as viscosity, pH, and subdenaturing concentrations of GdnHCl as well as of salts was studied. The rate of proton transfer was inversely proportional to solvent viscosity, suggesting that the rate-limiting step is the transfer of protons through the protein matrix. The rate is accelerated at lower pH values, and measurements of the fluorescence properties of tryptophan 57 suggest that the enhancement in rate is associated with an enhancement in protein dynamics. The rate of proton transfer is nearly independent of temperature, unlike the rate of the reverse process. When the stability of the protein was either decreased or increased by the addition of co-solutes, including the salts KCl, KNO(3), and K(2)SO(4), a significant decrease in the rate of proton transfer was observed in all cases. The lack of correlation between the rate of proton transfer and the stability of the protein suggests that the structure is tuned to ensure maximum efficiency of the dynamics that control the proton transfer function of the protein.

Proton conductance of influenza virus M2 protein in planar lipid bilayers

Purified M2 protein from the Udorn strain of influenza virus was reconstituted into planar lipid bilayers from liposomes. In 1 mM HCl, the single-channel conductance was measured as 6 pS with open probability of < or =0.03. The current voltage curve is linear over the achievable voltage range. The current amplitude is amantadine sensitive. In HCl solutions, the single-channel current was essentially invariant with changes in [Cl(-)], [Na(+)], and [tetraethylammonium] ([TEA(+)]), but dependent on [H(+)]. The reversal potential, determined with asymmetrical hydrogen chloride solution, is very close to the equilibrium potential of hydrogen. This appears to be the first report of single-channel proton currents with the full-length M2 protein.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Assembly and budding of influenza virus

Influenza viruses are causative agents of an acute febrile respiratory disease called influenza (commonly known as "flu") and belong to the Orthomyxoviridae family. These viruses possess segmented, negative stranded RNA genomes (vRNA) and are enveloped, usually spherical and bud from the plasma membrane (more specifically, the apical plasma membrane of polarized epithelial cells). Complete virus particles, therefore, are not found inside infected cells. Virus particles consist of three major subviral components, namely the viral envelope, matrix protein (M1), and core (viral ribonucleocapsid [vRNP]). The viral envelope surrounding the vRNP consists of a lipid bilayer containing spikes composed of viral glycoproteins (HA, NA, and M2) on the outer side and M1 on the inner side. Viral lipids, derived from the host plasma membrane, are selectively enriched in cholesterol and glycosphingolipids. M1 forms the bridge between the viral envelope and the core. The viral core consists of helical vRNP containing vRNA (minus strand) and NP along with minor amounts of NEP and polymerase complex (PA, PB1, and PB2). For viral morphogenesis to occur, all three viral components, namely the viral envelope (containing lipids and transmembrane proteins), M1, and the vRNP must be brought to the assembly site, i.e. the apical plasma membrane in polarized epithelial cells. Finally, buds must be formed at the assembly site and virus particles released with the closure of buds. Transmembrane viral proteins are transported to the assembly site on the plasma membrane via the exocytic pathway. Both HA and NA possess apical sorting signals and use lipid rafts for cell surface transport and apical sorting. These lipid rafts are enriched in cholesterol, glycosphingolipids and are relatively resistant to neutral detergent extraction at low temperature. M1 is synthesized on free cytosolic polyribosomes. vRNPs are made inside the host nucleus and are exported into the cytoplasm through the nuclear pore with the help of M1 and NEP. How M1 and vRNPs are directed to the assembly site on the plasma membrane remains unclear. The likely possibilities are that they use a piggy-back mechanism on viral glycoproteins or cytoskeletal elements. Alternatively, they may possess apical determinants or diffuse to the assembly site, or a combination of these pathways. Interactions of M1 with M1, M1 with vRNP, and M1 with HA and NA facilitate concentration of viral components and exclusion of host proteins from the budding site. M1 interacts with the cytoplasmic tail (CT) and transmembrane domain (TMD) of glycoproteins, and thereby functions as a bridge between the viral envelope and vRNP. Lipid rafts function as microdomains for concentrating viral glycoproteins and may serve as a platform for virus budding. Virus bud formation requires membrane bending at the budding site. A combination of factors including concentration of and interaction among viral components, increased viscosity and asymmetry of the lipid bilayer of the lipid raft as well as pulling and pushing forces of viral and host components are likely to cause outward curvature of the plasma membrane at the assembly site leading to bud formation. Eventually, virus release requires completion of the bud due to fusion of the apposing membranes, leading to the closure of the bud, separation of the virus particle from the host plasma membrane and release of the virus particle into the extracellular environment. Among the viral components, M1 contains an L domain motif and plays a critical role in budding. Bud completion requires not only viral components but also host components. However, how host components facilitate bud completion remains unclear. In addition to bud completion, influenza virus requires NA to release virus particles from sialic acid residues on the cell surface and spread from cell to cell. Elucidation of both viral and host factors involved in viral morphogenesis and budding may lead to the development of drugs interfering with the steps of viral morphogenesis and in disease progression.