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

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.

A Study of the Activity of Adamantyl Amines against Mutant Influenza A M2 Channels Identified a Polycyclic Cage Amine Triple Blocker, Explored by Molecular Dynamics Simulations and Solid-State NMR

We compared the anti-influenza potencies of 57 adamantyl amines and analogs against influenza A virus with serine-31 M2 proton channel, usually termed as WT M2 channel, which is amantadine sensitive. We also tested a subset of these compounds against viruses with the amantadine-resistant L26F, V27A, A30T, G34E M2 mutant channels. Four compounds inhibited WT M2 virus in vitro with mid-nanomolar potency, with 27 compounds showing sub-micromolar to low micromolar potency. Several compounds inhibited L26F M2 virus in vitro with sub-micromolar to low micromolar potency, but only three compounds blocked L26F M2-mediated proton current as determined by electrophysiology (EP). One compound was found to be a triple blocker of WT, L26F, V27A M2 channels by EP assays, but did not inhibit V27A M2 virus in vitro, and one compound inhibited WT, L26F, V27A M2 in vitro without blocking V27A M2 channel. One compound blocked only L26F M2 channel by EP, but did not inhibit virus replication. The triple blocker compound is as long as rimantadine, but could bind and block V27A M2 channel due to its larger girth as revealed by molecular dynamics simulations, while MAS NMR informed on the interaction of the compound with M2(18-60) WT or L26F or V27A.

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.

Effects of Cholesterol on the Partitioning of a Drug Molecule in Lipid Bilayers

Drug molecules either bind to membrane-bound targets or permeate through cell membranes to reach intracellular targets, and hence, their membrane partition and permeation are of great importance. Here, we studied the effects of cholesterol on the partition of amantadine, an antiflu drug molecule, into 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers using molecular dynamics simulations. The membrane partition of amantadine is sensitive to the cholesterol mole fraction (x_chol). In the absence of cholesterol, amantadine is stably bound in membranes, but at x_chol = 32%, it can escape to the aqueous phase, in agreement with recent experiments. The reduced membrane partition of amantadine at a high cholesterol content is mainly due to the perturbation of the bilayer structure and dynamics. Surrounding lipids stabilize amantadine by having their tails wrapped around the drug molecule, and this ability is compromised when cholesterol is present to increase the order in lipid tails. The atomic details on interactions with lipids and perturbations by cholesterol revealed here provide insight into membrane partition and delivery of drug molecules to their targets.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

In Silico Study of Membrane Lipid Composition Regulating Conformation and Hydration of Influenza Virus B M2 Channel

The proton conduction of transmembrane influenza virus B M2 (BM2) proton channel is possibly mediated by the membrane environment, but the detailed molecular mechanism is challenging to determine. In this work, how membrane lipid composition regulates the conformation and hydration of BM2 channel is elucidated in silico. The appearance of several important hydrogen-bond networks has been discovered, as the addition of negatively charged lipid palmitoyloleoyl phosphatidylglycerol (POPG) and cholesterol reduces membrane fluidity and augments membrane rigidity. A more rigid membrane environment is beneficial to expand the channel, allow more water to enter the channel, promote channel hydration, and then even affect the proton conduction facilitated by the hydrated channel. Thus, membrane environment could be identified as an important influence factor of conformation and hydration of BM2. These findings can provide a unique perspective for understanding the mechanism of membrane lipid composition regulating conformation and hydration of BM2 and have important significance to the further study of anti-influenza virus B drugs.

An effective molecular blocker of ion channel of M2 protein as anti-influenza a drug

