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

The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands

Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

Real-time tracking of drug binding to influenza A M2 reveals a high energy barrier

The drug Rimantadine binds to two different sites in the M2 protein from influenza A, a peripheral site and a pore site that is the primary site of efficacy. It remained enigmatic that pore binding did not occur in certain detergent micelles, and in particular incomplete binding was observed in a mixture of lipids selected to match the viral membrane. Here we show that two effects are responsible, namely changes in the protein upon pore binding that prevented detergent solubilization, and slow binding kinetics in the lipid samples. Using 55-100 kHz magic-angle spinning NMR, we characterize kinetics of drug binding in three different lipid environments: DPhPC, DPhPC with cholesterol and viral mimetic membrane lipid bilayers. Slow pharmacological binding kinetics allowed the characterization of spectral changes associated with non-specific binding to the protein periphery in the kinetically trapped pore-apo state. Resonance assignments were determined from a set of proton-detected 3D spectra. Chemical shift changes associated with functional binding in the pore of M2 were tracked in real time in order to estimate the activation energy. The binding kinetics are affected by pH and the lipid environment and in particular cholesterol. We found that the imidazole-imidazole hydrogen bond at residue histidine 37 is a stable feature of the protein across several lipid compositions. Pore binding breaks the imidazole-imidazole hydrogen bond and limits solubilization in DHPC detergent.

In Situ Spectroscopic Detection of Large-Scale Reorientations of Transmembrane Helices During Influenza A M2 Channel Opening

Viroporins are small ion channels in membranes of enveloped viruses that play key roles during viral life cycles. To use viroporins as drug targets against viral infection requires in-depth mechanistic understanding and, with that, methods that enable investigations under in situ conditions. Here, we apply surface-enhanced infrared absorption (SEIRA) spectroscopy to Influenza A M2 reconstituted within a solid-supported membrane, to shed light on the mechanics of its viroporin function. M2 is a paradigm of pH-activated proton channels and controls the proton flux into the viral interior during viral infection. We use SEIRA to track the large-scale reorientation of M2's transmembrane α-helices in situ during pH-activated channel opening. We quantify this event as a helical tilt from 26° to 40° by correlating the experimental results with solid-state nuclear magnetic resonance-informed computational spectroscopy. This mechanical motion is impeded upon addition of the inhibitor rimantadine, giving a direct spectroscopic marker to test antiviral activity. The presented approach provides a spectroscopic tool to quantify large-scale structural changes and to track the function and inhibition of the growing number of viroporins from pathogenic viruses in future studies.

pH-dependent endocytosis mechanisms for influenza A and SARS-coronavirus

The ongoing SARS-CoV-2 pandemic and the influenza epidemics have revived the interest in understanding how these highly contagious enveloped viruses respond to alterations in the physicochemical properties of their microenvironment. By understanding the mechanisms and conditions by which viruses exploit the pH environment of the host cell during endocytosis, we can gain a better understanding of how they respond to pH-regulated anti-viral therapies but also pH-induced changes in extracellular environments. This review provides a detailed explanation of the pH-dependent viral structural changes preceding and initiating viral disassembly during endocytosis for influenza A (IAV) and SARS coronaviruses. Drawing upon extensive literature from the last few decades and latest research, I analyze and compare the circumstances in which IAV and SARS-coronavirus can undertake endocytotic pathways that are pH-dependent. While there are similarities in the pH-regulated patterns leading to fusion, the mechanisms and pH activation differ. In terms of fusion activity, the measured activation pH values for IAV, across all subtypes and species, vary between approximately 5.0 to 6.0, while SARS-coronavirus necessitates a lower pH of 6.0 or less. The main difference between the pH-dependent endocytic pathways is that the SARS-coronavirus, unlike IAV, require the presence of specific pH-sensitive enzymes (cathepsin L) during endosomal transport. Conversely, the conformational changes in the IAV virus under acidic conditions in endosomes occur due to the specific envelope glycoprotein residues and envelope protein ion channels (viroporins) getting protonated by H+ ions. Despite extensive research over several decades, comprehending the pH-triggered conformational alterations of viruses still poses a significant challenge. The precise mechanisms of protonation mechanisms of certain during endosomal transport for both viruses remain incompletely understood. In absence of evidence, further research is needed.

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 μs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 μs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

Amino acid deprotonation rates from classical force fields

