Structure of the photosynthetic reaction centre from Rhodobacter sphaeroides at 2.65 A resolution: cofactors and protein-cofactor interactions

New development of photoinduced electron-transfer catalytic systems

As an alternative to conventional charge-separation functional molecular models based on multi-step long-range electron transfer (ET) within redox cascades, simple donor-acceptor dyads have been developed to attain a long-lived and high-energy charge-separated (CS) state without significant loss of excitation energy. In particular, a simple molecular electron donor-acceptor dyad, 9-mesityl-10-methylacridinium ion (Acr+-Mes), is capable of fast charge separation but extremely slow charge recombination. Such a simple molecular dyad has significant advantages with regard to synthetic feasibility, providing a variety of applications for photoinduced ET catalytic systems, including efficient photocatalytic systems for the solar energy conversion and construction of organic solar cells.

Trapped conformational states of semiquinone (D+*QB-*) formed by B-branch electron transfer at low temperature in Rhodobacter sphaeroides reaction centers

The reaction center (RC) from Rhodobacter sphaeroides captures light energy by electron transfer between quinones QA and QB, involving a conformational gating step. In this work, conformational states of D+*QB-* were trapped (80 K) and studied using EPR spectroscopy in native and mutant RCs that lack QA in which QB was reduced by the bacteriopheophytin along the B-branch. In mutant RCs frozen in the dark, a light induced EPR signal due to D+*QB-* formed in 30% of the sample with low quantum yield (0.2%-20%) and decayed in 6 s. A small signal with similar characteristics was also observed in native RCs. In contrast, the EPR signal due to D+*QB-* in mutant RCs illuminated while freezing formed in approximately 95% of the sample did not decay (tau >107 s) at 80 K (also observed in the native RC). In all samples, the observed g-values were the same (g = 2.0026), indicating that all active QB-*'s were located in a proximal conformation coupled with the nonheme Fe2+. We propose that before electron transfer at 80 K, the majority (approximately 70%) of QB, structurally located in the distal site, was not stably reducible, whereas the minority (approximately 30%) of active configurations was in the proximal site. The large difference in the lifetimes of the unrelaxed and relaxed D+*QB-* states is attributed to the relaxation of protein residues and internal water molecules that stabilize D+*QB-*. These results demonstrate energetically significant conformational changes involved in stabilizing the D+*QB-* state. The unrelaxed and relaxed states can be considered to be the initial and final states along the reaction coordinate for conformationally gated electron transfer.

Protein-cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: II. Geometry of the hydrogen bonds to the primary quinone formula by 1H and 2H ENDOR spectroscopy

The geometry of the hydrogen bonds to the two carbonyl oxygens of the semiquinone Q(A)(. -) in the reaction center (RC) from the photosynthetic purple bacterium Rhodobacter sphaeroides R-26 were determined by fitting a spin Hamiltonian to the data derived from (1)H and (2)H ENDOR spectroscopies at 35 GHz and 80 K. The experiments were performed on RCs in which the native Fe(2+) (high spin) was replaced by diamagnetic Zn(2+) to prevent spectral line broadening of the Q(A)(. -) due to magnetic coupling with the iron. The principal components of the hyperfine coupling and nuclear quadrupolar coupling tensors of the hydrogen-bonded protons (deuterons) and their principal directions with respect to the quinone axes were obtained by spectral simulations of ENDOR spectra at different magnetic fields on frozen solutions of deuterated Q(A)(. -) in H(2)O buffer and protonated Q(A)(. -) in D(2)O buffer. Hydrogen-bond lengths were obtained from the nuclear quadrupolar couplings. The two hydrogen bonds were found to be nonequivalent, having different directions and different bond lengths. The H-bond lengths r(OH) are 1.73 +/- 0.03 Angstrom and 1.60 +/- 0.04 Angstrom, from the carbonyl oxygens O(1) and O(4) to the NH group of Ala M260 and the imidazole nitrogen N(delta) of His M219, respectively. The asymmetric hydrogen bonds of Q(A)(. -) affect the spin density distribution in the quinone radical and its electronic structure. It is proposed that the H-bonds play an important role in defining the physical properties of the primary quinone, which affect the electron transfer processes in the RC.

Crystal structure of phosphatidylglycerophosphatase (PGPase), a putative membrane-bound lipid phosphatase, reveals a novel binuclear metal binding site and two "proton wires"

Phosphatidylglycerophosphatase (PGPase), an enzyme involved in lipid metabolism, catalyzes formation of phosphatidylglycerol from phosphatidylglycerophosphate. Phosphatidylglycerol is a multifunctional phospholipid, found in the biological membranes of many organisms. Here, we report the crystal structure of Listeria monocytogenes PGPase at 1.8 A resolution. PGPase, an all-helical molecule, forms a homotetramer. Each protomer contains an independent active site with two metal ions, Ca(2+) and Mg(2+), forming a hetero-binuclear center located in a hydrophilic cavity near the surface of the molecule. The binuclear center, conserved ligands, metal-bound water molecules, and an Asp-His dyad form the active site. The catalytic mechanism of this enzyme is likely to proceed via binuclear metal activated nucleophilic water. The binuclear metal-binding active-site environment of this structure should provide insights into substrate binding and metal-dependent catalysis. A long channel with inter-linked linear water chains, termed "proton wires," is observed at the tetramer interface. Comparison of similar water chain structures in photosynthetic reaction centers (RCs), Cytochrome f, gramicidin, and bacteriorhodopsin, suggests that PGPase may conduct protons via proton wires.

Charge Delocalization in the Special-Pair Radical Cation of Mutant Reaction Centers of Rhodobacter sphaeroides from Stark Spectra and Nonadiabatic Spectral Simulations

Stark and absorption spectra for the hole-transfer band of the bacteriochlorophyll special pair in the wild-type and L131LH, M160LH, and L131LH/M160LH mutants of the bacterial reaction center of Rhodobacter sphaeroides are presented, along with extensive analyses based on nonadiabatic spectral simulations. Dramatic changes in the Stark spectra are induced by the mutations, changes that are readily interpreted in terms of the redox-energy asymmetry and degree of charge localization in the special-pair radical cation. The effect of mutagenesis on key properties such as the electronic coupling within the special pair and the reorganization energy associated with intervalence hole transfer are determined for the first time. Results for the L131LH and M160LH/L131LH mutants indicate that these species can be considered as influencing the special pair primarily through modulation of the redox asymmetry, as is usually conceptualized, but M160LH is shown to develop a wide range of effects that can be interpreted in terms of significant mutation-induced structural changes in and around the special pair. The nonadiabatic spectra simulations are performed using both a simple two-state 1-mode and an extensive four-state 70-mode model, which includes the descriptions of additional electronic states and explicitly treats the major vibrational modes involved. Excellent agreement between the two simulation approaches is obtained. The simple model is shown to reproduce key features of the Stark effect of the main intervalence transition, while the extensive model quantitatively reproduces most features of the observed spectra for both the electronic and the phase-phonon regions, thus giving a more comprehensive description of the effect of the mutations on the properties of the special-pair radical cation. These results for a series of closely related mixed-valence complexes show that the Stark spectra provide a sensitive indicator for the properties of the mixed-valence complexes and should serve as an instructive example on the application of nonadiabatic simulations to the study of mixed-valence complexes in general as well as other chemical systems akin to the photosynthetic special pair.