Design of a drug compound that can effectively bind to the M2 ion channel and block the diffusion of hydrogen ions (H⁺) through and inhibit influenza A virus replication is an important task. Known anti-influenza drugs amantadine and rimantadine have a weak effect on influenza A virus. A new class of positively charged, +2 e.u., molecules is proposed here to block diffusion of H⁺ ion through the M2 channel. Several drug candidates, derivatives of a lead compound (diazabicyclooctane), were proposed and investigated. Molecular dynamics of thermal fluctuations of M2 protein structure and ionization-conformation coupling of all the ionizable residues were simulated at physiological pH. The influence of the most probable mutations of key drug-binding amino acid residues in the M2 ion channel were investigated too. It is shown that the suggested new blocker has high binding affinity for the M2 channel. There are two in-channel binding sites of high affinity, the first one has H-bonds with two of four serine residues Ser-31A (B) or Ser-31C(D), and the second one has H-bonds with two of four histidine residues His-37A (B), or His-37C(D). The main advantage of the new drug molecule is the positive charge, +2 e.u., which creates a positive electrostatic potential barrier (in addition to a steric one) for a transfer of H⁺ ion through M2 channel and may serve as an effective anti-influenza A virus drug.Communicated by Ramaswamy H. Sarma.

Membrane pore formation and ion selectivity of the Ebola virus delta peptide

The Ebola virus delta peptide homo-oligomerizes in the host cell membrane to form amphipathic pores that alter the membrane properties.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

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.

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.

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.

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

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

Conformational Activation of a Transmembrane Proton Channel from Constant pH Molecular Dynamics

Proton-coupled transmembrane proteins play important roles in human health and diseases; however, detailed mechanisms are often elusive. Experimentally resolving proton positions and structural details is challenging, and conventional molecular dynamics simulations are performed with preassigned and fixed protonation states. To address this challenge, here we illustrate the use of the state-of-the-art continuous constant pH molecular dynamics (CpHMD) to directly describe the activation of the M2 channel of influenza virus, for which abundant experimental data are available. Starting from the closed crystal structure, simulation reveals a pH-dependent conformational switch to an activated state that resembles the open crystal structure. Importantly, simulation affords the free energy of channel opening coupled to the titration of a histidine tetrad, thereby providing a thermodynamic mechanism for M2 activation, that is consistent with NMR data and resolves the controversy with crystal structures obtained at different pH values. This work illustrates the utility of CpHMD in offering previously unattainable conformational details and thermodynamic information for proton-coupled transmembrane channels and transporters.

Structure and Mechanism of the Influenza A M218-60 Dimer of Dimers

We report a magic angle spinning (MAS) NMR structure of the drug-resistant S31N mutation of M218-60 from Influenza A. The protein was dispersed in diphytanoyl-sn-glycero-3-phosphocholine lipid bilayers, and the spectra and an extensive set of constraints indicate that M218-60 consists of a dimer of dimers. In particular, ∼280 structural constraints were obtained using dipole recoupling experiments that yielded well-resolved (13)C-(15)N, (13)C-(13)C, and (1)H-(15)N 2D, 3D, and 4D MAS spectra, all of which show cross-peak doubling. Interhelical distances were measured using mixed (15)N/(13)C labeling and with deuterated protein, MAS at ωr/2π = 60 kHz, ω0H/2π = 1000 MHz, and (1)H detection of methyl-methyl contacts. The experiments reveal a compact structure consisting of a tetramer composed of four transmembrane helices, in which two opposing helices are displaced and rotated in the direction of the membrane normal relative to a four-fold symmetric arrangement, yielding a two-fold symmetric structure. Side chain conformations of the important gating and pH-sensing residues W41 and H37 are found to differ markedly from four-fold symmetry. The rmsd of the structure is 0.7 Å for backbone heavy atoms and 1.1 Å for all heavy atoms. This two-fold symmetric structure is different from all of the previous structures of M2, many of which were determined in detergent and/or with shorter constructs that are not fully active. The structure has implications for the mechanism of H(+) transport since the distance between His and Trp residues on different helices is found to be short. The structure also exhibits two-fold symmetry in the vicinity of the binding site of adamantyl inhibitors, and steric constraints may explain the mechanism of the drug-resistant S31N mutation.

Efficient Multistate Reactive Molecular Dynamics Approach Based on Short-Range Effective Potentials

