A fast and simple method to calculate protonation states in proteins.

A simple model for electrostatic interactions in proteins, based on a distance and position dependent screening of the electrostatic potential, is presented. It is applied in conjunction with a Monte Carlo algorithm to calculate pK(alpha) values of ionizable groups in proteins. The purpose is to furnish a simple, fast, and sufficiently accurate model to be incorporated into molecular dynamic simulations. This will allow for dynamic protonation calculations and for coupling between changes in structure and protonation state during the simulation. The best method of calculating protonation states available today is based on solving the linearized Poisson-Boltzmann equation on a finite difference grid. However, this model consumes far too much computer time to be a practical alternative. Tests are reported for fixed structures on bacteriorhodopsin, lysozyme, myoglobin, and calbindin. The studies include comparisons with Poisson-Boltzmann calculations with dielectric constants 4 and 20 inside the protein, a model with uniform dielectric constant 80 and distance-dependent dielectric models. The accuracy is comparable to that of Poisson-Boltzmann calculations with dielectric constant 20, and it is considerably better than that with epsilon = 4. The time to calculate the protonation at one pH value is at least 100 times less than that of a Poisson-Boltzmann calculation. Proteins 1999;36:474-483.     

Histidine residues regulate the transition of photoexcited rhodopsin to its active conformation, metarhodopsin II.

The biologically active photoproduct of rhodopsin, metarhodopsin II (M II), exists in a pH-sensitive equilibrium with its precursor, metarhodopsin I (M I). Increasing acidity favors M II, with the midpoint of the pH titration curve at pH 6.4. To test the long-standing proposal that histidine protonation regulates this conformational transition, we characterized mutant rhodopsins in which each of the 6 histidines was replaced by phenylalanine or cysteine. Only mutants substituted at the 3 conserved histidines showed abnormal M I-M II equilibria. Those in which His-211 was replaced by phenylalanine or cysteine formed little or no M II at either extreme of pH, whereas mutants substituted at His-65 or at His-152 showed enhanced sensitivity to protons. The simplest interpretation of these results is that His-211 is the site where protonation strongly stabilizes the M II conformation and that His-65 and His-152 are sites where protonation modestly destabilizes the M II conformation.     

Structural properties of human carbonic anhydrase II at pH 9.5.

The structure of human carbonic anhydrase II at pH 9.5 has been studied by X-ray crystallographic methods to 2.2 A resolution. These studies complement those performed under acidic conditions in which the catalytically-important proton-shuttle group, His-64, exhibits conformational mobility about side-chain torsion angle chi 1. However, no structural changes are observed in the conformation of His-64 at high pH. Therefore, we conclude that the protonation of His-64 (as well as zinc-bound hydroxide) may be a factor which contributes to the predominantly "out" conformation for His-64 observed at low pH.     

Crystal structures of myoglobin-ligand complexes at near-atomic resolution.

We have used x-ray crystallography to determine the structures of sperm whale myoglobin (Mb) in four different ligation states (unligated, ferric aquomet, oxygenated, and carbonmonoxygenated) to a resolution of better than 1.2 A. Data collection and analysis were performed in as much the same way as possible to reduce model bias in differences between structures. The structural differences among the ligation states are much smaller than previously estimated, with differences of <0.25 A root-mean-square deviation among all atoms. One structural parameter previously thought to vary among the ligation states, the proximal histidine (His-93) azimuthal angle, is nearly identical in all the ferrous complexes, although the tilt of the proximal histidine is different in the unligated form. There are significant differences, however, in the heme geometry, in the position of the heme in the pocket, and in the distal histidine (His-64) conformations. In the CO complex the majority conformation of ligand is at an angle of 18 +/- 3 degrees with respect to the heme plane, with a geometry similar to that seen in encumbered model compounds; this angle is significantly smaller than reported previously by crystallographic studies on monoclinic Mb crystals, but still significantly larger than observed by photoselection. The distal histidine in unligated Mb and in the dioxygenated complex is best described as having two conformations. Two similar conformations are observed in MbCO, in addition to another conformation that has been seen previously in low-pH structures where His-64 is doubly protonated. We suggest that these conformations of the distal histidine correspond to the different conformational substates of MbCO and MbO(2) seen in vibrational spectra. Full-matrix refinement provides uncertainty estimates of important structural parameters. Anisotropic refinement yields information about correlated disorder of atoms; we find that the proximal (F) helix and heme move approximately as rigid bodies, but that the distal (E) helix does not.     

