Bicarbonate rescues damaged proton-transfer pathway in photosystem II

The membrane-protein complex photosystem II (PSII) catalyzes photosynthetic water oxidation. Proton transfer plays an integral role in the catalytic cycle of water oxidation by maintaining charge balance to regulate and ensure the efficiency of the process. The hydrogen-bonded amino-acid residues that surround the oxygen-evolving complex (OEC) provide an efficient pathway for proton removal. Hence, it is crucial to identify these pathways to provide deeper insights into the proton-transfer mechanisms. In this study, we have used bicarbonate as a mobile exogenous proton-transfer reagent to recover the activity lost by site-directed mutations in order to identify amino-acid residues participating in the proton-transfer pathway. We find that bicarbonate restores efficient S-state cycling in D2-K317A PSII core complexes, but not in D1-D61A and CP43-R357K PSII core complexes, indicating that bicarbonate chemical rescue can be used to differentiate single-point mutations affecting the pathways of proton transfer from mutations that affect other aspects of the water-oxidation mechanism.     

Deuterium isotope effect of proton pumping in cytochrome c oxidase

In mitochondria and many aerobic bacteria cytochrome c oxidase is the terminal enzyme of the respiratory chain where it catalyses the reduction of oxygen to water. The free energy released in this process is used to translocate (pump) protons across the membrane such that each electron transfer to the catalytic site is accompanied by proton pumping. To investigate the mechanism of electron-proton coupling in cytochrome c oxidase we have studied the pH-dependence of the kinetic deuterium isotope effect of specific reaction steps associated with proton transfer in wild-type and structural variants of cytochrome c oxidases in which amino-acid residues in proton-transfer pathways have been modified. In addition, we have solved the structure of one of these mutant enzymes, where a key component of the proton-transfer machinery, Glu286, was modified to an Asp. The results indicate that the P3-->F3 transition rate is determined by a direct proton-transfer event to the catalytic site. In contrast, the rate of the F3-->O4 transition, which involves simultaneous electron transfer to the catalytic site and is characteristic of any transition during CytcO turnover, is determined by two events with similar rates and different kinetic isotope effects. These reaction steps involve transfer of protons, that are pumped, via a segment of the protein including Glu286 and Arg481.     

Suppression of the back proton-transfer from Asp85 to the retinal Schiff base in bacteriorhodopsin: a theoretical analysis of structural elements

The transfer of a proton from the retinal Schiff base to the nearby Asp85 protein group is an essential step in the directional proton-pumping by bacteriorhodopsin. To avoid the wasteful back reprotonation of the Schiff base from Asp85, the protein must ensure that, following Schiff base deprotonation, the energy barrier for back proton-transfer from Asp85 to the Schiff base is larger than that for proton-transfer from the Schiff base to Asp85. Here, three structural elements that may contribute to suppressing the back proton-transfer from Asp85 to the Schiff base are investigated: (i) retinal twisting; (ii) hydrogen-bonding distances in the active site; and (iii) the number and location of internal water molecules. The impact of the pattern of bond twisting on the retinal deprotonation energy is dissected by performing an extensive set of quantum-mechanical calculations. Structural rearrangements in the active site, such as changes of the Thr89:Asp85 distance and relocation of water molecules hydrogen-bonding to the Asp85 acceptor group, may participate in the mechanism which ensures that following the transfer of the Schiff base proton to Asp85 the protein proceeds with the subsequent photocycle steps, and not with back proton transfer from Asp85 to the Schiff base.     

Intramolecular proton-transfer reactions in a membrane-bound proton pump: the effect of pH on the peroxy to ferryl transition in cytochrome c oxidase

In the membrane-bound redox-driven proton pump cytochrome c oxidase, electron- and proton-transfer reactions must be coupled, which requires controlled modulation of the kinetic and/or thermodynamic properties of proton-transfer reactions through the membrane-spanning part of the protein. In this study we have investigated proton-transfer reactions through a pathway that is used for the transfer of both substrate and pumped protons in cytochrome c oxidase from Rhodobacter sphaeroides. Specifically, we focus on the formation of the so-called F intermediate, which is rate limited by an internal proton-transfer reaction from a possible branching point in the pathway, at a glutamic-acid residue (E(I-286)), to the binuclear center. We have also studied the reprotonation of E(I-286) from the bulk solution. Evaluation of the data in terms of a model presented in this work gives a rate of internal proton transfer from E(I-286) to the proton acceptor at the catalytic site of 1.1 x 10(4) s(-1). The apparent pK(a) of the donor (E(I-286)), determined from the pH dependence of the F-formation kinetics, was found to be 9.4, while the pK(a) of the proton acceptor at the catalytic site is likely to be > or = 2.5 pH units higher. In the pH range up to pH 10 the proton equilibrium between the bulk solution and E(I-286) was much faster than 10(4) s(-1), while in the pH range above pH 10 the proton uptake from solution is rate limiting for the overall reaction. The apparent second-order rate constant for proton transfer from the bulk solution to E(I-286) is >10(13) M(-1) s(-1), which indicates that the proton uptake is assisted by a local buffer consisting of protonatable residues at the protein surface.     