Nonbonded interactions between molecules usually include the van der Waals force and computationally expensive long-range electrostatic interactions. This article develops a more efficient approach: the effective-interaction multistate empirical-valence-bond (EI-MS-EVB) model. The EI-MS-EVB method relies on a mapping of all interactions onto a short-range and thus, computationally efficient effective potential. The effective potential is tabulated by matching its force to known trajectories obtained from the full-potential empirical multistate empirical-valence-bond (MS-EVB) model. The effective pairwise interaction depends on and is uniquely determined by the atomic configuration of the system, varying only with respect to the hydrogen-bonding topology. By comparing the EI-MS-EVB and full MS-EVB calculations of several equilibrium and dynamic properties important to hydrated excess proton solvation and transport, we show that the EI-MS-EVB model produces very accurate results for the specific system in which the tabulated potentials were generated. The EI-MS-EVB potential also transfers reasonably well to similar systems with different temperatures and box sizes. The EI-MS-EVB method also reduces the computational cost of the nonbonded interactions by about 1 order of magnitude in comparison with the full algorithm.

High-resolution structures of the M2 channel from influenza A virus reveal dynamic pathways for proton stabilization and transduction

The matrix 2 (M2) protein from influenza A virus is a proton channel that uses His37 as a selectivity filter. Here we report high-resolution (1.10 Å) cryogenic crystallographic structures of the transmembrane domain of M2 at low and high pH. These structures reveal that waters within the pore form hydrogen-bonded networks or "water wires" spanning 17 Å from the channel entrance to His37. Pore-lining carbonyl groups are well situated to stabilize hydronium via second-shell interactions involving bridging water molecules. In addition, room temperature crystallographic structures indicate that water becomes increasingly fluid with increasing temperature and decreasing pH, despite the higher electrostatic field. Complementary molecular dynamics simulations reveal a collective switch of hydrogen bond orientations that can contribute to the directionality of proton flux as His37 is dynamically protonated and deprotonated in the conduction cycle.

The influenza m2 cytoplasmic tail changes the proton-exchange equilibria and the backbone conformation of the transmembrane histidine residue to facilitate proton conduction

