Dynamics of the glutamic acid 242 side chain in cytochrome c oxidase

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a₃/Cu_B site, and for proton pumping. When present, the side chain of Glu-242 is orientated “down” towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity “above” Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an “up” position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the “up” position is preferred energetically when the cavity contains four water molecules, but the “down” position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the “resting” state of the enzyme, and that the “up/down” equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O₂ reduction.     

Towards the mechanism of proton pumping by the haem-copper oxidases

The haem-copper oxidases comprise a large family of enzymes that is widespread among aerobic organisms. These remarkable membrane-bound proteins catalyse the respiratory reduction of dioxygen to water, and conserve free energy from this reaction by operating as proton pumps. The mechanism of redox-dependent proton translocation has been elusive despite the availability of high resolution crystal structures from several oxidases. Here, we discuss some recent as well as some older results that may shed light on this mechanism. We conclude that proton-pumping is initiated by vectorial proton transfer from a conserved glutamic acid (Glu242 in the bovine enzyme) to a proton acceptor above the haem groups, and that this primary event is mechanistically coupled to electron transfer from haem a to the binuclear haem a₃/Cu_B centre. Subsequently, Glu242 is reprotonated from the negatively charged side of the membrane. Next this proton is transferred to the binuclear site to complete the chemistry, Glu242 is reprotonated once more, and the “prepumped” proton is ejected on the opposite side of the membrane. The different kinetics of electron-coupled proton transfer in different steps of the catalytic cycle may be related to differences in the driving force due to different E_m values of the electron acceptor in the binuclear site.     

Redox-coupled proton pumping in cytochrome c oxidase: Further insights from computer simulation

The membrane-bound enzyme cytochrome c oxidase, the terminal member in the respiratory chain, converts oxygen into water and generates an electrochemical gradient by coupling the electron transfer to proton pumping across the membrane. Here we have investigated the dynamics of an excess proton and the surrounding protein environment near the active sites. The multi-state empirical valence bond (MS-EVB) molecular dynamics method was used to simulate the explicit dynamics of proton transfer through the critically important hydrophobic channel between Glu242 (bovine notation) and the D-propionate of heme a₃ (PRDa₃) for the first time. The results from these molecular dynamics simulations indicate that the PRDa₃ can indeed re-orientate and dissociate from Arg438, despite the high stability of such an ion pair, and has the ability to accept protons via bound water molecules. Any large conformational change of the adjacent heme a D-propionate group is, however, sterically blocked directly by the protein. Free energy calculations of the PRDa₃ side chain isomerization and the proton translocation between Glu242 and the PRDa₃ site have also been performed. The results exhibit a redox state-dependent dynamical behavior and indicate that reduction of the low-spin heme a may initiate internal transfer of the pumped proton from Glu242 to the PRDa₃ site.     

Limited reversibility of transmembrane proton transfer assisting transmembrane electron transfer in a dihaem-containing succinate:quinone oxidoreductase

Membrane protein complexes can support both the generation and utilisation of a transmembrane electrochemical proton potential (Δp), either by supporting transmembrane electron transfer coupled to protolytic reactions on opposite sides of the membrane or by supporting transmembrane proton transfer. The first mechanism has been unequivocally demonstrated to be operational for Δp-dependent catalysis of succinate oxidation by quinone in the case of the dihaem-containing succinate:menaquinone reductase (SQR) from the Gram-positive bacterium Bacillus licheniformis. This is physiologically relevant in that it allows the transmembrane potential Δp to drive the endergonic oxidation of succinate by menaquinone by the dihaem-containing SQR of Gram-positive bacteria. In the case of a related but different respiratory membrane protein complex, the dihaem-containing quinol:fumarate reductase (QFR) of the ?-proteobacterium Wolinella succinogenes, evidence has been obtained that both mechanisms are combined, so as to facilitate transmembrane electron transfer by proton transfer via a both novel and essential compensatory transmembrane proton transfer pathway (“E-pathway”). Although the reduction of fumarate by menaquinol is exergonic, it is obviously not exergonic enough to support the generation of a Δp. This compensatory “E-pathway” appears to be required by all dihaem-containing QFR enzymes and results in the overall reaction being electroneutral. However, here we show that the reverse reaction, the oxidation of succinate by quinone, as catalysed by W. succinogenes QFR, is not electroneutral. The implications for transmembrane proton transfer via the E-pathway are discussed.     

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

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

Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump which converts atmospheric oxygen to water and couples the oxygen reduction reaction to the creation of a membrane proton gradient. The structure of the enzyme has been solved; however, the mechanism of proton pumping is still poorly understood. Recent calculations from this group indicate that one of the histidine ligands of enzyme's Cu_B center, His291, may play the role of the pumping element. In this paper, we report on the results of calculations that combined first principles DFT and continuum electrostatics to evaluate the energetics of the key energy generating step of the model—the transfer of the chemical proton to the binuclear center of the enzyme, where the hydroxyl group is converted to water, and the concerted expulsion of the proton from δ-nitrogen of His291 ligand of Cu_B center. We show that the energy generated in this step is sufficient to push a proton against an electrochemical membrane gradient of about 200 mV. We have also re-calculated the pK_a of His291 for an extended model in which the whole Fe_a3-Cu_B center with their ligands is treated by DFT. Two different DFT functionals (B3LYP and PBE0), and various dielectric models of the protein have been used in an attempt to estimate potential errors of the calculations. Although current methods of calculations do not allow unambiguous predictions of energetics in proteins within few pK_a units, as required in this case, the present calculation provides further support for the proposed His291 model of CcO pump and makes a specific prediction that could be targeted in the experimental test.     

