Differential effects of the Zn-His-Bkb vs Zn-His-[Asp/Glu] triad on Zn-core stability and reactivity

The most common partner of the Zn-bound His is the Asp/Glu carboxylate side chain in catalytic Zn sites and the backbone (Bkb) carbonyl group in structural Zn sites. To elucidate the factors governing the selection of the second-shell partner of the Zn-bound His in structural/catalytic Zn sites, systematic studies using density functional theory and continuum dielectric calculations were performed to determine the relative contributions of the second-shell Bkb carbonyl and the Asp/Glu carboxylate to the Zn-core stability and reactivity. The results show that the contributions of the second-shell Bkb carbonyl and Asp/Glu carboxylate to the Zn-core stability depend mainly on the solvent accessibility of the Zn-site and the composition of the Zn-core. They reveal the advantage of a second-shell Bkb carbonyl in anionic Zn cavities: it stabilizes anionic, buried Zn-cores more than the corresponding negatively charged Asp/Glu carboxylate, thus explaining the absence of the Zn-His-Asp/Glu triad in structural [Zn(Cys)3(His)]- cores. They also reveal the advantage of a second-shell Asp/Glu carboxylate in catalytic Zn-cores: relative to a Bkb carbonyl group, it increases (i) the HOMO energy of the cationic/neutral zinc core, (ii) the reactivity of the attacking Zn-bound OH-, (iii) electron transfer to the substrate, and (iv) the stability of the metal complex upon electron transfer. Furthermore, a second-shell Asp/Glu carboxylate could facilitate product release in the common cationic catalytic cores, by acting as a proton acceptor of the Zn-bound His creating an Asp...His- dyad that stabilizes the zinc dication more than the respective Bkb...His0 dyad.     

First-second shell interactions in metal binding sites in proteins: a PDB survey and DFT/CDM calculations

The role of the second shell in the process of metal binding and selectivity in metalloproteins has been elucidated by combining Protein Data Bank (PDB) surveys of Mg, Mn, Ca, and Zn binding sites with density functional theory/continuum dielectric methods (DFT/CDM). Peptide backbone groups were found to be the most common second-shell ligand in Mg, Mn, Ca, and Zn binding sites, followed (in decreasing order) by Asp/Glu, Lys/Arg, Asn/Gln, and Ser/Thr side chains. Aromatic oxygen- or nitrogen-containing side chains (Tyr, His, and Trp) and sulfur-containing side chains (Cys and Met) are seldom found in the second coordination layer. The backbone and Asn/Gln side chain are ubiquitous in the metal second coordination layer as their carbonyl oxygen and amide hydrogen can act as a hydrogen-bond acceptor and donor, respectively, and can therefore partner practically every first-shell ligand. The second most common outer-shell ligand, Asp/Glu, predominantly hydrogen bonds to a metal-bound water or Zn-bound histidine and polarizes the H-O or H-N bond. In certain cases, a second-shell Asp/Glu could affect the protonation state of the metal ligand. It could also energetically stabilize a positively charged metal complex more than a neutral ligand such as the backbone and Asn/Gln side chain. As for the first shell, the second shell is predicted to contribute to the metal selectivity of the binding site by discriminating between metal cations of different ionic radii and coordination geometries. The first-shell-second-shell interaction energies decay rapidly with increasing solvent exposure of the metal binding site. They are less favorable but are of the same order of magnitude as compared to the respective metal-first-shell interaction energies. Altogether, the results indicate that the structure and properties of the second shell are dictated by those of the first layer. The outer shell is apparently designed to stabilize/protect the inner-shell and complement/enhance its properties.     

Catalytic mechanism of guanine deaminase: an ONIOM and molecular dynamics study

