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.     

Dynamics of proton exchange in a model phosphonic acid‐functionalized polymer

The dynamics and mechanism of proton exchange in phosphonic acid‐functionalized polymers were studied using poly(vinyl‐phosphonic acid) (PVPA) as a model system along with quantum chemical calculations and Born–Oppenheimer molecular dynamics (BOMD) simulations at the B3LYP/TZVP level as model calculations. This theoretical study began with searching for the smallest, most active polymer segments and their intermediate conformations which could be involved in the local proton‐exchange process. The B3LYP/TZVP results confirmed that a low local dielectric environment and excess proton conditions are required to generate the intermediate conformations, and the shapes of the potential energy curves of the proton exchange between the two phosphonic acid functional groups are sensitive to the local conformational changes. In contrast, a high local dielectric environment increases the energy barriers, thereby preventing the proton from returning to the original functional group. Based on the static results, a mechanism for the proton exchange between the two functional groups involving fluctuations in the local dielectric environment and a local conformational change was proposed. The BOMD results confirmed the proposed mechanism by showing that the activation energies for the proton exchange in the hydrogen bond between two immobilized phosphonic acid moieties, obtained from the exponential relaxation behaviors of the envelopes of the velocity autocorrelation functions and the ¹H Nuclear Magnetic Resonance (NMR) line‐shape analyses, are too low to be the rate‐determining process. Instead, coupled librational motion in the backbone which leads to the interconversion between the two intermediate conformations possesses higher activation energy, and therefore represents one of the most important rate‐determining processes. These findings suggested that the rate of the proton exchange in the model phosphonic acid‐functionalized polymer is determined by the polymer mobility which, in this case, is the large‐amplitude librational motion of the vinyl backbone. © 2015 Wiley Periodicals, Inc.     

Electrostatic study of the proton pumping mechanism in bovine heart cytochrome C oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme of the cell respiratory chain in mitochondria and aerobic bacteria. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping across the inner-mitochondrial membrane, a process that results in a membrane electrochemical proton gradient. Although the structure of the enzyme has been solved for several organisms, the molecular mechanism of proton pumping remains unknown. In the present paper, continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine heart cytochrome c oxidase for different redox states of the enzyme along its catalytic cycle. Three different computational models of the enzyme were employed to test the stability of the results. The energetics and pH dependence of the P-->F, F-->O, and O-->E steps of the cycle have been investigated. On the basis of electrostatic calculations, two possible schemes of redox-linked proton pumping are discussed. The first scheme involves His291 as a pump element, whereas the second scheme involves a group linked to propionate D of heme a(3). In both schemes, loading of the pump site is coupled to ET between the two hemes of the enzyme, while transfer of a chemical proton is accompanied by ejection of the pumped H(+). The two models, as well as the energetics results are compared with recent experimental kinetic data. The proton pumping across the membrane is an endergonic process, which requires a sufficient amount of energy to be provided by the chemical reaction in the active site. In our calculations, the conversion of OH(-) to H(2)O provides 520 meV of energy to displace pump protons from a loading site and overall about 635 meV for each electron passing through the system. Assuming that the two charges are translocated per electron against the membrane potential of 200 meV, the model predicts an overall efficiency of 63%.     

Dielectric losses and proton conductivity of polyantimonic acid membranes

Composite membranes consisting of polyvinyl alcohol and disperse particles of polyantimonic acid with different ratios of components were studied. Processes related to the relaxation and proton transport in a heterogeneous system were considered. The proton transport in the membranes was shown to occur both within the particle (along the chains of alternating wcapacitorater molecules and oxonium ions), forming macrodipoles, and in the intercrystallite space in the polymer matrix. The relaxation was described by the two-layer Maxwell–Wagner capacitor with an effective relaxation time. The relaxation activation energy obtained from the dielectric loss data coincided with the proton conductivity activation energy in the frequency-independent region, which points to the same nature of ion transport.     

Intermolecular proton transport along the hydrogen bond and the dielectric properties of ice crystals