Conformational differences between various myoglobin ligated states as monitored by 1H NMR spectroscopy

In paramagnetic metmyoglobin, cyanomyoglobin (CNMb), and deoxymyoglobin, His-36 has a high pK (approximately 8), and the NMR titration behavior of the H-2 resonance is perturbed, due to the presence at low pH of a hydrogen bond with Glu-38, which is broken at high pH. The His-36 H-4 resonance shows no shift with pK approximately 8 because of two opposing chemical shift effects but monitors the titration of nearby Glu-36 (pK = 5.6). In diamagnetic derivatives [(carbon monoxy)myoglobin (COMb) and oxymyoglobin (oxyMb)], the titration behavior of His-36 H-2 and H-4 resonances is normalized (pK approximately 6.8). The very slight alkaline Bohr effect in sperm whale myoglobin (Mb) is interpreted in terms of the pK change of His-36 from deoxyMb to oxyMb and compensating pK changes in the opposite direction of other unspecified groups. In sperm whale COMb at 40 degrees C, the distal histidine (His-64) and His-97 have pK values of 5.0 and 5.9. The meso proton resonances remote from these groups do not show a titration shift, but the nearby gamma-meso proton (pK = 5.3) responds to titration of both histidines, and the upfield Val-68 methyl at -2.3 ppm (pK = 4.7) witnesses the titration of nearby His-64. At 20 degrees C, the latter resonance is reduced in size, and a second resonance occurs at -2.8 ppm, which is insensitive to pH and, hence, more remote from His-64. Both resonances arise from two conformations of Val-68 in slow equilibrium. In oxyMb at 20 degrees C, only the latter resonance is observed, presumably because of the steric restrictions imposed by the hydrogen bond between ligand and His-64 in oxyMb, which is not present in COMb. In oxyMb the pK of His-97 (5.6) is similar to that of the meso proton resonances (5.5) and to the pK of other pH-dependent processes, including the very small acid Bohr effect. It is likely that these processes are controlled by the titration of His-97.(ABSTRACT TRUNCATED AT 250 WORDS)     

Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation

For the structure and function of proteins, the pH of the solution is one of the determining parameters. Current molecular dynamics (MD) simulations account for the solution pH only in a limited way by keeping each titratable site in a chosen protonation state. We present an algorithm that generates trajectories at a Boltzmann distributed ensemble of protonation states by a combination of MD and Monte Carlo (MC) simulation. The algorithm is useful for pH-dependent structural studies and to investigate in detail the titration behavior of proteins. The method is tested on the acidic residues of the protein hen egg white lysozyme. It is shown that small structural changes may have a big effect on the pK(A) values of titratable residues.     

Assessing the acid-base and conformational properties of histidine residues in human prion protein (125-228) by means of pK(a) calculations and molecular dynamics simulations

A thorough study of the acid-base behavior of the four histidines and the other titratable residues of the structured domain of human prion protein (125-228) is presented. By using multi-tautomer electrostatic calculations, average titration curves have been built for all titratable residues, using the whole bundles of NMR structures determined at pH 4.5 and 7.0. According to our results, (1) only histidine residues are likely to be involved in the first steps of the pH-driven conformational transition of prion protein; (2) the pK(a)'s of His140 and His177 are approximately 7.0, whereas those of His155 and His187 are < 5.5. 10-ns long molecular dynamics simulations have been performed on five different models, corresponding to the most significant combinations of histidine protonation states. A critical comparison between the available NMR structures and our computational results (1) confirms that His155 and His187 are the residues whose protonation is involved in the conformational rearrangement of huPrP in mildly acidic condition, and (2) shows how their protonation leads to the destructuration of the C-terminal part of HB and to the loss of the last turn of HA that represent the crucial microscopic steps of the rearrangement.     

