Potentiation of proton transfer function by electrostatic interactions in photosynthetic reaction centers from Rhodobacter sphaeroides: First results from site-directed mutation of the H subunit.

The x-ray crystallographic structure of the photosynthetic reaction center (RC) has proven critical in understanding biological electron transfer processes. By contrast, understanding of intraprotein proton transfer is easily lost in the immense richness of the details. In the RC of Rhodobacter (Rb.) sphaeroides, the secondary quinone (QB) is surrounded by amino acid residues of the L subunit and some buried water molecules, with M- and H-subunit residues also close by. The effects of site-directed mutagenesis upon RC turnover and quinone function have implicated several L-subunit residues in proton delivery to QB, although some species differences exist. In wild-type Rb. sphaeroides, Glu L212 and Asp L213 represent an inner shell of residues of particular importance in proton transfer to QB. Asp L213 is crucial for delivery of the first proton, coupled to transfer of the second electron, while Glu L212, possibly together with Asp L213, is necessary for delivery of the second proton, after the second electron transfer. We report here the first study, by site-directed mutagenesis, of the role of the H subunit in QB function. Glu H173, one of a cluster of strongly interacting residues near QB, including Asp L213, was altered to Gln. In isolated mutant RCs, the kinetics of the first electron transfer, leading to formation of the semiquinone, QB-, and the proton-linked second electron transfer, leading to the formation of fully reduced quinol, were both greatly retarded, as observed previously in the Asp L213 --> Asn mutant. However, the first electron transfer equilibrium, QA-QB <==> QAQB-, was decreased, which is opposite to the effect of the Asp L213 --> Asn mutation. These major disruptions of events coupled to proton delivery to QB were largely reversed by the addition of azide (N3-). The results support a major role for electrostatic interactions between charged groups in determining the protonation state of certain entities, thereby controlling the rate of the second electron transfer. It is suggested that the essential electrostatic effect may be to "potentiate" proton transfer activity by raising the pK of functional entities that actually transfer protons in a coupled fashion with the second electron transfer. Candidates include buried water (H3O+) and Ser L223 (serine-OH2+), which is very close to the O5 carbonyl of the quinone.

Isotropic and anisotropic spin-spin interactions and a quantum phase transition in a dinuclear Cu(II) compound