Proton-transport mechanisms in cytochrome c oxidase revealed by studies of kinetic isotope effects

Cytochrome c oxidase (CytcO) is a membrane-bound enzyme, which catalyzes the reduction of di-oxygen to water and uses a major part of the free energy released in this reaction to pump protons across the membrane. In the Rhodobacter sphaeroides aa? CytcO all protons that are pumped across the membrane, as well as one half of the protons that are used for O? reduction, are transferred through one specific intraprotein proton pathway, which holds a highly conserved Glu286 residue. Key questions that need to be addressed in order to understand the function of CytcO at a molecular level are related to the timing of proton transfers from Glu286 to a "pump site" and the catalytic site, respectively. Here, we have investigated the temperature dependencies of the H/D kinetic-isotope effects of intramolecular proton-transfer reactions in the wild-type CytcO as well as in two structural CytcO variants, one in which proton uptake from solution is delayed and one in which proton pumping is uncoupled from O? reduction. These processes were studied for two specific reaction steps linked to transmembrane proton pumping, one that involves only proton transfer (peroxy-ferryl, P→F, transition) and one in which the same sequence of proton transfers is also linked to electron transfer to the catalytic site (ferryl-oxidized, F→O, transition). An analysis of these reactions in the framework of theory indicates that that the simpler, P→F reaction is rate-limited by proton transfer from Glu286 to the catalytic site. When the same proton-transfer events are also linked to electron transfer to the catalytic site (F→O), the proton-transfer reactions might well be gated by a protein structural change, which presumably ensures that the proton-pumping stoichiometry is maintained also in the presence of a transmembrane electrochemical gradient. Furthermore, the present study indicates that a careful analysis of the temperature dependence of the isotope effect should help us in gaining mechanistic insights about CytcO.     

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.     

Cytochrome oxidase: pathways for electron tunneling and proton transfer

 Electrons from cytochrome c, the substrate of cytochrome oxidase, a redox-linked proton pump, are accepted by Cu_A in subunit II. From there they are transferred to the proton pumping machinery in subunit I, cytochrome a and cytochrome a ₃–Cu_B. The reduction of the latter site, which is the dioxygen reducing unit, is coupled to proton uptake. Dioxygen reduction involves a peroxide and a ferryl ion intermediate, and it is the transition between these and back to the resting oxidized enzyme that are coupled to proton pumping. The X-ray structures suggest electron–transfer pathways that can account for the observed rates provided that the reorganization energies are small. They also reveal two proton-transfer pathways, and mutagenesis experiments have shown that one is used for proton uptake during the initial reduction of cytochrome a ₃–Cu_B, whereas the other mediates transfer of the pumped protons. Received: 23 March 1998 / Accepted: 11 May 1998     

Redesign of the proton-pumping machinery of cytochrome c oxidase: proton pumping does not require Glu(I-286)

One of the putative proton-transfer pathways leading from solution toward the binuclear center in many cytochrome c oxidases is the D-pathway, so-called because it starts with a highly conserved aspartate [D(I-132)] residue. Another highly conserved amino acid residue in this pathway, glutamate(I-286), has been indicated to play a central role in the proton-pumping machinery of mitochondrial-type enzymes, a role that requires a movement of the side chain between two distinct positions. In the present work we have relocated the glutamate to the opposite side of the proton-transfer pathway by constructing the double mutant EA(I-286)/IE(I-112). This places the side chain in about the same position in space as in the original enzyme, but does not allow for the same type of movement. The results show that the introduction of the second-site mutation, IE(I-112), in the EA(I-286) mutant enzyme results in an increase of the enzyme activity by a factor of >10. In addition, the double mutant enzyme pumps approximately 0.4 proton per electron. This observation restricts the number of possible mechanisms for the operation of the redox-driven proton pump. The proton-pumping machinery evidently does require the presence of a protonatable/polar residue at a specific location in space, presumably to stabilize an intact water chain. However, this residue does not necessarily have to be at a strictly conserved location in the amino acid sequence. In addition, the results indicate that E(I-286) is not the "proton gate" of cytochrome c oxidase controlling the flow of pumped protons from one to the other side of the membrane.     