The influenza M2 protein forms an acid-activated tetrameric proton channel important for the virus lifecycle. Residue His37 in the transmembrane domain is responsible for channel activation and proton selectivity. While the structure and dynamics of His37 have been well studied in TM peptide constructs, it has not been investigated in the presence of the full cytoplasmic domain, which increases the proton conductivity by 2-fold compared to the TM peptide. We report here (13)C and (15)N chemical shifts of His37 in the cytoplasmic-containing M2(21-97) and show that cationic histidines are already present at neutral pH, in contrast to the TM peptide, indicating that the cytoplasmic domain shifts the protonation equilibria. Quantification of the imidazole (15)N intensities yielded two resolved proton dissociation constants (pKa's) of 7.1 and 5.4, which differ from the TM result but resemble the M2(18-60) result, suggesting cooperative proton binding. The average His37 pKa is higher for M2(21-97) than for the shorter constructs. We attribute this higher pKa to direct and indirect effects of the cytoplasmic domain, which is rich in acidic residues. 2D (13)C-(13)C correlation spectra reveal seven His37 Cα-Cβ cross peaks at different pH, some of which are unique to the cytoplasmic-containing M2 and correspond to more ideal α-helical conformations. Based on the pH at which these chemical shifts appear and their side chain structures, we assign these conformations to His37 in differently charged tetramers. Thus, the cytoplasmic domain facilitates proton conduction through the transmembrane pore by modifying the His37-water proton exchange equilibria and the His37 backbone conformational distribution.

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.

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

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

Multiscale simulation reveals a multifaceted mechanism of proton permeation through the influenza A M2 proton channel

The influenza A virus M2 channel (AM2) is crucial in the viral life cycle. Despite many previous experimental and computational studies, the mechanism of the activating process in which proton permeation acidifies the virion to release the viral RNA and core proteins is not well understood. Herein the AM2 proton permeation process has been systematically characterized using multiscale computer simulations, including quantum, classical, and reactive molecular dynamics methods. We report, to our knowledge, the first complete free-energy profiles for proton transport through the entire AM2 transmembrane domain at various pH values, including explicit treatment of excess proton charge delocalization and shuttling through the His37 tetrad. The free-energy profiles reveal that the excess proton must overcome a large free-energy barrier to diffuse to the His37 tetrad, where it is stabilized in a deep minimum reflecting the delocalization of the excess charge among the histidines and the cost of shuttling the proton past them. At lower pH values the His37 tetrad has a larger total charge that increases the channel width, hydration, and solvent dynamics, in agreement with recent 2D-IR spectroscopic studies. The proton transport barrier becomes smaller, despite the increased charge repulsion, due to backbone expansion and the more dynamic pore water molecules. The calculated conductances are in quantitative agreement with recent experimental measurements. In addition, the free-energy profiles and conductances for proton transport in several mutants provide insights for explaining our findings and those of previous experimental mutagenesis studies.

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.

Theoretical frameworks for multiscale modeling and simulation

Biomolecular systems have been modeled at a variety of scales, ranging from explicit treatment of electrons and nuclei to continuum description of bulk deformation or velocity. Many challenges of interfacing between scales have been overcome. Multiple models at different scales have been used to study the same system or calculate the same property (e.g., channel conductance). Accurate modeling of biochemical processes under in vivo conditions and the bridging of molecular and subcellular scales will likely soon become reality.

Multiscale Multiphysics and Multidomain Models I: Basic Theory

This work extends our earlier two-domain formulation of a differential geometry based multiscale paradigm into a multidomain theory, which endows us the ability to simultaneously accommodate multiphysical descriptions of aqueous chemical, physical and biological systems, such as fuel cells, solar cells, nanofluidics, ion channels, viruses, RNA polymerases, molecular motors and large macromolecular complexes. The essential idea is to make use of the differential geometry theory of surfaces as a natural means to geometrically separate the macroscopic domain of solvent from the microscopic domain of solute, and dynamically couple continuum and discrete descriptions. Our main strategy is to construct energy functionals to put on an equal footing of multiphysics, including polar (i.e., electrostatic) solvation, nonpolar solvation, chemical potential, quantum mechanics, fluid mechanics, molecular mechanics, coarse grained dynamics and elastic dynamics. The variational principle is applied to the energy functionals to derive desirable governing equations, such as multidomain Laplace-Beltrami (LB) equations for macromolecular morphologies, multidomain Poisson-Boltzmann (PB) equation or Poisson equation for electrostatic potential, generalized Nernst-Planck (NP) equations for the dynamics of charged solvent species, generalized Navier-Stokes (NS) equation for fluid dynamics, generalized Newton's equations for molecular dynamics (MD) or coarse-grained dynamics and equation of motion for elastic dynamics. Unlike the classical PB equation, our PB equation is an integral-differential equation due to solvent-solute interactions. To illustrate the proposed formalism, we have explicitly constructed three models, a multidomain solvation model, a multidomain charge transport model and a multidomain chemo-electro-fluid-MD-elastic model. Each solute domain is equipped with distinct surface tension, pressure, dielectric function, and charge density distribution. In addition to long-range Coulombic interactions, various non-electrostatic solvent-solute interactions are considered in the present modeling. We demonstrate the consistency between the non-equilibrium charge transport model and the equilibrium solvation model by showing the systematical reduction of the former to the latter at equilibrium. This paper also offers a brief review of the field.

Multiscale multiphysics and multidomain models--flexibility and rigidity

The emerging complexity of large macromolecules has led to challenges in their full scale theoretical description and computer simulation. Multiscale multiphysics and multidomain models have been introduced to reduce the number of degrees of freedom while maintaining modeling accuracy and achieving computational efficiency. A total energy functional is constructed to put energies for polar and nonpolar solvation, chemical potential, fluid flow, molecular mechanics, and elastic dynamics on an equal footing. The variational principle is utilized to derive coupled governing equations for the above mentioned multiphysical descriptions. Among these governing equations is the Poisson-Boltzmann equation which describes continuum electrostatics with atomic charges. The present work introduces the theory of continuum elasticity with atomic rigidity (CEWAR). The essence of CEWAR is to formulate the shear modulus as a continuous function of atomic rigidity. As a result, the dynamics complexity of a macromolecular system is separated from its static complexity so that the more time-consuming dynamics is handled with continuum elasticity theory, while the less time-consuming static analysis is pursued with atomic approaches. We propose a simple method, flexibility-rigidity index (FRI), to analyze macromolecular flexibility and rigidity in atomic detail. The construction of FRI relies on the fundamental assumption that protein functions, such as flexibility, rigidity, and energy, are entirely determined by the structure of the protein and its environment, although the structure is in turn determined by all the interactions. As such, the FRI measures the topological connectivity of protein atoms or residues and characterizes the geometric compactness of the protein structure. As a consequence, the FRI does not resort to the interaction Hamiltonian and bypasses matrix diagonalization, which underpins most other flexibility analysis methods. FRI's computational complexity is of O(N(2)) at most, where N is the number of atoms or residues, in contrast to O(N(3)) for Hamiltonian based methods. We demonstrate that the proposed FRI gives rise to accurate prediction of protein B-Factor for a set of 263 proteins. We show that a parameter free FRI is able to achieve about 95% accuracy of the parameter optimized FRI. An interpolation algorithm is developed to construct continuous atomic flexibility functions for visualization and use with CEWAR.

Conformational analysis of the full-length M2 protein of the influenza A virus using solid-state NMR

The influenza A M2 protein forms a proton channel for virus infection and mediates virus assembly and budding. While extensive structural information is known about the transmembrane helix and an adjacent amphipathic helix, the conformation of the N-terminal ectodomain and the C-terminal cytoplasmic tail remains largely unknown. Using two-dimensional (2D) magic-angle-spinning solid-state NMR, we have investigated the secondary structure and dynamics of full-length M2 (M2FL) and found them to depend on the membrane composition. In 2D (13)C DARR correlation spectra, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)-bound M2FL exhibits several peaks at β-sheet chemical shifts, which result from water-exposed extramembrane residues. In contrast, M2FL bound to cholesterol-containing membranes gives predominantly α-helical chemical shifts. Two-dimensional J-INADEQUATE spectra and variable-temperature (13)C spectra indicate that DMPC-bound M2FL is highly dynamic while the cholesterol-containing membranes significantly immobilize the protein at physiological temperature. Chemical-shift prediction for various secondary-structure models suggests that the β-strand is located at the N-terminus of the DMPC-bound protein, while the cytoplasmic domain is unstructured. This prediction is confirmed by the 2D DARR spectrum of the ectodomain-truncated M2(21-97), which no longer exhibits β-sheet chemical shifts in the DMPC-bound state. We propose that the M2 conformational change results from the influence of cholesterol, and the increased helicity of M2FL in cholesterol-rich membranes may be relevant for M2 interaction with the matrix protein M1 during virus assembly and budding. The successful determination of the β-strand location suggests that chemical-shift prediction is a promising approach for obtaining structural information of disordered proteins before resonance assignment.