Lipids in photosynthetic reaction centres: structural roles and functional holes

Photosynthetic proteins power the biosphere. Reaction centres, light harvesting antenna proteins and cytochrome b(6)f (or bc(1)) complexes are expressed at high levels, have been subjected to an intensive spectroscopic, biochemical and mutagenic analysis, and several have been characterised to an informatively high resolution by X-ray crystallography. In addition to revealing the structural basis for the transduction of light energy, X-ray crystallography has brought molecular insights into the relationships between these multicomponent membrane proteins and their lipid environment. Lipids resolved in the X-ray crystal structures of photosynthetic proteins bind light harvesting cofactors, fill intra-protein cavities through which quinones can diffuse, form an important part of the monomer-monomer interface in multimeric structures and may facilitate structural flexibility in complexes that undergo partial disassembly and repair. It has been proposed that individual lipids influence the biophysical properties of reaction centre cofactors, and so affect the rate of electron transfer through the complex. Lipids have also been shown to be important for successful crystallisation of photosynthetic proteins. Comparison of the three types of reaction centre that have been structurally characterised reveals interesting similarities in the position of bound lipids that may point towards a generic requirement to reinforce the structure of the core electron transfer domain. The crystallographic data are also providing new opportunities to find molecular explanations for observed effects of different types of lipid on the structure, mechanism and organisation of reaction centres and other photosynthetic proteins.

Photoactivation of the photosynthetic reaction center of Blastochloris viridis in the crystalline state

Photoactivation in crystals of the bacterial reaction center of Blastochloris viridis was investigated by near-infrared spectroscopy. The bleaching of the special pair absorption at 970 nm and the simultaneous rise of the special pair cation absorption at 1300 nm were measured in response to transient irradiation by a HeNe laser over 5 orders of magnitude in laser power. The resulting power-saturation curve can be used to estimate the true extent of photoactivation achieved in a prior time-resolved crystallographic experiment (Baxter et al. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 5982-5987). The overall extent of photoactivation was 50%, which demonstrates that the time-resolved crystallographic method can be applied to the optically dense reaction center crystals. Measurement of the charge-recombination rate, however, suggests the presence of a long-lived P+ state within the crystal.

Electron-transfer dynamics of photosynthetic reaction centers in thermoresponsive soft materials

Poly(ethylene glycol)-grafted, lipid-based, thermoresponsive, soft nanostructures are shown to serve as scaffolding into which reconstituted integral membrane proteins, such as the bacterial photosynthetic reaction centers (RCs) can be stabilized, and their packing arrangement, and hence photophysical properties, can be controlled. The self-assembled nanostructures exist in two distinct states: a liquid-crystalline gel phase at temperatures above 21 degrees C and a non-birefringent, reduced viscosity state at lower temperatures. Characterization of the effect of protein introduction on the mesoscopic structure of the materials by 31P NMR and small-angle X-ray scattering shows that the expanded lamellar structure of the protein-free material is retained. At reduced temperatures, however, the aggregate structure is found to convert from a two-dimensional normal hexagonal structure to a three-dimensional cubic phase upon introduction of the RCs. Structural and functional characteristics of the RCs were determined by ground-state and femtosecond transient absorption spectroscopy. Time-resolved results indicate that the kinetics of primary electron transfer for the RCs in the low-viscosity cold phase of the self-assembled nanostructures are identical to those observed in a detergent-solubilized state in buffered aqueous solutions (approximately 4 ps) over a wide range of protein concentrations and experimental conditions. This is also true for RCs held within the lamellar gel phase at low protein concentrations and at short sample storage times. In contrast are kinetics from samples that are prepared with high RC concentrations and stored for several hours, which display additional kinetic components with extended electron-transfer times (approximately 10-12 ps). This observation is tentatively attributed to energy transfer between RCs that have laterally (in-plane) organized within the lipid bilayers of the lamellar gel phase prior to charge separation. These results not only demonstrate the use of soft nanostructures as a matrix in which to stabilize and organize membrane proteins but also suggest the possibility of using them to control the interactions between proteins and thus to tune their collective optical/electronic properties.

Mechanism of Carotenoid Singlet Excited State Energy Transfer in Modified Bacterial Reaction Centers

Ultrafast transient laser spectroscopy has been used to investigate carotenoid singlet excited state energy transfer in various Rhodobacter (Rb.) sphaeroides reaction centers (RCs) modified either genetically or chemically. The pathway and efficiency of energy transfer were examined as a function of the structures and energies of the donor and acceptor molecules. On the donor side, carotenoids with various extents of pi-electron conjugation were examined. RCs studied include those from the anaerobically grown wild-type strain containing the carotenoid spheroidene, which has 10 conjugated carbon-carbon double bonds; the GA strain containing neurosporene, which has nine conjugated double bonds; and aerobically grown wild-type cells, as well as aerobically grown H(M182)L mutant, both containing the carbonyl-containing carotenoid spheroidenone, which has 11 conjugated double bonds. By varying the structure of the carotenoid, we observed the effect of altering the energies of the carotenoid excited states on the rate of energy transfer. Both S(1)- and S(2)-mediated carotenoid-to-bacteriochlorophyll energy transfer processes were observed. The highest transfer efficiency, from both the S(1) and S(2) states, was observed using the carotenoid with the shortest chain. The S(1)-mediated carotenoid-to- bacteriochlorophyll energy transfer efficiencies were determined to be 96%, 84%, and 73% for neurosporene, spheroidene, and spheroidenone, respectively. The S(2)-mediated energy transfer efficiencies follow the same trend but could not be determined quantitatively because of limitations in the time resolution of the instrumentation. The dependence of the energy transfer rate on the energetics of the energy transfer acceptor was verified by performing measurements with RCs from the H(M182)L mutant. In this mutant, the bacteriochlorophyll (denoted B(B)) located between the carotenoid and the RC special pair (P) is replaced by a bacteriopheophytin (denoted phi(B)), where the Q(X) and Q(Y) bands of phi(B) are 1830 and 1290 cm(-1), respectively, higher in energy than those of B(B). These band shifts associated with phi(B) in the H(M182)L mutant significantly alter the spectral overlap between the carotenoid and phi(B), resulting in a significant decrease of the transfer efficiency from the carotenoid S(1) state to phi(B). This leaves energy transfer from the carotenoid S(2) state to phi(B) as the dominant channel. Largely because of this change in mechanism, the overall efficiency of energy transfer from the carotenoid to P decreases to less than 50% in this mutant. Because the spectral signature of phi(B) is different from that of B(A) in this mutant, we were able to demonstrate clearly that the carotenoid-to-P energy transfer is via phi(B). This finding supports the concept that, in wild-type RCs, the carotenoid-to-P energy transfer occurs through the cofactor located at the B(B) position.

How photosynthetic reaction centers control oxidation power in chlorophyll pairs P680, P700, and P870