A comparative analysis of the folding and misfolding pathways of the third PDZ domain of PSD95 investigated under different pH conditions

Equilibrium unfolding at neutral pH of the third PDZ domain of PSD95 is well described by the presence of a partly unfolded intermediate that presents association phenomena. After some days' incubation annular and fibrillar structures form from the oligomers. At pH values below 3, however, differential scanning calorimetry shows that PDZ3 seems to unfold under a two-state scheme. Kinetic measurements followed by dynamic light scattering, ThT and ANS fluorescence reveal that the misfolding pathway still exists despite the absence of any populated intermediates and shows an irreversible assembling of the supramacromolecular structures as well as an appreciable lag-phase, contrary to what is found in similar experiments at neutral pH. Moreover, as shown by transmission-electron-microscopy images, the annular structures seen at neutral pH completely disappear from incubated solutions. According to the structural information, this titration behavior appears to be the consequence of a conformational equilibrium that depends on the protonation of some Glu residues located at the C-terminal α3 helix and at the hairpin formed by strands β2 and β3. Our calculations suggest that the enthalpic contribution of these interactions may well be as much as 40kJ·mol(-1). The possible regulatory role of this equilibrium upon PDZ3 functionality and amyloid formation is briefly discussed.     

Structural dynamics of myoglobin: FTIR-TDS study of NO migration and binding

By using Fourier transform infrared photolysis difference spectroscopy combined with temperature derivative spectroscopy at cryogenic temperatures, we have measured infrared spectra of the stretching absorption on nitric oxide (NO) in the heme-bound and photodissociated states of ferrous and ferric nitrosyl myoglobin (MbNO) and a few site-specific Mb mutants. The NO absorption was utilized as a sensitive local probe of ligand interactions with active-site residues and movements within the protein. By comparison with results obtained in previous spectroscopic and structural studies of carbonmonoxy myoglobin (MbCO), the MbNO data were interpreted in structural terms. In the NO-bound state, conformational heterogeneity was inferred from the appearance of multiple bands, arising from different electrostatic interactions with active site residues, most importantly, His-64. In ferrous MbNO, a primary photoproduct site similar to site B of MbCO was found, as indicated by a characteristic NO stretching spectrum. In ferric MbNO, the His-64 side chain appears to interfere with trapping of NO in this site; only a very weak photoproduct spectrum was observed in Mb variants in which His-64 was present. Upon extended illumination, the photoproduct spectrum changed in a characteristic way, indicating that NO readily migrates to a secondary docking site C, the Xe4 cavity, in which the ligand performs librational motions on the picosecond time scale. This docking site may play a role in the physiological NO scavenging reaction. Surprisingly, NO cannot be trapped at all in secondary docking site D, the Xe1 cavity.     

Calculating pH-dependent free energy of proteins by using Monte Carlo protonation probabilities of ionizable residues

Protein folding, stability, and function are usually influenced by pH. And free energy plays a fundamental role in analysis of such pH-dependent properties. Electrostatics-based theoretical framework using dielectric solvent continuum model and solving Poisson-Boltzmann equation numerically has been shown to be very successful in understanding the pH-dependent properties. However, in this approach the exact computation of pH-dependent free energy becomes impractical for proteins possessing more than several tens of ionizable sites (e.g.>30), because exact evaluation of the partition function requires a summation over a vast number of possible protonation microstates. Here we present a method which computes the free energy using the average energy and the protonation probabilities of ionizable sites obtained by the well-established Monte Carlo sampling procedure. The key feature is to calculate the entropy by using the protonation probabilities. We used this method to examine a well-studied protein (lysozyme) and produced results which agree very well with the exact calculations. Applications to the optimum pH of maximal stability of proteins and protein–DNA interactions have also resulted in good agreement with experimental data. These examples recommend our method for application to the elucidation of the pH-dependent properties of proteins.     

