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.

Chemical Probes for Blocking of Influenza A M2 Wild-type and S31N Channels

We report on using the synthetic aminoadamantane-CH₂-aryl derivatives 1-6 as sensitive probes for blocking M2 S31N and influenza A virus (IAV) M2 wild-type (WT) channels as well as virus replication in cell culture. The binding kinetics measured using electrophysiology (EP) for M2 S31N channel are very dependent on the length between the adamantane moiety and the first ring of the aryl headgroup realized in 2 and 3 and the girth and length of the adamantane adduct realized in 4 and 5. Study of 1-6 shows that, according to molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations, all bind in the M2 S31N channel with the adamantyl group positioned between V27 and G34 and the aryl group projecting out of the channel with the phenyl (or isoxazole in 6) embedded in the V27 cluster. In this outward binding configuration, an elongation of the ligand by only one methylene in rimantadine 2 or using diamantane or triamantane instead of adamantane in 4 and 5, respectively, causes incomplete entry and facilitates exit, abolishing effective block compared to the amantadine derivatives 1 and 6. In the active M2 S31N blockers 1 and 6, the phenyl and isoxazolyl head groups achieve a deeper binding position and high k_on/low k_off and high k_on/high k_off rate constants, compared to inactive 2-5, which have much lower k_on and higher k_off. Compounds 1-5 block the M2 WT channel by binding in the longer area from V27-H37, in the inward orientation, with high k_on and low k_off rate constants. Infection of cell cultures by influenza virus containing M2 WT or M2 S31N is inhibited by 1-5 or 1-4 and 6, respectively. While 1 and 6 block infection through the M2 block mechanism in the S31N variant, 2-4 may block M2 S31N virus replication in cell culture through the lysosomotropic effect, just as chloroquine is thought to inhibit SARS-CoV-2 infection.

The boundary lipid around DMPC-spanning influenza A M2 transmembrane domain channels: Its structure and potential for drug accommodation