Conclusion We have attempted here to explain some of the dielectric properties of ice crystals by a concept involving the activated transport of a proton along the hydrogen bond. The choice of the almost symmetrical curve with two minima for a proton in the hydrogen bridge and the rejection of tunnel transport of a proton makes it possible to suggest a theory which explains adequately the dependence of the low frequency dielectric constant on temperature and the phenomenon of the dielectric relaxation of ice. The electrical polarization of an ice crystal according to our model is produced by the defect structure of the ice, the ion pairs, which arises when the proton moves between two neutral water molecules.This particular model does not involve orientation or Bjerrum defects in the ice structure and raises a doubt as to the necessity of introducing them into explanations of the dielectric properties of ice.Moscow Physicotechnical Institute. Translated from Zhurnal Strukturnoi Khimii, Vol. 11, No. 3, pp. 415–420, May–June, 1970.     

Energetics and dynamics of proton transfer reactions along short water wires

Proton transfer (pT) reactions in biochemical processes are often mediated by chains of hydrogen-bonded water molecules. We use hybrid density functional calculations to study pT along quasi one-dimensional water arrays that connect an imidazolium-imidazole proton donor-acceptor pair. We characterize the structures of intermediates and transition states, the energetics, and the dynamics of the pT reactions, including vibrational contributions to kinetic isotope effects. In molecular dynamics simulations of pT transition paths, we find that for short water chains with four water molecules, the pT reactions are semi-concerted. The formation of a high-energy hydronium intermediate next to the proton-donating group is avoided by a simultaneous transfer of a proton from the donor to the first water molecule, and from the first water molecule into the water chain. Lowering the dielectric constant of the environment and increasing the water chain length both reduce the barrier for pT. We study the effect of the driving force on the energetics of the pT reaction by changing the proton affinity of the donor and acceptor groups through halogen and methyl substitutions. We find that the barrier of the pT reaction depends linearly on the proton affinity of the donor but is nearly independent of the proton affinity of the acceptor, corresponding to Br?nsted slopes of one and zero, respectively.     

Functional dynamics of ion channels: modulation of proton movement by conformational switches

Detailed comparative studies of proton relay in native and chemically modified gramicidin channels provide a unique opportunity to uncover the structural basis of biological proton transport. The function of ion channels hinges on their ability to provide surrogate solvation in narrow pore filters so as to overcome the dielectric barrier presented by biological membranes. In the potassium channel KcsA and in the cation channel gramicidin, permeant selectivity and mobility are determined by the proteinaceous matrix via hydrogen bonding, charge-dipole, and dipole-dipole interactions. In particular, main-chain carbonyl groups in these pore interiors play an essential role in the solvation of alkali ions and of protons. In this study, molecular dynamics simulations reveal how the translocation of H(+) is controlled by nanosecond conformational transitions exchanging distorted states of the peptidic backbone in the single-file region of a dioxolane-linked analogue of the gramicidin dimer. These results underline the functional role of channel dynamics and provide a mechanism for the modulation of proton currents by fluctuating dipoles.     

"Proton holes" in long-range proton transfer reactions in solution and enzymes: A theoretical analysis

Proton transfers are fundamental to chemical processes in solution and biological systems. Often, the well-known Grotthuss mechanism is assumed where a series of sequential "proton hops" initiates from the donor and combines to produce the net transfer of a positive charge over a long distance. Although direct experimental evidence for the sequential proton hopping has been obtained recently, alternative mechanisms may be possible in complex molecular systems. To understand these events, all accessible protonation states of the mediating groups should be considered. This is exemplified by transfers through water where the individual water molecules can exist in three protonation states (water, hydronium, and hydroxide); as a result, an alternative to the Grotthuss mechanism for a proton transfer through water is to generate a hydroxide by first protonating the acceptor and then transfer the hydroxide toward the donor through water. The latter mechanism can be most generally described as the transfer of a "proton hole" from the acceptor to the donor where the "hole" characterizes the deprotonated state of any mediating molecule. This pathway is distinct and is rarely considered in the discussion of proton-transfer processes. Using a calibrated quantum mechanical/molecular mechanical (QM/MM) model and an effective sampling technique, we study proton transfers in two solution systems and in Carbonic Anhydrase II. Although the relative weight of the "proton hole" and Grotthuss mechanisms in a specific system is difficult to determine precisely using any computational approach, the current study establishes an energetics motivated framework that hinges on the donor/acceptor pKa values and electrostatics due to the environment to argue that the "proton hole" transfer is likely as important as the classical Grotthuss mechanism for proton transport in many complex molecular systems.     

Isotope effects associated with tunneling and double proton transfer in the hydrogen bonds of benzoic acid