Theoretical and computational analysis of the membrane potential generated by cytochrome c oxidase upon single electron injection into the enzyme

We have developed theory and the computational scheme for the analysis of the kinetics of the membrane potential generated by cytochrome c oxidase upon single electron injection into the enzyme. The theory allows one to connect the charge motions inside the enzyme to the membrane potential observed in the experiments by using data from the “dielectric topography” map of the enzyme that we have created. The developed theory is applied for the analysis of the potentiometric data recently reported by the Wikström group [I. Belevich, D.A. Bloch, N. Belevich, M. Wikström and M.I. Verkhovsky, Exploring the proton pump mechanism of cytochrome c oxidase in real time, Proc. Natl. Acad. Sci. U. S. A. 104 (2007) 2685-2690] on the O to E transition in Paracoccus denitrificans oxidase. Our analysis suggests, that the electron transfer to the binuclear center is coupled to a proton transfer (proton loading) to a group just “above” the binuclear center of the enzyme, from which the pumped proton is subsequently expelled by the chemical proton arriving to the binuclear center. The identity of the pump site could not be determined with certainty, but could be localized to the group of residues His326 (His291 in bovine), propionates of heme a₃, Arg 473/474, and Trp164. The analysis also suggests that the dielectric distance from the P-side to Fe a is 0.4 or larger. The difficulties and pitfalls of quantitative interpretation of potentiometric data are discussed.     

Effect of Mutations of Arginine 94 on Proton Pumping,Electron Transfer,and Superoxide Anion Generation in Cytochrome b of the bc1 Complex from Rhodobacter sphaeroides

Proton transfer involving internal water molecules that provide hydrogen bonds and facilitate proton diffusion has been identified in some membrane proteins. Arg-94 in cytochrome b of the Rhodobacter sphaeroides bc₁ complex is fully conserved and is hydrogen-bonded to the heme propionate and a chain of water molecules. To further elucidate the role of Arg-94, we generated the mutations R94A, R94D, and R94N. The wild-type and mutant bc₁ complexes were purified and then characterized. The results show that substitution of Arg-94 decreased electron transfer activity and proton pumping capability and increased O₂^˙̄ production, suggesting the importance of Arg-94 in the catalytic mechanism of the bc₁ complex in R. sphaeroides. This also suggests that the transport of H⁺, O₂, and O₂^˙̄ in the bc₁ complex may occur by the same pathway.     

Subunit a of the F1F0 ATP Synthase Requires YidC and SecYEG for Membrane Insertion

The insertion of inner membrane proteins in Escherichia coli occurs almost exclusively via the SecYEG pathway, while some membrane proteins require the membrane protein insertase YidC. In vitro analysis demonstrates that subunit a of the F₁F₀ ATP synthase (F₀a) is strictly dependent on Ffh, SecYEG and YidC for its membrane insertion but independent of the proton motive force. The insertion of the first transmembrane segment of F₀a also depends on Ffh and SecYEG but not on YidC, whereas the insertion is strongly dependent on the proton motive force, unlike the full-length F₀a protein. These data demonstrate an extensive role of YidC in the assembly of the F₀ sector of the F₁F₀ ATP synthase.     

Essential arginine in subunit a and aspartate in subunit c of FoF1 ATP synthase: Effect of repositioning within Helix 4 of subunit a and Helix 2 of subunit c

