A QM/MM study of proton transport pathways in a [NiFe] hydrogenase

A theoretical QM/MM study of the [NiFe] hydrogenase from Desulfovibrio fructosovorans has been performed to investigate possible routes of proton transfer between the active site and the protein surface. We obtained the minimum energy paths, with a modified version of the nudged elastic band method, for a set of proposed pathways. The calculations were carried out for the crystallographic structure and for several structures of the protein obtained from a molecular dynamics simulation. The results show one of the studied pathways to be preferred for transport from the active site to the surface, but the preference is not so strong when transport occurs in the opposite direction.     

Molecular dynamics study of the proposed proton transport pathways in [FeFe]-hydrogenase

Possible proton transport pathways in Clostridium pasteurianum (CpI) [FeFe]-hydrogenase were investigated with molecular dynamics simulations. This study was undertaken to evaluate the functional pathway and provide insight into the hydrogen bonding features defining an active proton transport pathway. Three pathways were evaluated, two of which consist of water wires and one of predominantly amino acid residues. Our simulations suggest that protons are not transported through water wires. Instead, the five-residue motif (Glu282, Ser319, Glu279, H₂O, Cys299) was found to be the likely pathway, consistent with previously made experimental observations. The pathway was found to have a persistent hydrogen bonded core (residues Cys299 to Ser319), with less persistent hydrogen bonds at the ends of the pathway for both H₂ release and H₂ uptake. Single site mutations of the four residues have been shown experimentally to deactivate the enzyme. The theoretical evaluation of these mutations demonstrates redistribution of the hydrogen bonds in the pathway, resulting in enzyme deactivation. Finally, coupling between the protein dynamics near the proton transport pathway and the redox partner binding regions was also found as a function of H₂ uptake and H₂ release states, which may be indicative of a correlation between proton and electron movement within the enzyme.     

Water dynamics simulation as a tool for probing proton transfer pathways in a heptahelical membrane protein

The proton transfer pathway in a heptahelical membrane protein, the light-driven proton pump bacteriorhodopsin (BR), is probed by a combined approach of structural analysis of recent X-ray models and molecular dynamics (MD) simulations that provide the diffusion pathways of internal and external water molecules. Analyzing the hydrogen-bond contact frequencies of the water molecules with protein groups, the complete proton pathway through the protein is probed. Beside the well-known proton binding sites in the protein interior-the protonated Schiff base, Asp85 and Asp96, and the H(5)O(2) (+) complex stabilized by Glu204 and Glu194-the proton release and uptake pathways to the protein surfaces are described in great detail. Further residues were identified, by mutation of which the proposed pathways can be verified. In addition the diffusion pathway of water 502 from Lys216 to Asp96 is shown to cover the positions of the intruding waters 503 and 504 in the N-intermediate. The transiently established water chain in the N-state provides a proton pathway from Asp96 to the Schiff base in the M- to N-transition in a Grotthus-like mechanism, as concluded earlier from time-resolved Fourier transform infrared experiments [le Coutre et al., Proc Nat Acad Sci USA 1995;92:4962-4966].     

The Nitric-oxide Reductase from Paracoccus denitrificans Uses a Single Specific Proton Pathway

The NO reductase from Paracoccus denitrificans reduces NO to N₂O (2NO + 2H⁺ + 2e⁻ → N₂O + H₂O) with electrons donated by periplasmic cytochrome c (cytochrome c-dependent NO reductase; cNOR). cNORs are members of the heme-copper oxidase superfamily of integral membrane proteins, comprising the O₂-reducing, proton-pumping respiratory enzymes. In contrast, although NO reduction is as exergonic as O₂ reduction, there are no protons pumped in cNOR, and in addition, protons needed for NO reduction are derived from the periplasmic solution (no contribution to the electrochemical gradient is made). cNOR thus only needs to transport protons from the periplasm into the active site without the requirement to control the timing of opening and closing (gating) of proton pathways as is needed in a proton pump. Based on the crystal structure of a closely related cNOR and molecular dynamics simulations, several proton transfer pathways were suggested, and in principle, these could all be functional. In this work, we show that residues in one of the suggested pathways (denoted pathway 1) are sensitive to site-directed mutation, whereas residues in the other proposed pathways (pathways 2 and 3) could be exchanged without severe effects on turnover activity with either NO or O₂. We further show that electron transfer during single-turnover reduction of O₂ is limited by proton transfer and can thus be used to study alterations in proton transfer rates. The exchange of residues along pathway 1 showed specific slowing of this proton-coupled electron transfer as well as changes in its pH dependence. Our results indicate that only pathway 1 is used to transfer protons in cNOR.     