The catalytic mechanism of Bacillus subtilis guanine deaminase (bGD), a Zn metalloenzyme, has been investigated by a combination of quantum mechanical calculations using the multilayered ONIOM method and molecular dynamics simulations. In contrast to a previously proposed catalytic mechanism, which requires the bound guanine to assume a rare tautomeric state, the ONIOM calculations showed that the active-site residues of the enzyme do not affect the tautomeric state of guanine, and consequently the bound guanine is a tautomer that is the most abundant in aqueous solution. Two residues, Glutamate 55 and Aspartate 114, were found to play important roles in proton shuttling in the reaction. The proposed reaction path is initiated by proton transfer from a Zn-bound water to protonate Asp114. This process may be quite complex and rather dynamic in nature, as revealed by the molecular dynamics (MD) simulations, whereby another water may bridge the Zn-bound water and Asp114, which then is eliminated by positioning of guanine in the active site. The binding of guanine stabilizes protonated Asp114 by hydrogen bond formation. Asp114 can then transfer its proton to the N3 of the bound guanine, facilitating the nucleophilic attack on C2 of the guanine by the Zn-bound hydroxide to form a tetrahedral intermediate. This occurs with a rather low barrier. Glu55 then transfers a proton from the Zn-hydroxide to the amino group of the reaction intermediate and, at this point, the C2-N2 bond has lengthened by 0.2 A compared to guanine, making C2-N2 bond cleavage more facile. The C2-N2 bond breaks forming ammonia, with an energy barrier of approximately 8.8 kcal/mol. Ammonia leaves the active site, and xanthine is freed by the cleavage of the Zn-O2 bond, with a barrier approximately 8.4 kcal/mol. Along this reaction path, the highest barrier comes from C2-N2 bond cleavage, while the barrier from the cleavage of the Zn-O2 bond is slightly smaller. The Zn-O2 bond can be broken without the assistance of water during the release of xanthine.     

Proton exit channels in bovine cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal transmembrane enzyme of the respiratory electron transport chain in aerobic cells. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping, a process which 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 and proton exit pathways remain unknown. In our previous work, the continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine cytochrome c oxidase for different redox states of the enzyme. A possible mechanism of oxygen reduction and proton pumping via His291 was proposed. In this paper, using electrostatic calculations, we examine the proton exit pathways in the enzyme. By monitoring the changes of the protonation states, proton affinities, and energies of electrostatic interactions between the titratable groups in different redox states of CcO, we identified the clusters of strongly interacting residues. Using these data, we detected four possible proton exit points on the periplasmic side of the membrane (Lys171B/Asp173B, His24B/Asp25B, Asp51, and Asp300). We then were able to trace the proton exit pathways and to evaluate the energy profiles along the paths. On the basis of energetic considerations and the conservation of the residues in a protein sequence, the most likely exit pathway is one via the Lys171B/Asp173B site. The obtained results are fully consistent with our His291 model of proton pumping, and provide a rationale for the absence of proton leaking in CcO between the pumping strokes.     

Theoretical study toward understanding the catalytic mechanism of pyruvate decarboxylase

Density functional calculations are employed to explore the mechanisms of all elementary reaction steps involved in the catalytic cycle of pyruvate decarboxylase (PDC). Different models are constructed for mimicking the involvement of some key residues in a certain step. The effect of the protein framework on the potential energy profiles of active site models is approximately modeled by fixing some freedoms, based on the crystal structure of the PDC enzyme from Saccharomyces cerevisiae (ScPDC). Our calculations confirm that Glu51 is the most important residue in the formation of the ylide and the release of acetaldehyde via the proton relay between Glu51, N1', and the 4'-amino group of thiamine diphosphate. The presence of Glu477 and Asp28 residues makes the decarboxylation of lactylthiamin diphosphate (LThDP) an endothermic process with a significant free energy barrier. The protonation of the alpha-carbanion to form 2-(1-hydroxyethyl)-thiamin diphosphate is found to go through a concerted double proton transfer transition state involving both Asp28 and His115 residues. The final step, acetaldehyde release, is likely to proceed through a concerted transition state involving carbon-carbon bond-breaking and the deprotonation of the alpha-hydroxyl group. The decarboxylation of LThDP and the protonation of the alpha-carbanion are two rate-limiting steps, relative to the facile occurrence of the ylide formation and acetaldehyde release. The catalytic roles of residues Glu51, Glu477, Asp28, and Gly417 in the active site of ScPDC in individual steps elucidated from the present study are in good agreement with those derived from site-directed mutagenesis.     