Resonance raman investigations of site-directed mutants of myoglobin: effects of distal histidine replacement

The resonance Raman spectra of met-, deoxy-, and (carbonmonoxy)myoglobin (MbCO) are studied as a function of amino acid replacement at the distal histidine-E7 position. The synthetic wild type is found to be spectroscopically identical with the native material. The methionine and glycine replacements do not affect the met or deoxy spectra but do lead to distinct changes in the nu Fe-CO region of the MbCO spectrum. The native MbCO displays a pH-dependent population redistribution of the nu Fe-CO modes, while the analogous population in the mutant systems is found to be pH independent. This indicates that histidine-E7 is the titratable group in native MbCO. Moreover, the pH dependence of the population dynamics is found to be inconsistent with a simple two-state Henderson-Hasselbalch analysis. Instead, we suggest a four-state model involving the coupling of histidine protonation and conformational change. Within this model, the pK of the distal histidine is found to be 6.0 in the "open" configuration and 3.8 in the "closed" conformation. This corresponds to a 3 kcal/mol destabilization of the positively charged distal histidine within the hydrophobic pocket and suggests how protonation can lead to a larger population of the "open" conformation. At pH 7, the pocket is found to be "open" approximately 3% of the time. Further work, involving both IR and Raman measurements, allows the electron-nuclear coupling strengths of the various nu Fe-CO and nu C-O Raman modes to be determined. The slowly rebinding conformational state, corresponding to nu Fe-CO = 518 cm-1 (nu C-O = 1932 cm-1), displays unusually weak coupling of the Fe-CO mode to the Soret transition. Studies of the nu Fe-CO region as a function of temperature reveal that the equilibria between the conformational states are quenched in both the native and glycine mutant below the freezing point of the solvent. Unusual line narrowing of the nu Fe-CO modes at the phase transition is also observed in all samples studied. This line narrowing stands in marked contrast to the other heme Raman modes and suggests that Fe-CO librational motion and/or distal pocket vibrational (or conformational) excitations are involved in the line broadening at room temperature.     

Distinct histidine residues control the acid-induced activation and inhibition of the cloned K(ATP) channel

The modulation of K(ATP) channels during acidosis has an impact on vascular tone, myocardial rhythmicity, insulin secretion, and neuronal excitability. Our previous studies have shown that the cloned Kir6.2 is activated with mild acidification but inhibited with high acidity. The activation relies on His-175, whereas the molecular basis for the inhibition remains unclear. To elucidate whether the His-175 is indeed the protonation site and what other structures are responsible for the pH-induced inhibition, we performed these studies. Our data showed that the His-175 is the only proton sensor whose protonation is required for the channel activation by acidic pH. In contrast, the channel inhibition at extremely low pH depended on several other histidine residues including His-186, His-193, and His-216. Thus, proton has both stimulatory and inhibitory effects on the Kir6.2 channels, which attribute to two sets of histidine residues in the C terminus.     

Chromophore-protein interactions in the anthozoan green fluorescent protein asFP499