Proton transfer in carbonic anhydrase is controlled by electrostatics rather than the orientation of the acceptor

Combined quantum mechanical/molecular mechanical (QM/MM) simulations are carried out to analyze factors that dictate the proton transfer in carbonic anhydrase II (CAII), an enzyme that has been used as a prototypical example of long-range proton transfers in biomolecules. In contrast to the long-held conjecture in the experimental literature, the computed potentials of mean force (PMF) suggest that the proton transfer in CAII is not very sensitive to the orientation of the acceptor group (His 64) and, therefore, the number of water molecules that bridge the donor (zinc-water) and acceptor groups. Perturbative analysis indicates that a series of polar and charged residues close to the transfer pathways make the dominant contribution to the barrier and exothermicity of the proton transfer reaction, thus supporting the proposal from previous studies of Warshel and co-workers using a somewhat simpler QM/MM model that electrostatic interactions play a major role in the proton transfer in CAII. The PMF results are in striking contrast to previous analysis using the same QM/MM method but an ensemble of minimum energy path (MEP) calculations, which found a steep dependence of the barrier height on the number of bridging water molecules. Analysis of the configurations sampled in the PMF and MEP simulations suggests that this difference arises because the PMF simulations sample a largely stepwise mechanism while the local MEP calculations artificially favored concerted transfers due to the specific protocol used to generate the initial configurations. Therefore, this study presents a compelling argument for carrying out proper conformational sampling in the study of long-range proton transfers. Finally, we illustrate that Phi analysis, which has been widely used in protein folding studies, can potentially generate new mechanistic information for long-range proton transfers regarding the sequence of events. The results of the perturbation analysis and the Phi analysis provide opportunities for experimentally testing the mechanistic proposals from this study and our recent work in which a stepwise "proton hole" transfer pathway has been proposed.     

Proton transfer pathways in the mutant His-64-Ala of human carbonic anhydrase II

We have investigated the possible proton transfer pathways from the surface of the protein to the zinc-bound water molecule in the mutant His-64-Ala of human carbonic anhydrase II. Starting with an input of known crystallographic structures of the mutant, we model the proton pathways as hydrogen-bonded networks of proton conducting groups and bound solvent molecules. No proton path is detected in the mutant, in close agreement with the experimental observation of a 20-fold decrease in its catalytic efficiency compared to the wild-type enzyme. We also investigate in detail changes in hydration structure at the active site of the mutant and the resulting proton paths in the presence of an exogenous proton donor 4-methylimidazole (4-MI). The proton transfer pathways thus detected are correlated to the observed chemical rescue of catalytic activity by 4-MI.     

Application of classical molecular dynamics for evaluation of proton transfer mechanism on a protein

Proton transfer reactions on surfaces are prevalent in biology, chemistry and physics. In the present study, we employed classical Molecular Dynamics simulations to search for the presence of transient configurations that enable proton transfer, or proton sharing, between adjacent carboxylate groups on the protein surface. The results demonstrate that, during random fluctuations of the residues on the surface, there are repeated situations in which nearby carboxylates either share a common proton through a hydrogen bond, or are connected by a few water molecules that form conducting networks. These networks do not extend out of the common Coulomb cage of the participating residues and the lifetimes of the bridged structures are sufficiently long to allow passage of a proton between the carboxylates. The detection of domains capable of supporting a rapid proton transfer on a protein supports the notion that clusters of carboxylates are the operative elements of proton collecting antennae, as in bacteriorhodopsin, cytochrome c oxidase or the photosynthetic reaction center.     

Analyses of the large subunit histidine-rich motif expose an alternative proton transfer pathway in [NiFe] hydrogenases