We have investigated the perturbation of influenza A M2TM in DMPC bilayers. We have shown that (a) DSC and SAXS detect changes in membrane organization caused by small changes (micromolar) in M2TM or aminoadamantane concentration and aminoadamantane structure, by comparison of amantadine and spiro[pyrrolidine-2,2'-adamantane] (AK13), (b) that WAXS and MD can suggest details of ligand topology. DSC and SAXS show that at a low M2TM micromolar concentration in DPMC bilayers, two lipid domains are observed, which likely correspond to M2TM boundary lipids and bulk-like lipids. At higher M2TM concentrations, one domain only is identified, which constitutes essentially all of the lipid molecules behaving as boundary lipids. According to SAXS, WAXS, and DSC in the absence of M2TM, both aminoadamantane drugs exert a similar perturbing effect on the bilayer at low concentrations. At the same concentrations of the drug when M2TM is present, amantadine and, to a lesser extent, AK13 cause, according to WAXS, a significant disordering of chain-stacking, which also leads to the formation of two lipid domains. This effect is likely due, according to MD simulations, to the preference of the more lipophilic AK13 to locate closer to the lateral surfaces of M2TM when compared to amantadine, which forms stronger ionic interactions with phosphate groups. The preference of AK13 to concentrate inside the lipid bilayer close to the exterior of the hydrophobic M2TM helices may contribute to its higher binding affinity compared to amantadine.

Novel influenza inhibitors designed to target PB1 interactions with host importin RanBP5

In search of novel targets for influenza inhibitors, a site on PB1 was selected for its high conservation and probable interaction with a host protein, RanBP5, that is key to nuclear import of PB1, where it complexes with PB2, PA, and NP to transcribe viral RNA. Docking with libraries of drug-like compounds led to a selection of five candidates that bound tightly and with a pose likely to inhibit protein binding. These were purchased and tested in vitro, found to be active, and then one was synthetically expanded to explore the structure-activity relationship. The top candidates had a carboxylic acid converted to an ester and electron-withdrawing substituents added to a phenyl group in the original structure. Resistance was slow to develop, but cytotoxicity was moderately high. Nuclear localization of PB1 and in vitro polymerase activity were both strongly inhibited.

Unraveling the Binding, Proton Blockage, and Inhibition of Influenza M2 WT and S31N by Rimantadine Variants

Recently, the binding kinetics of a ligand-target interaction, such as the residence time of a small molecule on its protein target, are seen as increasingly important for drug efficacy. Here, we investigate these concepts to explain binding and proton blockage of rimantadine variants bearing progressively larger alkyl groups to influenza A virus M2 wild type (WT) and M2 S31N protein proton channel. We showed that resistance of M2 S31N to rimantadine analogues compared to M2 WT resulted from their higher k_off rates compared to the k_on rates according to electrophysiology (EP) measurements. This is due to the fact that, in M2 S31N, the loss of the V27 pocket for the adamantyl cage resulted in low residence time inside the M2 pore. Both rimantadine enantiomers have similar channel blockage and binding k_on and k_off against M2 WT. To compare the potency between the rimantadine variants against M2, we applied approaches using different mimicry of M2, i.e., isothermal titration calorimetry and molecular dynamics simulation, EP, and antiviral assays. It was also shown that a small change in an amino acid at site 28 of M2 WT, which does not line the pore, seriously affects M2 WT blockage kinetics.

Divalent copper complexes as influenza A M2 inhibitors

New M2 blockers effective against the ubiquitous amantadine-resistant S31N M2 mutation in influenza A are needed. Six copper complexes, 2, 4, 6, 8, 9, and 10, were synthesized and found to block both wild type and S31N M2. Free Cu²⁺ also blocks M2 S31N but not S31N/H37A. The copper complexes do not block M2 H37A (either S31 or S31N). The complexes were effective against three influenza A strains in cell-culture assays, but less toxic to cells than CuCl₂. For example 4, Cu(cyclooctylamineiminodiacetate), which was stable at pH > 4 in the buffers used, had an EC₅₀ against A/Calif/07/2009 H1N1 of 0.7 ± 0.1 μM with a CC₅₀ of 147 μM (therapeutic index, averaged over three strains, 67.8). In contrast, CuCl₂ had an EC₅₀ of 3.8 ± 0.9 μM and CC₅₀ of 19 μM. Because M2 H37 is highly conserved, these complexes show promise for further testing as drugs against all strains of influenza A.

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

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

Investigation of a recent rise of dual amantadine-resistance mutations in the influenza A M2 sequence

Background

The S31N amantadine-resistance mutation in the influenza A M2 sequence currently occurs more frequently in nature than the S31 wild type. Overcoming the resistance of the S31N mutation is the primary focus of M2 researchers who aim to develop novel antiviral therapies. Recent studies have noted a possible rise in frequency of the V27A/S31N double amantadine-resistance mutation in recent years. The purpose of this study is to investigate this recent rise in frequency of the double mutation and any possible bias of the other mutations toward co-occurrence with S31N or S31 strains.

Results

The primary dataset used for this study was comprised of 24,152 influenza A M2 channel sequences which were downloaded from UniProt. There is an increased frequency for the S31N/V27A dual AR mutation in recent years, especially in swine. A test for difference in two proportions indicates that the V27A mutation is co-occurring with S31N more often than expected (p-value < 0.001) when considering individual amino acid frequencies. At the same time, the different propensities for the V27A as compared to the V27T dual mutant may reflect differences in viral fitness or protein energetics, and this information could be exploited to focus drug development so as to reduce further drug insensitivity.

Conclusions

The development of the S31N/V27A variant in the Midwestern US swine may be a harbinger of novel human strain development. V27A/S31N is a possible path forward for the evolution of M2 which may convey a new level of drug resistance and should receive attention in drug design.

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

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

Aminoadamantanes with persistent in vitro efficacy against H1N1 (2009) influenza A

A series of 2-adamantanamines with alkyl adducts of various lengths were examined for efficacy against strains of influenza A including those having an S31N mutation in M2 proton channel that confer resistance to amantadine and rimantadine. The addition of as little as one CH₂ group to the methyl adduct of the amantadine/rimantadine analogue, 2-methyl-2-aminoadamantane, led to activity in vitro against two M2 S31N viruses A/Calif/07/2009 (H1N1) and A/PR/8/34 (H1N1) but not to a third A/WS/33 (H1N1). Solid state NMR of the transmembrane domain (TMD) with a site mutation corresponding to S31N shows evidence of drug binding. But electrophysiology using the full length S31N M2 protein in HEK cells showed no blockade. A wild type strain, A/Hong Kong/1/68 (H3N2) developed resistance to representative drugs within one passage with mutations in M2 TMD, but A/Calif/07/2009 S31N was slow (>8 passages) to develop resistance in vitro, and the resistant virus had no mutations in M2 TMD. The results indicate that 2-alkyl-2-aminoadamantane derivatives with sufficient adducts can persistently block p2009 influenza A in vitro through an alternative mechanism. The observations of an HA1 mutation, N160D, near the sialic acid binding site in both 6-resistant A/Calif/07/2009(H1N1) and the broadly resistant A/WS/33(H1N1) and of an HA1 mutation, I325S, in the 6-resistant virus at a cell-culture stable site suggest that the drugs tested here may block infection by direct binding near these critical sites for virus entry to the host cell.

Fluorescence anisotropy of diphenylhexatriene and its cationic Trimethylamino derivative in liquid dipalmitoylphosphatidylcholine liposomes: opposing responses to isoflurane

Background

The mechanism of action of volatile general anesthetics has not yet been resolved. In order to identify the effects of isoflurane on the membrane, we measured the steady-state anisotropy of two fluorescent probes that reside at different depths. Incorporation of anesthetic was confirmed by shifting of the main phase transition temperature.

Results

In liquid crystalline dipalmitoylphosphatidylcholine liposomes, isoflurane (7-25 mM in the bath) increases trimethylammonium-diphenylhexatriene fluorescence anisotropy by ~0.02 units and decreases diphenylhexatriene anisotropy by the same amount.

Conclusions

The anisotropy data suggest that isoflurane decreases non-axial dye mobility in the headgroup region, while increasing it in the tail region. We propose that these results reflect changes in the lateral pressure profile of the membrane.

Functional reconstitution of influenza A M2(22-62)

Amantadine-sensitive proton uptake by liposomes is currently the preferred method of demonstrating M2 functionality after reconstitution, to validate structural determination with techniques such as solid-state NMR. With strong driving forces (two decades each of both [K(+)] gradient-induced membrane potential and [H(+)] gradient), M2(22-62) showed a transport rate of 78 H(+)/tetramer-s (pH(o) 6.0, pH(i) 8.0, nominal V(m)=-114 mV), higher than previously measured for similar, shorter, and full-length constructs. Amantadine sensitivity of the conductance domain at pH 6.8 was also comparable to other published reports. Proton flux rate was optimal at protein densities of 0.05-1.0% (peptide wt.% in lipid). Rundown of total proton uptake after addition of valinomycin and CCCP, as detected by delayed addition of valinomycin, indicated M2-induced K(+) flux of 0.1K(+)/tetramer-s, and also demonstrated that the K(+) permeability, relative to H(+), was 2.8 × 10(-6). Transport rate, amantadine and cyclooctylamine sensitivity, acid activation, and H(+) selectivity were all consistent with full functionality of the reconstituted conductance domain. Decreased external pH increased proton uptake with an apparent pK(a) of 6.

Insight into the mechanism of the influenza A proton channel from a structure in a lipid bilayer

The M2 protein from the influenza A virus, an acid-activated proton-selective channel, has been the subject of numerous conductance, structural, and computational studies. However, little is known at the atomic level about the heart of the functional mechanism for this tetrameric protein, a His(37)-Trp(41) cluster. We report the structure of the M2 conductance domain (residues 22 to 62) in a lipid bilayer, which displays the defining features of the native protein that have not been attainable from structures solubilized by detergents. We propose that the tetrameric His(37)-Trp(41) cluster guides protons through the channel by forming and breaking hydrogen bonds between adjacent pairs of histidines and through specific interactions of the histidines with the tryptophan gate. This mechanism explains the main observations on M2 proton conductance.

Drug sensitivity, drug-resistant mutations, and structures of three conductance domains of viral porins

Recent controversies associated with the structure of the M2 protein from influenza A virus and the binding site of drug molecules amantadine and rimantadine motivated the comparison here of the drug binding to three viral porins including the M2 proteins from influenza A and B as well as the viral protein 'u' from HIV-1. While the M2 protein from influenza B does not normally bind amantadine, chimeras with the M2 protein from influenza A show blockage by amantadine. Similarly, Vpu does not normally bind rimantadine, but the single site mutation A18H converts a non-specific channel to a selective proton channel that is sensitive to rimantadine. The comparison of structures and amino acid sequences shows that the membrane protein sample environment can have a significant influence on the structural result. While a bilayer surface bound amphipathic helix has been characterized for AM2, such a helix may be possible for BM2 although it has evaded structural characterization in detergent micelles. A similar amphipathic helix seems less likely for Vpu. Even though the A18H Vpu mutant forms rimantadine sensitive proton channels, the binding of drug and its influence on the protein structure appears to be very different from that for the M2 proteins. Indeed, drug binding and drug resistance in these viral porins appears to result from a complex set of factors.

Gramicidin channels are internally gated

Gramicidin channels are archetypal molecular subjects for solid-state NMR studies and investigations of single-channel or cation conductance. Until now, the transitions between on and off conductance states have been thought, based on multichannel studies, to represent monomer <--> dimer reactions. Here we use a single-molecule deposition method (vesicle fusion to a planar bilayer) to show that gramicidin dimer channels do not normally dissociate when conductance terminates. Furthermore, the observation of two 13C peaks in solid-state NMR indicates very stable dichotomous conformations for both the first and second peptide bonds in the monomers, and a two-dimensional chemical exchange spectrum with a 12-s mixing time demonstrates that the Val1 carbonyl conformations exchange slowly, with lifetimes of several seconds. It is proposed that gramicidin channels are gated by small conformational changes in the channel near the permeation pathway. These studies demonstrate how regulation of conformations governing closed <--> open transitions may be achieved and studied at the molecular level.

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

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

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

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

Noncontact dipole effects on channel permeation. III. Anomalous proton conductance effects in gramicidin

Proton transport on water wires, of interest for many problems in membrane biology, is analyzed in side-chain analogs of gramicidin A channels. In symmetrical 0.1N HCl solutions, fluorination of channel Trp(11), Trp-(13), or Trp(15) side chains is found to inhibit proton transport, and replacement of one or more Trps with Phe enhances proton transport, the opposite of the effects on K(+) transport in lecithin bilayers. The current-voltage relations are superlinear, indicating that some membrane field-dependent process is rate limiting. The interfacial dipole effects are usually assumed to affect the rate of cation translocation across the channel. For proton conductance, however, water reorientation after proton translocation is anticipated to be rate limiting. We propose that the findings reported here are most readily interpreted as the result of dipole-dipole interactions between channel waters and polar side chains or lipid headgroups. In particular, if reorientation of the water column begins with the water nearest the channel exit, this hypothesis explains the negative impact of fluorination and the positive impact of headgroup dipole on proton conductance.

Proton transport through influenza A virus M2 protein reconstituted in vesicles

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

Tryptophan contributions to the empirical free-energy profile in gramicidin A/M heterodimer channels

Gramicidin A/gramicidin M heterodimer conductances were measured in planar lipid bilayers and found to form two distinguishable populations about halfway between the gramicidin A and gramicidin M homodimer conductances. This implies that the principle difference in the gramicidin A and gramicidin M transport free-energy profiles occurs at the channel center, where it would produce similar effects on the rate-limiting barrier for the two heterodimers. Kinetic analysis based on this and nearly all previously published homodimer conductance data for both gramicidin A and gramicidin M channels confirms this conclusion, indicating that the translocation step is approximately 100-fold slower in gramicidin M homodimers than in gramicidin A homodimers and that first- and second-ion exit-rate constants are higher by factors of 24 and 10, respectively. Assuming that the ratios of rate constants are related to the free-energy difference between gramicidin A and gramicidin M, we construct an effective ion-Trp free-energy interaction profile that has a minimum at the channel center.

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

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

Proton conductance of influenza virus M2 protein in planar lipid bilayers

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

Applied field nonequilibrium molecular dynamics simulations of ion exit from a β-barrel model of the L-type calcium channel

We present results of applied field nonequilibrium molecular dynamics simulations (AF NEMD) of a minimal beta-barrel model channel intended to represent an L-type calcium channel that suggests a possible relationship between glutamate side chain conformational changes and ion flux in calcium channels. The beta-barrel is used to provide a scaffolding for glutamate side chains and a confinement for electrolyte of dimensions similar to the expected channel structure. It was preloaded with ions to explore relative rates of ion exit for different occupancy configurations. Our simulations with an asymmetrical flexible selectivity filter represented by four glutamate side chains (EEEE), one of which differs in initial dihedrals from the other three, indicate a plausible mechanism for the observed anomalous mole fraction effect seen in calcium channels. Apparent rates of electric field-induced exit from channels preloaded with three Na+ ions are much higher than for channels with one Ca2+ followed by two Na+ ions, consistent with the common notion that Ca2+ block of Na+ current is due to competition between the Ca2+ and Na+ ions for the negatively charged (EEEE) locus. In our model, the Ca2+ ion ligates simultaneously to the four negatively charged glutamate side chains and sterically blocks the permeation pathway. Ca2+-relief of Ca2+-block is suggested by a much higher rate of exit for channels preloaded with three Ca2+ ions than for channels with two Ca2+ ions.

AFM visualization of mobile influenza A M2 molecules in planar bilayers

We report the observation of influenza A M2 (M2) incorporated in a dipalmitoylphosphatidylcholine (DPPC) supported planar bilayer on mica, formed by use of a modified vesicle fusion method from proteoliposomes and visualized with contact mode atomic force microscopy. Incubation of proteoliposomes in a hyperosmotic solution and increased DPPC/M2 weight ratios improved supported planar bilayer formation by M2/DPPC proteoliposomes. M2's extra-bilayer domains were observed as particles estimated to protrude 1-1.5 nm above the bilayer surface and <4 nm in diameter. Particle density was 5-18% of the nominal tetramer density. Movement of observable M2 particles was independent of the probe tip. The mean lateral diffusion coefficient (D) of M2 was 4.4 +/- 1.0 x 10(-14) cm(2)/s. Eighty-two percent of observable particles were mobile on the observable timescale (D > 6 x 10(-15) cm(2)/s). Protein-protein interactions were also observed directly.

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).