The isotope effects associated with double proton transfer in the hydrogen bonds of benzoic acid (BA) dimers have been measured using field-cycling (1)H NMR relaxometry and quasielastic neutron scattering. By studying mixed isotope (hydrogen and deuterium) samples, the dynamics of three isotopologues, BA-HH, BA-HD, and BA-DD, have been investigated. Low temperature measurements provide accurate measurements of the incoherent tunneling rate, k(0). This parameter scales accurately with the mass number, m, according to the formula k(0)=(E/m)e(-Fm) providing conclusive evidence that the proton transfer process is a strongly correlated motion of two hydrons. Furthermore, we conclude that the tunneling pathway is the same for the three isotopologue species. Measurements at higher temperatures illuminate the through barrier processes that are mediated via intermediate or excited vibrational states. In parallel with the investigation of proton transfer dynamics, the theoretical and experimental aspects of studying spin-lattice relaxation in single crystals of mixed isotope samples are investigated in depth. Heteronuclear dipolar interactions between (1)H and (2)H isotopes contribute significantly to the overall proton spin-lattice relaxation and it is shown that these must be modeled correctly to obtain accurate values for the proton transfer rates. Since the sample used in the NMR measurements was a single crystal, full account of the orientation dependence of the spin-lattice relaxation with respect to the applied B field was incorporated into the data analysis.     

Modeling energy landscapes of proton motion in nonaqueous, tethered proton wires

We have modeled structures and energetics of anhydrous proton-conducting wires: tethered hydrogen-bonded chains of the form ···HX···HX···HX···, with functional groups HX = imidazole, triazole, and formamidine; formic, sulfonic, and phosphonic acids. We have applied density functional theory (DFT) to model proton wires up to 19 units long, where each proton carrier is linked to an effective backbone to mimic polymer tethering. This approach allows the direct calculation of hydrogen bond strengths. The proton wires were found to be stabilized by strong hydrogen bonds (up to 50 kJ/mol) whose strength correlates with the proton affinity of HX [related to pK(b)(HX)] and not to pK(a)(HX) as is often assumed. Geometry optimizations and ab initio molecular dynamics near 400 K on imidazole-based proton wires both predict that adding a proton to the end of such wires causes the excess charge to embed into the interior segments of these wires. Proton translocation energy landscapes for imidazole-based wires are sensitive to the imidazole attachment point (head or feet) and to wire architecture (linear or interdigitated). Linear imidazole wires with head-attachment exhibit low barriers for intrawire proton motion, rivaling proton diffusion in liquid imidazole. Excess charge relaxation from the edge of wires is found to be dominated by long-range Grotthuss shuttling for distances as long as 42 ?, especially for interdigitated wires. For imidazole, we predict that proton translocation is controlled by the energetics of desorption from the proton wire, even for relatively long wires (600 imidazole units). Proton desorption energies show no correlation with functional group properties, suggesting that proton desorption is a collective process in proton wires.     

Long-lived states to monitor protein unfolding by proton NMR

The relaxation of long-lived states (LLS) corresponds to the slow return to statistical thermal equilibrium between symmetric and antisymmetric proton spin states. This process is remarkably sensitive to the presence of external spins and can be used to obtain information about partial unfolding of proteins. We detected the appearance of a destabilized conformer of ubiquitin when urea is added to the protein in its native state. This conformer shows increased mobility in the C-terminus, which significantly extends the lifetimes of proton LLS magnetisation in Ser-65. These changes could not be detected by conventional measurements of T(1) and T(2) relaxation times of protons, and would hardly be sensed by carbon-13 or nitrogen-15 relaxation measurements. Conformers with similar dynamic and structural features, as revealed by LLS relaxation times, could be observed, in the absence of urea, in two ubiquitin mutants, L67S and L69S.     

Effect of noncovalent interactions on conformers of the n-butylbenzene monomer studied by mass analyzed threshold ionization spectroscopy and basis-set convergent ab initio computations

Two conformational isomers of the aromatic hydrocarbon n-butylbenzene have been studied using two-color MATI (mass analyzed threshold ionization) spectroscopy to explore the effect of conformation on ionization dynamics. Cationic states of g auche-conformer III and anti- conformers IV were selectively produced by two-color excitation via the respective S 1 origins. Adiabatic ionization potentials of the gauche- and anti-conformations were determined to be 70146 and 69872 +/- 5 cm (-1) respectively. Spectral features and vibrational modes are interpreted with the aid of MP2/cc-pVDZ ab initio calculations, and ionization-induced changes in the molecular conformations are discussed. Complete basis set (CBS) ab initio studies at MP2 level reveal reliable energetics for all four n-butylbenzene conformers observed in earlier two-color REMPI (resonance enhanced multiphoton ionization) spectra. For the S 0 state, the energies of conformer III, IV and V are above conformer I by 130, 289, 73 cm (-1), respectively. Furthermore, the combination of the CBS calculations with the measured REMPI, MATI spectra allowed the determination of the energetics of all four conformers in the S 1 and D 0 states.     