Competition between Li+ and Mg2+ in metalloproteins. Implications for lithium therapy

Lithium is used (in the form of soluble salts) to treat bipolar disorder and has been considered as a possible drug in treating chronic neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases. One of the proposed mechanisms of Li(+) action involves a competition between the alien Li(+) and native Mg(2+) for metal-binding sites and subsequent inhibition of key enzymes involved in specific neurotransmission pathways, but not vital Mg(2+) proteins in the cell. This raises the following intriguing questions: Why does Li(+) replace Mg(2+) only in enzymes involved in bipolar disorder, but not in Mg(2+) proteins essential to cells? In general, what factors allow monovalent Li(+) to displace divalent Mg(2+) in proteins? Specifically, how do the composition, overall charge, and solvent exposure of the metal-binding site as well as a metal-bound phosphate affect the selectivity of Li(+) over Mg(2+)? Among the many possible factors, we show that the competition between Mg(2+) and Li(+) depends on the net charge of the metal complex, which is determined by the numbers of metal cations and negatively charged ligands, as well as the relative solvent exposure of the metal cavity. The protein itself is found to select Mg(2+) over the monovalent Li(+) by providing a solvent-inaccessible Mg(2+)-binding site lined by negatively charged Asp/Glu, whereas the cell machinery was found to select Mg(2+) among other competing divalent cations in the cellular fluids such as Ca(2+) and Zn(2+) by maintaining a high concentration ratio of Mg(2+) to its biogenic competitor in various biological compartments. The calculations reveal why Li(+) replaces Mg(2+) only in enzymes that are known targets of Li(+) therapy, but not in Mg(2+) enzymes essential to cells, and also reveal features common to the former that differ from those in the latter proteins.     

A model study of the efficiency of the Asp–His–Ser triad

The observation of the Asp–His–Ser triad (Asp: aspartate, His: histidine, Ser: serine) triad both in mammalian and bacterial proteases suggests a special efficiency. A series of B3LYP/D95*(d,p) calculations on various [X–H^βY]^? dyads (as part of the [X–H^βY–H^αAc]^? model triad, HAc: acetic acid) made from eight different anions X^? and 15 different coupling elements H^βY was done to analyze the molecular origin of this efficiency. The X^? anion acts merely as an electron density donor independent of its chemical nature, and the evolutionary selection of Asp for the catalytic triad therefore seems to be caused by the pH of the triads environment. As the linking proton H^β moves from Y^? to X^?, electron density is effectively moved from X^? to Y^? thereby increasing the proton affinity (PA) of the [X–HY]^? dyad, which finally leads to the deprotonization of the HAc molecule. The degree to which the position of H^α controls the PA is dominatly determined by the coupling element HY. The model calculations indicate that 4‐methyl‐1H‐imidazole (HMim) is a very efficient coupling element, which suggest that the evolutionary convergence to the Asp–His–Ser is not only controlled by the ready availability of the imidazole motive in His but also by its high efficiency. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

A molecular dynamics exploration of the catalytic mechanism of yeast cytosine deaminase