At the heart of photosynthetic reaction centers (RCs) are pairs of chlorophyll a (Chla), P700 in photosystem I (PSI) and P680 in photosystem II (PSII) of cyanobacteria, algae, or plants, and a pair of bacteriochlorophyll a (BChla), P870 in purple bacterial RCs (PbRCs). These pairs differ greatly in their redox potentials for one-electron oxidation, E(m). For P680, E(m) is 1,100-1,200 mV, but for P700 and P870, E(m) is only 500 mV. Calculations with the linearized Poisson-Boltzmann equation reproduce these measured E(m) differences successfully. Analyzing the origin for these differences, we found as major factors in PSII the unique Mn(4)Ca cluster (relative to PSI and PbRC), the position of P680 close to the luminal edge of transmembrane alpha-helix d (relative to PSI), local variations in the cd loop (relative to PbRC), and the intrinsically higher E(m) of Chla compared with BChla (relative to PbRC).

Analysis of absorption spectra of purple bacterial reaction centers in the near infrared region by higher order derivative spectroscopy

Reaction centers (RCs) of purple bacteria are uniquely suited objects to study the mechanisms of the photosynthetic conversion of light energy into chemical energy. A recently introduced method of higher order derivative spectroscopy [I.K. Mikhailyuk, H. Lokstein, A.P. Razjivin, A method of spectral subband decomposition by simultaneous fitting the initial spectrum and a set of its derivatives, J. Biochem. Biophys. Methods 63 (2005) 10-23] was used to analyze the NIR absorption spectra of RC preparations from Rhodobacter (R.) sphaeroides strain 2R and Blastochloris (B.) viridis strain KH, containing bacteriochlorophyll (BChl) a and b, respectively. Q(y) bands of individual RC porphyrin components (BChls and bacteriopheophytins, BPheo) were identified. The results indicate that the upper exciton level P(y+) of the photo-active BChl dimer in RCs of R. sphaeroides has an absorption maximum of 810nm. The blue shift of a complex integral band at approximately 800nm upon oxidation of the RC is caused primarily by bleaching of P(y+), rather than by an electrochromic shift of the absorption band(s) of the monomeric BChls. Likewise, the disappearance of a band peaking at 842nm upon oxidation of RCs from B. viridis indicates that this band has to be assigned to P(y+). A blue shift of an absorption band at approximately 830nm upon oxidation of RCs of B. viridis is also essentially caused by the disappearance of P(y+), rather than by an electrochromic shift of the absorption bands of monomeric BChls. Absorption maxima of the monomeric BChls, B(B) and B(A) are at 802 and 797nm, respectively, in RCs of R. sphaeroides at room temperature. BPheo co-factors H(B) and H(A) peak at 748 and 758nm, respectively, at room temperature. For B. viridis RCs the spectral positions of H(B) and H(A) were found to be 796 and 816nm, respectively, at room temperature.

Effects of electric fields on proton transport through water chains

Molecular dynamics simulations on quantum energy surfaces are carried out to study the effects of perturbing electric fields on proton transport (PT) in protonated water chains. As an idealized model of a hydrophobic cavity in the interior of a protein the water molecules are confined into a carbon nanotube (CNT). The water chain connects a hydrated hydronium ion (H3O+) at one end of the CNT and an imidazole molecule at the other end. Without perturbing electric fields PT from the hydronium proton donor to the imidazole acceptor occurs on a picosecond time scale. External perturbations to PT are created by electric fields of varying intensities, normal to the CNT axis, generated by a neutral pair of charges on the nanotube wall. For fields above approximately 0.5 VA, the hydronium ion is effectively trapped at the CNT center, and PT blocked. Fields of comparable strength are generated inside proteins by nearby polar/charged amino acids. At lower fields the system displays a rich dynamic behavior, where the excess charge shuttles back and forth along the water chain before reaching the acceptor group on the picosecond time scale. The effects of the perturbing field on the proton movement are analyzed in terms of structural and dynamic properties of the water chain. The implications of these observations on PT in biomolecular systems and its control by external perturbing fields are discussed.

Reconstruction of a kinetic model of the chromatophore vesicles from Rhodobacter sphaeroides

We present a molecular model of a chromatophore vesicle from Rhodobacter sphaeroides. These vesicles are ideal benchmark systems for molecular and systemic simulations, because they have been well studied, they are small, and they are naturally separated from their cellular environment. To set up a photosynthetic chain working under steady-state conditions, we compiled from the experimental literature the specific activities and geometries that have been determined for their constituents. This data then allowed defining the stoichiometries for all membrane proteins. This article contains the kinetic part of the reconstructed model, while the spatial reconstruction is presented in a companion article. By considering the transport properties of the Cytochrome c(2) and ubiquinone pools, we show that their size and oxidation states allow for an efficient buffering of the statistical fluctuations that arise from the small size of the vesicles. Stoichiometric and kinetic considerations indicate that a typical chromatophore vesicle of Rb. sphaeroides with a diameter of 45 nm should contain approximately five bc(1) monomers.

Ubiquinone (coenzyme Q10) binding sites: Low dielectric constant of the gate prevents the escape of the semiquinone

The photosynthetic reaction center (RC) from purple bacteria is frequently used as a model for the interaction of ubiquinones (coenzyme Q) with membrane proteins. Single-turnover flash activation of RC leads to formation of the semiquinone (SQ) of the secondary acceptor quinone after odd flashes and quinol after even flashes. The ubiquinol escapes the binding site in 1 ms, while the SQ does not leave the binding site for at least 5 min. Observed difference between these times suggests a large energetic barrier for the SQ. However, high apparent dielectric constant in the vicinity of the quinone ring (>or=25) results in a relatively small electrostatic energy of SQ stabilization. To resolve this apparent contradiction I suggest that a significant part of the kinetic stabilization of the SQ is achieved by the special topology of the binding site in which quinone can exit the binding site only by moving its headgroup toward the center of the membrane. The large energetic penalty of transferring the charged headgroup to the membrane dielectric can explain the observed kinetic stability of the SQ.

Triplet state spectra and dynamics of geometric isomers of carotenoids

The observation of preferential binding of cis-carotenoids in purple bacterial photosynthetic reaction centers versus trans-isomers in antenna pigment protein complexes has led to the hypothesis that the natural selection of stereoisomers has physiological significance. In order to test this hypothesis, we have undertaken a systematic series of investigations comparing the optical spectroscopic properties and excited state dynamics of cis and trans isomers of carotenoids. The present work compares the triplet state spectra, lifetimes, and energy transfer rates of all-trans-spheroidene and 13,14-locked-cis-spheroidene, the latter of which is incapable of isomerizing to the all-trans configuration, and therefore provides a unique opportunity to examine the triplet state properties of a structurally stable cis molecule. The data reveal only small differences in spectra, decay dynamics, and transfer times and suggest there is little intrinsic advantage in either triplet energy transfer or triplet state decay arising from the inherently different isomeric forms of cis compared to trans carotenoids.

Low-frequency resonance Raman studies of the H(M202)G cavity mutant of bacterial photosynthetic reaction centers