Simulating the proton transfer in gramicidin A by a sequential dynamical Monte Carlo method

The large interest in long-range proton transfer in biomolecules is triggered by its importance for many biochemical processes such as biological energy transduction and drug detoxification. Since long-range proton transfer occurs on a microsecond time scale, simulating this process on a molecular level is still a challenging task and not possible with standard simulation methods. In general, the dynamics of a reactive system can be described by a master equation. A natural way to describe long-range charge transfer in biomolecules is to decompose the process into elementary steps which are transitions between microstates. Each microstate has a defined protonation pattern. Although such a master equation can in principle be solved analytically, it is often too demanding to solve this equation because of the large number of microstates. In this paper, we describe a new method which solves the master equation by a sequential dynamical Monte Carlo algorithm. Starting from one microstate, the evolution of the system is simulated as a stochastic process. The energetic parameters required for these simulations are determined by continuum electrostatic calculations. We apply this method to simulate the proton transfer through gramicidin A, a transmembrane proton channel, in dependence on the applied membrane potential and the pH value of the solution. As elementary steps in our reaction, we consider proton uptake and release, proton transfer along a hydrogen bond, and rotations of water molecules that constitute a proton wire through the channel. A simulation of 8 mus length took about 5 min on an Intel Pentium 4 CPU with 3.2 GHz. We obtained good agreement with experimental data for the proton flux through gramicidin A over a wide range of pH values and membrane potentials. We find that proton desolvation as well as water rotations are equally important for the proton transfer through gramicidin A at physiological membrane potentials. Our method allows to simulate long-range charge transfer in biological systems at time scales, which are not accessible by other methods.     

470—The proton pump at work: Part I. The basic concept of excess proton/defect proton translocation

The Nagle-Morowitz proton pump, which is based on proton transport in water and ice, is shown to be inapplicable to weakly H-bonded proton transport systems. It is demonstrated that in weakly or non H-bonded systems protons can migrate either as excess protons H using a high-lying proton conduction band or as defect protons H′ using intermediate energy levels. As H? the protons act as positive charge carriers endowed with a very high mobility. As H′ they act as negative charge carriers. The H′ transport mechanism involves thermally activated proton tunneling. Owing to the very large mobility differences between H? and H′, > 10⁶, very high potentials can be self-generated in any situation which creates a concentration gradient. Proton conductivity data on inorganic model compounds are presented. Applying these results to proton transport across biomembranes, transmembrane potentials U_m, acidification Δ pH and transient phenomena can be explained as result of H? and H′ translocation.     

A polarizable multistate empirical valence bond model for proton transport in aqueous solution

A multistate empirical valence bond model for proton transport in water, which explicitly includes solvent polarization, is presented. Polarization is included for each valence-bond state via induced point dipoles, and the model is parametrized to be used with an effective path integral derived potential surface, so as to include quantum effects of the transferring proton. The new model is shown to reproduce ab initio geometries and energetics for small protonated clusters. It is also shown that the new model gives a diffusion constant for the excess proton in water, which is in good agreement with experiment, and that the qualitative features of ab initio path integral simulations [D. Marx, M. E. Tuckerman, J. Hutter, and M. Parrinello, Nature (London) 397, 601 (1999)] are well reproduced.     

Excited-state double proton transfer dynamics of model DNA base pairs: 7-hydroxyquinoline dimers

The excited-state double proton transfer of model DNA base pairs, 7-hydroxyquinoline dimers, in benzene has been investigated using picosecond time-resolved fluorescence spectroscopy. Upon excitation, whereas singly hydrogen-bonded noncyclic dimers do not go through tautomerization within the relaxation time of 1400 ps, doubly hydrogen-bonded cyclic dimers undergo excited-state double proton transfer on the time scale of 25 ps to form tautomeric dimers, which subsequently undergo a conformational change in 180 ps to produce singly hydrogen-bonded tautomers. The rate constant of the double proton transfer reaction is temperature-independent, showing a large kinetic isotope effect of 5.2, suggesting that the rate is governed mostly by tunneling.     