Exploring Histidine Conformations in the M2 Channel Lumen of the Influenza A Virus at Neutral pH via Molecular Simulations

The pH-regulated M2 proton channel from the influenza A virus has a His-tetrad in its transmembrane (TM) domain that is essential for proton conduction. At neutral pH, the tetrad has been observed in two distinct configurations, the "His-box" and "dimer-of-dimers", with similar backbone structures suggesting competing models for proton conduction. Here, we propose that both conformations can play a role. In support of this hypothesis, we used molecular dynamics simulations based on density functional theory to simulate the M2-TM domain and force-field-based simulations to estimate the relevant free-energy barriers. Both configurations are stable on accessible simulation time scales, and transitions between them occur faster than the millisecond time scale of proton conduction. Moreover, the deprotonation energy is too high for spontaneous conduction, consistent with their occurrence in the low-current regime. Our computations support a multiconfiguration model with different population levels, thereby connecting experimental data obtained under different conditions.

Structural and energetic analysis of drug inhibition of the influenza A M2 proton channel

The type A influenza virus matrix protein 2 (M2) is a highly selective proton channel in the viral envelope. Because of its crucial role in viral infection and replication, the M2 channel has been a target of anti-influenza drugs. Due to the occurrence of drug-resistant mutations in the M2 channel, existing anti-influenza drugs that block the M2 channel, such as amantadine and rimantadine, have lost their efficacy against these mutant channels. Recent experimental and computational efforts have made great progress in understanding the drug resistance mechanisms of these mutations as well as designing novel drug candidates to block the mutant M2 channels. In this review, we briefly summarize the structural characteristics of the M2 channel, and then we discuss these recent studies on drug resistance and drug design of the mutant channels, focusing on the structures and energetics. We show that structural biology experiments and molecular modeling have led to the successful design of novel drugs targeting mutant M2 channels.

