Proton transfer along the external proton channel of bacteriorhodopsin

The process of proton transfer along a proton channel is considered using bacteriorhodopsin as a model system, for which a large body of experimental data is available. The possible amino acid composition of the external proton half-channel of bacteriorhodopsin and the stepwise scheme of proton transfer consistent with experimental data are proposed. The rate of proton transfer between fixed centers is assessed for certain regions of this channel for which spectroscopic data are available.     

Quantum-mechanical descriptions of proton transfer in biosystems containing hydrogen-bonded chains

A mechanism of proton transfer along the proton channel of (F0) ATP-synthetase of a membrane is suggested. In the small polaron model the charged fault (an excess of protons or a proton hole) transfer is considered in a longitudinal electric field along an assumed chain which is formed by hydroxyl groups connected by strong H-bonded chains. A number of kinetic parameters are estimated. The theoretical data are compared with the experimental results.     

Ab initio molecular dynamics study of proton transfer in a polyglycine analog of the ion channel gramicidin A.

Proton transfer in biological systems is thought to often proceed through hydrogen-bonded chains of water molecules. The ion channel, gramicidin A (gA), houses within its helical structure just such a chain. Using the density functional theory based ab initio molecular dynamics Car-Parrinello method, the structure and dynamics of proton diffusion through a polyglycine analog of the gA ion channel has been investigated. In the channel, a proton, which is initially present as hydronium (H3O+), rapidly forms a strong hydrogen bond with a nearest neighbor water, yielding a transient H5O2+ complex. As in bulk water, strong hydrogen bonding of this complex to a second neighbor solvation shell is required for proton transfer to occur. Within gA, this second neighbor shell included not only a channel water molecule but also a carbonyl of the channel backbone. The present calculations suggest a transport mechanism in which a priori carbonyl solvation is a requirement for proton transfer.     

The water channel of cytochrome c oxidase: inferences from inhibitor studies.

Cytochrome c oxidase couples electron transfer to proton transfer from inside the mitochondrion to the cytosol. Protons pass through a channel; it is closed except when protons are pumped. Electron transfer is also coupled to a water cycle. Water moves into and out of the oxidase during electron transfer, presumably through a channel. The three processes are coupled because of the common dependence on electron transfer. If water and protons had to pass through the same channel for the proton to pass, it might be possible to block the pore by entraining small molecules in the flow. The data in this report indicate that there is a correlation between the ability of a compound to inhibit the oxidase and its size. Formamide and formaldehyde are potent inhibitors. Larger and smaller molecules are poor inhibitors. Formamide introduces an internal block in electron transfer. It is a slow-onset, reversible inhibitor, dependent on turnover to manifest its effects. Vesicular oxidase is less influenced by formamide than is soluble oxidase; formamide must pass a permeability barrier to act. The data are consistent with a proton channel with constrictions at both ends that open to yield a pore of approximately 4 A.     

A role for the protein in internal electron transfer to the catalytic center of cytochrome c oxidase

Internal electron transfer (ET) to heme a(3) during anaerobic reduction of oxidized bovine heart cytochrome c oxidase (CcO) was studied under conditions where heme a and Cu(A) were fully reduced by excess hexaamineruthenium. The data show that ET to heme a(3) is controlled by the state of ionization of a single protolytic residue with a pK(a) of 6.5 +/- 0.2. On the basis of the view that ET to the catalytic site is limited by coupled proton transfer, this pK(a) was attributed to Glu60 which is located at the entrance of the proton-conducting K channel on the matrix side of CcO. It is proposed that Glu60 controls proton entry into the channel. However, even with this channel open, there is the second factor that regulates ET, and this is ascribed to the rate of proton diffusion in the channel. In addition, it is concluded that proton transfer in the K channel is reversibly inhibited by the detergent Triton X-100. It is also found that the rate of ET to heme a(3) in the as-isolated resting enzyme and in CcO "activated" by reaction of fully reduced enzyme with O(2) is the same, implying that the catalytic sites of these two forms of oxidized enzyme are essentially identical.     