Yeast cytosine deaminase (yCD), a zinc metalloenzyme of significant biomedical interest, is investigated by a series of molecular dynamics simulations in its free form and complexed with its reactant (cytosine), product (uracil), several reaction intermediates, and an intermediate analogue. Quantum chemical calculations, used to construct a model for the catalytic Zn ion with its ligands (two cysteines, a histidine, and one water) show, by comparison with crystal structure data, that the cysteines are deprotonated and the histidine is monoprotonated. The simulations suggest that Glu64 plays a critical role in the catalysis by yCD. The rotation of the Glu64 side-chain carboxyl group that can be protonated or deprotonated permits it to act as a proton shuttle between the Zn-bound water and cytosine and subsequent reaction intermediates. Free energy methods are used to obtain the barriers for these rotations, and they are sufficiently small to permit rotation on a nanosecond time scale. In the course of the reaction, cytosine reorients to a geometry to favor nucleophilic attack by a Zn-bound hydroxide. A stable position for a reaction product, ammonia, was located in the active site, and the free energy of exchange with a water molecule was evaluated. The simulations also reveal small motions of the C-terminus and the loop that contains Phe114 that may be important for reactant binding and product release.     

A QM/MM study on the catalytic mechanism of pyruvate decarboxylase

Pyruvate decarboxylase (PDC) is a typical thiamin diphosphate (ThDP)-dependent enzyme with widespread applications in industry. Though studies regarding the reaction mechanism of PDC have been reported, they are mainly focused on the formation of ThDP ylide and some elementary steps in the catalytic cycle, studies about the whole catalytic cycle of PDC are still not completed. In these previous studies, a major controversy is whether the key active residues (Glu473, Glu50′, Asp27′, His113′, His114′) are protonated or ionized during the reaction. To explore the catalytic mechanism and the role of key residues in the active site, three whole-enzyme models were considered, and the combined QM/MM calculations on the nonoxidative decarboxylation of pyruvate to acetaldehyde catalyzed by PDC were performed. According to our computational results, the fundamental reaction pathways, the complete energy profiles of the whole catalytic cycle, and the specific role of key residues in the common steps were obtained. It is also found that the same residue with different protonation states will lead to different reaction pathways and energy profiles. The mechanism derived from the model in which the residues (Glu473, Glu50′, Asp27′, His113′, His114′) are in their protonated states is most consistent with experimental observations. Therefore, extreme care must be taken when assigning the protonation states in the mechanism study. Because the experimental determination of protonation state is currently difficult, the combined QM/MM method provides an indirect means for determining the active-site protonation state.     

Substituent Effect on Proton Affinity of Imidazole in Cu,Zn-Superoxide Dismutase

To investigate whether the proton-accepting ability of imidazole in Cu,Zn-superoxide dismutase （SOD） was possibly modulated by Zn（Ⅱ） or not, the proton affinity （Ap） of N^3 in imidazole group was calculated by density functional theory （DFT） with B3LYP functional. It was found that Zn（Ⅱ） attenuates the Ap, because of its electron-withdrawing effect, while the three ligands connected with Zn（Ⅱ） （residues of two His and one Asp） exert an opposite effect, owing to their electron-donating ability. This finding suggested that the three ligands should play a role in the normal function of Cu,Zn-SOD and should be taken into consideration in the future study.     

Computational study on mechanistic details of the aminoethanol rearrangement catalyzed by the vitamin B12-dependent ethanolamine ammonia lyase: His and Asp/Glu acting simultaneously as catalytic auxiliaries