A highly conserved histidine-rich region with unknown function was recognized in the large subunit of [NiFe] hydrogenases. The HxHxxHxxHxH sequence occurs in most membrane-bound hydrogenases, but only two of these histidines are present in the cytoplasmic ones. Site-directed mutagenesis of the His-rich region of the T. roseopersicina membrane-attached Hyn hydrogenase disclosed that the enzyme activity was significantly affected only by the replacement of the His104 residue. Computational analysis of the hydrogen bond network in the large subunits indicated that the second histidine of this motif might be a component of a proton transfer pathway including Arg487, Asp103, His104 and Glu436. Substitutions of the conserved amino acids of the presumed transfer route impaired the activity of the Hyn hydrogenase. Western hybridization was applied to demonstrate that the cellular level of the mutant hydrogenases was similar to that of the wild type. Mostly based on theoretical modeling, few proton transfer pathways have already been suggested for [NiFe] hydrogenases. Our results propose an alternative route for proton transfer between the [NiFe] active center and the surface of the protein. A novel feature of this model is that this proton pathway is located on the opposite side of the large subunit relative to the position of the small subunit. This is the first study presenting a systematic analysis of an in silico predicted proton translocation pathway in [NiFe] hydrogenases by site-directed mutagenesis.     

Three ways in, one way out: water dynamics in the trans-membrane domains of the inner membrane translocase AcrB

Powered by proton-motive force, the inner membrane translocase AcrB is the engine of the AcrAB-TolC efflux pump in Escherichia coli. As proton conduction in proteins occurs along hydrogen-bonded networks of polar residues and water molecules, knowledge of the protein-internal water distribution and water-interacting residues allows drawing conclusions to possible pathways of proton conduction. Here, we report a series of 6× 50 ns independent molecular dynamics simulations of asymmetric AcrB embedded in a phospholipid/water environment. Simulating each monomer in its proposed protonation state, we calculated for each trans-membrane domain the average water distribution, identified residues interacting with these waters and quantified each residue's frequency of water hydrogen bond contact. Combining this information we find three possible routes of proton transfer connecting a continuously hydrated region of known key residues in the TMD interior to bulk water by one cytoplasmic and up to three periplasm water channels in monomer B and A. We find that water access of the trans-membrane domains is regulated by four groups of residues in a combination of side chain re-orientations and shifts of trans-membrane helices. Our findings support a proton release event via Arg971 during the C intermediate or in the transition to A, and proton uptake occurring in the A or B state or during a so far unknown intermediate in between B and C where cytoplasmic water access is still possible. Our simulations suggest experimentally testable hypotheses, which have not been investigated so far.     

Structural features of [NiFeSe] and [NiFe] hydrogenases determining their different properties: a computational approach

Hydrogenases are metalloenzymes that catalyze the reversible reaction \textH₂ \leftrightarrows 2\textH ⁺ + 2\texte ^- {\text{H}}_{2} \leftrightarrows 2{\text{H}}^{ + } + 2{\text{e}}^{ - } , being potentially useful in H₂ production or oxidation. [NiFeSe] hydrogenases are a particularly interesting subgroup of the [NiFe] class that exhibit tolerance to O₂ inhibition and produce more H₂ than standard [NiFe] hydrogenases. However, the molecular determinants responsible for these properties remain unknown. Hydrophobic pathways for H₂ diffusion have been identified in [NiFe] hydrogenases, as have proton transfer pathways, but they have never been studied in [NiFeSe] hydrogenases. Our aim was, for the first time, to characterize the H₂ and proton pathways in a [NiFeSe] hydrogenase and compare them with those in a standard [NiFe] hydrogenase. We performed molecular dynamics simulations of H₂ diffusion in the [NiFeSe] hydrogenase from Desulfomicrobium baculatum and extended previous simulations of the [NiFe] hydrogenase from Desulfovibrio gigas (Teixeira et al. in Biophys J 91:2035–2045, 2006). The comparison showed that H₂ density near the active site is much higher in [NiFeSe] hydrogenase, which appears to have an alternative route for the access of H₂ to the active site. We have also determined a possible proton transfer pathway in the [NiFeSe] hydrogenase from D. baculatum using continuum electrostatics and Monte Carlo simulation and compared it with the proton pathway we found in the [NiFe] hydrogenase from D. gigas (Teixeira et al. in Proteins 70:1010–1022, 2008). The residues constituting both proton transfer pathways are considerably different, although in the same region of the protein. These results support the hypothesis that some of the special properties of [NiFeSe] hydrogenases could be related to differences in the H₂ and proton pathways.     

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N(2)O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666-70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb(3) oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.     

Molecular dynamics of a protein surface: ion-residues interactions