Despite their similar fold topologies, anthozoan fluorescent proteins (FPs) can exhibit widely different optical properties, arising either from chemical modification of the chromophore itself or from specific interactions of the chromophore with the surrounding protein moiety. Here we present a structural and spectroscopic investigation of the green FP asFP499 from the sea anemone Anemonia sulcata var. rufescens to explore the effects of the protein environment on the chromophore. The optical absorption and fluorescence spectra reveal two discrete species populated in significant proportions over a wide pH range. Moreover, multiple protonation reactions are evident from the observed pH-dependent spectral changes. The x-ray structure of asFP499, determined by molecular replacement at a resolution of 1.85 A, shows the typical beta-barrel fold of the green FP from Aequorea victoria (avGFP). In its center, the chromophore, formed from the tripeptide Gln(63)-Tyr(64)-Gly(65), is tightly held by multiple hydrogen bonds in a polar cage that is structurally quite dissimilar to that of avGFP. The x-ray structure provides interesting clues as to how the spectroscopic properties are fine tuned by the chromophore environment.     

1H nuclear magnetic resonance titration curves and microenvironments of aromatic residues in bovine pancreatic ribonuclease A

The aromatic region of the NMR spectrum of bovine pancreatic ribonuclease A was analyzed in order to clarify the nature of the microenvironments surrounding the individual histidine, tyrosine, and phenylalanine residues and the interactions with inhibitors. The NMR titration curves of ring protons of six tyrosine and three phenylalanine residues as well as four histidine residues were determined at 37 degrees C between pH 1.5 and pH 11.5 under various conditions. The titration curves were analyzed on the basis of a scheme of a simple proton dissociation sequence and the most probable values were obtained for the macroscopic pK values and intrinsic chemical shifts. The microenvironments surrounding the residues and the effects of inhibitors are discussed on the basis of these results. Based on the titration curves of ring protons, the six tyrosine residues were classified into the following four groups: (1) titratable and different chemical shifts for C(delta) and C(epsilon) protons (two tyrosine residues), (2) titratable but similar chemical shifts for C(delta) and C(epsilon) protons (two tyrosine residues), (3) not titratable and different chemical shifts for C(delta) and C(epsilon) protons (one tyrosine residues), and (4) not titratable and similar chemical shifts for C(delta) and C(epsilon) protons (one tyrosine residue). The resonance signals of ring protons were tentatively assigned to tyrosine and phenylalanine residues. The NMR titration curves of His-48 ring protons were continuous in solution containing 0.2 M sodium acetate but were discontinuous in solution containing 0.3 M NaCl because the NMR signals disappeared at pH values between 5 and 6.5. The effects of addition of formate, acetate, propionate, and ethanol were investigated in order to elucidate the mechanism of the continuity of the titration curves of His-48 in the presence of acetate ion. The NMR signal of His-48 C(2) protons was observed at pH 6 in the presence of acetate and propionate ions but was not observed in the presence of formate ion or ethanol. This indicated that both the alkyl chain and the anionic carboxylate group are necessary for the continuity of the titration curves of His-48 ring protons. Based on the results, the mechanism of the effects of acetate ion is discussed.     

A fast and accurate computational approach to protein ionization

We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.     

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF

Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.     

Protonation changes upon ligand binding to trypsin and thrombin: structural interpretation based on pK(a) calculations and ITC experiments

The protonation states of a protein and a ligand can be altered upon complex formation. Such changes can be detected experimentally by isothermal titration calorimetry (ITC). For a series of ligands binding to the serine proteases trypsin and thrombin, we previously performed an extensive ITC and crystallographic study and were able to identify protonation changes for four complexes. However, since ITC measures only the overall proton exchange, it does not provide structural insights into the functional groups involved in the proton transfer. Using Poisson-Boltzmann calculations based on our recently developed PEOE_PB charges, we compute pK(a) values for all complexes of our former study in order to reveal the residues with altered protonation states. The results indicate that His57, a member of the catalytic triad, is responsible for the most relevant pK(a) shifts leading to the experimentally detected protonation changes. This finding is in contrast to our previous assumption that the observed protonation changes occur at the carboxylic group of the ligands. The newly detected proton acceptor is used for a revised factorization of the ITC data, which is necessary whenever the protonation inventory changes upon complexation. The pK(a) values of complexes showing no protonation change in the ITC experiment are reliably predicted in most cases, whereas predictions of strongly coupled systems remain problematic.     