Noncontact dipole effects on channel permeation. VI. 5F- and 6F-Trp gramicidin channel currents

Fluorination of peptide side chains has been shown to perturb gramicidin channel conductance without significantly changing the average side chain structure, which, it is hoped, will allow detailed analysis of electrostatic modulation of current flow. Here we report a 1312-point potassium current-voltage-concentration data set for homodimeric channels formed from gramicidin A (gA) or any of eight fluorinated Trp analogs in both lecithin and monoglyceride bilayers. We fit the data with a three-barrier, two-site, two-ion (3B2S) kinetic model. The fluorination-induced changes in the rate constants were constrained by the same factor in both lipids. The rate constant changes were converted to transition-state free-energy differences for comparison with previous electrostatic potential energy differences based on an ab initio force field. The model allowed a reasonably good fit (chi = 8.29 with 1271 degrees of freedom). The measured changes were subtle. Nevertheless, the fitted energy perturbations agree well with electrostatic predictions for five of the eight peptides. For the other three analogs, the fitted changes suggested a reduced translocation barrier rather than the reduced exit barrier as predicted by electrostatics.

Spectroscopic analysis of low pH and lipid-induced structural changes in type A botulinum neurotoxin relevant to membrane channel formation and translocation

Botulinum neurotoxin (BoNT) is an extremely toxic protein to animals and humans. In its mode of action, one of its subunits mediates its translocation by integrating itself into the membrane bilayer. We have examined the membrane channel activity of type A BoNT (BoNT/A) and its heavy (H) chain in planar lipid membrane under various pH conditions to understand the possible role of the channel activity in the translocation of the BoNT/A light (L) chain under physiological conditions. Only BoNT/A H chain, and not the BoNT/A, exhibited membrane channel activity for translocation of ions. The H chain-induced increase in conductance did not require a pH gradient across the lipid membrane, although it was enhanced by a pH gradient. To understand the molecular basis of the membrane channel activity and the translocation of the L chain, the secondary structure of BoNT/A and its H and L chains were analyzed using circular dichroism (CD) and Fourier-transform infrared (FT-IR) spectroscopy at different pH values. BoNT/A showed no structural alternation upon acidifying the buffer pH. However, an increase in beta-sheet content of BoNT/A H chain at low pH was noted when examined by FT-IR. The L chain structure significantly changed with decrease in pH, and the change was mostly reversible. In addition, the neurotoxin and its subunit chains induced a partially reversible aggregation of liposomes at low pH, which indicated their integration into the lipid bilayer. Temperature-induced denaturation studies of BoNT/A H chain indicated major structural reorganization upon its interaction with membrane, especially at low pH.