Identification of the proton pathway in bacterial reaction centers: decrease of proton transfer rate by mutation of surface histidines at H126 and H128 and chemical rescue by imidazole identifies the initial proton donors.

The pathway for proton transfer to Q(B) was studied in the reaction center (RC) from Rhodobacter sphaeroides. The binding of Zn(2+) or Cd(2+) to the RC surface at His-H126, His-H128, and Asp-H124 inhibits the rate of proton transfer to Q(B), suggesting that the His may be important for proton transfer [Paddock, M. L., Graige, M. S., Feher, G. and Okamura, M. Y. (1999) Proc. Natl. Acad. Sci. U.S.A. 96, 6183-6188]. To assess directly the role of the histidines, mutant RCs were constructed in which either one or both His were replaced with Ala. In the single His mutant RCs, no significant effects were observed. In contrast, in the double mutant RC at pH 8.5, the observed rates of proton uptake associated with both the first and the second proton-coupled electron-transfer reactions k(AB)(()(1)()) [Q(A)(-)(*)Q(B)-Glu(-) + H(+) --> Q(A)(-)(*)Q(B)-GluH --> Q(A)Q(B)(-)(*)-GluH] and k(AB)(()(2)()) [Q(A)(-)(*)Q(B)(-)(*) + H(+) --> Q(A)(-)(*)(Q(B)H)(*) --> Q(A)(Q(B)H)(-)], were found to be slowed by factors of approximately 10 and approximately 4, respectively. Evidence that the observed changes in the double mutant RC are due to a reduction in the proton-transfer rate constants are provided by the observations: (i) k(AB)(1) at pH approximately pK(a) of GluH became biphasic, indicating that proton transfer is slower than electron transfer and (ii) k(AB)(2) became independent of the driving force for electron transfer, indicating that proton transfer is the rate-limiting step. These changes were overcome by the addition of exogenous imidazole which acts as a proton donor in place of the imidazole groups of His that were removed in the double mutant RC. Thus, we conclude that His-H126 and His-H128 facilitate proton transfer into the RC, acting as RC-bound proton donors at the entrance of the proton-transfer pathways.     

Kinetics of proton-linked flavin conformational changes in p-hydroxybenzoate hydroxylase

p-Hydroxybenzoate hydroxylase (PHBH) is an FAD-dependent monooxygenase that catalyzes the hydroxylation of p-hydroxybenzoate (pOHB) to 3,4-dihydroxybenzoate in an NADPH-dependent reaction. Two structural features are coupled to control the reactivity of PHBH with NADPH: a proton-transfer network that allows protons to be passed between the sequestered active site and solvent and a flavin that adopts two positions: "in", where the flavin is near pOHB, and "out", where the flavin is near NADPH. PHBH uses the proton-transfer network to test for the presence of a suitable aromatic substrate before allowing the flavin to adopt the NADPH-accessible conformation. In this work, kinetic analysis of the His72Asn mutant, with a disrupted proton-transfer network, showed that flavin movement could occur in the presence or absence of NADPH but that NADPH stimulated movement to the reactive conformation required for hydride transfer. Substrate and solvent isotope effects on the transient kinetics of reduction of the His72Asn mutant showed that proton transfer was linked to flavin movement and that the conformational change occurred in a step separate from that of hydride transfer. Proton transfers during the reductive half-reaction were observed directly in the wild-type enzyme by performing experiments in the presence of a fluorescent pH-indicator dye in unbuffered solutions. NADPH binding caused rapid proton release from the enzyme, followed by proton uptake after flavin reduction. Solvent and substrate kinetic isotope effects showed that proton-coupled flavin movement and reduction also occurred in different steps in wild-type PHBH. These results allow a detailed kinetic scheme to be proposed for the reductive half-reaction of the wild-type enzyme. Three kinetic models considered for substrate-induced isomerization are analyzed in the Appendix.     

Cytochrome bc1 complexes of microorganisms.