Molecular basis for redox-Bohr and cooperative effects in cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies of oxidized and reduced high-resolution structures at pH 7.6

The tetraheme cytochrome c3 is a small metalloprotein with ca. 13,000 Da found in sulfate-reducing bacteria, which is believed to act as a partner of hydrogenase. The three-dimensional structure of the oxidized and reduced forms of cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774 at pH 7.6 were determined using high-resolution X-ray crystallography and were compared with the previously determined oxidized form at pH 4.0. Theoretical calculations were performed with both structures, using continuum electrostatic calculations and Monte Carlo sampling of protonation and redox states, in order to understand the molecular basis of the redox-Bohr and cooperativity effects related to the coupled transfer of electrons and protons. We were able to identify groups that showed redox-linked conformational changes. In particular, Glu61, His76, and propionate D of heme II showed important contributions to the redox-cooperativity, whereas His76, propionate A of heme I, and propionate D of heme IV were the key residues for the redox-Bohr effect. Upon reduction, an important movement of the backbone region surrounding hemes I and II was also identified, that, together with a few redox-linked conformational changes in side-chain residues, results in a significant decrease in the solvent accessibility of hemes I and II.     

The role of His-18 in amyloid formation by human islet amyloid polypeptide

The 37-residue islet amyloid polypeptide (IAPP) is the major protein component of the amyloid deposits found in type-II diabetes. IAPP is stored in a relatively low pH environment in the pancreatic secretory granules prior to its release to the extracellular environment. Human IAPP contains a single histidine at position 18. Aggregation of IAPP is considerably faster at a lower pH (4.0 +/- 0.3) than at high pH (8.8 +/- 0.3), as judged by turbidity and thioflavine-T fluorescence studies. The rate of aggregation at low pH increases drastically in the presence of salt. CD experiments show that the conversion of largely unstructured monomers to beta-sheet-rich structures is faster at high pH. TEM studies show that fibrils are formed at both pH values but are more prevalent at pH 8.8 (+/-0.3). Both the free N terminus of IAPP and His-18 will titrate over the pH range studied. An N-terminal acetylated fragment consisting of residues 8-37 of human IAPP was also studied to isolate contributions from the protonation of His-18. Previous studies have shown that this fragment forms protofibrils that are very similar to those formed by intact IAPP. The effects of varying the protonation state of His-18 in the 8-37 analogue indicate that the rate of aggregation and fibril formation is noticeably faster when His-18 is deprotonated, similar to the wild type. However, the pH-dependent effects are larger for full-length IAPP than for the disulfide-truncated, acetylated analogue. TEM studies indicate differences in the morphology of the deposits formed at high and low pH. These results are discussed in light of recent structural models of IAPP fibrils.     

pH modulates the quinone position in the photosynthetic reaction center from Rhodobacter sphaeroides in the neutral and charge separated states

The structure of the photosynthetic reaction-center from Rhodobacter sphaeroides has been determined at four different pH values (6.5, 8.0, 9.0, 10.0) in the neutral and in charge separated states. At pH 8.0, in the neutral state, we obtain a resolution of 1.87 A, which is the best ever reported for the bacterial reaction center protein. Our crystallographic data confirm the existence of two different binding positions of the secondary quinone (QB). We observe a new orientation of QB in its distal position, which shows no ring-flip compared to the orientation in the proximal position. Datasets collected for the different pH values show a pH-dependence of the population of the proximal position. The new orientation of QB in the distal position and the pH-dependence could be confirmed by continuum electrostatics calculations. Our calculations are in agreement with the experimentally observed proton uptake upon charge separation. The high resolution of our crystallographic data allows us to identify new water molecules and external residues being involved in two previously described hydrogen bond proton channels. These extended proton-transfer pathways, ending at either of the two oxo-groups of QB in its proximal position, provide additional evidence that ring-flipping is not required for complete protonation of QB upon reduction.