The role of Trp side chains in tuning single proton conduction through gramicidin channels

We present an extensive set of measurements of proton conduction through gramicidin A (gA), B (gB), and M (gM) homodimer channels which have 4, 3, or 0 Trp residues at each end of the channel, respectively. In gA we find a shoulder separating two domains of conductance increasing with concentration, confirming the results of Eisenman, G., B. Enos, J. Hagglund, and J. Sandblom. 1980. Ann. NY. Acad. Sci. 339:8-20. In gB, the shoulder is shifted by approximately 1/2 pH unit to higher H(+) concentrations and is very sharply defined. No shoulder appears in the gM data, but an associated transition from sublinear to superlinear I-V values occurs at a 100-fold higher [H(+)] in gM than in gA. The data in the low concentration domain are analyzed using a configuration space model of single-proton conduction, assuming that the difference in the proton potential of mean force (PMF) between gA and its analogs is constant, similar to the results of Anderson, D., R. B. Shirts, T. A. Cross, and D. D. Busath. 2001. Biophys. J. 81:1255-1264. Our results suggest that the average amplitudes of the calculated proton PMFs are nearly correct, but that the water reorientation barrier calculated for gA by molecular dynamics using the PM6 water model (Pomès, R., and B. Roux. 1997. Biophys. J. 72:246a) must be reduced in amplitude by 1.5 kcal/mol or more, and is not rate-limiting for gA.