We report electron-paramagnetic resonance (EPR) studies at similar to 9.5 GHz (X band) and similar to 34 GHz (Q band) of powder and single-crystal samples of the compound Cu(2)[TzTs](4) [N-thiazol-2-yl-toluenesulfonamidatecopper(II)], C(40)H(36)Cu(2)N(8)O(8)S(8), having copper(II) ions in dinuclear units. Our data allow determining an antiferromagnetic interaction J(0)=(-113 +/- 1) cm(-1) (H(ex)=-J(0)S(1)center dot S(2)) between Cu(II) ions in the dinuclear unit and the anisotropic contributions to the spin-spin coupling matrix D (H(ani)=S(1)center dot D center dot S(2)), a traceless symmetric matrix with principal values D/4=(0.198 +/- 0.003) cm(-1) and E/4=(0.001 +/- 0.003) cm(-1) arising from magnetic dipole-dipole and anisotropic exchange couplings within the units. In addition, the single-crystal EPR measurements allow detecting and estimating very weak exchange couplings between neighbor dinuclear units, with an estimated magnitude parallel to J(')parallel to=(0.060 +/- 0.015) cm(-1). The interactions between a dinuclear unit and the ""environment"" of similar units in the structure of the compound produce a spin dynamics that averages out the intradinuclear dipolar interactions. This coupling with the environment leads to decoherence, a quantum phase transition that collapses the dipolar interaction when the isotropic exchange coupling with neighbor dinuclear units equals the magnitude of the intradinuclear dipolar coupling. Our EPR experiments provide a new procedure to follow the classical exchange-narrowing process as a shift and collapse of the line structure (not only as a change of the resonance width), which is described with general (but otherwise simple) theories of magnetic resonance. Using complementary procedures, our EPR measurements in powder and single-crystal samples allow measuring simultaneously three types of interactions differing by more than three orders of magnitude (between 113 cm(-1) and 0.060 cm(-1)).

The phylogenetic relationships of Caulobacter, Asticcacaulis and Brevundimonas species and their taxonomic implications

The phylogenetic relationships among the species of Caulobacter, Asticcacaulis and Brevundimonas were studied by comparison of their 16S rDNA sequences. The analysis of almost complete sequences confirmed the early evolutionary divergence of the freshwater and marine species of Caulobacter reported previously [Stahl, D. A., Key, R,, Flesher, B, & Smit, J. (1992), J Bacteriol 174, 2193-2198]. The freshwater species formed two distinct clusters. One cluster contained the species Caulobacter bacteroides, Caulobacter crescentus, Caulobacter fusiformis and Caulobacter henricii. C, bacteroides and C, fusiformis are very closely related (sequence identity 99.8%). The second cluster was not exclusive and contained the species Caulobacter intermedius, Caulobacter subvibrioides and Caulobacter variabilis, as well as Brevundimonas diminuta and Brevundimonas vesicularis, The marine species Caulobacter halobacteroides and Caulobacter maris were very closely related, with a sequence identity of 99.7%, These two species were most closely but distantly related to the marine hyphal/budding bacteria Hyphomonas jannaschiana and Hirschia baltica, which formed a deep phylogenetic line with Rhodobacter sphaeroides and Rhodobacter capsulatus, Caulobacter leidyia is unrelated to the other species of Caulobacter and belongs to the alpha-4 subclass of the Proteobacteria, forming a distinct cluster with Asticcacaulis excentricus and Asticcacaulis biprosthecium, The taxonomic implications of the polyphyletic nature of the genus Caulobacter and the absence of a type culture for the type species of the genus, Caulobacter vibrioides, are discussed.

A Sco homologue plays a role in defence against oxidative stress in pathogenic Neisseria

Sco proteins are found in mitochondria and in a variety of oxidase positive bacteria. Although Sco is required for the formation of the Cu-A centre in a cytochrome oxidase of the aa(3) type, it was observed that oxidases with a Cu-A centre are not present in many bacteria that contain a Sco homologue. Two bacteria of this type are the pathogens Neisseria meningitidis and Neisseria gonorrhoeae. The sco genes of N. gonorrhoeae strain 1291 and N. meningitidis strain MC58 were cloned, inactivated by inserting a kanamycin resistance cassette and used to make knockout mutants by allelic exchange. Both N. gonorrhoeae and N. meningitidis sco mutants were highly sensitive to oxidative killing by paraquat, indicating that Sco is involved in protection against oxidative stress in these bacteria. (C) 2003 Federation of European Biochemical Societies. Published by Elsevier Science B.V. All rights reserved.

Periplasmic nitrate reductases: structural and spectroscopic studies

Dissertação apresentada para obtenção do grau de Doutor em Bioquímica, especialidade Bioquímica-Física, pela Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa

Modeling chemotaxis reveals the role of reversed phosphotransfer and a bi-functional kinase-phosphatase

Understanding how multiple signals are integrated in living cells to produce a balanced response is a major challenge in biology. Two-component signal transduction pathways, such as bacterial chemotaxis, comprise histidine protein kinases (HPKs) and response regulators (RRs). These are used to sense and respond to changes in the environment. Rhodobacter sphaeroides has a complex chemosensory network with two signaling clusters, each containing a HPK, CheA. Here we demonstrate, using a mathematical model, how the outputs of the two signaling clusters may be integrated. We use our mathematical model supported by experimental data to predict that: (1) the main RR controlling flagellar rotation, CheY6, aided by its specific phosphatase, the bifunctional kinase CheA3, acts as a phosphate sink for the other RRs; and (2) a phosphorelay pathway involving CheB2 connects the cytoplasmic cluster kinase CheA3 with the polar localised kinase CheA2, and allows CheA3-P to phosphorylate non-cognate chemotaxis RRs. These two mechanisms enable the bifunctional kinase/phosphatase activity of CheA3 to integrate and tune the sensory output of each signaling cluster to produce a balanced response. The signal integration mechanisms identified here may be widely used by other bacteria, since like R. sphaeroides, over 50% of chemotactic bacteria have multiple cheA homologues and need to integrate signals from different sources.

Modelo de poços quânticos para a transferência de elétrons

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Hydrogen production and consumption of organic acids by a phototrophic microbial consortium

Alternative fuel sources have been extensively studied. Hydrogen gas has gained attention because its combustion releases only water, and it can be produced by microorganisms using organic acids as substrates. The aim of this study was to enrich a microbial consortium of photosynthetic purple non-sulfur bacteria from an Upflow Anaerobic Sludge Blanket reactor (UASB) using malate as carbon source. After the enrichment phase, other carbon sources were tested, such as acetate (30 mmol l(-1)), butyrate (17 mmol l(-1)), citrate (11 mmol l(-1)), lactate (23 mmol l(-1)) and malate (14.5 mmol l(-1)). The reactors were incubated at 30 degrees C under constant illumination by 3 fluorescent lamps (81 mu mol m(-2) s(-1)). The cumulative hydrogen production was 7.8, 9.0, 7.9, 5.6 and 13.9 mmol H-2 l(-1) culture for acetate, butyrate, citrate, lactate and malate, respectively. The maximum hydrogen yield was 0.6, 1.4, 0.7, 0.5 and 0.9 mmol H-2 mmol(-1) substrate for acetate, butyrate, citrate, lactate and malate, respectively. The consumption of substrates was 43% for acetate, 37% for butyrate, 100% for citrate, 49% for lactate and 100% for malate. Approximately 26% of the clones obtained from the Phototrophic Hydrogen-Producing Bacterial Consortium (PHPBC) were similar to Rhodobacter, Rhodospirillum and Rhodopseudomonas, which have been widely cited in studies of photobiological hydrogen production. Clones similar to the genus Sulfurospirillum (29% of the total) were also found in the microbial consortium. Copyright (C) 2012, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

The photosynthetic bacterial reaction center in native and artificial envirnoments: effects on light-induced electron transfer

In this thesis we focussed on the characterization of the reaction center (RC) protein purified from the photosynthetic bacterium Rhodobacter sphaeroides. In particular, we discussed the effects of native and artificial environment on the light-induced electron transfer processes. The native environment consist of the inner antenna LH1 complex that copurifies with the RC forming the so called core complex, and the lipid phase tightly associated with it. In parallel, we analyzed the role of saccharidic glassy matrices on the interplay between electron transfer processes and internal protein dynamics. As a different artificial matrix, we incorporated the RC protein in a layer-by-layer structure with a twofold aim: to check the behaviour of the protein in such an unusual environment and to test the response of the system to herbicides. By examining the RC in its native environment, we found that the light-induced charge separated state P+QB - is markedly stabilized (by about 40 meV) in the core complex as compared to the RC-only system over a physiological pH range. We also verified that, as compared to the average composition of the membrane, the core complex copurifies with a tightly bound lipid complement of about 90 phospholipid molecules per RC, which is strongly enriched in cardiolipin. In parallel, a large ubiquinone pool was found in association with the core complex, giving rise to a quinone concentration about ten times larger than the average one in the membrane. Moreover, this quinone pool is fully functional, i.e. it is promptly available at the QB site during multiple turnover excitation of the RC. The latter two observations suggest important heterogeneities and anisotropies in the native membranes which can in principle account for the stabilization of the charge separated state in the core complex. The thermodynamic and kinetic parameters obtained in the RC-LH1 complex are very close to those measured in intact membranes, indicating that the electron transfer properties of the RC in vivo are essentially determined by its local environment. The studies performed by incorporating the RC into saccharidic matrices evidenced the relevance of solvent-protein interactions and dynamical coupling in determining the kinetics of electron transfer processes. The usual approach when studying the interplay between internal motions and protein function consists in freezing the degrees of freedom of the protein at cryogenic temperature. We proved that the “trehalose approach” offers distinct advantages with respect to this traditional methodology. We showed, in fact, that the RC conformational dynamics, coupled to specific electron transfer processes, can be modulated by varying the hydration level of the trehalose matrix at room temperature, thus allowing to disentangle solvent from temperature effects. The comparison between different saccharidic matrices has revealed that the structural and dynamical protein-matrix coupling depends strongly upon the sugar. The analyses performed in RCs embedded in polyelectrolyte multilayers (PEM) structures have shown that the electron transfer from QA - to QB, a conformationally gated process extremely sensitive to the RC environment, can be strongly modulated by the hydration level of the matrix, confirming analogous results obtained for this electron transfer reaction in sugar matrices. We found that PEM-RCs are a very stable system, particularly suitable to study the thermodynamics and kinetics of herbicide binding to the QB site. These features make PEM-RC structures quite promising in the development of herbicide biosensors. The studies discussed in the present thesis have shown that, although the effects on electron transfer induced by the native and artificial environments tested are markedly different, they can be described on the basis of a common kinetic model which takes into account the static conformational heterogeneity of the RC and the interconversion between conformational substates. Interestingly, the same distribution of rate constants (i.e. a Gamma distribution function) can describe charge recombination processes in solutions of purified RC, in RC-LH1 complexes, in wet and dry RC-PEM structures and in glassy saccharidic matrices over a wide range of hydration levels. In conclusion, the results obtained for RCs in different physico-chemical environments emphasize the relevance of the structure/dynamics solvent/protein coupling in determining the energetics and the kinetics of electron transfer processes in a membrane protein complex.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

Direct electron transfer to cytochrome c oxidase investigated by electrochemistry and time-resolved surface-enhanced infrared absorption spectroscopy

Cytochrom c Oxidase (CcO), der Komplex IV der Atmungskette, ist eine der Häm-Kupfer enthaltenden Oxidasen und hat eine wichtige Funktion im Zellmetabolismus. Das Enzym enthält vier prosthetische Gruppen und befindet sich in der inneren Membran von Mitochondrien und in der Zellmembran einiger aerober Bakterien. Die CcO katalysiert den Elektronentransfer (ET) von Cytochrom c zu O2, wobei die eigentliche Reaktion am binuklearen Zentrum (CuB-Häm a3) erfolgt. Bei der Reduktion von O2 zu zwei H2O werden vier Protonen verbraucht. Zudem werden vier Protonen über die Membran transportiert, wodurch eine elektrochemische Potentialdifferenz dieser Ionen zwischen Matrix und Intermembranphase entsteht. Trotz ihrer Wichtigkeit sind Membranproteine wie die CcO noch wenig untersucht, weshalb auch der Mechanismus der Atmungskette noch nicht vollständig aufgeklärt ist. Das Ziel dieser Arbeit ist, einen Beitrag zum Verständnis der Funktion der CcO zu leisten. Hierzu wurde die CcO aus Rhodobacter sphaeroides über einen His-Anker, der am C-Terminus der Untereinheit II angebracht wurde, an eine funktionalisierte Metallelektrode in definierter Orientierung gebunden. Der erste Elektronenakzeptor, das CuA, liegt dabei am nächsten zur Metalloberfläche. Dann wurde eine Doppelschicht aus Lipiden insitu zwischen die gebundenen Proteine eingefügt, was zur sog. proteingebundenen Lipid-Doppelschicht Membran (ptBLM) führt. Dabei musste die optimale Oberflächenkonzentration der gebundenen Proteine herausgefunden werden. Elektrochemische Impedanzspektroskopie(EIS), Oberflächenplasmonenresonanzspektroskopie (SPR) und zyklische Voltammetrie (CV) wurden angewandt um die Aktivität der CcO als Funktion der Packungsdichte zu charakterisieren. Der Hauptteil der Arbeit betrifft die Untersuchung des direkten ET zur CcO unter anaeroben Bedingungen. Die Kombination aus zeitaufgelöster oberflächenverstärkter Infrarot-Absorptionsspektroskopie (tr-SEIRAS) und Elektrochemie hat sich dafür als besonders geeignet erwiesen. In einer ersten Studie wurde der ET mit Hilfe von fast scan CV untersucht, wobei CVs von nicht-aktivierter sowie aktivierter CcO mit verschiedenen Vorschubgeschwindigkeiten gemessen wurden. Die aktivierte Form wurde nach dem katalytischen Umsatz des Proteins in Anwesenheit von O2 erhalten. Ein vier-ET-modell wurde entwickelt um die CVs zu analysieren. Die Methode erlaubt zwischen dem Mechanismus des sequentiellen und des unabhängigen ET zu den vier Zentren CuA, Häm a, Häm a3 und CuB zu unterscheiden. Zudem lassen sich die Standardredoxpotentiale und die kinetischen Koeffizienten des ET bestimmen. In einer zweiten Studie wurde tr-SEIRAS im step scan Modus angewandt. Dafür wurden Rechteckpulse an die CcO angelegt und SEIRAS im ART-Modus verwendet um Spektren bei definierten Zeitscheiben aufzunehmen. Aus diesen Spektren wurden einzelne Banden isoliert, die Veränderungen von Vibrationsmoden der Aminosäuren und Peptidgruppen in Abhängigkeit des Redoxzustands der Zentren zeigen. Aufgrund von Zuordnungen aus der Literatur, die durch potentiometrische Titration der CcO ermittelt wurden, konnten die Banden versuchsweise den Redoxzentren zugeordnet werden. Die Bandenflächen gegen die Zeit aufgetragen geben dann die Redox-Kinetik der Zentren wieder und wurden wiederum mit dem vier-ET-Modell ausgewertet. Die Ergebnisse beider Studien erlauben die Schlussfolgerung, dass der ET zur CcO in einer ptBLM mit größter Wahrscheinlichkeit dem sequentiellen Mechanismus folgt, was dem natürlichen ET von Cytochrom c zur CcO entspricht.

Stable transduction of bovine TLR4 and bovine MD-2 into LPS-nonresponsive cells and soluble CD14 promote the ability to respond to LPS

The interaction of bovine cells with lipopolysaccharide (LPS) was explored using human embryo kidney (HEK) 293 cell line stably transduced with bovine toll-like receptor-4 (TLR4) alone or in combination with bovine MD-2. These lines and mock-transduced HEK293 cells were tested by flow cytometry for LPS-fluorescein isothiocyanate (LPS-FITC) binding, nuclear factor kappa B (NFkappaB) activation, interleukin-8 (IL-8) production and interferon-beta mRNA expression/interferon (IFN) type I production. Whereas bovine TLR4 was sufficient to promote binding of high concentrations of LPS-FITC, both bovine TLR4 and MD-2 were required for activation by LPS, as assessed by NFkappaB activation and IL-8 production. Induction of IFN bioactivity was not observed in doubly transduced HEK293 cells, and no evidence for IFN-beta mRNA induction in response to LPS was obtained, although cells responded by IFN-beta mRNA expression to stimulation by Sendai virus and poly-inosinic acid-poly-cytidylic acid (poly(I:C)). Cells stably transduced with both bovine TLR4 and bovine MD-2 responded to LPS by IL-8 production, in decreasing order, in the presence of fetal bovine serum (FCS), of human serum, and of human serum albumin (HSA). The reduced activity in the presence of HSA could be restored by the addition of soluble CD14 (sCD14) but not of LPS binding protein (LBP). This is in contrast to macrophages which show a superior response to LPS in the presence of HSA when compared with macrophages stimulated by LPS in the presence of FCS. This suggests that macrophages but not HEK293 cells express factors rendering LPS stimulation serum-independent. Stably double-transduced cells reacted, in decreasing order, to LPS from Rhodobacter sphaeroides, to LPS from Escherichia coli, to synthetic lipd-IVa (compound 406), to diphosphoryl-lipid-A (S. minnesota) and to monophosphoryl-lipid-A (S. minnesota). They failed to react to the murine MD-2/TLR4 ligand taxol. This resembles the reactivity of bovine macrophages with regard to sensitivity (ED(50)) and order of potency but is distinct from the reactivity pattern of other species. This formally establishes that in order to react to LPS, cattle cells require serum factors (e.g. sCD14) and cell-expressed factors such as MD-2 and TLR4. The cell lines described are the first of a series expressing defined pattern recognition receptors (PRR) of bovine origin. They will be useful in the study of the interaction of the bovine TLR4-MD-2 complex and Gram-negative bovine pathogens, e.g. the agents causing Gram-negative bovine mastitis.

Mimicking respiratory phosphorylation using purified enzymes

The enzymes of oxidative phosphorylation are a striking example of the functional association of multiple enzyme complexes, working together to form ATP from cellular reducing equivalents. These complexes, such as cytochrome c oxidase or the ATP synthase, are typically investigated individually and therefore, their functional interplay is not well understood. Here, we present methodology that allows the co-reconstitution of purified terminal oxidases and ATP synthases in synthetic liposomes. The enzymes are functionally coupled via proton translocation where upon addition of reducing equivalents the oxidase creates and maintains a transmembrane electrochemical proton gradient that energizes the synthesis of ATP by the F1F0 ATP synthase. The method has been tested with the ATP synthases from Escherichia coli and spinach chloroplasts, and with the quinol and cytochrome c oxidases from E. coli and Rhodobacter sphaeroides, respectively. Unlike in experiments with the ATP synthase reconstituted alone, the setup allows in vitro ATP synthesis under steady state conditions, with rates up to 90 ATP×s(-1)×enzyme(-1). We have also used the novel system to study the phenomenon of "mild uncoupling" as observed in mitochondria upon addition of low concentrations of ionophores (e.g. FCCP, SF6847) and the recoupling effect of 6-ketocholestanol. While we could reproduce the described effects, our data with the in vitro system does not support the idea of a direct interaction between a mitochondrial protein and the uncoupling agents as proposed earlier.

Implementación de algoritmos de ensamblaje de genomas en sistemas de memoria compartida y memoria distribuida

El desarrollo de algoritmos ensambladores de genes y la utilización de estos está viviendo un aumento muy espectacular en los últimos años. Debido a las mejoras ofrecidas en los dispositivos hardware de los numerosos supercomputadores que existen hoy en día se pueden realizar experimentos científicos de una manera más asequible que hace unos años. Este proyecto servirá como introducción en el complejo mundo de algoritmos científicos, más concretamente en algoritmos ensambladores de genomas. Veremos de primera mano cómo utilizar estas nuevas tecnologías, con ejemplos sencillos, pero con un desarrollo lo bastante importante para darnos una idea del funcionamiento de todas las fases de experimentación que engloban los algoritmos ensambladores y la utilización de la programación paralela en supercomputadores. Concretamente en este proyecto se van a analizar exhaustivamente una serie de algoritmos ensambladores que serán probados en uno de los supercomputadores más potentes de España, el Magerit 2. En estas pruebas vamos a proceder al ensamblado de genomas de tres tipos de organismos como bacterias (Staphylococcus Aureus, y Rhodobacter Sphaeroides) y una prueba gran escala con el genoma del Cromosoma 14 del Homo Sapiens Sapiens (Ser humano). Después procederemos a la comparación de todos los resultados obtenidos para poder comprobar que algoritmos realizan mejor su trabajo y ajustar dicha decisión a las necesidades que tenemos actualmente para buscar un algoritmo eficaz.

Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.