The rearrangement of aminoethanol catalyzed by ethanolamine ammonia lyase is investigated by computational means employing DFT (B3LYP/6-31G) and ab initio molecular orbital theory (QCISD/cc-pVDZ). The study aims at providing a detailed account on various crucial aspects, in particular a distinction between a direct intramolecular migration of the partially protonated NH(2) group vs elimination of NH(4)(+). Three mechanistic scenarios were explored: (i) According to the calculations, irrespective of the nature of the protonating species, intramolecular migration of the NH(3) group is energetically less demanding than elimination of NH(4)(+). However, all computed activation enthalpies exceed the experimentally derived activation enthalpy (15 kcal/mol) associated with the rate-determining step, i.e., the hydrogen abstraction from the 5'-deoxyadenosine by the product radical. For example, when imidazole is used as a model system for His interacting with the NH(3) group of the substrate, the activation enthalpy for the migration process amounts to 27.4 kcal/mol. If acetic acid is employed to mimic Asp or Glu, the activation enthalpy is somewhat lower, being equal to 24.2 kcal/mol. (ii) For a partial deprotonation of the substrate 2 at the OH group, the rearrangement mechanism consists of the dissociation of an NH(2) radical from C(2) and its association at C(1) atom. For all investigated proton acceptors (i.e., OH(-), HCOO(-), CH(3)COO(-), CH(2)NH, imidazole), the activation enthalpy for the dissociation step also exceeds 15 kcal/mol. Typical data are 20.2 kcal/mol for Ac(-) and 23.8 kcal/mol for imidazole. (iii) However, in a synergistic action of partial protonation of the NH(2) group and partial deprotonation of the OH group by the two conceivable catalytic auxiliaries Asp/Glu and His, the activation enthalpy computed is compatible with the experimental data. For imidazole and acetate as model systems, the activation enthalpy is equal to 13.7 kcal/mol. This synergistic action of the two catalytic groups is expected to take place in a physiologically realistic pH range of 6-9.5, and the present computational findings may help to further characterize the yet unknown structural details of the ethanolamine ammonia lyase's active site.     

Copper-zinc superoxide dismutase: theoretical insights into the catalytic mechanism

The mechanism for the toxic superoxide radical disproportionation to molecular oxygen and hydrogen peroxide by copper-zinc superoxide dismutase (CuZnSOD) has been studied using the B3LYP hybrid density functional. On the basis of the X-ray structure of the enzyme, the molecular system investigated includes the first-shell protein ligands of the two metal centers as well as the second-shell ligand Asp122. The substrates of the model reaction are two superoxide radical anions, approaching the copper center at the beginning of two half-reactions: the first part of the catalytic cycle involving Cu+ oxidation and the second part reducing Cu2+ back to its initial state. The quantitative free energy profile of the reaction is obtained and discussed in connection with the experimental data on the reduction potentials and CuZnSOD kinetics. The optimized structures are analyzed and compared to the experimental ones. The two transition states alternate the protonation state of His61 and correspond to histidine Cu-His61-Zn bridge rupture/reformation. Modifications applied to the initial model allow the importance of Asp122 for catalysis to be estimated.     

IDENTIFICATION OF THE PROTON ACCEPTOR OF SCHIFF BASE DEPROTONATION IN BACTERIORHODOPSIN: A FOURIER-TRANSFORM-INFRARED STUDY OF THE MUTANT ASP85 → GLU IN ITS NATURAL LIPID ENVIRONMENT

Abstract— In order to assign the proton acceptor for Schiff base deprotonation in bacteriorhodopsin to a specific Asp residue, the photoreaction of the Asp85 → Glu mutant, as expressed in Halobacterium sp . GRB, was investigated by static low-temperature and time-resolved infrared difference spec-troscopy. Measurements were also performed on the mutant protein labeled with [4-¹³C]Asp which allowed discrimination between Asp and Glu residues. 14,15-di¹³C-retinal was incorporated to distinguish amide-II absorbance changes from changes of the ethylenic mode of the chromophore. In agreement with earlier UV-VIS measurements, our data show that from both the 540 and 610 nm species present in a pH-dependent equilibrium, intermediates similar to K and L can be formed. The 14 ms time-resolved spectrum of the 540 nm species shows that a glutamic acid becomes protonated in the M-like intermediate, whereas the comparable difference spectrum of the 610 nm species demonstrates that in the initial state a glutamic acid is already protonated. In conjunction with earlier observations of protonation of an Asp residue in wild-type M, the data provide direct evidence that the proton acceptor in the deprotonation reaction of the Schiff base is Asp85.     

A theoretical study on the catalytic mechanism of Mus musculus adenosine deaminase