Low-frequency (90-435 cm(-1)) NIR-excitation (875-900 nm) resonance Raman (RR) studies are reported for the H(M202)G cavity mutant of bacterial photosynthetic reaction centers (RCs) from Rb. sphaeroides that was first described by Goldsmith et al. [(1996) Biochemistry 35: 2421-2428]. In this mutant, the His residue that axially ligates the Mg ion of the M-side bacteriochlorophyll (BChl) of the special pair primary donor (P) is replaced by a non-ligating Gly residue. Regardless, the Mg ion of P(M) in the H(M202)G RCs remains pentacoordinates and is presumably ligated by a water molecule, although this axial ligand has not been definitively identified. The low-frequency RR studies of the H(M202)G RCs are accompanied by studies of RCs exchanged with D(2)O and incubated with imidazole (Im). The RR studies of the cavity mutant RCs reveal the following: (1) The structure of P(M) in the H(M202)G RCs is different from that of the wild-type, consistent with an altered BChl core. (2) A water ligand for P(M) in the H(M202)G RCs is generally consistent with the low-frequency RR spectra. The Mg-OH(2) stretching vibration is tentatively assigned to a band at 318 cm(-1), a frequency higher than that of the Mg-His stretch of the native pigment ( approximately approximately 235 cm(-1)). (3) The BChl core structure of P(M) in the cavity mutant is rendered similar (but not identical) to that of the wild-type when the adventitious water axial ligand is replaced by Im. (4) Exchange with D(2)O results in more global structural changes, likely involving the protein, which in turn affect the structure of the BChls in P. (5) Assignment of the low-frequency vibrational spectrum of P is generally more complex than originally suggested.

Kinetic analysis of the thermal stability of the photosynthetic reaction center from Rhodobacter sphaeroides

The temperature-induced denaturation of the photosynthetic reaction center from Rhodobacter sphaeroides has been studied through the changes that occur in the absorption spectrum of the bound chromophores on heating. At elevated temperatures, the characteristic absorbance bands of the bacteriochlorins bound to the polypeptides within the reaction center are lost, and are replaced by features typical of unbound bacteriochlorophyll and bacteriopheophytin. The kinetics of the spectral changes cannot be explained by a direct conversion from the functional to the denatured form of the protein, and require the presence of at least one intermediate. Possible mechanisms for the transformation via an intermediate are examined using a global analysis of the kinetic data, and the most likely mechanism is shown to involve a reversible transformation between the native state and an off-pathway intermediate, coupled to an irreversible transformation to the denatured state. The activation energies for the transformations between the three components are calculated from the effect of temperature on the individual rate constants, and the likely structural changes of the protein during the temperature-induced transformation are discussed.

Electronic Transitions of the Soret Band of Reaction Centers from Rhodobacter sphaeroides Studied by Femtosecond Transient Absorbance Spectroscopy

The Soret band of reaction centers from Rhodobacter sphaeroides has been systematically studied using femtosecond transient absorption spectroscopy. When the excitation wavelength was scanned over the entire Soret band, the approximate absorption spectra of the bacteriochlorophyll dimer, the monomer bacteriochlorophylls, and the bacteriopheophytins within the Soret band were determined by analyzing the ground state bleaching with about 100 fs resolution. The main contribution of H is on the blue end of the spectrum, peaking near 350 nm, P absorbs mostly on the red side of the spectrum, but probably has multiple bands, and the main absorbance of B likely lies between H and P, overlapping with P on the red side (particularly near 390 nm). The energy transfer from B to P in the QY band takes about 300 fs when Soret-band excitation is used and the time constant of overall energy transfer from H to B to P in the QY band when H is specifically excited near 350 nm is about 500 fs. Internal conversion after Soret-band excitation is the rate-limiting step for the energy-transfer process. The time constant of internal conversion for B and P is less than 300 fs, and for H it is about 500 fs.

Investigation of the Effects of Different Carotenoids on the Absorption and CD Signals of Light Harvesting 1 Complexes

Absorption and circular dichroism (CD) spectra of light-harvesting (LH)1 complexes from the purple bacteria Rhodobacter (Rba.) sphaeroides and Rhodospirillum (Rsp.) rubrum are presented. The complexes exhibit very low intensity, highly nonconservative, near-infrared (NIR) CD spectra. Absorption and CD spectra from several mutant and reconstituted LH1 complexes, with the carotenoid neurosporene and the precursor phytoene replacing the wild-type (WT) carotenoids, are also examined. The experiments show that the position of the carotenoid bands as well as the bacteriochlorophyll (BChl)/carotenoid ratio affect the NIR CD spectra: bluer bands and larger ratios make the NIR CD signal more conservative. Modeling results that support this finding are presented. This study, combined with the theoretical approach of the companion paper, where modeling of such complexes is presented and discussed in detail, provide a complete explanation of the origin of the nonconservative NIR CD spectra of LH1 and B820.

Protein–cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: Effect of hydrogen bonding on the electronic and geometric structure of the primary quinone. A density functional theory study

The effect of hydrogen bonding to the primary quinone (Q(A) and Q(*)(-)(A)) in bacterial reaction centers was studied using density functional theory (DFT) calculations. The charge neutral state Q(A) was investigated by optimizing the hydrogen atom positions of model systems extracted from 15 different X-ray structures. From this analysis, mean values of the H-bond lengths and directions were derived. It was found that the N(delta)-H of His M219 forms a shorter H-bond to Q(A) than the N-H of Ala M260. The H-bond of His M219 is linear and more twisted out of the quinone plane. The radical anion Q(*)(-)(A) in the protein environment was investigated by using a mixed quantum mechanics/molecular mechanics (QM/MM) approach. Two geometry optimizations with a different number of flexible atoms were performed. H-bond lengths were obtained and spectroscopic parameters calculated, i.e. the hyperfine and nuclear quadrupole couplings of magnetic nuclei coupled to the radical. Good agreement was found with the results provided by EPR/ENDOR spectroscopy. This implies that the calculated lengths and directions of the H-bonds to Q(*)(-)(A) are reliable values. From a comparison of the neutral and reduced state of Q(A) it was concluded that the H-bond distances are shortened by approximately 0.17 Angstroms (His M219) and approximately 0.13 Angstroms (Ala M260) upon single reduction of the quinone. It is shown that the point-dipole approximation can not be used for an estimation of H-bond lengths from measured hyperfine couplings in a system with out-of-plane H-bonding. In contrast, the evaluation of the nuclear quadrupole couplings of (2)H nuclei substituted in the hydrogen bonds yields H-bond lengths close to the values that were deduced from DFT geometry optimizations. The significance of hydrogen bonding to the quinone cofactors in biological systems is discussed.

Probing the Contribution of Electronic Coupling to the Directionality of Electron Transfer in Photosynthetic Reaction Centers

Subpicosecond transient absorption studies are reported for a set of Rhodobacter (R.) capsulatus bacterial photosynthetic reaction centers (RCs) designed to probe the origins of the unidirectionality of charge separation via one of two electron transport chains in the native pigment-protein complex. All of the RCs have been engineered to contain a heterodimeric primary electron donor (D) consisting of a bacteriochlorophyll (BChl) and a bacteriopheophytin (BPh). The BPh component of the M heterodimer (Mhd) or L heterodimer (Lhd) is introduced by substituting a Leu for His M200 or His L173, respectively. Previous work on primary charge separation in heterodimer mutants has not included the Lhd RC from R. capsulatus, which we report for the first time. The Lhd and Mhd RCs are used as controls against which we assess RCs that combine the heterodimer mutations with a second mutation (His substituted for Leu at M212) that results in replacement of the native L-side BPh acceptor with a BChl (beta). The transient absorption spectra reveal clear evidence for charge separation to the normally inactive M-side BPh acceptor (H(M)) in Lhd-beta RCs to form D+H(M)- with a yield of approximately 6%. This state also forms in Mhd-beta RCs but with about one-quarter the yield. In both RCs, deactivation to the ground state is the predominant pathway of D decay, as it is in the Mhd and Lhd single mutants. Analysis of the results indicates an upper limit ofV2L/V2m < or = 4 for the contribution of the electronic coupling elements to the relative rates of electron transfer to the L versus M sides of the wild-type RC. In comparison to the L/M rate ratio (kL/kM) approximately 30 for wild-type RCs, our findings indicate that electronic factors contribute approximately 35% at most to directionality with the other 65% deriving from energetic considerations, which includes differences in free energies, reorganization energies, and contributions of one- and two-step mechanisms on the two sides of the RC.