Collective diffusion model for ion conduction through microscopic channels

Ion conduction through microscopic channels is of central importance in both biology and nanotechnology. To better understand the current-voltage (I-V) dependence of ion channels, here we describe and prove a collective diffusion model that quantitatively relates the spontaneous ion permeation at equilibrium to the stationary ionic fluxes driven by small voltages. The model makes it possible to determine the channel conductance in the linear I-V range from equilibrium simulations without the application of a voltage. To validate the theory, we perform molecular-dynamics simulations on two channels-a conical-shaped nanopore and the transmembrane pore of an α-hemolysin-under both equilibrium and nonequilibrium conditions. The simulations reveal substantial couplings between the motions of cations and anions, which are effectively captured by the collective coordinate in the model. Although the two channels exhibit very different linear ranges in the I-V curves, in both cases the channel conductance at small voltages is in reasonable agreement with the prediction from the equilibrium simulation. The simulations also suggest that channel charges, rather than geometric asymmetry, play a more prominent role in current rectification.

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

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

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

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

Multiscale molecular dynamics simulations of membrane proteins

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

Computationally Efficient Multiconfigurational Reactive Molecular Dynamics

It is a computationally demanding task to explicitly simulate the electronic degrees of freedom in a system to observe the chemical transformations of interest, while at the same time sampling the time and length scales required to converge statistical properties and thus reduce artifacts due to initial conditions, finite-size effects, and limited sampling. One solution that significantly reduces the computational expense consists of molecular models in which effective interactions between particles govern the dynamics of the system. If the interaction potentials in these models are developed to reproduce calculated properties from electronic structure calculations and/or ab initio molecular dynamics simulations, then one can calculate accurate properties at a fraction of the computational cost. Multiconfigurational algorithms model the system as a linear combination of several chemical bonding topologies to simulate chemical reactions, also sometimes referred to as "multistate". These algorithms typically utilize energy and force calculations already found in popular molecular dynamics software packages, thus facilitating their implementation without significant changes to the structure of the code. However, the evaluation of energies and forces for several bonding topologies per simulation step can lead to poor computational efficiency if redundancy is not efficiently removed, particularly with respect to the calculation of long-ranged Coulombic interactions. This paper presents accurate approximations (effective long-range interaction and resulting hybrid methods) and multiple-program parallelization strategies for the efficient calculation of electrostatic interactions in reactive molecular simulations.

Expanding the view of proton pumping in cytochrome c oxidase through computer simulation

In cytochrome c oxidase (CcO), a redox-driven proton pump, protons are transported by the Grotthuss shuttling via hydrogen-bonded water molecules and protonatable residues. Proton transport through the D-pathway is a complicated process that is highly sensitive to alterations in the amino acids or the solvation structure in the channel, both of which can inhibit proton pumping and enzymatic activity. Simulations of proton transport in the hydrophobic cavity showed a clear redox state dependence. To study the mechanism of proton pumping in CcO, multi-state empirical valence bond (MS-EVB) simulations have been conducted, focusing on the proton transport through the D-pathway and the hydrophobic cavity next to the binuclear center. The hydration structures, transport pathways, effects of residues, and free energy surfaces of proton transport were revealed in these MS-EVB simulations. The mechanistic insight gained from them is herein reviewed and placed in context for future studies.