F_oF₁ ATP synthase couples proton flow through the integral membrane portion F_o (ab₂c₁₀) to ATP-synthesis in the extrinsic F₁-part ((αβ)₃γδε) (Escherichia coli nomenclature and stoichiometry). Coupling occurs by mechanical rotation of subunits c₁₀γε relative to (αβ)₃δab₂. Two residues were found to be essential for proton flow through ab₂c_10, namely Arg₂₁₀ in subunit a (aR210) and Asp₆₁ in subunits c (cD61). Their deletion abolishes proton flow, but “horizontal” repositioning, by anchoring them in adjacent transmembrane helices, restores function. Here, we investigated the effects of “vertical” repositioning aR210, cD61, or both by one helical turn towards the N- or C-termini of their original helices. Other than in the horizontal the vertical displacement changes the positions of the side chains within the depth of the membrane. Mutant aR210A/aN214R appeared to be short-circuited in that it supported proton conduction only through EF₁-depleted EF_o, but not in EF_oEF_1, nor ATP-driven proton pumping. Mutant cD61N/cM65D grew on succinate, retained the ability to synthesize ATP and supported passive proton conduction but apparently not ATP hydrolysis-driven proton pumping.     

pK(a) values for the unfolded state under native conditions explain the pH-dependent stability of PGB1

Understanding the role of electrostatics in protein stability requires knowledge of these interactions in both the folded and unfolded states. Electrostatic interactions can be probed experimentally by characterizing ionization equilibria of titrating groups, parameterized as pK_a values. However, pK_a values of the unfolded state are rarely accessible under native conditions, where the unfolded state has a very low population. Here, we report pK_a values under nondenaturing conditions for two unfolded fragments of the protein G B1 domain that mimic the unfolded state of the intact protein. pK_a values were determined for carboxyl groups by monitoring their pH-dependent ¹³C chemical shifts. Monte Carlo simulations using a Gaussian chain model provide corrections for changes in electrostatic interactions that arise from fragmentation of the protein. Most pK_a values for the unfolded state agree well with model values, but some residues show significant perturbations that can be rationalized by local electrostatic interactions. The pH-dependent stability was calculated from the experimental pK_a values of the folded and unfolded states and compared to experimental stability data. The use of experimental pK_a values for the unfolded state results in significantly improved agreement with experimental data, as compared to calculations based on model data alone.     

Electron transfer patterns of the di-heme protein cytochrome c4 from Pseudomonas stutzeri

We report kinetic data for the two-step electron transfer (ET) oxidation and reduction of the two-domain di-heme redox protein Pseudomonas stutzeri cytochrome (cyt) c₄ by [Co(bipy)₃]^2+/3+ (bipy = 2,2′-bipyridine). Following earlier reports, the data accord with both bi- and tri-exponential kinetics. A complete kinetic scheme includes both “cooperative” intermolecular ET between each heme group and the external reaction partner, and intramolecular ET between the two heme groups. A new data analysis scheme shows unequivocally that two-ET oxidation and reduction of P. stutzeri cyt c₄ is entirely dominated by intermolecular ET between the heme groups and the external reaction partner in the ms time range, with virtually no contribution from intramolecular interheme ET in this time range. This is in striking contrast to two-ET electrochemical oxidation or reduction of P. stutzeri cyt c₄ for which fast, ms to sub-ms intramolecular interheme ET is a crucial step. The rate constant dependence on the solvent viscosity has disclosed strong coupling to both a (set of) frictionally damped solvent/protein nuclear modes and intramolecular friction-less “ballistic” modes, indicative of notable protein structural mobility in the overall two-ET process. We suggest that conformational protein mobility blocks intramolecular interheme ET in bulk homogeneous solution but triggers opening of this gated ET channel in the electrochemical environment or in the membrane environment of natural respiratory cyt c₄ function.     

Effects of bleaching on the pharmacological and toxicological activities elicited by the aqueous extracts prepared from two “fire corals” collected in the Mexican Caribbean

Hydrocorals of the genus Millepora are abundant skeleton-forming inhabitants of coral reefs around the world. These species are popularly known as “fire corals” since contact with them causes severe pain, skin eruptions and blisters as a result of the release of unidentified toxins. Millepora species associate with photosynthetic dinoflagellates of the genus Symbiodinium (“zooxanthellae”), and up to now the role of these symbionts in the toxic effects induced by the “fire corals” is unknown. In this study, we compared the hemolytic, vasoconstrictor, and phospholipase A₂ (PLA₂) activities of the crude aqueous extracts prepared from normal and bleached specimens of two hydrocorals collected in the Mexican Caribbean, Milleporaalcicornis and Millepora complanata. Electrophoretic analysis revealed some differences between the protein profiles of the extracts prepared from normal and bleached specimens. Bleaching decreased, but not abolished, the hemolytic effect induced by the hydrocorals extracts and the phospholipase A₂ activity of M. complanata extract. Furthermore, it did not modify the enzymatic activity of M. alcicornis extract and vasoconstriction elicited by both extracts. Our results suggest that the presence of the symbionts does not importantly influence the pharmacological and toxic effects induced by Millepora ssp. extracts, and indicate that cnidarians are the main source of the bioactive compounds.     