On the role of basic residues in adapting the reaction centre-LH1 complex for growth at elevated temperatures in purple bacteria

The purple photosynthetic bacterium Thermochromatium tepidum is a moderate thermophile, with a growth optimum of 48-50 degrees C. The X-ray crystal structure of the reaction centre from this organism has been determined, and compared with that from mesophilic bacteria such as Blastochloris viridis and Rhodobacter sphaeroides (Nogi T et al. (2000) Proc Natl Acad Sci USA 97: 13561-13566). Structural features that could contribute to the enhanced thermal stability of the Thermochromatium tepidum reaction centre were discussed, including three arginine residues exposed at the periplasmic side of the membrane that are not present in reaction centres from mesophilic organisms, and potentially could increase the affinity of the complex for the surrounding membrane. In the present report these arginine residues, plus a histidine identified from an extensive sequence alignment, were engineered into structurally homologous positions in the Rhodobacter sphaeroides reaction centre, and the effect on the thermal stability of the Rhodobacter sphaeroides complex was examined. We find that these residues do not enhance the thermal stability of the reaction centre, as assessed by absorbance spectroscopy of the bacteriochlorin cofactors in membrane-bound reaction centres. Possible roles of these residues in the Thermochromatium tepidum reaction centre are discussed, and it is proposed that they facilitate stronger binding of the reaction centre to the encircling LH1 antenna complex, through ionic interactions with acidic residues at the C-terminal end of the LH1 alpha-polypeptide. Such an interaction could enhance the stability of the so-called 'RC-LH1 core' complex that is formed between the reaction centre and the LH1 antenna, and which represents the minimal functional photosynthetic unit in all known purple photosynthetic bacteria. Stronger bonding interactions between the two complexes could also contribute to an increase in the rigidity of the photosynthetic membrane in Thermochromatium tepidum, in accord with the general finding that the cytoplasmic membrane from thermophilic eubacteria is less fluid than its counterpart in mesophilic bacteria.

Replacement or exclusion of the B-branch bacteriopheophytin in the purple bacterial reaction centre: The HB cofactor is not required for assembly or core function of the Rhodobacter sphaeroides complex

All of the membrane-embedded cofactors of the purple bacterial reaction centre have well-defined functional or structural roles, with the exception of the bacteriopheophytin (H(B)) located approximately half-way across the membrane on the so-called inactive- or B-branch of cofactors. Sequence alignments indicate that this bacteriochlorin cofactor is a conserved feature of purple bacterial reaction centres, and a pheophytin is also found at this position in the Photosystem-II reaction centre. Possible structural or functional consequences of replacing the H(B) bacteriopheophytin by bacteriochlorophyll were investigated in the Rhodobacter sphaeroides reaction centre through mutagenesis of residue Leu L185 to His (LL185H). Results from absorbance spectroscopy indicated that the LL185H mutant assembled with a bacteriochlorophyll at the H(B) position, but this did not affect the capacity of the reaction centre to support photosynthetic growth, or change the kinetics of charge separation along the A-branch of cofactors. It was also found that mutation of residue Ala M149 to Trp (AM149W) caused the reaction centre to assemble without an H(B) bacteriochlorin, demonstrating that this cofactor is not required for correct assembly of the reaction centre. The absence of a cofactor at this position did not affect the capacity of the reaction centre to support photosynthetic growth, or the kinetics of A-branch electron transfer. A combination of X-ray crystallography and FTIR difference spectroscopy confirmed that the H(B) cofactor was absent in the AM149W mutant, and that this had not produced any significant disturbance of the adjacent ubiquinol reductase (Q(B)) site. The data are discussed with respect to possible functional roles of the H(B) bacteriopheophytin, and we conclude that the reason(s) for conservation of a bacteriopheophytin cofactor at this position in purple bacterial reaction centres are likely to be different from those underlying conservation of a pheophytin at the analogous position in Photosystem-II.

Lipids in and around photosynthetic reaction centres

Reaction centres are membrane-embedded pigment–protein complexes that transduce the energy of sunlight into a biologically useful form. The most heavily studied reaction centres are the PS-I (Photosystem I) and PS-II complexes from oxygenic phototrophs, and the reaction centre from purple photosynthetic bacteria. A great deal is known about the compositions and structures of these reaction centres, and the mechanism of light-activated transmembrane electron transfer, but less is known about how they interact with other components of the photosynthetic membrane, including the membrane lipids. X-ray crystallography has provided high-resolution structures for PS-I and the purple bacterial reaction centre, and revealed binding sites for a number of lipids, either embedded in the protein interior or attached to the protein surface. These lipids play a variety of roles, including the binding of cofactors and the provision of structural support. The challenges of modelling surface-associated electron density features such as lipids, detergents, small amphiphiles and ions are discussed.

Control of Quinone Redox Potentials in Photosystem II: Electron Transfer and Photoprotection

In O(2)-evolving complex Photosystem II (PSII), an unimpeded transfer of electrons from the primary quinone (Q(A)) to the secondary quinone (Q(B)) is essential for the efficiency of photosynthesis. Recent PSII crystal structures revealed the protein environment of the Q(A/B) binding sites. We calculated the plastoquinone (Q(A/B)) redox potentials (E(m)) for one-electron reduction with a full account of the PSII protein environment. We found two different H-bond patterns involving Q(A) and D2-Thr217, resulting in an upshift of E(m)(Q(A)) by 100 mV if the H bond between Q(A) and Thr is present. The formation of this H bond to Q(A) may be the origin of a photoprotection mechanism, which is under debate. At the Q(B) side, the formation of a H bond between D2-Ser264 and Q(B) depends on the protonation state of D1-His252. Q(B) adopts the high-potential form if the H bond to Ser is present. Conservation of this residue and H-bond pattern for Q(B) sites among bacterial photosynthetic reaction centers (bRC) and PSII strongly indicates their essential requirement for electron transfer function.

Ubiquinone reduction in the photosynthetic reaction centre of Rhodobacter sphaeroides: interplay between electron transfer, proton binding and flips of the quinone ring

This review is focused on reactions that gate (control) the electron transfer between the primary quinone QA and secondary quinone QB in the photosynthetic reaction centre of Rhodobacter sphaeroides. The results on electron and proton transfer are discussed in relation to structural information and to the steered molecular dynamics simulations of the QB ring flip in its binding pocket. Depending on the initial position of QB in the pocket and on certain conditions, the rate of electron transfer is suggested to be limited either by the quinone ring flip or by the charge-compensating proton equilibration between the surface and the buried QB site.