Kinetics of proton transport into influenza virions by the viral M2 channel

M2 protein of influenza A viruses is a tetrameric transmembrane proton channel, which has essential functions both early and late in the virus infectious cycle. Previous studies of proton transport by M2 have been limited to measurements outside the context of the virus particle. We have developed an in vitro fluorescence-based assay to monitor internal acidification of individual virions triggered to undergo membrane fusion. We show that rimantadine, an inhibitor of M2 proton conductance, blocks the acidification-dependent dissipation of fluorescence from a pH-sensitive virus-content probe. Fusion-pore formation usually follows internal acidification but does not require it. The rate of internal virion acidification increases with external proton concentration and saturates with a pK(m) of ∼4.7. The rate of proton transport through a single, fully protonated M2 channel is approximately 100 to 400 protons per second. The saturating proton-concentration dependence and the low rate of internal virion acidification derived from authentic virions support a transporter model for the mechanism of proton transfer.

Mechanism of function of viral channel proteins and implications for drug development

Viral channel-forming proteins comprise a class of viral proteins which, similar to their host companions, are made to alter electrochemical or substrate gradients across lipid membranes. These proteins are active during all stages of the cellular life cycle of viruses. An increasing number of proteins are identified as channel proteins, but the precise role in the viral life cycle is yet unknown for the majority of them. This review presents an overview about these proteins with an emphasis on those with available structural information. A concept is introduced which aligns the transmembrane domains of viral channel proteins with those of host channels and toxins to give insights into the mechanism of function of the viral proteins from potential sequence identities. A summary of to date investigations on drugs targeting these proteins is given and discussed in respect of their mode of action in vivo.

Membrane protein structure and dynamics from NMR spectroscopy

We review the current state of membrane protein structure determination using solid-state nuclear magnetic resonance (NMR) spectroscopy. Multidimensional magic-angle-spinning correlation NMR combined with oriented-sample experiments has made it possible to measure a full panel of structural constraints of membrane proteins directly in lipid bilayers. These constraints include torsion angles, interatomic distances, oligomeric structure, protein dynamics, ligand structure and dynamics, and protein orientation and depth of insertion in the lipid bilayer. Using solid-state NMR, researchers have studied potassium channels, proton channels, Ca(2+) pumps, G protein-coupled receptors, bacterial outer membrane proteins, and viral fusion proteins to elucidate their mechanisms of action. Many of these membrane proteins have also been investigated in detergent micelles using solution NMR. Comparison of the solid-state and solution NMR structures provides important insights into the effects of the solubilizing environment on membrane protein structure and dynamics.

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.

NMR detection of pH-dependent histidine-water proton exchange reveals the conduction mechanism of a transmembrane proton channel