Time-resolved measurements indicated that protons could propagate on the surface of a protein or a membrane by a special mechanism that enhanced the shuttle of the proton toward a specific site. It was proposed that a suitable location of residues on the surface contributes to the proton shuttling function. In this study, this notion was further investigated by the use of molecular dynamics simulations, where Na(+) and Cl(-) are the ions under study, thus avoiding the necessity for quantum mechanical calculations. Molecular dynamics simulations were carried out using as a model a few Na(+) and Cl(-) ions enclosed in a fully hydrated simulation box with a small globular protein (the S6 of the bacterial ribosome). Three independent 10-ns-long simulations indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. However, it was noted that close to some domains the ions extended their duration near the surface, thus suggesting that the local electrostatic potential hindered their diffusion to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of the molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.     

Simulating proton translocations in proteins: probing proton transfer pathways in the Rhodobacter sphaeroides reaction center.

A general method for simulating proton translocations in proteins and for exploring the role of different proton transfer pathways is developed and examined. The method evaluates the rate constants for proton transfer processes using the energetics of the relevant proton configurations. The energies (DeltaG((m))) of the different protonation states are evaluated in two steps. First, the semimicroscopic version of the protein dipole Langevin dipole (PDLD/S) method is used to evaluate the intrinsic energy of bringing the protons to their protein sites, when the charges of all protein ionized residues are set to zero. Second, the interactions between the charged groups are evaluated by using a Coulomb's Law with an effective dielectric constant. This approach, which was introduced in an earlier study by one of the authors of the current report, allows for a very fast determination of any DeltaG((m)) and for practical evaluation of the time-dependent proton population: That is, the rate constants for proton transfer processes are evaluated by using the corresponding DeltaG((m)) values and a Marcus type relationship. These rate constants are then used to construct a master equation, the integration of which by a fourth-order Runge-Kutta method yields the proton population as a function of time. The integration evaluates, 'on the fly,' the changes of the rate constants as a result of the time-dependent changes in charge-charge interaction, and this feature benefits from the fast determination of DeltaG((m)). The method is demonstrated in a preliminary study of proton translocation processes in the reaction center of Rhodobacter sphaeroides. It is found that proton transfer across water chains involves significant activation barriers and that ionized protein residues probably are involved in the proton transfer pathways. The potential of the present method in analyzing mutation experiments is discussed briefly and illustrated. The present study also examines different views of the nature of proton translocations in proteins. It is shown that such processes are controlled mainly by the electrostatic interaction between the proton site and its surroundings rather than by the local bond rearrangements of water molecules that are involved in the proton pathways. Thus, the overall rate of proton transport frequently is controlled by the highest barrier along the conduction pathway. Proteins 1999;36:484-500.     

Dynamical study,hydrogen bond analysis,and constant pH simulations of the beta carbonic anhydrase of Methanobacterium thermoautotrophicum

Within the five classes (α, β, γ, δ, and ζ) of carbonic anhydrases (CAs) the first two, containing mammal and plant representatives, are the most studied among all CAs. In this study, we have focused our investigation on the beta-class carbonic anhydrase of Methanobacterium thermoautotrophicum. We investigated both the importance of the Asp-Arg dyad near the catalytic zinc-bound water and the possible roles that water molecules within the active site and residues near the entrance of the catalytic cleft have on the first step of the enzyme’s reaction mechanism. Hydrogen-bonding analysis of selected residues within the active site and constant pH replica exchange molecular dynamics constant pH replica exchange simulations were performed. The latter was done in order to evaluate the pK_a values of possible proton acceptors. We found an intricate hydrogen-bonding network involving two acidic residues within the active site, Asp16 and Asp34, and the catalytic water molecule. We also observed a very strong interaction between the zinc-bound water and residues Asp34 and Arg36. This interaction was not significantly affected by the change in the protonation state of both the catalytic water and aspartate residue 34. The pK_a analysis show that the effect of the R36A mutation affects not only the possible proton acceptors, but also the catalytic water itself.     

Proton transfer within the active-site cavity of carbonic anhydrase III

The maximal turnover rate of CO2 hydration catalyzed by the carbonic anhydrases is limited by proton transfer steps from the zinc-bound water to solution, steps that regenerate the catalytically active zinc-bound hydroxide. Catalysis of CO2 hydration by wild-type human carbonic anhydrase III (HCA III) (k(cat) = 2 ms (-1)) is the least efficient among the carbonic anhydrases in its class, in part because it lacks an efficient proton shuttle residue. We have used site-directed mutagenesis to test positions within the active-site cavity of HCA III for their ability to carry out proton transfer by replacing various residues with histidine. Catalysis by wild-type HCA III and these six variants was determined from the initial velocity of hydration of CO2 measured by stopped-flow spectrophotometry and from the exchange of 18O between CO2 and H2O at chemical equilibrium by mass spectrometry. The results show that histidine at three positions (Lys64His, Arg67His and Phe131His) have the capacity to transfer protons during catalysis, enhancing maximal velocity of CO2 hydration and 18O exchange from 4- to 15-fold compared with wild-type HCA III. Histidine residues at the other three positions (Trp5His, Tyr7His, Phe20His) showed no firm evidence for proton transfer. These results are discussed in terms of the stereochemistry of the active-site cavity and possible proton transfer pathways.     