An inverting basket model for AE1 transport

We describe an "inverting basket" model for transport in the erythrocyte anion exchanger, AE1. The inverting basket is formed by the side chains of three putative key residues, two positively (Lys 826 and Arg 730) and one negatively (Glu 681) charged residue. We have tentatively chosen seven transmembrane helices, TM1, TM2, TM4, TM8, TM10, TM12 and TM13 to form a conical channel using the well-established Glu 681 of TM8 and candidates Lys 826 and Arg 730 of TM12-13 and TM10, respectively, to form the inverting basket. We assume that these residues bind to an anion and shift from outward facing (C(o)) to inward facing (C(i)) conformation without significant backbone movements to transport an anion across the membrane. The transition of the complex (residues and ion) from outward facing (C(o)) to inward facing (C(i)) constitutes one "basket" inversion. The barrier to inversion is composed of two major components: that of the anhydrous complex, which we refer to as a steric energy barrier and a dehydration effect due to the removal of charges in the complex from water in the channel. The steric barrier is dependent on the side chain charge and configuration and on the ion charge and size. The dehydration effect, for our model, ameliorates the steric barrier, in the case of the empty complex but less so for the monovalent and divalent ions. We conclude, that it is possible for a seven-helix bundle to have a steric barrier to basket inversion, but that hydration effects in thin hydrophobic barrier models may be more complex than usually envisioned.