The acid-activated proton channel formed by the influenza M2 protein is important for the life cycle of the virus. A single histidine, His37, in the M2 transmembrane domain (M2TM) is responsible for pH activation and proton selectivity of the channel. Recent studies suggested three models for how His37 mediates proton transport: a shuttle mechanism involving His37 protonation and deprotonation, a H-bonded imidazole-imidazolium dimer model, and a transporter model involving large protein conformational changes in synchrony with proton conduction. Using magic-angle-spinning (MAS) solid-state NMR spectroscopy, we examined the proton exchange and backbone conformational dynamics of M2TM in a virus-envelope-mimetic membrane. At physiological temperature and pH, (15)N NMR spectra show fast exchange of the imidazole (15)N between protonated and unprotonated states. To quantify the proton exchange rates, we measured the (15)N T(2) relaxation times and simulated them for chemical-shift exchange and fluctuating N-H dipolar fields under (1)H decoupling and MAS. The exchange rate is 4.5 × 10(5) s(-1) for Nδ1 and 1.0 × 10(5) s(-1) for Nε2, which are approximately synchronized with the recently reported imidazole reorientation. Binding of the antiviral drug amantadine suppressed both proton exchange and ring motion, thus interfering with the proton transfer mechanism. By measuring the relative concentrations of neutral and cationic His as a function of pH, we determined the four pK(a) values of the His37 tetrad in the viral membrane. Fitting the proton current curve using the charge-state populations from these pK(a)'s, we obtained the relative conductance of the five charge states, which showed that the +3 channel has the highest time-averaged unitary conductance. At physiologically relevant pH, 2D correlation spectra indicated that the neutral and cationic histidines do not have close contacts, ruling out the H-bonded dimer model. Moreover, a narrowly distributed nonideal helical structure coexists with a broadly distributed ideal helical conformation without interchange on the sub-10 ms time scale, thus excluding the transporter model in the viral membrane. These data support the shuttle mechanism of proton conduction, whose essential steps involve His-water proton exchange facilitated by imidazole ring reorientations.

Transmembrane communication: general principles and lessons from the structure and function of the M2 proton channel, K⁺ channels, and integrin receptors

Signal transduction across biological membranes is central to life. This process generally happens through communication between different domains and hierarchical coupling of information. Here, we review structural and thermodynamic principles behind transmembrane (TM) signal transduction and discuss common themes. Communication between signaling domains can be understood in terms of thermodynamic and kinetic principles, and complex signaling patterns can arise from simple wiring of thermodynamically coupled domains. We relate this to functions of several signal transduction systems: the M2 proton channel from influenza A virus, potassium channels, integrin receptors, and bacterial kinases. We also discuss key features in the structural rearrangements responsible for signal transduction in these systems.

Study on phylogenetic relationships, variability, and correlated mutations in M2 proteins of influenza virus A

M2 channel, an influenza virus transmembrane protein, serves as an important target for antiviral drug design. There are still discordances concerning the role of some residues involved in proton transfer as well as the mechanism of inhibition by commercial drugs. The viral M2 proteins show high conservativity; about 3/4 of the positions are occupied by one residue in over 95%. Nine M2 proteins from the H3N2 strain and possibly two proteins from H2N2 strains make a phylogenic cluster closely related to 2RLF. The variability range is limited to 4 residues/position with one exception. The 2RLF protein stands out by the presence of 2 serines at the positions 19 and 50, which are in most other M2 proteins occupied by cysteines. The study of correlated mutations shows that there are several positions with significant mutational correlation that have not been described so far as functionally important. That there are 5 more residues potentially involved in the M2 mechanism of action. The original software used in this work (Consensus Constructor, SSSSg, Corm, Talana) is freely accessible as stand-alone offline applications upon request to the authors. The other software used in this work is freely available online for noncommercial purposes at public services on bioinformatics such as ExPASy or NCBI. The study on mutational variability, evolutionary relationship, and correlated mutation presented in this paper is a potential way to explain more completely the role of significant factors in proton channel action and to clarify the inhibition mechanism by specific drugs.

A Combined Metadynamics and Umbrella Sampling Method for the Calculation of Ion Permeation Free Energy Profiles

Free energy calculations are one of the most useful methods for the study of ion transport mechanisms through confined spaces such as protein ion channels. Their reliability depends on a correctly defined reaction coordinate (RC). A straight line is usually not a proper RC for such complicated processes so in this work a combined metadynamics/umbrella sampling (MTD/US) method is proposed. In the combined method, the ion transport pathway is first identified by the MTD method and then the free energy profile or potential of mean force (PMF) along the path is calculated using umbrella sampling. This combined method avoids the discontinuity problem often associated with normal umbrella sampling calculations that assume a straight line RC and it provides a more physically accurate PMF for such processes. The method is demonstrated for the proton transport process through the protein channel of aquaporin-1.

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.