Nucleophile activation in PD…(D/E)xK metallonucleases: An experimental and computational pKa study

Metallonucleases conduct metal-dependent nucleic acid hydrolysis. While metal ions serve in multiple mechanistic capacities in these enzymes, precisely how the attacking water is activated remains unclear for those lacking an obvious general base. All arguments hinge on appropriate pK_as for active site moieties very close to this species, and measurement of the pK_a of a specific water molecule is difficult to access experimentally. Here we describe a computational approach for exploring the local electrostatic influences on the water-derived nucleophile in metallonucleases featuring the common PD…(D/E)xK motif. We utilized UHBD to predict the pK_as of active site groups, including that of a water molecule positioned to act as a nucleophile. The pK_a of a Mg(II)-ligated water molecule hydrogen bonded to the conserved Lys70 in a Mg(II)-PvuII enzyme complex was calculated to be 6.5. The metal and the charge on the Lys group were removed in separate experiments; both resulted in the elevation of the pK_a of this water molecule, consistent with contributions from both moieties to lowering this pK_a. This behavior is preserved among other PD…(D/E)xK metallonucleases. pK_as extracted from the pH dependence of the single turnover rate constant are compared to previous experimental data and the above predicted pK_as.     

Prediction of Secondary Ionization of the Phosphate Group in Phosphotyrosine Peptides

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK_a values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK_a calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK_a values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK_a of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein—phospholipase. The agreement between the computed and the experimental pK_a values is better than ±0.3 pH units for the optimized solute dielectric constant of 11–13. This approach is promising and its extension to other phospho-amino acids is in progress.     

Crystal Structure of Glyceraldehyde-3-Phosphate Dehydrogenase 1 from Methicillin-Resistant Staphylococcus aureus MRSA252 Provides Novel Insights into Substrate Binding and Catalytic Mechanism

The dreaded pathogen Staphylococcus aureus is one of the causes of morbidity and mortality worldwide. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), one of the key glycolytic enzymes, is irreversibly oxidized under oxidative stress and is responsible for sustenance of the pathogen inside the host. With an aim to elucidate the catalytic mechanism and identification of intermediates involved, we describe in this study different crystal structures of GAPDH1 from methicillin-resistant S. aureus MRSA252 (SaGAPDH1) in apo and holo forms of wild type, thioacyl intermediate, and ternary complexes of active-site mutants with physiological substrate d-glyceraldehyde-3-phosphate (G3P) and coenzyme NAD⁺. A new phosphate recognition site, “new P_i” site, similar to that observed in GAPDH from Thermotoga maritima, is reported here, which is 3.40 Å away from the “classical P_i” site. Ternary complexes discussed are representatives of noncovalent Michaelis complexes in the ground state. d-G3P is bound to all the four subunits of C151S.NAD and C151G.NAD in more reactive hydrate (gem-di-ol) form. However, in C151S + H178N.NAD, the substrate is bound to two chains in aldehyde form and in gem-di-ol form to the other two. This work reports binding of d-G3P to the C151G mutant in an inverted manner for the very first time. The structure of the thiaocyl complex presented here is formed after the hydride transfer. The C3 phosphate of d-G3P is positioned at the “P_s” site in the ternary complexes but at the “new P_i” site in the thioacyl complex and C1-O1 bond points opposite to His178 disrupting the alignment between itself and NE2 of His178. A new conformation (Conformation I) of the 209-215 loop has also been identified, where the interaction between phosphate ion at the “new P_i” site and conserved Gly212 is lost. Altogether, inferences drawn from the kinetic analyses and crystal structures suggest the “flip-flop” model proposed for the enzyme mechanism.     

Cytochrome c1 exhibits two binding sites for cytochrome c in plants