Model channel ion currents in NaCl-extended simple point charge water solution with applied-field molecular dynamics

Using periodic boundary conditions and a constant applied field, we have simulated current flow through an 8.125-A internal diameter, rigid, atomistic channel with polar walls in a rigid membrane using explicit ions and extended simple point charge water. Channel and bath currents were computed from 10 10-ns trajectories for each of 10 different conditions of concentration and applied voltage. An electric field was applied uniformly throughout the system to all mobile atoms. On average, the resultant net electric field falls primarily across the membrane channel, as expected for two conductive baths separated by a membrane capacitance. The channel is rarely occupied by more than one ion. Current-voltage relations are concentration dependent and superlinear at high concentrations.

Noncontact dipole effects on channel permeation. IV. Kinetic model of 5F-Trp(13) gramicidin A currents

Nonlinear least squares fitting was used to assign rate constants for the three-barrier, two-site, double-occupancy, single-filing kinetic model for previously reported current-voltage relations of (5F-Indole)Trp(13) gramicidin A and gramicidin A channels (, 75:2830-2844). By judicious coupling of parameters, it was possible to reduce the parameter space from 64 parameters to 24, and a reasonable fit consistent with other experimental data was obtained. The main features of the fit were that fluorination increased the rate constant for translocation by a factor of 2.33, consistent with a free energy change in the translocation barrier of -0.50 kcal/mol, and increased first-ion binding affinity by a factor of 1.13, primarily by decreasing the first-ion exit rate constant. The translocation rate constant was 5.62 times slower in diphytanoyl phosphatidylcholine (DPhPC) bilayers than in monoolein (GMO) bilayers (coupled for the four combinations of peptide and salt), suggesting a 44.2-mV difference in the projection of the interfacial dipole into the channel. Thus fluorination caused increased currents in DPhPC bilayers, where a high interfacial dipole potential makes translocation more rate limiting because the translocation barrier was reduced, and decreased currents in GMO bilayers, where ion exit or entry is rate limiting because these barriers were increased.