Acid ionization constants (pKa's) of titratable amino acid side chains have received a large amount of experimental and theoretical attention. In many situations, however, the rates of protonation and deprotonation, kon & koff, may also be important, for example in understanding the mechanism of action of proton channels or membrane proteins that couple proton transport to other processes. Protonation and deprotonation involve making and breaking of covalent bonds, which cannot be studied by classical force fields. However, environment effects on the rates should be captured by such methods. Here we present an approach for estimating deprotonation rates based on Warshel's extension of Marcus's theory of electron transfer, with input from molecular simulations. The missing bond dissociation energy is represented by a constant term determined by fitting the pKa value in solution. The statistics of the energy gap between protonated and deprotonated states is used to compute free energy curves of the two states and thus free energy barriers, from which the rate can be deduced. The method is applied to Glu, Asp, and His in bulk solution and select membrane proteins: the M2 proton channel, bacteriorhodopsin, and cytochrome c oxidase.

Clustering of tetrameric influenza M2 peptides in lipid bilayers investigated by 19F solid-state NMR

The influenza M2 protein forms a drug-targeted tetrameric proton channel to mediate virus uncoating, and carries out membrane scission to enable virus release. While the proton channel function of M2 has been extensively studied, the mechanism by which M2 catalyzes membrane scission is still not well understood. Previous fluorescence and electron microscopy studies indicated that M2 tetramers concentrate at the neck of the budding virus in the host plasma membrane. However, molecular evidence for this clustering is scarce. Here, we use ¹⁹F solid-state NMR to investigate M2 clustering in phospholipid bilayers. By mixing equimolar amounts of 4F-Phe47 labeled M2 peptide and CF₃-Phe47 labeled M2 peptide and measuring F-CF₃ cross peaks in 2D ¹⁹F¹⁹F correlation spectra, we show that M2 tetramers form nanometer-scale clusters in lipid bilayers. This clustering is stronger in cholesterol-containing membranes and phosphatidylethanolamine (PE) membranes than in cholesterol-free phosphatidylcholine and phosphatidylglycerol membranes. The observed correlation peaks indicate that Phe47 sidechains from different tetramers are less than ~2 nm apart. ¹H¹⁹F correlation peaks between lipid chain protons and fluorinated Phe47 indicate that Phe47 is more deeply inserted into the lipid bilayer in the presence of cholesterol than in its absence, suggesting that Phe47 preferentially interacts with cholesterol. Static ³¹P NMR spectra indicate that M2 induces negative Gaussian curvature in the PE membrane. These results suggest that M2 tetramers cluster at cholesterol- and PE-rich regions of cell membranes to cause membrane curvature, which in turn can facilitate membrane scission in the last step of virus budding and release.

Emulating Membrane Protein Environments─How Much Lipid Is Required for a Native Structure: Influenza S31N M2

This report investigates the homotetrameric membrane protein structure of the S31N M2 protein from Influenza A virus in the presence of a high molar ratio of lipid. The structured regions of this protein include a single transmembrane helix and an amphipathic helix. Two structures of the S31N M2 conductance domain from Influenza A virus have been deposited in the Protein Data Bank (PDB). These structures present different symmetries about the channel main axis. We present new magic angle spinning and oriented sample solid-state NMR spectroscopic data for S31N M2 in liquid crystalline lipid bilayers using protein tetramer:lipid molar ratios ranging from 1:120 to 1:240. The data is consistent with an essentially 4-fold-symmetric structure very similar to the M2 WT structure that also has a single conformation for the four monomers, except at the His37 and Trp41 functional sites when characterized in samples with a high molar ratio of lipid. While detergent solubilization is well recognized today as a nonideal environment for small membrane proteins, here we discuss the influence of a high lipid to protein ratio for samples of the S31N M2 protein to stabilize an essentially 4-fold-symmetric conformation of the M2 membrane protein. While it is generally accepted that the chemical and physical properties of the native environment of membrane proteins needs to be reproduced judiciously to achieve the native protein structure, here we show that not only the character of the emulated membrane environment is important but also the abundance of the environment is important for achieving the native structure. This is a critical finding as a membrane protein spectroscopist's goal is always to generate a sample with the highest possible protein sensitivity while obtaining spectra of the native-like structure.

Multiscale Simulation of an Influenza A M2 Channel Mutant Reveals Key Features of Its Markedly Different Proton Transport Behavior

The influenza A M2 channel, a prototype for viroporins, is an acid-activated viroporin that conducts protons across the viral membrane, a critical step in the viral life cycle. Four central His37 residues control channel activation by binding subsequent protons from the viral exterior, which opens the Trp41 gate and allows proton flux to the interior. Asp44 is essential for maintaining the Trp41 gate in a closed state at high pH, resulting in asymmetric conduction. The prevalent D44N mutant disrupts this gate and opens the C-terminal end of the channel, resulting in increased conduction and a loss of this asymmetric conduction. Here, we use extensive Multiscale Reactive Molecular Dynamics (MS-RMD) and quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations with an explicit, reactive excess proton to calculate the free energy of proton transport in this M2 mutant and to study the dynamic molecular-level behavior of D44N M2. We find that this mutation significantly lowers the barrier of His37 deprotonation in the activated state and shifts the barrier for entry to the Val27 tetrad. These free energy changes are reflected in structural shifts. Additionally, we show that the increased hydration around the His37 tetrad diminishes the effect of the His37 charge on the channel's water structure, facilitating proton transport and enabling activation from the viral interior. Altogether, this work provides key insight into the fundamental characteristics of PT in WT M2 and how the D44N mutation alters this PT mechanism, and it expands understanding of the role of emergent mutations in viroporins.

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

With an estimated 1 billion people affected across the globe, influenza is one of the most serious health concerns worldwide. Therapeutic treatments have encompassed a number of key functional viral proteins, mainly focused on the M2 proton channel and neuraminidase. This review highlights the efforts spent in targeting the M2 proton channel, which mediates the proton transport toward the interior of the viral particle as a preliminary step leading to the release of the fusion peptide in hemagglutinin and the fusion of the viral and endosomal membranes. Besides the structural and mechanistic aspects of the M2 proton channel, attention is paid to the challenges posed by the development of efficient small molecule inhibitors and the evolution toward novel ligands and scaffolds motivated by the emergence of resistant strains.

Determination of Histidine Protonation States in Proteins by Fast Magic Angle Spinning NMR

Histidine residues play important structural and functional roles in proteins, such as serving as metal-binding ligands, mediating enzyme catalysis, and modulating proton channel activity. Many of these activities are modulated by the ionization state of the imidazole ring. Here we present a fast MAS NMR approach for the determination of protonation and tautomeric states of His at frequencies of 40-62 kHz. The experiments combine 1H detection with selective magnetization inversion techniques and transferred echo double resonance (TEDOR)-based filters, in 2D heteronuclear correlation experiments. We illustrate this approach using microcrystalline assemblies of HIV-1 CACTD-SP1 protein.

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.

Protonation-Deprotonation Switch Controls the Amyloid-like Misfolding of Nucleic-Acid-Binding Domains of TDP-43

Nutrient starvation stress acidifies the cytosol and leads to the formation of large protein assemblies and misfolded aggregates. However, how starvation stress is sensed at the molecular level and leads to protein misfolding is poorly understood. TDP-43 is a vital protein, which, under stress-like conditions, associates with stress granule proteins via its functional nucleic-acid-binding domains (TDP-43^tRRM) and misfolds to form aberrant aggregates. Here, we show that the monomeric N form of TDP-43^tRRM forms a misfolded amyloid-like protein assembly, β form, in a pH-dependent manner and identified the critical protein side-chain residue whose protonation triggers its misfolding. We systematically mutated the three buried ionizable residues, D105, H166, and H256, to neutral amino acids to block the pH-dependent protonation-deprotonation titration of their side chain and studied their effect on the N-to-β transition. We observed that D105A and H256Q resembled TDP-43^tRRM in their pH-dependent misfolding behavior. However, H166Q retains the N-like secondary structure under low-pH conditions and does not show pH-dependent misfolding to the β form. These results indicate that H166 is the critical side-chain residue whose protonation triggers the misfolding of TDP-43^tRRM and shed light on how stress-induced misfolding of proteins during neurodegeneration could begin from site-specific triggers.

Water orientation and dynamics in the closed and open influenza B virus M2 proton channels

The influenza B M2 protein forms a water-filled tetrameric channel to conduct protons across the lipid membrane. To understand how channel water mediates proton transport, we have investigated the water orientation and dynamics using solid-state NMR spectroscopy and molecular dynamics (MD) simulations. 13C-detected water 1H NMR relaxation times indicate that water has faster rotational motion in the low-pH open channel than in the high-pH closed channel. Despite this faster dynamics, the open-channel water shows higher orientational order, as manifested by larger motionally-averaged 1H chemical shift anisotropies. MD simulations indicate that this order is induced by the cationic proton-selective histidine at low pH. Furthermore, the water network has fewer hydrogen-bonding bottlenecks in the open state than in the closed state. Thus, faster dynamics and higher orientational order of water molecules in the open channel establish the water network structure that is necessary for proton hopping.

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.

Gut bacteria are essential for normal cuticle development in herbivorous turtle ants

Across the evolutionary history of insects, the shift from nitrogen-rich carnivore/omnivore diets to nitrogen-poor herbivorous diets was made possible through symbiosis with microbes. The herbivorous turtle ants Cephalotes possess a conserved gut microbiome which enriches the nutrient composition by recycling nitrogen-rich metabolic waste to increase the production of amino acids. This enrichment is assumed to benefit the host, but we do not know to what extent. To gain insights into nitrogen assimilation in the ant cuticle we use gut bacterial manipulation, ¹⁵N isotopic enrichment, isotope-ratio mass spectrometry, and ¹⁵N nuclear magnetic resonance spectroscopy to demonstrate that gut bacteria contribute to the formation of proteins, catecholamine cross-linkers, and chitin in the cuticle. This study identifies the cuticular components which are nitrogen-enriched by gut bacteria, highlighting the role of symbionts in insect evolution, and provides a framework for understanding the nitrogen flow from nutrients through bacteria into the insect cuticle.

Highlighting membrane protein structure and function: A celebration of the Protein Data Bank

Biological membranes define the boundaries of cells and compartmentalize the chemical and physical processes required for life. Many biological processes are carried out by proteins embedded in or associated with such membranes. Determination of membrane protein (MP) structures at atomic or near-atomic resolution plays a vital role in elucidating their structural and functional impact in biology. This endeavor has determined 1198 unique MP structures as of early 2021. The value of these structures is expanded greatly by deposition of their three-dimensional (3D) coordinates into the Protein Data Bank (PDB) after the first atomic MP structure was elucidated in 1985. Since then, free access to MP structures facilitates broader and deeper understanding of MPs, which provides crucial new insights into their biological functions. Here we highlight the structural and functional biology of representative MPs and landmarks in the evolution of new technologies, with insights into key developments influenced by the PDB in magnifying their impact.

Influenza A M2 Inhibitor Binding Understood through Mechanisms of Excess Proton Stabilization and Channel Dynamics

Prevalent resistance to inhibitors that target the influenza A M2 proton channel has necessitated a continued drug design effort, supported by a sustained study of the mechanism of channel function and inhibition. Recent high-resolution X-ray crystal structures present the first opportunity to see how the adamantyl amine class of inhibitors bind to M2 and disrupt and interact with the channel’s water network, providing insight into the critical properties that enable their effective inhibition in wild-type M2. In this work, we examine the hypothesis that these drugs act primarily as mechanism-based inhibitors by comparing hydrated excess proton stabilization during proton transport in M2 with the interactions revealed in the crystal structures, using the Multiscale Reactive Molecular Dynamics (MS-RMD) methodology. MS-RMD, unlike classical molecular dynamics, models the hydrated proton (hydronium-like cation) as a dynamic excess charge defect and allows bonds to break and form, capturing the intricate interactions between the hydrated excess proton, protein atoms, and water. Through this, we show that the ammonium group of the inhibitors is effectively positioned to take advantage of the channel’s natural ability to stabilize an excess protonic charge and act as a hydronium mimic. Additionally, we show that the channel is especially stable in the drug binding region, highlighting the importance of this property for binding the adamantane group. Finally, we characterize an additional hinge point near Val27, which dynamically responds to charge and inhibitor binding. Altogether, this work further illuminates a dynamic understanding of the mechanism of drug inhibition in M2, grounded in the fundamental properties that enable the channel to transport and stabilize excess protons, with critical implications for future drug design efforts.

Increased Dissociation of Adamantanamines in Influenza A M2 S31N with Partial Block by Rimantadine

The ubiquitous mutation from serine (WT) to asparagine at residue 31 (S31N) in the influenza A M2 channel renders it insensitive to amantadine (AMT) and rimantadine (RMT) block, but it is unknown whether the inhibition results from weak binding or incomplete block. Two-electrode voltage clamp (TEVC) of transfected Xenopus oocytes revealed that the M2 S31N channel is essentially fully blocked by AMT at 10 mM, demonstrating that, albeit weak, AMT binding in a channel results in complete block of its proton current. In contrast, RMT achieves only a modest degree of block in the M2 S31N channel at 1 mM, with very little increase in block at 10 mM, indicating that the RMT binding site in the channel saturates with only modest block. From exponential curve fits to families of proton current wash-in and wash-out traces, the association rate constant (k₁) is somewhat decreased for both AMT and RMT in the S31N, but the dissociation rate constant (k₂) is dramatically increased compared with WT. The potentials of mean force (PMF) from adaptive biasing force (ABF) molecular dynamics simulations predict that rate constants should be exquisitely sensitive to the charge state of the His37 selectivity filter of M2. With one exception out of eight cases, predictions from the simulations with one and three charged side chains bracket the experimental rate constants, as expected for the acidic bath used in the TEVC assay. From simulations, the weak binding can be accounted for by changes in the potentials of mean force, but the partial block by RMT remains unexplained.

Revealing weak histidine 15N homonuclear scalar couplings using Solid-State Magic-Angle-Spinning NMR spectroscopy

The tautomeric structure and chemistry of the histidine imidazole ring play active roles in many structurally and functionally important proteins and polypeptides. While in NMR spectroscopy histidine chemical shifts (e.g. ¹⁵N, ¹³C, and ¹H) have been commonly used to characterize the tautomeric structure, hydrogen bonding, and torsion angles, homonuclear ¹⁵N scalar couplings in histidine have rarely been reported. Here, we propose double spin-echo sequences to compare the observed signals with and without a 90° pulse between the two spin-echo periods, such that their signal ratio as a function of the echo time solely depends on homonuclear scalar couplings, allowing for measuring weak homonuclear scalar couplings without influence from transverse dephasing effects, thus capable of revealing hydrogen-bond mediated ¹⁵N-¹⁵N J-couplings that can provide direct and definitive evidence for the formation of N^…H^…N hydrogen-bonding associated with the imidazole ring. We used two ¹³C,¹⁵N labeled histidine samples recrystallized from solutions at pH 6.3 and pH 11.0 to demonstrate the feasibility of this method and reveal the existence of a weak two-bond scalar coupling between the N_δ1 and N_ε2 sites in the histidine imidazole ring in three tautomeric states and the presence of a hydrogen-bond mediated scalar coupling between the N_δ1 site in the imidazole ring and the backbone N_α site in the histidine neutral τ and π states. Our results demonstrate that weak ¹⁵N homonuclear scalar couplings can be measured even when their values are less than their corresponding intrinsic natural linewidths, thus providing direct and definitive evidence for the formation of N^…H^…N hydrogen bonding that is associated with the histidine imidazole ring.

Site-specific resolution of anionic residues in proteins using solid-state NMR spectroscopy

NMR spectroscopy is commonly used to infer site-specific acid dissociation constants (pK_a) since the chemical shift is sensitive to the protonation state. Methods that probe atoms nearest to the functional groups involved in acid/base chemistry are the most sensitive for determining the protonation state. In this work, we describe a magic-angle-spinning (MAS) solid-state NMR approach to measure chemical shifts on the side chain of the anionic residues aspartate and glutamate. This method involves a combination of double quantum spectroscopy in the indirect dimension and REDOR dephasing to provide a sensitive and resolved view of these amino acid residues that are commonly involved in enzyme catalysis and membrane protein transport. To demonstrate the applicability of the approach, we carried out measurements using a microcrystalline soluble protein (ubiquitin) and a membrane protein embedded in lipid bilayers (EmrE). Overall, the resolution available from the double quantum dimension and confidence in identification of aspartate and glutamate residues from the REDOR filter make this method the most convenient for characterizing protonation states and deriving pK_a values using MAS solid-state NMR.

Imidazole-Imidazole Hydrogen Bonding in the pH-Sensing Histidine Side Chains of Influenza A M2

The arrangement of histidine side chains in influenza A M2 tetramer determines their pKa values, which define pH-controlled proton conduction critical to the virus lifecycle. Both water-associated and hydrogen-bonded imidazole-imidazolium histidine quaternary structures have been proposed, based on crystal structures and NMR chemical shifts, respectively. Here we show, using the conduction domain construct of M2 in lipid bilayers, that the imidazole rings are hydrogen bonded even at a pH of 7.8 in the neutral charge state. An intermolecular 8.9 ± 0.3 Hz 2hJNN hydrogen bond is observed between H37 Nε and Nδ recorded in a fully protonated sample with 100 kHz magic-angle spinning. This interaction could not be detected in the drug-bound sample.

M2 amphipathic helices facilitate pH-dependent conformational transition in influenza A virus

The matrix-2 (M2) protein from influenza A virus is a tetrameric, integral transmembrane (TM) protein that plays a vital role in viral replication by proton flux into the virus. The His37 tetrad is a pH sensor in the center of the M2 TM helix that activates the channel in response to the low endosomal pH. M2 consists of different regions that are believed to be involved in membrane targeting, packaging, nucleocapsid binding, and proton transport. Although M2 has been the target of many experimental and theoretical studies that have led to significant insights into its structure and function under differing conditions, the main mechanism of proton transport, its conformational dynamics, and the role of the amphipathic helices (AHs) on proton conductance remain elusive. To this end, we have applied explicit solvent constant pH molecular dynamics using the multisite λ-dynamics approach (CpHMD^MSλD) to investigate the buried ionizable residues comprehensively and to elucidate their effect on the conformational transition. Our model recapitulates the pH-dependent conformational transition of M2 from closed to open state when the AH domain is included in the M2 construct, revealing the role of the amphipathic helices on this transition and shedding light on the proton-transport mechanism. This work demonstrates the importance of including the amphipathic helices in future experimental and theoretical studies of ion channels. Finally, our work shows that explicit solvent CpHMD^MSλD provides a realistic pH-dependent model for membrane proteins.

Observation of the Imidazole-Imidazolium Hydrogen Bonds Responsible for Selective Proton Conductance in the Influenza A M2 Channel

The integral membrane M2 protein is a 97-residue membrane protein that assembles as a tetramer to conduct protons at a slow rate (10²-10³/s) when activated by low pH. The proton conductance mechanism has been extensively debated in the literature, but it is accepted that the proton conductance is facilitated by hydrogen bonds involving the His37 residues. However, the hydrogen bonding partnership remains unresolved. Here, we report on the measurement of ¹⁵N-¹⁵N J-couplings of ¹⁵N His37-labeled full length M2 (M2FL) protein from Influenza A virus embedded in synthetic liquid crystalline lipid bilayers using two-dimensional J-resolved NMR spectroscopy. We experimentally observed the hydrogen-bond mediated J-couplings between N_δ1 and N_ε2 of adjacent His37 imidazole rings, providing direct evidence for the existence of various imidazolium-imidazole hydrogen-bonding geometries in the histidine tetrad at low pH, thus validating the proton conduction mechanism in the M2FL protein by which the proton is transferred through the breaking and reforming of the hydrogen bonds between pairs of His37 residues.

Proton-Induced Conformational and Hydration Dynamics in the Influenza A M2 Channel

The influenza A M2 protein is an acid-activated proton channel responsible for acidification of the inside of the virus, a critical step in the viral life cycle. This channel has four central histidine residues that form an acid-activated gate, binding protons from the outside until an activated state allows proton transport to the inside. While previous work has focused on proton transport through the channel, the structural and dynamic changes that accompany proton flux and enable activation have yet to be resolved. In this study, extensive Multiscale Reactive Molecular Dynamics simulations with explicit Grotthuss-shuttling hydrated excess protons are used to explore detailed molecular-level interactions that accompany proton transport in the +0, + 1, and +2 histidine charge states. The results demonstrate how the hydrated excess proton strongly influences both the protein and water hydrogen-bonding network throughout the channel, providing further insight into the channel's acid-activation mechanism and rectification behavior. We find that the excess proton dynamically, as a function of location, shifts the protein structure away from its equilibrium distributions uniquely for different pH conditions consistent with acid-activation. The proton distribution in the xy-plane is also shown to be asymmetric about the channel's main axis, which has potentially important implications for the mechanism of proton conduction and future drug design efforts.

X-ray Crystal Structure of the Influenza A M2 Proton Channel S31N Mutant in Two Conformational States: An Open and Shut Case

The amantadine-resistant S31N mutant of the influenza A M2 proton channel has become prevalent in currently circulating viruses. Here, we have solved an X-ray crystal structure of M2(22-46) S31N that contains two distinct conformational states within its asymmetric unit. This structure reveals the mechanism of adamantane resistance in both conformational states of the M2 channel. In the Inward_open conformation, the mutant Asn31 side chain faces the channel pore and sterically blocks the adamantane binding site. In the Inward_closed conformation, Asn31 forms hydrogen bonds with carbonyls at the monomer-monomer interface, which twists the monomer helices and constricts the channel pore at the drug binding site. We also examine M2(19-49) WT and S31N using solution NMR spectroscopy and show that distribution of the two conformational states is dependent on both detergent choice and experimental pH.

Viroporins in the Influenza Virus

Influenza is a highly contagious virus that causes seasonal epidemics and unpredictable pandemics. Four influenza virus types have been identified to date: A, B, C and D, with only A–C known to infect humans. Influenza A and B viruses are responsible for seasonal influenza epidemics in humans and are responsible for up to a billion flu infections annually. The M2 protein is present in all influenza types and belongs to the class of viroporins, i.e., small proteins that form ion channels that increase membrane permeability in virus-infected cells. In influenza A and B, AM2 and BM2 are predominantly proton channels, although they also show some permeability to monovalent cations. By contrast, M2 proteins in influenza C and D, CM2 and DM2, appear to be especially selective for chloride ions, with possibly some permeability to protons. These differences point to different biological roles for M2 in types A and B versus C and D, which is also reflected in their sequences. AM2 is by far the best characterized viroporin, where mechanistic details and rationale of its acid activation, proton selectivity, unidirectionality, and relative low conductance are beginning to be understood. The present review summarizes the biochemical and structural aspects of influenza viroporins and discusses the most relevant aspects of function, inhibition, and interaction with the host.

Elucidating Relayed Proton Transfer through a His-Trp-His Triad of a Transmembrane Proton Channel by Solid-State NMR

Proton transfer through membrane-bound ion channels is mediated by both water and polar residues of proteins, but the detailed molecular mechanism is challenging to determine. The tetrameric influenza A and B virus M2 proteins form canonical proton channels that use an HxxxW motif for proton selectivity and gating. The BM2 channel also contains a second histidine (His), H27, equidistant from the gating tryptophan, which leads to a symmetric H¹⁹xxxW²³xxxH²⁷ motif. The proton-dissociation constants (pK_a's) of H19 in BM2 were found to be much lower than the pK_a's of H37 in AM2. To determine if the lower pK_a's result from H27-facilitated proton dissociation of H19, we have now investigated a H27A mutant of BM2 using solid-state NMR. ¹⁵N NMR spectra indicate that removal of the second histidine converted the protonation and tautomeric equilibria of H19 to be similar to the H37 behavior in AM2, indicating that the peripheral H27 is indeed the origin of the low pK_a's of H19 in wild-type BM2. Measured interhelical distances between W23 sidechains indicate that the pore constriction at W23 increases with the H19 tetrad charge but is independent of the H27A mutation. These results indicate that H27 both accelerates proton dissociation from H19 to increase the inward proton conductance and causes the small reverse conductance of BM2. The proton relay between H19 and H27 is likely mediated by the intervening gating tryptophan through cation-π interactions. This relayed proton transfer may exist in other ion channels and has implications for the design of imidazole-based synthetic proton channels.

Design of a light-gated proton channel based on the crystal structure of Coccomyxa rhodopsin

The light-driven proton pump Coccomyxa subellipsoidea rhodopsin (CsR) provides-because of its high expression in heterologous host cells-an opportunity to study active proton transport under controlled electrochemical conditions. In this study, solving crystal structure of CsR at 2.0-Å resolution enabled us to identify distinct features of the membrane protein that determine ion transport directivity and voltage sensitivity. A specific hydrogen bond between the highly conserved Arg⁸³ and the nearby nonconserved tyrosine (Tyr¹⁴) guided our structure-based transformation of CsR into an operational light-gated proton channel (CySeR) that could potentially be used in optogenetic assays. Time-resolved electrophysiological and spectroscopic measurements distinguished pump currents from channel currents in a single protein and emphasized the necessity of Arg⁸³ mobility in CsR as a dynamic extracellular barrier to prevent passive conductance. Our findings reveal that molecular constraints that distinguish pump from channel currents are structurally more confined than was generally expected. This knowledge might enable the structure-based design of novel optogenetic tools, which derive from microbial pumps and are therefore ion specific.

BEST and SOFAST experiments for resonance assignment of histidine and tyrosine side chains in 13C/15N labeled proteins

Aromatic amino-acid side chains are essential components for the structure and function of proteins. We present herein a set of NMR experiments for time-efficient resonance assignment of histidine and tyrosine side chains in uniformly ¹³C/¹⁵N-labeled proteins. The use of band-selective ¹³C pulses allows to deal with linear chains of coupled spins, thus avoiding signal loss that occurs in branched spin systems during coherence transfer. Furthermore, our pulse schemes make use of longitudinal ¹H relaxation enhancement, Ernst-angle excitation, and simultaneous detection of ¹H and ¹³C steady-state polarization to achieve significant signal enhancements.

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na⁺ and K⁺ ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the K_v1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine’s NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of −70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

Inhibitors of the M2 Proton Channel Engage and Disrupt Transmembrane Networks of Hydrogen-Bonded Waters

Water-mediated interactions play key roles in drug binding. In protein sites with sparse polar functionality, a small-molecule approach is often viewed as insufficient to achieve high affinity and specificity. Here we show that small molecules can enable potent inhibition by targeting key waters. The M2 proton channel of influenza A is the target of the antiviral drugs amantadine and rimantadine. Structural studies of drug binding to the channel using X-ray crystallography have been limited because of the challenging nature of the target, with the one previously solved crystal structure limited to 3.5 Å resolution. Here we describe crystal structures of amantadine bound to M2 in the Inward_closed conformation (2.00 Å), rimantadine bound to M2 in both the Inward_closed (2.00 Å) and Inward_open (2.25 Å) conformations, and a spiro-adamantyl amine inhibitor bound to M2 in the Inward_closed conformation (2.63 Å). These X-ray crystal structures of the M2 proton channel with bound inhibitors reveal that ammonium groups bind to water-lined sites that are hypothesized to stabilize transient hydronium ions formed in the proton-conduction mechanism. Furthermore, the ammonium and adamantyl groups of the adamantyl-amine class of drugs are free to rotate in the channel, minimizing the entropic cost of binding. These drug-bound complexes provide the first high-resolution structures of drugs that interact with and disrupt networks of hydrogen-bonded waters that are widely utilized throughout nature to facilitate proton diffusion within proteins.

Structural and Functional Properties of Viral Membrane Proteins

Viruses have developed a large variety of transmembrane proteins to carry out their infectious cycles. Some of these proteins are simply anchored to membrane via transmembrane helices. Others, however, adopt more interesting structures to perform tasks such as mediating membrane fusion and forming ion-permeating channels. Due to the dynamic or plastic nature shown by many of the viral membrane proteins, structural and mechanistic understanding of these proteins has lagged behind their counterparts in prokaryotes and eukaryotes. This chapter provides an overview of the use of NMR spectroscopy to unveil the transmembrane and membrane-proximal regions of viral membrane proteins, as well as their interactions with potential therapeutics.

Hidden motions and motion-induced invisibility: Dynamics-based spectral editing in solid-state NMR

Solid-state nuclear magnetic resonance (ssNMR) spectroscopy enables the structural characterization of a diverse array of biological assemblies that include amyloid fibrils, non-amyloid aggregates, membrane-associated proteins and viral capsids. Such biological samples feature functionally relevant molecular dynamics, which often affect different parts of the sample in different ways. Solid-state NMR experiments' sensitivity to dynamics represents a double-edged sword. On the one hand, it offers a chance to measure dynamics in great detail. On the other hand, certain types of motion lead to signal loss and experimental inefficiencies that at first glance interfere with the application of ssNMR to overly dynamic proteins. Dynamics-based spectral editing (DYSE) ssNMR methods leverage motion-dependent signal losses to simplify spectra and enable the study of sub-structures with particular motional properties.

Structure and Dynamics of Membrane Proteins from Solid-State NMR

Solid-state nuclear magnetic resonance (SSNMR) spectroscopy elucidates membrane protein structures and dynamics in atomic detail to yield mechanistic insights. By interrogating membrane proteins in phospholipid bilayers that closely resemble biological membranes, SSNMR spectroscopists have revealed ion conduction mechanisms, substrate transport dynamics, and oligomeric interfaces of seven-transmembrane helix proteins. Research has also identified conformational plasticity underlying virus-cell membrane fusions by complex protein machineries, and β-sheet folding and assembly by amyloidogenic proteins bound to lipid membranes. These studies collectively show that membrane proteins exhibit extensive structural plasticity to carry out their functions. Because of the inherent dependence of NMR frequencies on molecular orientations and the sensitivity of NMR frequencies to dynamical processes on timescales from nanoseconds to seconds, SSNMR spectroscopy is ideally suited to elucidate such structural plasticity, local and global conformational dynamics, protein-lipid and protein-ligand interactions, and protonation states of polar residues. New sensitivity-enhancement techniques, resolution enhancement by ultrahigh magnetic fields, and the advent of 3D and 4D correlation NMR techniques are increasingly aiding these mechanistically important structural studies.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Transport-Relevant Protein Conformational Dynamics and Water Dynamics on Multiple Time Scales in an Archetypal Proton Channel: Insights from Solid-State NMR

The influenza M2 protein forms a tetrameric proton channel that conducts protons from the acidic endosome into the virion by shuttling protons between water and a transmembrane histidine. Previous NMR studies have shown that this histidine protonates and deprotonates on the microsecond time scale. However, M2's proton conduction rate is 10-1000 s^-1, more than 2 orders of magnitude slower than the histidine-water proton-exchange rate. M2 is also known to be conformationally plastic. To address the disparity between the functional time scale and the time scales of protein conformational dynamics and water dynamics, we have now investigated a W41F mutant of the M2 transmembrane domain using solid-state NMR. ¹³C chemical shifts of the membrane-bound peptide indicate the presence of two distinct tetramer conformations, whose concentrations depend exclusively on pH and hence the charge-state distribution of the tetramers. High-temperature 2D correlation spectra indicate that these two conformations interconvert at a rate of ∼400 s^-1 when the +2 and +3 charge states dominate, which gives the first experimental evidence of protein conformational motion on the transport time scale. Protein ¹³C-detected water ¹H T₂ relaxation measurements show that channel water relaxes an order of magnitude faster than bulk water and membrane-associated water, indicating that channel water undergoes nanosecond motion in a pH-independent fashion. These results connect motions on three time scales to explain M2's proton-conduction mechanism: picosecond-to-nanosecond motions of water molecules facilitate proton Grotthuss hopping, microsecond motions of the histidine side chain allow water-histidine proton transfer, while millisecond motions of the entire four-helix bundle constitute the rate-limiting step, dictating the number of protons released into the virion.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Entropy and Polarity Control the Partition and Transportation of Drug-like Molecules in Biological Membrane

Partition and transportation of drug in the plasma membrane of a mammalian cell are the prerequisite for its function on target protein. Therefore, comprehensive understanding of the physicochemical properties and mechanism behind these complex phenomena is crucial in pharmaceutical research. By using the state-of-art molecular simulations with polarization effect implicitly or explicitly included, we studied the permeation behavior of 2-aminoethoxydiphenyl borate (2-APB), a broad-spectrum modulator for a number of membrane proteins. We showed that the protonation state and therefore the polarity of the drug is critical for its partition, and that the drug is likely to switch between different protonation states along its permeation pathway. By changing the degrees of freedom, protonation further affects the thermodynamic of the permeation pathway of 2-APB, leading to different entropic contributions. A survey on 54 analog structures with similar backbone to 2-APB showed that delicate balance between entropy and polarity plays an important role in drugs’ potency.

Functional studies reveal the similarities and differences between AM2 and BM2 proton channels from influenza viruses

AM2 and BM2 proton channels are attractive antiviral drug targets due to their essential roles during influenza virus replication. Although both AM2 and BM2 are proton-selective ion channels, they share little sequence similarity except for the HXXXW sequence, which suggests that their proton conductance properties might differ. To test this hypothesis, we applied two-electrode voltage clamp electrophysiological assays to study the specific conductance, leakage current, channel activation, and inhibition of AM2 and BM2 proton channels. It was found that BM2 channel has a higher specific conductance than AM2 channel at pH5.5. Unlike AM2 channel, whose proton conductance is asymmetric (from viral exterior to interior), BM2 channel is capable of conducting proton in both directions. Moreover, BM2 requires a more acidic pH for channel activation than AM2, as revealed by its lower pK_a values. Finally, both AM2 and BM2 can be inhibited by Cu(II) and Cu(I). Overall, the results from this side-by-side comparison of AM2 and BM2 channels reveal the structure-function relationships of these two viroporins, and such information might be important for the designing of novel ion channels.

XFEL structures of the influenza M2 proton channel: Room temperature water networks and insights into proton conduction

The M2 proton channel of influenza A is a drug target that is essential for the reproduction of the flu virus. It is also a model system for the study of selective, unidirectional proton transport across a membrane. Ordered water molecules arranged in "wires" inside the channel pore have been proposed to play a role in both the conduction of protons to the four gating His37 residues and the stabilization of multiple positive charges within the channel. To visualize the solvent in the pore of the channel at room temperature while minimizing the effects of radiation damage, data were collected to a resolution of 1.4 Å using an X-ray free-electron laser (XFEL) at three different pH conditions: pH 5.5, pH 6.5, and pH 8.0. Data were collected on the Inward_open state, which is an intermediate that accumulates at high protonation of the His37 tetrad. At pH 5.5, a continuous hydrogen-bonded network of water molecules spans the vertical length of the channel, consistent with a Grotthuss mechanism model for proton transport to the His37 tetrad. This ordered solvent at pH 5.5 could act to stabilize the positive charges that build up on the gating His37 tetrad during the proton conduction cycle. The number of ordered pore waters decreases at pH 6.5 and 8.0, where the Inward_open state is less stable. These studies provide a graphical view of the response of water to a change in charge within a restricted channel environment.

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

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

Structural Basis for Asymmetric Conductance of the Influenza M2 Proton Channel Investigated by Solid-State NMR Spectroscopy

The influenza M2 protein forms an acid-activated proton channel that is essential for virus replication. The transmembrane H37 selects for protons under low external pH while W41 ensures proton conduction only from the N terminus to the C terminus and prevents reverse current under low internal pH. Here, we address the molecular basis for this asymmetric conduction by investigating the structure and dynamics of a mutant channel, W41F, which permits reverse current under low internal pH. Solid-state NMR experiments show that W41F M2 retains the pH-dependent α-helical conformations and tetrameric structure of the wild-type (WT) channel but has significantly altered protonation and tautomeric equilibria at H37. At high pH, the H37 structure is shifted toward the π tautomer and less cationic tetrads, consistent with faster forward deprotonation to the C terminus. At low pH, the mutant channel contains more cationic tetrads than the WT channel, consistent with faster reverse protonation from the C terminus. ¹⁵N NMR spectra allow the extraction of four H37 pK_as and show that the pK_as are more clustered in the mutant channel compared to WT M2. Moreover, binding of the antiviral drug, amantadine, at the N-terminal pore at low pH did not convert all histidines to the neutral state, as seen in WT M2, but left half of all histidines cationic, unambiguously demonstrating C-terminal protonation of H37 in the mutant. These results indicate that asymmetric conduction in WT M2 is due to W41 inhibition of C-terminal acid activation by H37. When Trp is replaced by Phe, protons can be transferred to H37 bidirectionally with distinct rate constants.

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.

Beyond Structural Biology to Functional Biology: Solid-State NMR Experiments and Strategies for Understanding the M2 Proton Channel Conductance

In terms of structural biology, solid-state NMR experiments and strategies have been well established for resonance assignments, leading to the determination of three-dimensional structures of insoluble membrane proteins in their native-like environment. It is also known that NMR has the unique capabilities to characterize structure-function relationships of membrane-bound biological systems beyond structural biology. Here, we report on solid-state NMR experiments and strategies for extracting functional activities on a sub-millisecond time scale. Specifically, we use the His37-labeled full length M2 (M2FL) protein of the Influenza A virus embedded in synthetic lipid bilayers as an example to characterize the proton conduction mechanism and kinetics. The integral membrane M2 protein assembles as a tetrameric bundle to form a proton-conducting channel that is activated by low pH and is essential for the viral lifecycle. Our results present convincing evidence for the formation of imidazolium-imidazole hydrogen bonds in the His37 tetrad at low pH and that these hydrogen bonds have a low barrier that facilitates the proton conduction mechanism in the M2FL protein. Moreover, it has been possible to measure hydronium ion exchange between water and the protons in the His37 NH bonds based on chemical exchange spectroscopy with minimized spin diffusion. The results identify an exchange rate constant of ∼4000 s^-1 for pH 5.8 at -10 °C.