In plants, channeling of cytochrome c molecules between complexes III and IV has been purported to shuttle electrons within the supercomplexes instead of carrying electrons by random diffusion across the intermembrane bulk phase. However, the mode plant cytochrome c behaves inside a supercomplex such as the respirasome, formed by complexes I, III and IV, remains obscure from a structural point of view. Here, we report ab-initio Brownian dynamics calculations and nuclear magnetic resonance-driven docking computations showing two binding sites for plant cytochrome c at the head soluble domain of plant cytochrome c₁, namely a non-productive (or distal) site with a long heme-to-heme distance and a functional (or proximal) site with the two heme groups close enough as to allow electron transfer. As inferred from isothermal titration calorimetry experiments, the two binding sites exhibit different equilibrium dissociation constants, for both reduced and oxidized species, that are all within the micromolar range, thus revealing the transient nature of such a respiratory complex. Although the docking of cytochrome c at the distal site occurs at the interface between cytochrome c₁ and the Rieske subunit, it is fully compatible with the complex III structure. In our model, the extra distal site in complex III could indeed facilitate the functional cytochrome c channeling towards complex IV by building a “floating boat bridge” of cytochrome c molecules (between complexes III and IV) in plant respirasome.     

Formamide probes a role for water in the catalytic cycle of cytochrome c oxidase

Formamide is a slow-onset inhibitor of mitochondrial cytochrome c oxidase that is proposed to act by blocking water movement through the protein. In the presence of formamide the redox level of mitochondrial cytochrome c oxidase evolves over the steady state as the apparent electron transfer rate from cytochrome a to cytochrome a₃ slows. At maximal inhibition cytochrome a and cytochrome c are fully reduced, whereas cytochrome a₃ and Cu_B remain fully oxidized consistent with the idea that formamide interferes with electron transfer between cytochrome a and the oxygen reaction site. However, transient kinetic studies show that intrinsic rates of electron transfer are unchanged in the formamide-inhibited enzyme. Formamide inhibition is demonstrated for another member of the heme-oxidase family, cytochrome c oxidase from Bacillus subtilis, but the onset of inhibition is much quicker than for mitochondrial oxidase. If formamide inhibition arises from a steric blockade of water exchange during catalysis then water exchange in the smaller bacterial oxidase is more open. Subunit III removal from the mitochondrial oxidase hastens the onset of formamide inhibition suggesting a role for subunit III in controlling water exchange during the cytochrome c oxidase reaction.     

Energy and electron transfer in the photosynthetic reaction center complex of Acidiphilium rubrum containing Zn-bacteriochlorophyll a studied by femtosecond up-conversion spectroscopy

A photosynthetic reaction center (RC) complex was isolated from a purple bacterium, Acidiphilium rubrum. The RC contains bacteriochlorophyll a containing Zn as a central metal (Zn-BChl a) and bacteriopheophytin a (BPhe a) but no Mg-BChl a. The absorption peaks of the Zn-BChl a dimer (P_Zn), the accessory Zn-BChl a (B_Zn), and BPhe a (H) at 4 K in the RC showed peaks at 875, 792, and 753 nm, respectively. These peaks were shorter than the corresponding peaks in Rhodobacter sphaeroides RC that has Mg-BChl a. The kinetics of fluorescence from P_Zn^*, measured by fluorescence up-conversion, showed the rise and the major decay with time constants of 0.16 and 3.3 ps, respectively. The former represents the energy transfer from B_Zn^* to P_Zn, and the latter, the electron transfer from P_Zn to H. The angle between the transition dipoles of B_Zn and P_Zn was estimated to be 36° based on the fluorescence anisotropy. The time constants and the angle are almost equal to those in the Rb. sphaeroides RC. The high efficiency of A. rubrum RC seems to be enabled by the chemical property of Zn-BChl a and by the L168HE modification of the RC protein that modifies P_Zn.     

Sequential assembly of photosynthetic units in Rhodobacter sphaeroides as revealed by fast repetition rate analysis of variable bacteriochlorophyll a fluorescence

The development of functional photosynthetic units in Rhodobacter sphaeroides was followed by near infra-red fast repetition rate (IRFRR) fluorescence measurements that were correlated to absorption spectroscopy, electron microscopy and pigment analyses. To induce the formation of intracytoplasmic membranes (ICM) (greening), cells grown aerobically both in batch culture and in a carbon-limited chemostat were transferred to semiaerobic conditions. In both aerobic cultures, a low level of photosynthetic complexes was observed, which were composed of the reaction center and the LH1 core antenna. Interestingly, in the batch cultures the reaction centers were essentially inactive in forward electron transfer and exhibited low photochemical yields F_V/F_M, whereas the chemostat culture displayed functional reaction centers with a rather rapid (1-2 ms) electron transfer turnover, as well as a high F_V/F_M of ∼0.8. In both cases, the transfer to semiaerobiosis resulted in rapid induction of bacteriochlorophyll a synthesis that was reflected by both an increase in the number of LH1-reaction center and peripheral LH2 antenna complexes. These studies establish that photosynthetic units are assembled in a sequential manner, where the appearance of the LH1-reaction center cores is followed by the activation of functional electron transfer, and finally by the accumulation of the LH2 complexes.