Water as a structural element in a channel: gating in the Kcsa channel, and implications for voltage-gated ion channels

Water is becoming understood as a structural element in proteins. Here we are concerned with one particular type of protein, ion channels. The S. Lividans KcsA K(+) channel, the X-ray structure of which is known, is gated by protons (i.e, by a drop in pH). Ab initio calculations suggest that an H(5)O(2) group, partially charged, connects the E118 residues in the gating region, when the four residues have a -2 net charge, but that the hydrogen bonding is not strong enough to do this when the charge becomes -1. The H(5)O(2) group would block the channel, in the -2 state, and prevent motion of the four transmembrane (TM) segments of the protein, by binding them. With the weaker bond in the -1 state, the TM segments would be able to separate (as they have been found to do experimentally, opening the channel. Voltage gated channels have four additional TM segments for each of the four domains of the channel protein. These appear to allow motion of protons; in fact there is evidence that the initial step in gating must be the transfer of a proton. We have earlier shown that the transfer of a single proton between two methylamines under the influence of a field is possible, as proton tunneling. Subsequent steps are hypothesized to result from four proton transfer cascades of about three protons each, triggered by the initial proton transfer. We suggest that the extra 4 TM segments of the voltage gated channel act as a voltage to proton-current transducer. Water, held by hydrogen bonds, is also suggested as the source of the accessibility data found with MTS reagents, based largely on simulations, our earlier Monte Carlo simulations as well as molecular dynamics studies reported by others. These waters may also play a structural role in the protein.     

Coupling of light-induced electron transfer to proton uptake in photosynthesis

Light energy is transformed into chemical energy in photosynthesis by coupling a light-induced electron transfer to proton uptake. The resulting proton gradient drives ATP synthesis. In this study, we monitored the light-induced reactions in a 100-kDa photosynthetic protein from 30 ns to 35 s by FTIR difference spectroscopy. The results provide detailed mechanistic insights into the electron and proton transfer reactions of the QA to QB transition: reduction of QA in picoseconds induces protonation of histidines, probably of His126 and His128 in the H subunit at the entrance of the proton uptake channel, and of Asp210 in the L subunit inside the channel at 12 micros and 150 micros. This seems to be a prerequisite for the reduction of QB, mainly at 150 micros. QA- is reoxidized at 1.1 ms, and a proton is transferred from Asp210 to Glu212 in the L subunit, the proton donor to QB-. Notably, our data indicate that QB is not reduced directly by QA- but presumably through an intermediary electron donor.     

Dynamic mechanism of proton transfer in mannitol 2-dehydrogenase from Pseudomonas fluorescens: mobile GLU292 controls proton relay through a water channel that connects the active site with bulk solvent

The active site of mannitol 2-dehydrogenase from Pseudomonas fluorescens (PfM2DH) is connected with bulk solvent through a narrow protein channel that shows structural resemblance to proton channels utilized by redox-driven proton pumps. A key element of the PfM2DH channel is the "mobile" Glu(292), which was seen crystallographically to adopt distinct positions up and down the channel. It was suggested that the "down → up" conformational change of Glu(292) could play a proton relay function in enzymatic catalysis, through direct proton shuttling by the Glu or because the channel is opened for water molecules forming a chain along which the protons flow. We report evidence from site-directed mutagenesis (Glu(292) → Ala) substantiated by data from molecular dynamics simulations that support a role for Glu(292) as a gate in a water chain (von Grotthuss-type) mechanism of proton translocation. Occupancy of the up and down position of Glu(292) is influenced by the bonding and charge state of the catalytic acid base Lys(295), suggesting that channel opening/closing motions of the Glu are synchronized to the reaction progress. Removal of gatekeeper control in the E292A mutant resulted in a selective, up to 120-fold slowing down of microscopic steps immediately preceding catalytic oxidation of mannitol, consistent with the notion that formation of the productive enzyme-NAD(+)-mannitol complex is promoted by a corresponding position change of Glu(292), which at physiological pH is associated with obligatory deprotonation of Lys(295) to solvent. These results underscore the important role of conformational dynamics in the proton transfer steps of alcohol dehydrogenase catalysis.     

Replacing Asn207 by aspartate at the neck of the D channel in the aa3-type cytochrome c oxidase from Rhodobacter sphaeroides results in decoupling the proton pump

Cytochrome oxidase catalyzes the reduction of O2 to water and conserves the considerable free energy available from this reaction in the form of a proton motive force. For each electron, one proton is electrogenically pumped across the membrane. Of particular interest is the mechanism by which the proton pump operates. Previous studies of the oxidase from Rhodobacter sphaeroides have shown that all of the pumped protons enter the enzyme through the D channel and that a point mutant, N139D, in the D channel completely eliminates proton pumping without reducing oxidase activity. N139 is one of three asparagines near the entrance of the D channel, where there is a narrowing or neck, through which a single file of water molecules pass. In the current work, it is shown that replacement of a second asparagine in this region by an asparate, N207D, also decouples the proton pump without altering the oxidase activity of the enzyme. Previous studies demonstrated that the N139D mutant results in an increase in the apparent pKa of E286, a functionally critical residue that is located 20 A away from N139 at the opposite end of the D channel. In the current work, it is shown that the N207 mutation also increases the apparent pKa of E286. This finding reinforces the proposal that the elimination of proton pumping is the result of an increase of the apparent proton affinity of E286, which, in turn, prevents the timely proton transfer to a proton accepter group within the exit channel of the proton pump.     

Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study.

The proton transfer activity of the light-driven proton pump, bacteriorhodopsin (bR) in the photochemical cycle might imply internal water molecules. The free energy of inserting water molecules in specific sites along the bR transmembrane channel has been calculated using molecular dynamics simulations based on a microscopic model. The existence of internal hydration is related to the free energy change on transfer of a water molecule from bulk solvent into a specific binding site. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each hydrated model from molecular dynamics simulations of the creation of water molecules into specific protein-binding sites. A rigorous statistical mechanical formulation allowing the calculation of the free energy of transfer of water molecules from the bulk to a protein cavity is used to estimate the probabilities of occupancy in the putative bR proton channel. The channel contains a region lined primarily by nonpolar side-chains. Nevertheless, the results indicate that the transfer of four water molecules from bulk water to this apparently hydrophobic region is thermodynamically permitted. The column forms a continuous hydrogen-bonded chain over 12 A between a proton donor, Asp 96, and the retinal Schiff base acceptor. The presence of two water molecules in direct hydrogen-bonding association with the Schiff base is found to be strongly favorable thermodynamically. The implications of these results for the mechanism of proton transfer in bR are discussed.     

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.     

Solution NMR structure of the V27A drug resistant mutant of influenza A M2 channel

The M2 protein of influenza A virus forms a proton-selective channel that is required for viral replication. It is the target of the anti-influenza drugs, amantadine and rimantadine. Widespread drug resistant mutants, however, has greatly compromised the effectiveness of these drugs. Here, we report the solution NMR structure of the highly pathogenic, drug resistant mutant V27A. The structure reveals subtle structural differences from wildtype that maybe linked to drug resistance. The V27A mutation significantly decreases hydrophobic packing between the N-terminal ends of the transmembrane helices, which explains the looser, more dynamic tetrameric assembly. The weakened channel assembly can resist drug binding either by destabilizing the rimantadine-binding pocket at Asp44, in the case of the allosteric inhibition model, or by reducing hydrophobic contacts with amantadine in the pore, in the case of the pore-blocking model. Moreover, the V27A structure shows a substantially increased channel opening at the N-terminal end, which may explain the faster proton conduction observed for this mutant. Furthermore, due to the high quality NMR data recorded for the V27A mutant, we were able to determine the structured region connecting the channel domain to the C-terminal amphipathic helices that was not determined in the wildtype structure. The new structural data show that the amphipathic helices are packed much more closely to the channel domain and provide new insights into the proton transfer pathway.     

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.     

Functional properties of threefold and fourfold channels in ferritin deduced from electrostatic calculations

The iron storage protein ferritin contains threefold and fourfold symmetric channels that are thought to provide pathways for the transfer of Fe(2+) ions in and out of the protein. Using the known crystal structure of the ferritin protein, we perform electrostatic potential energy calculations to elucidate the functional properties of these channels. The threefold channel is shown to be responsible for the transit of Fe(2+) ions. Monovalent ions can also diffuse through the threefold channel but presence of divalent ions in the pore retards this process leading to a selectivity mechanism similar to the one observed in calcium channels. The fourfold channel is found to be impermeant to all cations with the possible exception of protons. Because proton transfer is essential to maintain the electroneutrality of the protein during iron deposition, we suggest that the function of the fourfold channel is to form a "proton wire" that facilitates their transfer in and out of ferritin.     

Time-resolved detection of transient movement of helices F and G in doubly spin-labeled bacteriorhodopsin

Photo-excited structural changes of the light-driven proton pump bacteriorhodopsin were monitored using double-site-directed spin labeling combined with electron paramagnetic resonance (EPR) spectroscopy. The inter-spin distances between nitroxides attached at residue positions 100 and 226, 101 and 160, and 101 and 168 were determined for the BR initial state and the trapped M photo-intermediate. Distance changes that occur during the photocycle were followed with millisecond time resolution under physiological conditions at 293 K. The kinetic analysis of the EPR data and comparison with the absorbance changes in the visible spectrum reveal an outward movement of helix F during the late M intermediate and a subsequent approach of helix G toward the proton channel. The displacements of the cytoplasmic moieties of these helices amount to 0.1-0.2 nm. We propose that the resulting opening of the proton channel decreases the pK of the proton donor D96 and facilitates proton transfer to the Schiff base during the M-to-N transition.     

Crystal structure of transhydrogenase domain III at 1.2 A resolution

The nicotinamide nucleotide transhydrogenases (TH) of mitochondria and bacteria are membrane-intercalated proton pumps that transduce substrate binding energy and protonmotive force via protein conformational changes. In mitochondria, TH utilizes protonmotive force to promote direct hydride ion transfer from NADH to NADP, which are bound at the distinct extramembranous domains I and III, respectively. Domain II is the membrane-intercalated domain and contains the enzyme's proton channel. This paper describes the crystal structure of the NADP(H) binding domain III of bovine TH at 1.2 A resolution. The structure reveals that NADP is bound in a manner inverted from that previously observed for nucleotide binding folds. The non-classical binding mode exposes the NADP(H) nicotinamide ring for direct contact with NAD(H) in domain I, in accord with biochemical data. The surface of domain III surrounding the exposed nicotinamide is comprised of conserved residues presumed to form the interface with domain I during hydride ion transfer. Further, an adjacent region contains a number of acidic residues, forming a surface with negative electrostatic potential which may interact with extramembranous loops of domain II. Together, the distinctive surface features allow mechanistic considerations regarding the NADP(H)-promoted conformation changes that are involved in the interactions of domain III with domains I and II for hydride ion transfer and proton translocation.     

Proton Transfer in the K-Channel Analog of B-Type Cytochrome c Oxidase from Thermus thermophilus

A key enzyme in aerobic metabolism is cytochrome c oxidase (CcO), which catalyzes the reduction of molecular oxygen to water in the mitochondrial and bacterial membranes. Substrate electrons and protons are taken up from different sides of the membrane and protons are pumped across the membrane, thereby generating an electrochemical gradient. The well-studied A-type CcO uses two different entry channels for protons: the D-channel for all pumped and two consumed protons, and the K-channel for the other two consumed protons. In contrast, the B-type CcO uses only a single proton input channel for all consumed and pumped protons. It has the same location as the A-type K-channel (and thus is named the K-channel analog) without sharing any significant sequence homology. In this study, we performed molecular-dynamics simulations and electrostatic calculations to characterize the K-channel analog in terms of its energetic requirements and functionalities. The function of Glu-15B as a proton sink at the channel entrance is demonstrated by its rotational movement out of the channel when it is deprotonated and by its high pK_A value when it points inside the channel. Tyr-244 in the middle of the channel is identified as the valve that ensures unidirectional proton transfer, as it moves inside the hydrogen-bond gap of the K-channel analog only while being deprotonated. The electrostatic energy landscape was calculated for all proton-transfer steps in the K-channel analog, which functions via proton-hole transfer. Overall, the K-channel analog has a very stable geometry without large energy barriers.     

What really prevents proton transport through aquaporin? Charge self-energy versus proton wire proposals

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.     

Proton Transfer in the K-Channel Analog of B-Type Cytochrome c Oxidase from Thermus thermophilus

A key enzyme in aerobic metabolism is cytochrome c oxidase (CcO), which catalyzes the reduction of molecular oxygen to water in the mitochondrial and bacterial membranes. Substrate electrons and protons are taken up from different sides of the membrane and protons are pumped across the membrane, thereby generating an electrochemical gradient. The well-studied A-type CcO uses two different entry channels for protons: the D-channel for all pumped and two consumed protons, and the K-channel for the other two consumed protons. In contrast, the B-type CcO uses only a single proton input channel for all consumed and pumped protons. It has the same location as the A-type K-channel (and thus is named the K-channel analog) without sharing any significant sequence homology. In this study, we performed molecular-dynamics simulations and electrostatic calculations to characterize the K-channel analog in terms of its energetic requirements and functionalities. The function of Glu-15B as a proton sink at the channel entrance is demonstrated by its rotational movement out of the channel when it is deprotonated and by its high pK_A value when it points inside the channel. Tyr-244 in the middle of the channel is identified as the valve that ensures unidirectional proton transfer, as it moves inside the hydrogen-bond gap of the K-channel analog only while being deprotonated. The electrostatic energy landscape was calculated for all proton-transfer steps in the K-channel analog, which functions via proton-hole transfer. Overall, the K-channel analog has a very stable geometry without large energy barriers.     

Structural elements involved in proton translocation by cytochrome c oxidase as revealed by backbone amide hydrogen-deuterium exchange of the E286H mutant

Cytochrome c oxidase is the terminal electron acceptor in the respiratory chains of aerobic organisms and energetically couples the reduction of oxygen to water to proton pumping across the membrane. The mechanisms of proton uptake, gating, and pumping have yet to be completely elucidated at the molecular level for these enzymes. For Rhodobacter sphaeroides CytcO (cytochrome aa3), it appears as though the E286 side chain of subunit I is a branching point from which protons are shuttled either to the catalytic site for O2 reduction or to the acceptor site for pumped protons. Amide hydrogen-deuterium exchange mass spectrometry was used to investigate how mutation of this key branching residue to histidine (E286H) affects the structures and dynamics of four redox intermediate states. A functional characterization of this mutant reveals that E286H CytcO retains approximately 1% steady-state activity that is uncoupled from proton pumping and that proton transfer from H286 is significantly slowed. Backbone amide H-D exchange kinetics indicates that specific regions of CytcO, perturbed by the E286H mutation, are likely to be involved in proton gating and in the exit pathway for pumped protons. The results indicate that redox-dependent conformational changes around E286 are essential for internal proton transfer. E286H CytcO, however, is incapable of these specific conformational changes and therefore is insensitive to the redox state of the enzyme. These data support a model where the side chain conformation of E286 controls proton translocation in CytcO through its interactions with the proton gate, which directs the flow of protons either to the active site or to the exit pathway. In the E286H mutant, the proton gate does not function properly and the exit channel is unresponsive. These results provide new insight into the structure and mechanism of proton translocation by CytcO.