The catalytic mechanism of Mus musculus adenosine deaminase (ADA) has been studied by quantum mechanics and two‐layered ONIOM calculations. Our calculations show that the previously proposed mechanism, involving His238 as the general base to activate the Zn‐bound water, has a high activation barrier of about 28 kcal/mol at the proposed rate‐determining nucleophilic addition step, and the corresponding calculated kinetic isotope effects are significantly different from the recent experimental observations. We propose a revised mechanism based on calculations, in which Glu217 serves as the general base to abstract the proton of the Zn‐bound water, and the protonated Glu217 then activates the substrate for the subsequent nucleophilic addition. The rate‐determining step is the proton transfer from Zn‐OH to 6‐NH₂ of the tetrahedral intermediate, in which His238 serves as a proton shuttle for the proton transfer. The calculated kinetic isotope effects agree well with the experimental data, and calculated activation energy is also consistent with the experimental reaction rate. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Distinguishing Protonation States of Histidine Ligands to the Oxidized Rieske Iron–Sulfur Cluster through 15N Vibrational Frequency Shifts

The Rieske [2Fe–2S] cluster is a vital component of many oxidoreductases, including mitochondrial cytochrome bc1; its chloroplast equivalent, cytochrome b6f; one class of dioxygenases; and arsenite oxidase. The Rieske cluster acts as an electron shuttle and its reduction is believed to couple with protonation of one of the cluster′s His ligands. In cytochromes bc1 and b6f, for example, the Rieske cluster acts as the first electron acceptor in a modified Q cycle. The protonation states of the cluster′s His ligands determine its ability to accept a proton and possibly an electron through a hydrogen bond to the electron carrier, ubiquinol. Experimental determination of the protonation states of a Rieske cluster′s two His ligands by NMR spectroscopy is difficult, due to the close proximity of the two paramagnetic iron atoms of the cluster. Therefore, this work reports density functional calculations and proposes that difference vibrational spectroscopy with ¹⁵N isotopic substitution may be used to assign the protonation states of the His ligands of the oxidized Rieske [2Fe–2S] complex.     

Molecular orbital study of the effects of ionic amino acid residues on proton transfer energetics in the active site of carboxypeptidase a

As part of a study of the catalytic mechanism of carboxypeptidase A, two proton transfers in the system Glu 270 H₂OZnHis 69-Asp 142 are treated with an ab initio MO method. Results show that the proton transfers are strongly affected by the environment. It is predicted that the proton between His 69 and Asp 142 is covalently bonded to His 69.     

All-electron calculations of the nucleation structures in metal-induced zinc-finger folding: role of the Peptide backbone

Although the folding of individual protein domains has been extensively studied both experimentally and theoretically, protein folding induced by a metal cation has been relatively understudied. Almost all folding mechanisms emphasize the role of the side-chain interactions rather than the peptide backbone in the protein folding process. Herein, we focus on the thermodynamics of the coupled metal binding and protein folding of a classical zinc-finger (ZF) peptide, using all-electron calculations to obtain the structures of possible nucleation centers and free energy calculations to determine their relative stability in aqueous solution. The calculations indicate that a neutral Cys first binds to hexahydrated Zn2+ via its ionized sulfhydryl group and neutral backbone oxygen, with the release of four water molecules and a proton. Another nearby Cys then binds in the same manner as the first one, yielding a fully dehydrated Zn2+. Subsequently, two His ligands from the C-terminal part of the peptide successively dislodge the Zn-bound backbone oxygen atoms to form the native-like Zn-(Cys)2(His)2 complex. Each successive Zn complex accumulates increasingly favorable and native interactions, lowering the energy of the ZF polypeptide, which concomitantly becomes more compact, reducing the search volume, thus guiding folding to the native state. In the protein folding process, not only the side chains but also the backbone peptide groups play a critical role in stabilizing the nucleation structures and promoting the hydrophobic core formation.     

Factors controlling the reactivity of zinc finger cores