The cytochrome bc1 complex is the most widely occurring electron transfer complex capable of energy transduction. Cytochrome bc1 complexes are found in the plasma membranes of phylogenetically diverse photosynthetic and respiring bacteria, and in the inner mitochondrial membrane of all eucaryotic cells. In all of these species the bc1 complex transfers electrons from a low-potential quinol to a higher-potential c-type cytochrome and links this electron transfer to proton translocation. Most bacteria also possess alternative pathways of quinol oxidation capable of circumventing the bc1 complex, but these pathways generally lack the energy-transducing, protontranslocating activity of the bc1 complex. All cytochrome bc1 complexes contain three electron transfer proteins which contain four redox prosthetic groups. These are cytochrome b, which contains two b heme groups that differ in their optical and thermodynamic properties; cytochrome c1, which contains a covalently bound c-type heme; and a 2Fe-2S iron-sulfur protein. The mechanism which links proton translocation to electron transfer through these proteins is the proton motive Q cycle, and this mechanism appears to be universal to all bc1 complexes. Experimentation is currently focused on understanding selected structure-function relationships prerequisite for these redox proteins to participate in the Q-cycle mechanism. The cytochrome bc1 complexes of mitochondria differ from those of bacteria, in that the former contain six to eight supernumerary polypeptides, in addition to the three redox proteins common to bacteria and mitochondria. These extra polypeptides are encoded in the nucleus and do not contain redox prosthetic groups. The functions of the supernumerary polypeptides of the mitochondrial bc1 complexes are generally not known and are being actively explored by genetically manipulating these proteins in Saccharomyces cerevisiae.     

Molecular modeling of the bacterial outer membrane receptor energizer, ExbBD/TonB, based on homology with the flagellar motor, MotAB

The MotA/MotB proteins serve as the motor that drives bacterial flagellar rotation in response to the proton motive force (pmf). They have been shown to comprise a transmembrane proton pathway. The ExbB/ExbD/TonB protein complex serves to energize transport of iron siderophores and vitamin B₁₂ across the outer membrane of the Gram-negative bacterial cell using the pmf. These two protein complexes have the same topology and are homologous. Based on molecular data for the MotA/MotB proteins, we propose simple three-dimensional channel structures for both MotA/MotB and ExbB/ExbD/TonB using modeling methods. Features of the derived channels are discussed, and two possible proton transfer pathways for the ExbBD/TonB system are proposed. These analyses provide a guide for molecular studies aimed at elucidating the mechanism by which chemiosmotic energy can be transferred either between two adjacent membranes to energize outer membrane transport or to the bacterial flagellum to generate torque.     

Molecular modeling of the bacterial outer membrane receptor energizer,ExbBD/TonB,based on homology with the flagellar motor,MotAB

The MotA/MotB proteins serve as the motor that drives bacterial flagellar rotation in response to the proton motive force (pmf). They have been shown to comprise a transmembrane proton pathway. The ExbB/ExbD/TonB protein complex serves to energize transport of iron siderophores and vitamin B12 across the outer membrane of the Gram-negative bacterial cell using the pmf. These two protein complexes have the same topology and are homologous. Based on molecular data for the MotA/MotB proteins, we propose simple three-dimensional channel structures for both MotA/MotB and ExbB/ExbD/TonB using modeling methods. Features of the derived channels are discussed, and two possible proton transfer pathways for the ExbBD/TonB system are proposed. These analyses provide a guide for molecular studies aimed at elucidating the mechanism by which chemiosmotic energy can be transferred either between two adjacent membranes to energize outer membrane transport or to the bacterial flagellum to generate torque.     

Uncoupling protein-1 (UCP1) contributes to the basal proton conductance of brown adipose tissue mitochondria

Proton leak pathways uncouple substrate oxidation from ATP synthesis in mitochondria. These pathways are classified as basal (not regulated) or inducible (activated and inhibited). Previously it was found that over half of the basal proton conductance of muscle mitochondria was catalyzed by the adenine nucleotide translocase (ANT), an abundant mitochondrial anion carrier protein. To determine whether ANT is the unique protein catalyst, or one of many proteins that catalyze basal proton conductance, we measured proton leak kinetics in mitochondria isolated from brown adipose tissue (BAT). BAT can express another mitochondrial anion carrier, UCP1, at concentrations similar to ANT. Basal proton conductance was measured under conditions where UCP1 and ANT were catalytically inactive and was found to be lower in mitochondria from UCP1 knockout mice compared to wild-type. Ablation of another abundant inner membrane protein, nicotinamide nucleotide transhydrogenase, had no effect on proton leak kinetics in mitochondria from liver, kidney or muscle, showing that basal proton conductance is not catalyzed by all membrane proteins. We identify UCP1 as a second protein propagating basal proton leak, lending support to the hypothesis that basal leak pathways are perpetrated by members of the mitochondrial anion carrier family but not by other mitochondrial inner membrane proteins.     

Localized control of proton transfer through the D-pathway in cytochrome c oxidase: application of the proton-inventory technique