Proton paths in the sarcoplasmic reticulum Ca(2+) -ATPase

The sarcoplasmic reticulum Ca(2+)-ATPase (SERCA1a) pumps Ca(2+) and countertransport protons. Proton pathways in the Ca(2+) bound and Ca(2+)-free states are suggested based on an analysis of crystal structures to which water molecules were added. The pathways are indicated by chains of water molecules that interact favorably with the protein. In the Ca(2+) bound state Ca(2)E1, one of the proposed Ca(2+) entry paths is suggested to operate additionally or alternatively as proton pathway. In analogs of the ADP-insensitive phosphoenzyme E2P and in the Ca(2+)-free state E2, the proton path leads between transmembrane helices M5 to M8 from the lumenal side of the protein to the Ca(2+) binding residues Glu-771, Asp-800 and Glu-908. The proton path is different from suggested Ca(2+) dissociation pathways. We suggest that separate proton and Ca(2+) pathways enable rapid (partial) neutralization of the empty cation binding sites. For this reason, transient protonation of empty cation binding sites and separate pathways for different ions are advantageous for P-type ATPases in general.     

Slow proton transfer through the pathways for pumped protons in cytochrome c oxidase induces suicide inactivation of the enzyme

In the absence of subunit III the aa(3)-type cytochrome c oxidase exhibits a shortened catalytic life span (total number of turnovers) due to an increased probability of undergoing irreversible inactivation during steady-state turnover. Inactivation results from structural alteration of the heme a(3)-Cu(B) active site in subunit I [Hosler (2004) Biochim. Biophys. Acta 1655, 332-339]. The absence of subunit III also dramatically slows proton uptake to the active site via the D proton pathway, as well as inhibiting the proton backflow/exit pathway that connects the active site/proton pump with the outer surface of the oxidase complex. Here we demonstrate that these phenomena are linked: slow proton delivery to the active site through these pathways induces suicide inactivation, thus shortening the catalytic life span of the enzyme. Mutations that inhibit the D pathway, but not the K pathway, increase the probability of suicide inactivation. Strong inhibition of the D pathway allows suicide inactivation to occur even in the presence of subunit III. Arachidonic acid, which stimulates proton uptake by the D pathway, retards suicide inactivation. Steady-state turnover in the presence of DeltaPsi and DeltapH, which inhibits proton uptake from the inner surface of the protein, enhances suicide inactivation. Simultaneous inhibition of proton uptake from both sides of the protein by a double mutation affecting the D pathway and the proton backflow/exit pathway greatly shortens the catalytic life span of the oxidase even in the presence of subunit III. Thus, maintenance of rapid proton transfer through the D pathway and the backflow/exit pathway is one mechanism by which subunit III normally functions to prevent suicide inactivation of cytochrome c oxidase. The experiments suggest that increased lifetimes of the heme a(3) oxoferryl intermediates as well as the anionic form of Glu286 of the D pathway cause suicide inactivation in the active site.     

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.     

Proton transport in carbonic anhydrase: Insights from molecular simulation

This article reviews the insights gained from molecular simulations of human carbonic anhydrase II (HCA II) utilizing non-reactive and reactive force fields. The simulations with a reactive force field explore protein transfer and transport via Grotthuss shuttling, while the non-reactive simulations probe the larger conformational dynamics that underpin the various contributions to the rate-limiting proton transfer event. Specific attention is given to the orientational stability of the His64 group and the characteristics of the active site water cluster, in an effort to determine both of their impact on the maximal catalytic rate. The explicit proton transfer and transport events are described by the multistate empirical valence bond (MS-EVB) method, as are alternative pathways for the excess proton charge defect to enter/leave the active site. The simulation results are interpreted in light of experimental results on the wild-type enzyme and various site-specific mutations of HCA II in order to better elucidate the key factors that contribute to its exceptional efficiency.