Although the Zn(2+) cation in Zn·Cys(4), Zn·Cys(3)His, Zn·Cys(2)His(2), and Zn(2)Cys(6) cores of zinc finger (Zf) proteins typically plays a structural role, the Zn-bound thiolates in some Zf cores are reactive. Such labile Zf cores can serve as drug targets for retroviral or cancer therapies. Previous studies showed that the reactivity of a Zn-bound thiolate toward electrophiles is significantly reduced if it forms S---NH hydrogen bonds with the backbone amide. However, we found several well-known inactive Zf cores containing Cys ligands with no H-bonding interactions. Here, we show that H bonds from the peptide backbone or bonds from a second Zn cation to Zn-bound S atoms suppress the reactivity not only of these S atoms, but also of Zn-bound S* atoms with no interactions. Indeed, two or more indirect NH---S hydrogen bonds raise the free energy barrier for methylation of a Zn-bound S* in a Cys(4) core more than a direct NH---S* hydrogen bond. These findings help to elucidate why several well-known Zf cores have Cys ligands with no H bonds, but are unreactive. They also help to provide guidelines for distinguishing labile Cys-rich Zn sites from structural ones, which in turn help to identify novel potential Zf drug targets.     

DFT/electrostatic calculations of pK(a) values in cytochrome c oxidase

Using classical electrostatic calculations, earlier we examined the dependence of the protonation state of bovine cytochrome c oxidase (CcO) on its redox state. Based on these calculations, we have proposed a model of CcO proton pumping that involves His291, one of the Cu(B) histidine ligands, which was found to respond to redox changes of the enzyme Fe(a)(3)-Cu(B) catalytic center. In this work, we employ combined density functional and continuum electrostatic calculations to evaluate the pK(a)() values of His291 and Glu242, two key residues of the model. The pK(a) values are calculated for different redox states of the enzyme, and the influence of different factors on the pK(a)'s is analyzed in detail. The calculated pK(a)() values of Glu242 are between 9.4 and 12.0, depending on the redox state of the protein, which is in excellent agreement with recent experimental measurements. Assuming the reduced state of heme a(3), His291 of the oxidized Cu(B) center possesses a pK(a)() between 2.1 and 4.0, while His291 of the reduced Cu(B) center has a pK(a) above 17. The obtained results support the proposal that the His291 ligand of the Cu(B) center in CcO is a proton pump element.     

Influence of the membrane potential on the protonation of bacteriorhodopsin: insights from electrostatic calculations into the regulation of proton pumping

Proton binding and release are elementary steps for the transfer of protons within proteins, which is a process that is crucial in biochemical catalysis and biological energy transduction. Local electric fields in proteins affect the proton binding energy compared to aqueous solution. In membrane proteins, also the membrane potential affects the local electrostatics and can thus be crucial for protein function. In this paper, we introduce a procedure to calculate the protonation probability of titratable sites of a membrane protein in the presence of a membrane potential. In the framework of continuum electrostatics, we use a modified Poisson-Boltzmann equation to include the influence of the membrane potential. Our method considers that in a transmembrane protein each titratable site is accessible for protons from only one side of the membrane depending on the hydrogen bond pattern of the protein. We show that the protonation of sites receiving their protons from different sides of the membrane is differently influenced by the membrane potential. In addition, the effect of the membrane potential is combined with the effect of the pH gradient resulting from proton pumping. Our method is applied to bacteriorhodopsin, a light-activated proton pump. We find that the proton pumping of this protein might be regulated by Asp115, a conserved residue for which no function has been identified yet. According to our calculations, the interaction of Asp115 with Asp85 leads to the protonation of the latter if the pH gradient or the membrane potential becomes too large. Since Asp85 is the primary proton acceptor in the photocycle, bacteriorhodopsin molecules in which Asp85 is protonated cannot pump protons. Furthermore, we estimate how the membrane potential affects the energetics of the individual proton-transfer reactions of the photocycle. Most reactions, except the initial proton transfer from the Schiff base to Asp85, are influenced. Our calculations give new insights into the mechanism with which bacteriorhodopsin senses the membrane potential and the pH gradient and how the proton pumping is regulated by these parameters.