Rewiring photosynthesis: engineering wrong-way electron transfer in the purple bacterial reaction centre

The purple bacterial reaction centre uses light energy to separate charge across the cytoplasmic membrane, reducing ubiquinone and oxidizing a c-type cytochrome. The protein possesses a macroscopic structural two-fold symmetry but displays a strong functional asymmetry, with only one of two available membrane-spanning branches of cofactors (the so-called A-branch) being used to catalyse photochemical charge separation. The factors underlying this functional asymmetry have been the subject of study for many years but are still not fully understood. Site-directed mutagenesis has been partially successful in rerouting electron transfer along the normally inactive B-branch, allowing comparison of the kinetics of equivalent electron transfer reactions on the two branches. Both the primary and secondary electron transfer steps on the B-branch appear to be considerably slower than their A-branch counterparts. The effectiveness of different mutations in rerouting electron transfer along the B-branch of cofactors is discussed.

Ubiquinol oxidation in the cytochrome bc1 complex: Reaction mechanism and prevention of short-circuiting

This review is focused on the mechanism of ubiquinol oxidation by the cytochrome bc1 complex (bc1). This integral membrane complex serves as a "hub" in the vast majority of electron transfer chains. The bc1 oxidizes a ubiquinol molecule to ubiquinone by a unique "bifurcated" reaction where the two released electrons go to different acceptors: one is accepted by the mobile redox active domain of the [2Fe-2S] iron-sulfur Rieske protein (FeS protein) and the other goes to cytochrome b. The nature of intermediates in this reaction remains unclear. It is also debatable how the enzyme prevents short-circuiting that could happen if both electrons escape to the FeS protein. Here, I consider a reaction mechanism that (i) agrees with the available experimental data, (ii) entails three traits preventing the short-circuiting in bc1, and (iii) exploits the evident structural similarity of the ubiquinone binding sites in the bc1 and the bacterial photosynthetic reaction center (RC). Based on the latter congruence, it is suggested that the reaction route of ubiquinol oxidation by bc1 is a reversal of that leading to the ubiquinol formation in the RC. The rate-limiting step of ubiquinol oxidation is then the re-location of a ubiquinol molecule from its stand-by site within cytochrome b into a catalytic site, which is formed only transiently, after docking of the mobile redox domain of the FeS protein to cytochrome b. In the catalytic site, the quinone ring is stabilized by Glu-272 of cytochrome b and His-161 of the FeS protein. The short circuiting is prevented as long as: (i) the formed semiquinone anion remains bound to the reduced FeS domain and impedes its undocking, so that the second electron is forced to go to cytochrome b; (ii) even after ubiquinol is fully oxidized, the reduced FeS domain remains docked to cytochrome b until electron(s) pass through cytochrome b; (iii) if cytochrome b becomes (over)reduced, the binding and oxidation of further ubiquinol molecules is hampered; the reason is that the Glu-272 residue is turned towards the reduced hemes of cytochrome b and is protonated to stabilize the surplus negative charge; in this state, this residue cannot participate in the binding/stabilization of a ubiquinol molecule.

Strong Effects of an Individual Water Molecule on the Rate of Light-driven Charge Separation in the Rhodobacter sphaeroides Reaction Center

The role of a water molecule (water A) located between the primary electron donor (P) and first electron acceptor bacteriochlorophyll (B(A)) in the purple bacterial reaction center was investigated by mutation of glycine M203 to leucine (GM203L). The x-ray crystal structure of the GM203L reaction center shows that the new leucine residue packs in such a way that water A is sterically excluded from the complex, but the structure of the protein-cofactor system around the mutation site is largely undisturbed. The results of absorbance and resonance Raman spectroscopy were consistent with either the removal of a hydrogen bond interaction between water A and the keto carbonyl group of B(A) or a change in the local electrostatic environment of this carbonyl group. Similarities in the spectroscopic properties and x-ray crystal structures of reaction centers with leucine and aspartic acid mutations at the M203 position suggested that the effects of a glycine to aspartic acid substitution at the M203 position can also be explained by steric exclusion of water A. In the GM203L mutant, loss of water A was accompanied by an approximately 8-fold slowing of the rate of decay of the primary donor excited state, indicating that the presence of water A is important for optimization of the rate of primary electron transfer. Possible functions of this water molecule are discussed, including a switching role in which the redox potential of the B(A) acceptor is rapidly modulated in response to oxidation of the primary electron donor.

Excitation-Energy Migration in Self-Assembled Cyclic Zinc(II)-Porphyrin Arrays: A Close Mimicry of a Natural Light-Harvesting System

The excitation-energy-hopping (EEH) times within two-dimensional cyclic zinc(II)-porphyrin arrays 5 and 6, which were prepared by intermolecular coordination and ring-closing metathesis reaction of olefins, were deduced by modeling the EEH process based on the anisotropy depolarization as well as the exciton-exciton annihilation dynamics. Assuming the number of energy-hopping sites N = 5 and 6, the two different experimental observables, that is, anisotropy depolarization and exciton-excition annihilation times, consistently give the EEH times of 8.0 +/- 0.5 and 5.3 +/- 0.6 ps through the 1,3-phenylene linkages of 5 and 6, respectively. Accordingly, the self-assembled cyclic porphyrin arrays have proven to be well-defined two-dimensional models for natural light-harvesting complexes.

Conformational regulation of charge recombination reactions in a photosynthetic bacterial reaction center

In bright light the photosynthetic reaction center (RC) of Rhodobacter sphaeroides stabilizes the P(+)(870).Q(-)(A) charge-separated state and thereby minimizes the potentially harmful effects of light saturation. Using X-ray diffraction we report a conformational change that occurs within the cytoplasmic domain of this RC in response to prolonged illumination with bright light. Our observations suggest a novel structural mechanism for the regulation of electron transfer reactions in photosynthesis.

Proteins, chlorophylls and lipids: X-ray analysis of a three-way relationship

Photosynthetic reaction centres and light harvesting complexes have been at the forefront of crystallographic studies of integral membrane proteins. In recent years, there have been spectacular advances in our understanding of the structure of (bacterio)chlorophyll-containing membrane proteins from oxygenic and anoxygenic phototrophs. In these complex structures, the protein scaffold encases different combinations of cofactors and interacts with several tightly bound lipid species that play a variety of hitherto unrecognized structural roles. Some of these lipids have relevance to the physiological function of the protein, whereas others are important for the formation of highly ordered crystals. The first site-directed mutagenesis studies of individual lipid binding sites have now underlined the importance of the lipid component for the structural stability of protein-cofactor-lipid complexes.

Electron transfer pathways and protein response to charge separation in photosynthetic reaction centers: time-resolved high-field ENDOR of the spin-correlated radical pair P865(+)QA(-)