In the reaction cycle of cytochrome c oxidase from Rhodobacter sphaeroides, one of the steps that are coupled to proton pumping, the oxo-ferryl-to-oxidized transition (F --> O), displays a large kinetic deuterium isotope effect of about 7. In this study we have investigated in detail the dependence of the kinetics of this reaction step ?k(FO)(chi) on the fraction (chi) D(2)O in the enzyme solution (proton-inventory technique). According to a simplified version of the Gross-Butler equation, from the shape of the graph describing k(FO)(chi)/k(FO)(0), conclusions can be drawn concerning the number of protonatable sites involved in the rate-limiting proton-transfer reaction step. Even though the proton-transfer reaction during the F --> O transition takes place over a distance of at least 30 A and involves a large number of protonatable sites, the proton-inventory analysis displayed a linear dependence, which indicates that the entire deuterium isotope effect of 7 is associated with a single protonatable site. On the basis of experiments with site-directed mutants of cytochrome c oxidase, this localized proton-transfer rate control is proposed to be associated with glutamate (I-286) in the D-pathway. Consequently, the results indicate that proton transfer from the glutamate controls the rate of all events during the F --> O reaction step. The proton-inventory analysis of the overall enzyme turnover reveals a nonlinear plot characteristic of at least two protonatable sites involved in the rate-limiting step in the transition state, which indicates that this step does not involve proton transfer through the same pathway (or through the same mechanism) as during the F --> O transition.     

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].     

Common themes and problems of bioenergetics and voltage-gated proton channels

The existence of a proton-selective pathway through a protein is a common feature of voltage-gated proton channels and a number of molecules that play pivotal roles in bioenergetics. Although the functions and structures of these molecules are quite diverse, the proton conducting pathways share a number of fundamental properties. Conceptual parallels include the translocation by hydrogen-bonded chain mechanisms, problems of supply and demand, equivalence of chemical and electrical proton gradients, proton wells, alternating access sites, pK(a) changes induced by protein conformational change, and heavy metal participation in proton transfer processes. An archetypal mechanism involves input and output proton pathways (hydrogen-bonded chains) joined by a regulatory site that switches the accessibility of the bound proton from one 'channel' to the other, by means of a pK(a) change, molecular movement, or both. Although little is known about the structure of voltage-gated proton channels, they appear to share many of these features. Evidently, nature has devised a limited number of mechanisms to accomplish various design strategies, and these fundamental mechanisms are repeated with minor variation in many superficially disparate molecules.     

Glutathione reductase: solvent equilibrium and kinetic isotope effects

Glutathione reductase catalyzes the NADPH-dependent reduction of oxidized glutathione (GSSG). The kinetic mechanism is ping-pong, and we have investigated the rate-limiting nature of proton-transfer steps in the reactions catalyzed by the spinach, yeast, and human erythrocyte glutathione reductases using a combination of alternate substrate and solvent kinetic isotope effects. With NADPH or GSSG as the variable substrate, at a fixed, saturating concentration of the other substrate, solvent kinetic isotope effects were observed on V but not V/K. Plots of Vm vs mole fraction of D2O (proton inventories) were linear in both cases for the yeast, spinach, and human erythrocyte enzymes. When solvent kinetic isotope effect studies were performed with DTNB instead of GSSG as an alternate substrate, a solvent kinetic isotope effect of 1.0 was observed. Solvent kinetic isotope effect measurements were also performed on the asymmetric disulfides GSSNB and GSSNP by using human erythrocyte glutathione reductase. The Km values for GSSNB and GSSNP were 70 microM and 13 microM, respectively, and V values were 62 and 57% of the one calculated for GSSG, respectively. Both of these substrates yield solvent kinetic isotope effects greater than 1.0 on both V and V/K and linear proton inventories, indicating that a single proton-transfer step is still rate limiting. These data are discussed in relationship to the chemical mechanism of GSSG reduction and the identity of the proton-transfer step whose rate is sensitive to solvent isotopic composition. Finally, the solvent equilibrium isotope effect measured with yeast glutathione reductase is 4.98, which allows us to calculate a fractionation factor for the thiol moiety of GSH of 0.456.     

Structural elements involved in electron-coupled proton transfer in cytochrome c oxidase

Haem-copper oxidases are the last components of the respiratory chains in aerobic organisms. These membrane-bound enzymes energetically couple the electron transfer (eT) reactions associated with reduction of dioxygen to water, to proton pumping across the membrane. Even though the mechanism of proton pumping at the molecular level still remains to be uncovered, recent progress has presented us with the structural features of the pumping machinery and detailed information about the eT and proton-transfer reactions associated with the pumping process.     

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.