On the maximum charge state and proton transfer reactivity of peptide and protein ions formed by electrospray ionization

A relatively simple model for calculation of the energetics of gas-phase proton transfer reactions and the maximum charge state of multiply protonated ions formed by electrospray ionization is presented. This model is based on estimates of the intrinsic proton transfer reactivity of sites of protonation and point charge Coulomb interactions. From this model, apparent gas-phase basicities (GB^app) of multiply protonated ions are calculated. Comparison of this value to the gas-phase basicity of the solvent from which an ion is formed enables a maximum charge state to be calculated. For 13 commonly electrosprayed proteins, our calculated maximum charge states are within an average of 6% of the experimental values reported in the literature. This indicates that the maximum charge state for proteins is determined by their gas-phase reactivity. Similar results are observed for peptides with many basic residues. For peptides with few basic residues, we find that the maximum charge state is better correlated to the charge state in solution. For low charge state ions, we find that the most basic sites Arg, Lys, and His are preferentially protonated. A significant fraction of the less basic residues Pro, Trp, and Gln are protonated in high charge state ions. The calculated GB^app of individual protonation sites varies dramatically in the high charge state ions. From these values, we calculate a reduced cross section for proton transfer reactivity that is significantly lower than the Langevin collision frequency when the GB^app of the ion is approximately equal to the GB of the neutral base.     

Factors influencing proton positions in biomolecules

Results of quantum mechanical calculations are presented that suggest a number of mechanisms whereby protons may be shifted from one group to another along an H bond. The first factor to be considered is a stretching of the bond that drastically raises the energy barrier to transfer. It is possible to predict barriers for an arbitrary system based only on results for a simple system and knowledge of the relevant bond length in the isolated subsystems. Factors that increase the intrinsic basicity of the B group in A-H-B lead not only to a lowering of the energy of the A-HB state relative to AH-B but also to a reduction in the barrier to transfer of the proton from A to B. Ions in the vicinity of the H bond exert a powerful influence and can shift the proton to the less basic group across a gradient of several pK units. Rather than shielding the proton from the external ion, the H bond acts instead to amplify the effects of the electric field. Reorientation of the A and B groups relative to one another, i.e., bends of the H bond, also produce surprisingly large changes in the relative energies of the AH-B and A-HB states. Such bends are capable of pushing the proton across to the normally less basic group, providing a mechanism of coupling conformational changes to proton ‘pumping’ activity. It is found that the high and low pH states of a given H bond can have dramatically differnt relative populations of the AH-B and A-HB configurations. These observations are explained in terms of fundamental concepts involving electrostatic interaction energies.     

Structural character and energetics of tyrosyl radical formation by electron/proton transfers of a covalently linked histidine-tyrosine: a model for cytochrome C oxidase

The structural, energetic, and electronic and IR spectroscopic properties for a model of the cross-linked histidine-tyrosine (His-Tyr) residues as found in cytochrome c oxidase (CcO) are investigated by ab initio methods. The formation of a His-Tyr radical is studied by two paths: proton release followed by electron release and vice versa. The energetics for the proton/electron releases of the Tyr depend modestly on the cross-linked His substituent and, more sensitively, on the charge of the cation attached to the imino N site of the His residue. Protonation of the imino N site significantly increases the electron ionization potential and decreases the proton dissociation energy, making them competitive processes. A positive charge placed at the imino N site, whose value is scanned from zero to one, shows a continuous increase in ionization potential and a decrease in proton dissociation energy, with the +1 limit agreeing well with the protonated imino N site result, indicating a dominant electrostatic effect. The charge populations and the spin density distributions of the His-Tyr model, the radical cation formed by electron ionization, the anion formed by proton dissociation, and the final His-Tyr radical depend sensitively on the substituents, implying a modulation role on the charge transfer between the phenol and imidazole rings, especially for the charged species. His-Tyr and protonated His-Tyr exhibit differences among their respective structural isomers with consequences on their IR absorptions. Small barriers between their pseudo-cis and pseudo-trans rotamers demonstrate the relative flexibility between the two rings, and these may facilitate proton release and charge transfer. The cation effect demonstrates that the cationized cross-linked His-Tyr should be the best candidate to mimic the covalently ring-linked histidine-tyrosine structure in CcO.