Recently we reported the first observation of time-resolved (TR) high-frequency (HF) electron nuclear double resonance (ENDOR) of the transient charge separated state P865(+)Q(-)A in purple photosynthetic bacterial reaction centers (RC) (Poluektov, O. G., et al. J. Am. Chem. Soc. 2004, 126, 1644-1645). The high resolution and orientational selectivity of HF ENDOR allows us to directly probe protein environments by spectrally selecting specific nuclei in isotopically labeled samples. A new phenomenon associated with the spin correlated radical pair (SCRP) nature of P865(+)Q(-)A was observed. The TR-HF ENDOR spectra of protein nuclei (protons) surrounding deuterated QA(-) exhibit a derivative-like, complicated line shape, which differs considerably from the HF ENDOR spectrum of the protein nuclei surrounding thermally equilibrated QA(-). Here, a theoretical analysis of these observations is presented that shows that the positions and amplitudes of ENDOR lines contain information on hyperfine interactions (HFI) of a particular nucleus (a proton of the protein) with both correlated electron spins. Thus, spin density delocalization in the protein environment between the SCRP donor and acceptor molecules can be revealed via HF ENDOR. Novel approaches for acquiring and analyzing SCRP ENDOR that simplify interpretation of the spectra are discussed. Furthermore, we report here that the positions of the ENDOR lines of the SCRP shift with an increase in the time after laser flash, which initiates electron transfer. These shifts provide direct spectroscopic evidence of reorganization of the protein environment to accommodate the donor-acceptor charge-separated state P865(+)QA(-).

Investigation of ubiquinol formation in isolated photosynthetic reaction centers by rapid-scan Fourier transform IR spectroscopy

Light-induced formation of ubiquinol-10 in Rhodobacter sphaeroides reaction centers was followed by rapid-scan Fourier transform IR difference spectroscopy, a technique that allows the course of the reaction to be monitored, providing simultaneously information on the redox states of cofactors and on protein response. The spectrum recorded between 4 and 29 ms after the second flash showed bands at 1,470 and 1,707 cm(-1), possibly due to a QH(-) intermediate state. Spectra recorded at longer delay times showed a different shape, with bands at 1,388 (+) and 1,433 (+) cm(-1) characteristic of ubiquinol. These spectra reflect the location of the ubiquinol molecule outside the Q(B) binding site. This was confirmed by Fourier transform IR difference spectra recorded during and after continuous illumination in the presence of an excess of exogenous ubiquinone molecules, which revealed the process of ubiquinol formation, of ubiquinone/ubiquinol exchange at the Q(B) site and between detergent micelles, and of Q(B)(-) and QH(2) reoxidation by external redox mediators. Kinetics analysis of the IR bands allowed us to estimate the ubiquinone/ubiquinol exchange rate between detergent micelles to approximately 1 s. The reoxidation rate of Q(B)(-) by external donors was found to be much lower than that of QH(2), most probably reflecting a stabilizing/protecting effect of the protein for the semiquinone form. A transient band at 1,707 cm(-1) observed in the first scan (4-29 ms) after both the first and the second flash possibly reflects transient protonation of the side chain of a carboxylic amino acid involved in proton transfer from the cytoplasm towards the Q(B) site.

The First Picoseconds in Bacterial Photosynthesis?Ultrafast Electron Transfer for the Efficient Conversion of Light Energy

In this Minireview, we describe the function of the bacterial reaction centre (RC) as the central photosynthetic energy-conversion unit by ultrafast spectroscopy combined with structural analysis, site-directed mutagenesis, pigment exchange and theoretical modelling. We show that primary energy conversion is a stepwise process in which an electron is transferred via neighbouring chromophores of the RC. A well-defined chromophore arrangement in a rigid protein matrix, combined with optimised energetics of the different electron carriers, allows a highly efficient charge-separation process. The individual molecular reactions at room temperature are well described by conventional electron-transfer theory.

Excitation Energy Migration in A Dodecameric Porphyrin Wheel

Intramolecular excitation energy hopping (EEH) time within a dodecameric porphyrin wheel C6ZA, in which six meso-meso linked zinc(II) diporphyrin (Z2) subunits are bridged by 1,3-phenylene spacers, is deduced by a Förster energy hopping model based on S(1)-S(1) exciton-exciton annihilation and anisotropy depolarization. Under the assumption that the energy hopping sites are six Z2 subunits, two different observables (e.g., exciton-exciton annihilation and anisotropy depolarization times) consistently give the EEH time of 4.0 +/- 0.4 ps via 1,3-phenylene spacer of C6ZA, which is faster than 9.4 ps of linear 2Z2 (1,3-phenylene-linked zinc(II) tetraporphyrin). As a consequence, C6ZA serves as a well-defined two-dimensional model for a light-harvesting complex.

Investigation of B-branch electron transfer by femtosecond time resolved spectroscopy in a Rhodobacter sphaeroides reaction centre that lacks the QA ubiquinone

The dynamics of electron transfer in a membrane-bound Rhodobacter sphaeroides reaction centre containing a combination of four mutations were investigated by transient absorption spectroscopy. The reaction centre, named WAAH, has a mutation that causes the reaction centre to assemble without a Q(A) ubiquinone (Ala M260 to Trp), a mutation that causes the replacement of the H(A) bacteriopheophytin with a bacteriochlorophyll (Leu M214 to His) and two mutations that remove acidic groups close to the Q(B) ubiquinone (Glu L212 to Ala and Asp L213 to Ala). Previous work has shown that the Q(B) ubiquinone is reduced by electron transfer along the so-called inactive cofactor branch (B-branch) in the WAAH reaction centre (M.C. Wakeham, M.G. Goodwin, C. McKibbin, M.R. Jones, Photo-accumulation of the P(+)Q(B)(-) radical pair state in purple bacterial reaction centres that lack the Q(A) ubiquinone, FEBS Letters 540 (2003) 234-240). In the present study the dynamics of electron transfer in the membrane-bound WAAH reaction centre were studied by femtosecond transient absorption spectroscopy, and the data analysed using a compartmental model. The analysis indicates that the yield of Q(B) reduction via the B-branch is approximately 8% in the WAAH reaction centre, consistent with results from millisecond time-scale kinetic spectroscopy. Possible contributions to this yield of the constituent mutations in the WAAH reaction centre and the membrane environment of the complex are discussed.

Studies of proton translocations in biological systems: simulating proton transport in carbonic anhydrase by EVB-based models

Proton transport (PTR) processes play a major role in bioenergetics and thus it is important to gain a molecular understanding of these processes. At present the detailed description of PTR in proteins is somewhat unclear and it is important to examine different models by using well-defined experimental systems. One of the best benchmarks is provided by carbonic anhydrase III (CA III), because this is one of the few systems where we have a clear molecular knowledge of the rate constant of the PTR process and its variation upon mutations. Furthermore, this system transfers a proton between several water molecules, thus making it highly relevant to a careful examination of the "proton wire" concept. Obtaining a correlation between the structure of this protein and the rate of the PTR process should help to discriminate between alternative models and to give useful clues about PTR processes in other systems. Obviously, obtaining such a correlation requires a correct representation of the "chemistry" of PTR between different donors and acceptors, as well as the ability to evaluate the free energy barriers of charge transfer in proteins, and to simulate long-time kinetic processes. The microscopic empirical valence bond (Warshel, A., and R. M. Weiss. 1980. J. Am. Chem. Soc. 102:6218-6226; and Aqvist, J., and A. Warshel. 1993. Chem. Rev. 93:2523-2544) provides a powerful way for representing the chemistry and evaluating the free energy barriers, but it cannot be used with the currently available computer times in direct simulation of PTR with significant activation barriers. Alternatively, one can reduce the empirical valence bond (EVB) to the modified Marcus' relationship and use semimacroscopic electrostatic calculations plus a master equation to determine the PTR kinetics (Sham, Y., I. Muegge, and A. Warshel. 1999. Proteins. 36:484-500). However, such an approximation does not provide a rigorous multisite kinetic treatment. Here we combine the useful ingredients of both approaches and develop a simplified EVB effective potential that treats explicitly the chain of donors and acceptors while considering implicitly the rest of the protein/solvent system. This approach can be used in Langevin dynamics simulations of long-time PTR processes. The validity of our new simplified approach is demonstrated first by comparing its Langevin dynamics results for a PTR along a chain of water molecules in water to the corresponding molecular dynamics simulations of the fully microscopic EVB model. This study examines dynamics of both models in cases of low activation barriers and the dependence of the rate on the energetics for cases with moderate barriers. The study of the dependence on the activation barrier is next extended to the range of higher barriers, demonstrating a clear correlation between the barrier height and the rate constant. The simplified EVB model is then examined in studies of the PTR in carbonic anhydrase III, where it reproduces the relevant experimental results without the use of any parameter that is specifically adjusted to fit the energetics or dynamics of the reaction in the protein. We also validate the conclusions obtained previously from the EVB-based modified Marcus' relationship. It is concluded that this approach and the EVB-based model provide a reliable, effective, and general tool for studies of PTR in proteins. Finally in view of the behavior of the simulated result, in both water and the CA III, we conclude that the rate of PTR in proteins is determined by the electrostatic energy of the transferred proton as long as this energy is higher than a few kcal/mol.

Resonance Raman characterization of Rhodobacter capsulatus reaction centers with lysine mutations near the accessory bacteriochlorophylls

Lysine residues have been introduced into Rhodobacter capsulatus reaction centers at M-polypeptide position 201 and at L-polypeptide position 178. These positions are in the proximity of ring V of the accessory bacterochlorophylls BA and BB, respectively. Resonance Raman studies indicate that the introduction of a Lys residue at either position M201 or L178 results in structural perturbations to the BChl cofactors. Lys at L178 directly interacts with BB, most likely via a hydrogen bond. The hydrogen bonding interaction is consistent with enhanced B branch electron transfer that is observed in RCs from the S(L178)K/G(M201)D/L(M212)H triple mutant versus the G(M201)D/L(M212)H double mutant. In contrast, the introduction of a Lys at M201 does not result in hydrogen bonding to the BA cofactor, in contrast to the introduction of a His at M201. Accordingly, the alkyl ammonium head group of the side chain of the Lys at M201 residue appears to be distant from BA.

Realistic Simulations of Proton Transport along the Gramicidin Channel: Demonstrating the Importance of Solvation Effects

The nature of proton transduction (PTR) through a file of water molecules, along the gramicidin A (gA) channel, has long been considered as being highly relevant to PTR in biological systems. Previous attempts to model this process implied that the so-called Grotthuss mechanism and the corresponding orientation of the water file plays a major role. The present work reexamines the PTR in gA by combining a fully microscopic empirical valence bond (EVB) model and a recently developed simplified EVB-based model with Langevin dynamics (LD) simulations. The full model is used first to evaluate the free energy profile for a stepwise PTR process. The corresponding results are then used to construct the effective potential of the simplified EVB. This later model is then used in Langevin dynamics simulations, taking into account the correct physics of possible concerted motions and the effect of the solvent reorganization. The simulations reproduce the observed experimental trend and lead to a picture that is quite different from that assumed previously. It is found that the PTR in gA is controlled by the change in solvation energy of the transferred proton along the channel axis. Although the time dependent electrostatic fluctuations of the channel and water dipoles play their usual role in modulating the proton-transfer process (Proc. Natl. Acad. Sci. U.S.A. 1984, 81, 444), the PTR rate is mainly determined by the free energy profile. Furthermore, the energetics of the reorientation of the unprotonated water file do not appear to provide a consistent way of assessing the activation barrier for the PTR process. It seems to us that in the case of gA, and probably other systems with significant electrostatic barriers for the transfer of the proton charge, the PTR rate is controlled by the electrostatic barrier. This finding has clear consequences with regards to PTR processes in biological systems.

Structural and functional studies on the tetraheme cytochrome subunit and its electron donor proteins: the possible docking mechanisms during the electron transfer reaction

The photosynthetic reaction centers (RCs) classified as the group II possess a peripheral cytochrome (Cyt) subunit, which serves as the electron mediator to the special-pair. In the cycle of the photosynthetic electron transfer reactions, the Cyt subunit accepts electrons from soluble electron carrier proteins, and re-reduces the photo-oxidized special-pair of the bacteriochlorophyll. Physiologically, high-potential cytochromes such as the cytochrome c2 and the high-potential iron-sulfur protein (HiPIP) function as the electron donors to the Cyt subunit. Most of the Cyt subunits possess four heme c groups, and it was unclear which heme group first accepts the electron from the electron donor. The most distal heme to the special-pair, the heme-1, has a lower redox potential than the electron donors, which makes it difficult to understand the electron transfer mechanism mediated by the Cyt subunit. Extensive mutagenesis combined with kinetic studies has made a great contribution to our understanding of the molecular interaction mechanisms, and has demonstrated the importance of the region close to the heme-1 in the electron transfer. Moreover, crystallographic studies have elucidated two high-resolution three-dimensional structures for the RCs containing the Cyt subunit, the Blastochloris viridis and Thermochromatium tepidum RCs, as well as the structures of their electron donors. An examination of the structural data also suggested that the binding sites for both the cytochrome c2 and the HiPIP are located adjacent to the solvent-accessible edge of the heme-1. In addition, it is also indicated by the structural and biochemical data that the cytochrome c2 and the HiPIP dock with the Cyt subunit by c2 is recognized through electrostatic interactions while hydrophobic interactions are important in the HiPIP docking.

The structure and function of the cytochrome c2: reaction center electron transfer complex from Rhodobacter sphaeroides

In the photosynthetic bacterium, Rhodobacter sphaeroides, the mobile electron carrier, cytochrome c2 (cyt c2) transfers an electron from reduced heme to the photooxidized bacteriochlorophyll dimer in the membrane bound reaction center (RC) as part of the light induced cyclic electron transfer chain. A complex between these two proteins that is active in electron transfer has been crystallized and its structure determined by X-ray diffraction. The structure of the cyt:RC complex shows the cyt c2 (cyt c2) positioned at the center of the periplasmic surface of the RC. The exposed heme edge from cyt c2 is in close tunneling contact with the electron acceptor through an intervening bridging residue, Tyr L162 located on the RC surface directly above the bacteriochlorophyll dimer. The binding interface between the two proteins can be divided into two regions: a short-range interaction domain and a long-range interaction domain. The short-range domain includes residues immediately surrounding the tunneling contact region around the heme and Tyr L162 that display close intermolecular contacts optimized for electron transfer. These include a small number of hydrophobic interactions, hydrogen bonds and a pi-cation interaction. The long-range interaction domain consists of solvated complementary charged residues; positively charged residues from the cyt and negatively charged residues from the RC that provide long range electrostatic interactions that can steer the two proteins into position for rapid association.