The Rieske Protein: A Case Study on the Pitfalls of Multiple Sequence Alignments and Phylogenetic Reconstruction

Something old, something new, something borrowed, something blue: the anaerobic microbial ancestry of aerobic respiration

Aerobic respiration evolved by bricolage, with modules cobbled together as microbial biochemistry coevolved with Earth's geochemistry. The mitochondrial electron transport chain represents a patchwork of respiratory modules inherited from microbial methanogenesis, iron oxidation, anoxygenic photosynthesis, and denitrification pathways, and preserves a biochemical record of Earth's redox environment over its four-billion-year history. Imprints of the anoxic early Earth are recognizable in Complex I's numerous iron-sulfur cofactors and vestigial binding sites for ferredoxin, nickel-iron, and molybdopterin, whereas the more recent advent of oxygen as a terminal electron acceptor necessitated use of heme and copper cofactors by Complex IV. Bricolage of respiratory complexes resulted in supercomplexes for improved electron transfer efficiency in some bacteria and archaea, and in many eukaryotes. Accessory subunits evolved to wrap mitochondrial supercomplexes for improved assembly and stability. Environmental microbes with 'fossil' proteins that are similar to ancestral forms of the respiratory complexes deserve further scrutiny and may reveal new insights on the evolution of aerobic respiration.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

The controversy on the ancestral arsenite oxidizing enzyme; deducing evolutionary histories with phylogeny and thermodynamics

The three presently known enzymes responsible for arsenic-using bioenergetic processes are arsenite oxidase (Aio), arsenate reductase (Arr) and alternative arsenite oxidase (Arx), all of which are molybdoenzymes from the vast group referred to as the Mo/W-bisPGD enzyme superfamily. Since arsenite is present in substantial amounts in hydrothermal environments, frequently considered as vestiges of primordial biochemistry, arsenite-based bioenergetics has long been predicted to be ancient. Conflicting scenarios, however, have been put forward proposing either Arr/Arx or Aio as operating in the ancestral metabolism. Phylogenetic data argue in favor of Aio whereas biochemical and physiological data led several authors to propose Arx/Arr as the most ancient anaerobic arsenite metabolizing enzymes. Here we combine phylogenetic approaches with physiological and biochemical experiments to demonstrate that the Arx/Arr enzymes could not have been functional in the Archaean geological eon. We propose that Arr reacts with menaquinones to reduce arsenate whereas Arx reacts with ubiquinone to oxidize arsenite, in line with thermodynamic considerations. The distribution of the quinone biosynthesis pathways, however, clearly indicates that the ubiquinone pathway is recent. An updated phylogeny of Arx furthermore reinforces the hypothesis of a recent emergence of this enzyme. We therefore conclude that anaerobic arsenite redox conversion in the Archaean must have been performed in a metabolism involving Aio.

Hybrid cluster proteins in a photosynthetic microalga

Hybrid cluster proteins (HCPs) are metalloproteins characterized by the presence of an iron-sulfur-oxygen cluster. These proteins occur in all three domains of life. In eukaryotes, HCPs have so far been found only in a few anaerobic parasites and photosynthetic microalgae. With respect to all species harboring an HCP, the green microalga Chlamydomonas reinhardtii stands out by the presence of four HCP genes. The study of the gene and protein structures as well as the phylogenetic analyses strongly support a model in which the HCP family in the alga has emerged from a single gene of alpha proteobacterial origin and then expanded by several rounds of duplications. The spectra and redox properties of HCP1 and HCP3, produced heterologously in Escherichia coli, were analyzed by electron paramagnetic resonance spectroscopy on redox-titrated samples. Both proteins contain a [4Fe-4S]-cluster as well as a [4Fe-2O-2S]-hybrid cluster with paramagnetic properties related to those of HCPs from Desulfovibrio species. Immunoblotting experiments combined with mass spectrometry-based proteomics showed that both nitrate and darkness contribute to the strong upregulation of the HCP levels in C. reinhardtii growing under oxic conditions. The link to the nitrate metabolism is discussed in the light of recent data on the potential role of HCP in S-nitrosylation in bacteria.

On the Natural History of Flavin-Based Electron Bifurcation

Electron bifurcation is here described as a special case of the continuum of electron transfer reactions accessible to two-electron redox compounds with redox cooperativity. We argue that electron bifurcation is foremost an electrochemical phenomenon based on (a) strongly inverted redox potentials of the individual redox transitions, (b) a high endergonicity of the first redox transition, and (c) an escapement-type mechanism rendering completion of the first electron transfer contingent on occurrence of the second one. This mechanism is proposed to govern both the traditional quinone-based and the newly discovered flavin-based versions of electron bifurcation. Conserved and variable aspects of the spatial arrangement of electron transfer partners in flavoenzymes are assayed by comparing the presently available 3D structures. A wide sample of flavoenzymes is analyzed with respect to conserved structural modules and three major structural groups are identified which serve as basic frames for the evolutionary construction of a plethora of flavin-containing redox enzymes. We argue that flavin-based and other types of electron bifurcation are of primordial importance to free energy conversion, the quintessential foundation of life, and discuss a plausible evolutionary ancestry of the mechanism.

Emergence of cytochrome bc complexes in the context of photosynthesis

The cytochrome bc (cyt bc) complexes are involved in Q-cycling; they oxidize membrane quinols by high-potential electron acceptors, such as cytochromes or plastocyanin, and generate transmembrane proton gradient. In several prokaryotic lineages, and also in plant chloroplasts, the catalytic core of the cyt bc complexes is built of a four-helical cytochrome b (cyt b) that contains three hemes, a three-helical subunit IV, and an iron-sulfur Rieske protein (cytochrome b₆ f-type complexes). In other prokaryotic lineages, and also in mitochondria, the cyt b subunit is fused with subunit IV, yielding a seven- or eight-helical cyt b with only two hemes (cyt bc₁ -type complexes). Here we present an updated phylogenomic analysis of the cyt b subunits of cyt bc complexes. This analysis provides further support to our earlier suggestion that (1) the ancestral version of cyt bc complex contained a small four-helical cyt b with three hemes similar to the plant cytochrome b₆ and (2) independent fusion events led to the formation of large cyts b in several lineages. In the search for a primordial function for the ancestral cyt bc complex, we address the intimate connection between the cyt bc complexes and photosynthesis. Indeed, the Q-cycle turnover in the cyt bc complexes demands high-potential electron acceptors. Before the Great Oxygenation Event, the biosphere had been highly reduced, so high-potential electron acceptors could only be generated upon light-driven charge separation. It appears that an ancestral cyt bc complex capable of Q-cycling has emerged in conjunction with the (bacterio)chlorophyll-based photosynthetic systems that continuously generated electron vacancies at the oxidized (bacterio)chlorophyll molecules.

The molecular evolution of the Qo motif

Quinol oxidation in the catalytic quinol oxidation site (Q_o site) of cytochrome (cyt) bc₁ complexes is the key step of the Q cycle mechanism, which laid the ground for Mitchell’s chemiosmotic theory of energy conversion. Bifurcated electron transfer upon quinol oxidation enables proton uptake and release on opposite membrane sides, thus generating a proton gradient that fuels ATP synthesis in cellular respiration and photosynthesis. The Q_o site architecture formed by cyt b and Rieske iron–sulfur protein (ISP) impedes harmful bypass reactions. Catalytic importance is assigned to four residues of cyt b formerly described as PEWY motif in the context of mitochondrial complexes, which we now denominate Q_o motif as comprehensive evolutionary sequence analysis of cyt b shows substantial natural variance of the motif with phylogenetically specific patterns. In particular, the Q_o motif is identified as PEWY in mitochondria, α- and ε-Proteobacteria, Aquificae, Chlorobi, Cyanobacteria, and chloroplasts. PDWY is present in Gram-positive bacteria, Deinococcus–Thermus and haloarchaea, and PVWY in β- and γ-Proteobacteria. PPWF only exists in Archaea. Distinct patterns for acidophilic organisms indicate environment-specific adaptations. Importantly, the presence of PDWY and PEWY is correlated with the redox potential of Rieske ISP and quinone species. We propose that during evolution from low to high potential electron-transfer systems in the emerging oxygenic atmosphere, cyt bc₁ complexes with PEWY as Q_o motif prevailed to efficiently use high potential ubiquinone as substrate, whereas cyt b with PDWY operate best with low potential Rieske ISP and menaquinone, with the latter being the likely composition of the ancestral cyt bc₁ complex.

A new class of marine Euryarchaeota group II from the Mediterranean deep chlorophyll maximum

We have analyzed metagenomic fosmid clones from the deep chlorophyll maximum (DCM), which, by genomic parameters, correspond to the 16S ribosomal RNA (rRNA)-defined marine Euryarchaeota group IIB (MGIIB). The fosmid collections associated with this group add up to 4 Mb and correspond to at least two species within this group. From the proposed essential genes contained in the collections, we infer that large sections of the conserved regions of the genomes of these microbes have been recovered. The genomes indicate a photoheterotrophic lifestyle, similar to that of the available genome of MGIIA (assembled from an estuarine metagenome in Puget Sound, Washington Pacific coast), with a proton-pumping rhodopsin of the same kind. Several genomic features support an aerobic metabolism with diversified substrate degradation capabilities that include xenobiotics and agar. On the other hand, these MGIIB representatives are non-motile and possess similar genome size to the MGIIA-assembled genome, but with a lower GC content. The large phylogenomic gap with other known archaea indicates that this is a new class of marine Euryarchaeota for which we suggest the name Thalassoarchaea. The analysis of recruitment from available metagenomes indicates that the representatives of group IIB described here are largely found at the DCM (ca. 50 m deep), in which they are abundant (up to 0.5% of the reads), and at the surface mostly during the winter mixing, which explains formerly described 16S rRNA distribution patterns. Their uneven representation in environmental samples that are close in space and time might indicate sporadic blooms.

Bioenergetic evolution in proteobacteria and mitochondria

Mitochondria are the energy-producing organelles of our cells and derive from bacterial ancestors that became endosymbionts of microorganisms from a different lineage, together with which they formed eukaryotic cells. For a long time it has remained unclear from which bacteria mitochondria actually evolved, even if these organisms in all likelihood originated from the α lineage of proteobacteria. A recent article (Degli Esposti M, et al. 2014. Evolution of mitochondria reconstructed from the energy metabolism of living bacteria. PLoS One 9:e96566) has presented novel evidence indicating that methylotrophic bacteria could be among the closest living relatives of mitochondrial ancestors. Methylotrophs are ubiquitous bacteria that live on single carbon sources such as methanol and methane; in the latter case they are called methanotrophs. In this review, I examine their possible ancestry to mitochondria within a survey of the common features that can be found in the central and terminal bioenergetic systems of proteobacteria and mitochondria. I also discuss previously overlooked information on methanotrophic bacteria, in particular their intracytoplasmic membranes resembling mitochondrial cristae and their capacity of establishing endosymbiotic relationships with invertebrate animals and archaic plants. This information appears to sustain the new idea that mitochondrial ancestors could be related to extant methanotrophic proteobacteria, a possibility that the genomes of methanotrophic endosymbionts will hopefully clarify.

Diversity of arsenite oxidizing bacterial communities in arsenic-rich deltaic aquifers in West Bengal, India

High arsenic (As) concentration in groundwater has affected human health, particularly in South-East Asia putting millions of people at risk. Biogeochemical cycling of As carried out by different bacterial groups are suggested to control the As fluxes in aquifers. A functional diversity approach in link with As precipitation was adopted to study bacterial community structures and their variation within the As contaminated Bengal Delta Plain (BDP) aquifers of India. Groundwater samples collected from two shallow aquifers in Karimpur II (West Bengal, India), during years 2010 and 2011, were investigated to trace the effects immediately after monsoon period (precipitation) on community structure and diversity of bacterial assemblages with a focus on arsenite oxidizing bacterial phyla for two successive years. The study focused on amplification, clone library generation and sequencing of the arsenite oxidase large sub-unit gene aioA and 16S rRNA marker, with respect to changes in elemental concentrations. New set of primers were designed to amplify the aioA gene as a phylogenetic marker to study taxonomically diverse arsenite oxidizing bacterial groups in these aquifers. The overall narrow distribution of bacterial communities based on aioA and 16S rRNA sequences observed was due to poor nutrient status and anoxic conditions in these As contaminated aquifers. Proteobacteria was the dominant phylum detected, within which Acidovorax, Hydrogenophaga, Albidiferax, Bosea, and Polymorphum were the major arsenite oxidizing bacterial genera based on the number of clones sequenced. The structure of bacterial assemblages including those of arsenite oxidizing bacteria seems to have been affected by increase in major elemental concentrations (e.g., As, Fe, S, and Si) within two sampling sessions, which was supported by statistical analyses. One of the significant findings of this study is detection of novel lineages of 16S rRNA-like bacterial sequences indicating presence of indigenous bacterial communities BDP wells that can play important role in biogeochemical cycling of elements including As.

Evolution of mitochondria reconstructed from the energy metabolism of living bacteria

The ancestors of mitochondria, or proto-mitochondria, played a crucial role in the evolution of eukaryotic cells and derived from symbiotic α-proteobacteria which merged with other microorganisms - the basis of the widely accepted endosymbiotic theory. However, the identity and relatives of proto-mitochondria remain elusive. Here we show that methylotrophic α-proteobacteria could be the closest living models for mitochondrial ancestors. We reached this conclusion after reconstructing the possible evolutionary pathways of the bioenergy systems of proto-mitochondria with a genomic survey of extant α-proteobacteria. Results obtained with complementary molecular and genetic analyses of diverse bioenergetic proteins converge in indicating the pathway stemming from methylotrophic bacteria as the most probable route of mitochondrial evolution. Contrary to other α-proteobacteria, methylotrophs show transition forms for the bioenergetic systems analysed. Our approach of focusing on these bioenergetic systems overcomes the phylogenetic impasse that has previously complicated the search for mitochondrial ancestors. Moreover, our results provide a new perspective for experimentally re-evolving mitochondria from extant bacteria and in the future produce synthetic mitochondria.

Free energy conversion in the LUCA: Quo vadis?

Living entities are unimaginable without means to harvest free energy from the environment, that is, without bioenergetics. The quest to understand the bioenergetic ways of early life therefore is one of the crucial elements to understand the emergence of life on our planet. Over the last few years, several mutually exclusive scenarios for primordial bioenergetics have been put forward, all of which are based on some sort of empirical observation, a remarkable step forward from the previous, essentially untestable, ab initio models. We here try to present and compare these scenarios while at the same time discuss their respective empirical weaknesses. The goal of this article is to harness crucial new expertise from the entire field by stimulating a larger part of the bioenergetics community to become involved in "origin-of-energy-metabolism" research. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Evolution of cytochrome bc complexes: from membrane-anchored dehydrogenases of ancient bacteria to triggers of apoptosis in vertebrates

This review traces the evolution of the cytochrome bc complexes from their early spread among prokaryotic lineages and up to the mitochondrial cytochrome bc1 complex (complex III) and its role in apoptosis. The results of phylogenomic analysis suggest that the bacterial cytochrome b6f-type complexes with short cytochromes b were the ancient form that preceded in evolution the cytochrome bc1-type complexes with long cytochromes b. The common ancestor of the b6f-type and the bc1-type complexes probably resembled the b6f-type complexes found in Heliobacteriaceae and in some Planctomycetes. Lateral transfers of cytochrome bc operons could account for the several instances of acquisition of different types of bacterial cytochrome bc complexes by archaea. The gradual oxygenation of the atmosphere could be the key evolutionary factor that has driven further divergence and spread of the cytochrome bc complexes. On the one hand, oxygen could be used as a very efficient terminal electron acceptor. On the other hand, auto-oxidation of the components of the bc complex results in the generation of reactive oxygen species (ROS), which necessitated diverse adaptations of the b6f-type and bc1-type complexes, as well as other, functionally coupled proteins. A detailed scenario of the gradual involvement of the cardiolipin-containing mitochondrial cytochrome bc1 complex into the intrinsic apoptotic pathway is proposed, where the functioning of the complex as an apoptotic trigger is viewed as a way to accelerate the elimination of the cells with irreparably damaged, ROS-producing mitochondria. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.

Multiple Rieske/cytb complexes in a single organism

Most organisms contain a single Rieske/cytb complex. This enzyme can be integrated in any respiratory or photosynthetic electron transfer chain that is quinone-based and sufficiently energy rich to allow for the turnover of three enzymes - a quinol reductase, a Rieske/cytb complex and a terminal oxidase. Despite this universal usability of the enzyme a variety of phylogenetically distant organisms have multiple copies thereof and no reason for this redundancy is obvious. In this review we present an overview of the distribution of multiple copies among species and describe their properties from the scarce experimental results, analysis of their amino acid sequences and genomic context. We discuss the predicted redox properties of the Rieske cluster in relation to the nature of the pool quinone. It appears that acidophilic iron-oxidizing bacteria specialized one of their two copies for reverse electron transfer, archaeal Thermoprotei adapted their three copies to the interaction with different oxidases and several, phylogenetically unrelated species imported a second complex with a putative heme ci that may confer some yet to be determined properties to the complex. These hypothesis and all the more the so far completely unexplained cases call for further studies and we put forward a number of suggestions for future research that we hope to be stimulating for the field. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.

The prokaryotic Mo/W-bisPGD enzymes family: a catalytic workhorse in bioenergetic

Over the past two decades, prominent importance of molybdenum-containing enzymes in prokaryotes has been put forward by studies originating from different fields. Proteomic or bioinformatic studies underpinned that the list of molybdenum-containing enzymes is far from being complete with to date, more than fifty different enzymes involved in the biogeochemical nitrogen, carbon and sulfur cycles. In particular, the vast majority of prokaryotic molybdenum-containing enzymes belong to the so-called dimethylsulfoxide reductase family. Despite its extraordinary diversity, this family is characterized by the presence of a Mo/W-bis(pyranopterin guanosine dinucleotide) cofactor at the active site. This review highlights what has been learned about the properties of the catalytic site, the modular variation of the structural organization of these enzymes, and their interplay with the isoprenoid quinones. In the last part, this review provides an integrated view of how these enzymes contribute to the bioenergetics of prokaryotes. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Insight into the evolution of the iron oxidation pathways

Iron is a ubiquitous element in the universe. Ferrous iron (Fe(II)) was abundant in the primordial ocean until the oxygenation of the Earth's atmosphere led to its widespread oxidation and precipitation. This change of iron bioavailability likely put selective pressure on the evolution of life. This element is essential to most extant life forms and is an important cofactor in many redox-active proteins involved in a number of vital pathways. In addition, iron plays a central role in many environments as an energy source for some microorganisms. This review is focused on Fe(II) oxidation. The fact that the ability to oxidize Fe(II) is widely distributed in Bacteria and Archaea and in a number of quite different biotopes suggests that the dissimilatory Fe(II) oxidation is an ancient energy metabolism. Based on what is known today about Fe(II) oxidation pathways, we propose that they arose independently more than once in evolution and evolved convergently. The iron paleochemistry, the phylogeny, the physiology of the iron oxidizers, and the nature of the cofactors of the redox proteins involved in these pathways suggest a possible scenario for the timescale in which each type of Fe(II) oxidation pathways evolved. The nitrate dependent anoxic iron oxidizers are likely the most ancient iron oxidizers. We suggest that the phototrophic anoxic iron oxidizers arose in surface waters after the Archaea/Bacteria-split but before the Great Oxidation Event. The neutrophilic oxic iron oxidizers possibly appeared in microaerobic marine environments prior to the Great Oxidation Event while the acidophilic ones emerged likely after the advent of atmospheric O(2). This article is part of a Special Issue entitled: The evolutionary aspects of bioenergetic systems.

On the universal core of bioenergetics

Living cells are able to harvest energy by coupling exergonic electron transfer between reducing and oxidising substrates to the generation of chemiosmotic potential. Whereas a wide variety of redox substrates is exploited by prokaryotes resulting in very diverse layouts of electron transfer chains, the ensemble of molecular architectures of enzymes and redox cofactors employed to construct these systems is stunningly small and uniform. An overview of prominent types of electron transfer chains and of their characteristic electrochemical parameters is presented. We propose that basic thermodynamic considerations are able to rationalise the global molecular make-up and functioning of these chemiosmotic systems. Arguments from palaeogeochemistry and molecular phylogeny are employed to discuss the evolutionary history leading from putative energy metabolisms in early life to the chemiosmotic diversity of extant organisms. Following the Occam's razor principle, we only considered for this purpose origin of life scenarios which are contiguous with extant life. This article is part of a Special Issue entitled: The evolutionary aspects of bioenergetic systems.

Arsenics as bioenergetic substrates

Although at low concentrations, arsenic commonly occurs naturally as a local geological constituent. Whereas both arsenate and arsenite are strongly toxic to life, a number of prokaryotes use these compounds as electron acceptors or donors, respectively, for bioenergetic purposes via respiratory arsenate reductase, arsenite oxidase and alternative arsenite oxidase. The recent burst in discovered arsenite oxidizing and arsenate respiring microbes suggests the arsenic bioenergetic metabolisms to be anything but exotic. The first goal of the present review is to bring to light the widespread distribution and diversity of these metabolizing pathways. The second goal is to present an evolutionary analysis of these diverse energetic pathways. Taking into account not only the available data on the arsenic metabolizing enzymes and their phylogenetical relatives but also the palaeogeochemical records, we propose a crucial role of arsenite oxidation via arsenite oxidase in primordial life. This article is part of a Special Issue entitled: The evolutionary aspects of bioenergetic systems.

Phylogeny of Rieske/cytb complexes with a special focus on the Haloarchaeal enzymes

Rieske/cytochrome b (Rieske/cytb) complexes are proton pumping quinol oxidases that are present in most bacteria and Archaea. The phylogeny of their subunits follows closely the 16S-rRNA phylogeny, indicating that chemiosmotic coupling was already present in the last universal common ancestor of Archaea and bacteria. Haloarchaea are the only organisms found so far that acquired Rieske/cytb complexes via interdomain lateral gene transfer. They encode two Rieske/cytb complexes in their genomes; one of them is found in genetic context with nitrate reductase genes and has its closest relatives among Actinobacteria and the Thermus/Deinococcus group. It is likely to function in nitrate respiration. The second Rieske/cytb complex of Haloarchaea features a split cytochrome b sequence as do Cyanobacteria, chloroplasts, Heliobacteria, and Bacilli. It seems that Haloarchaea acquired this complex from an ancestor of the above-mentioned phyla. Its involvement in the bioenergetic reaction chains of Haloarchaea is unknown. We present arguments in favor of the hypothesis that the ancestor of Haloarchaea, which relied on a highly specialized bioenergetic metabolism, that is, methanogenesis, and was devoid of quinones and most enzymes of anaerobic or aerobic bioenergetic reaction chains, integrated laterally transferred genes into its genome to respond to a change in environmental conditions that made methanogenesis unfavorable.

The ineluctable requirement for the trans-iron elements molybdenum and/or tungsten in the origin of life

An evolutionary tree of key enzymes from the Complex-Iron-Sulfur-Molybdoenzyme (CISM) superfamily distinguishes "ancient" members, i.e. enzymes present already in the last universal common ancestor (LUCA) of prokaryotes, from more recently evolved subfamilies. The majority of the presented subfamilies and, as a consequence, the Molybdo-enzyme superfamily as a whole, appear to have existed in LUCA. The results are discussed with respect to the nature of bioenergetic substrates available to early life and to problems arising from the low solubility of molybdenum under conditions of the primordial Earth.

Phylogenetic analysis reveals the surprising diversity of an oxygenase class

As metalloenzymes capable of transforming a broad range of substrates with high stereo- and regio-specificity, the multicomponent Rieske oxygenases (ROs) have been studied in bacterial systems for applications in bioremediation and industrial biocatalysis. These studies include genetic and biochemical investigations, determination of enzyme structure, phylogenetic analysis, and enzyme classification. Although RO terminal oxygenase components (RO-Os) share a conserved domain structure, their sequences are highly divergent and present significant challenges for identification and classification. Herein, we present the first global phylogenetic analysis of a broad range of RO-Os from diverse taxonomic groups. We employed objective, structure-based criteria to significantly reduce the inclusion of erroneously aligned sequences in the analysis. Our findings reveal that RO biochemical studies to date have been largely concentrated in an unexpectedly narrow portion of the RO-O sequence landscape. Additionally, our analysis demonstrates the existence two distinct groups of RO-O sequences. Finally, the sequence diversity recognized in this study necessitates a new RO-O classification scheme. We therefore propose a P450-like naming system. Our results reveal a diversity of sequence and potential catalytic functionality that has been wholly unappreciated in the RO literature. This study also demonstrates that many commonly used bioinformatic tools may not be sufficient to analyze the vast amount of data available in current databases. These findings facilitate the expanded exploration of RO catalytic capabilities in both biological and technological contexts and increase the potential for practical exploitation of their activities.

The Q cycle of cytochrome bc complexes: a structure perspective

Aspects of the crystal structures of the hetero-oligomeric cytochrome bc(1) and b(6)f ("bc") complexes relevant to their electron/proton transfer function and the associated redox reactions of the lipophilic quinones are discussed. Differences between the b(6)f and bc(1) complexes are emphasized. The cytochrome bc(1) and b(6)f dimeric complexes diverge in structure from a core of subunits that coordinate redox groups consisting of two bis-histidine coordinated hemes, a heme b(n) and b(p) on the electrochemically negative (n) and positive (p) sides of the complex, the high potential [2Fe-2S] cluster and c-type heme at the p-side aqueous interface and aqueous phase, respectively, and quinone/quinol binding sites on the n- and p-sides of the complex. The bc(1) and b(6)f complexes diverge in subunit composition and structure away from this core. b(6)f Also contains additional prosthetic groups including a c-type heme c(n) on the n-side, and a chlorophyll a and β-carotene. Common structure aspects; functions of the symmetric dimer. (I) Quinone exchange with the bilayer. An inter-monomer protein-free cavity of approximately 30Å along the membrane normal×25Å (central inter-monomer distance)×15Å (depth in the center), is common to both bc(1) and b(6)f complexes, providing a niche in which the lipophilic quinone/quinol (Q/QH(2)) can be exchanged with the membrane bilayer. (II) Electron transfer. The dimeric structure and the proximity of the two hemes b(p) on the electrochemically positive side of the complex in the two monomer units allow the possibility of two alternate routes of electron transfer across the complex from heme b(p) to b(n): intra-monomer and inter-monomer involving electron cross-over between the two hemes b(p). A structure-based summary of inter-heme distances in seven bc complexes, representing mitochondrial, chromatophore, cyanobacterial, and algal sources, indicates that, based on the distance parameter, the intra-monomer pathway would be favored kinetically. (III) Separation of quinone binding sites. A consequence of the dimer structure and the position of the Q/QH(2) binding sites is that the p-side QH(2) oxidation and n-side Q reduction sites are each well separated. Therefore, in the event of an overlap in residence time by QH(2) or Q molecules at the two oxidation or reduction sites, their spatial separation would result in minimal steric interference between extended Q or QH(2) isoprenoid chains. (IV) Trans-membrane QH(2)/Q transfer. (i) n/p-side QH(2)/Q transfer may be hindered by lipid acyl chains; (ii) the shorter less hindered inter-monomer pathway across the complex would not pass through the center of the cavity, as inferred from the n-side antimycin site on one monomer and the p-side stigmatellin site on the other residing on the same surface of the complex. (V) Narrow p-side portal for QH(2)/Q passage. The [2Fe-2S] cluster that serves as oxidant, and whose histidine ligand serves as a H(+) acceptor in the oxidation of QH(2), is connected to the inter-monomer cavity by a narrow extended portal, which is also occupied in the b(6)f complex by the 20 carbon phytyl chain of the bound chlorophyll.

Evolution of photosynthesis

Energy conversion of sunlight by photosynthetic organisms has changed Earth and life on it. Photosynthesis arose early in Earth's history, and the earliest forms of photosynthetic life were almost certainly anoxygenic (non-oxygen evolving). The invention of oxygenic photosynthesis and the subsequent rise of atmospheric oxygen approximately 2.4 billion years ago revolutionized the energetic and enzymatic fundamentals of life. The repercussions of this revolution are manifested in novel biosynthetic pathways of photosynthetic cofactors and the modification of electron carriers, pigments, and existing and alternative modes of photosynthetic carbon fixation. The evolutionary history of photosynthetic organisms is further complicated by lateral gene transfer that involved photosynthetic components as well as by endosymbiotic events. An expanding wealth of genetic information, together with biochemical, biophysical, and physiological data, reveals a mosaic of photosynthetic features. In combination, these data provide an increasingly robust framework to formulate and evaluate hypotheses concerning the origin and evolution of photosynthesis.

Heliobacterial Rieske/cytb complex

Data on structure and function of the Rieske/cytb complex from Heliobacteria are scarce. They indicate that the complex is related to the b (6) f complex in agreement with the phylogenetic position of the organism. It is composed of a diheme cytochrome c, and a Rieske iron-sulfur protein, together with transmembrane cytochrome b (6) and subunit IV. Additional small subunits may be part of the complex. The cofactor content comprises heme c (i), first discovered in the Q(i) binding pocket of b (6) f complexes. The redox midpoint potentials are more negative than in b (6) f complex in agreement with the lower redox midpoint potentials (by about 150 mV) of its reaction partners, menaquinone, and cytochrome c (553). The enzyme is implicated in cyclic electron transfer around the RCI. Functional studies are favored by the absence of antennae and the simple photosynthetic reaction chain but are hampered by the high oxygen sensitivity of the organism, its chlorophyll, and lipids.

The "green" phylogenetic clade of Rieske/cytb complexes

More than a decade ago, Heliobacteria were recognised to contain a Rieske/cytb complex in which the cytochrome b subunit is split into two separate proteins, a peculiar feature characteristic of the cyanobacterial and plastidic b (6) f complex. The common presence of RCI-type reaction centres further emphasise possible evolutionary links between Heliobacteria, Chlorobiaceae and Cyanobacteria. In this contribution, we further explore the evolutionary relationships among these three phototrophic lineages by both molecular phylogeny and consideration of phylogenetic marker traits of the superfamily of Rieske/cytb complexes. The combination of these two methods suggests the existence of a "green" clade involving many non-phototrophs in addition to the mentioned RCI-type photosynthetic organisms. Structural and functional idiosyncrasies are (re-)interpreted in the framework of evolutionary biology and more specifically evolutionary bioenergetics.

Arsenite oxidase from Ralstonia sp. 22: characterization of the enzyme and its interaction with soluble cytochromes

We characterized the aro arsenite oxidation system in the novel strain Ralstonia sp. 22, a beta-proteobacterium isolated from soil samples of the Salsigne mine in southern France. The inducible aro system consists of a heterodimeric membrane-associated enzyme reacting with a dedicated soluble cytochrome c(554). Our biochemical results suggest that the weak association of the enzyme to the membrane probably arises from a still unknown interaction partner. Analysis of the phylogeny of the aro gene cluster revealed that it results from a lateral gene transfer from a species closely related to Achromobacter sp. SY8. This constitutes the first clear cut case of such a transfer in the Aro phylogeny. The biochemical study of the enzyme demonstrates that it can accommodate in vitro various cytochromes, two of which, c(552) and c(554,) are from the parent species. Cytochrome c(552) belongs to the sox and not the aro system. Kinetic studies furthermore established that sulfite and sulfide, substrates of the sox system, are both inhibitors of Aro activity. These results reinforce the idea that sulfur and arsenic metabolism are linked.

The small subunit AroB of arsenite oxidase: lessons on the [2Fe-2S] Rieske protein superfamily

Here, we describe the characterization of the [2Fe-2S] clusters of arsenite oxidases from Rhizobium sp. NT-26 and Ralstonia sp. 22. Both reduced Rieske proteins feature EPR signals similar to their homologs from Rieske-cyt b complexes, with g values at 2.027, 1.88, and 1.77. Redox titrations in a range of pH values showed that both [2Fe-2S] centers have constant E(m) values up to pH 8 at approximately +210 mV. Above this pH value, the E(m) values of both centers are pH-dependent, similar to what is observed for the Rieske-cyt b complexes. The redox properties of these two proteins, together with the low E(m) value (+160 mV) of the Alcaligenes faecalis arsenite oxidase Rieske (confirmed herein), are in line with the structural determinants observed in the primary sequences, which have previously been deduced from the study of Rieske-cyt b complexes. Since the published E(m) value of the Chloroflexus aurantiacus Rieske (+100 mV) is in conflict with this sequence analysis, we re-analyzed membrane samples of this organism and obtain a new value (+200 mV). Arsenite oxidase activity was affected by quinols and quinol analogs, which is similar to what is found with the Rieske-cyt b complexes. Together, these results show that the Rieske protein of arsenite oxidase shares numerous properties with its counterpart in the Rieske-cyt b complex. However, two cysteine residues, strictly conserved in the Rieske-cyt b-Rieske and considered to be crucial for its function, are not conserved in the arsenite oxidase counterpart. We discuss the role of these residues.

Menaquinone as pool quinone in a purple bacterium

Purple bacteria have thus far been considered to operate light-driven cyclic electron transfer chains containing ubiquinone (UQ) as liposoluble electron and proton carrier. We show that in the purple gamma-proteobacterium Halorhodospira halophila, menaquinone-8 (MK-8) is the dominant quinone component and that it operates in the Q(B)-site of the photosynthetic reaction center (RC). The redox potentials of the photooxidized pigment in the RC and of the Rieske center of the bc(1) complex are significantly lower (E(m) = +270 mV and +110 mV, respectively) than those determined in other purple bacteria but resemble those determined for species containing MK as pool quinone. These results demonstrate that the photosynthetic cycle in H. halophila is based on MK and not on UQ. This finding together with the unusual organization of genes coding for the bc(1) complex in H. halophila suggests a specific scenario for the evolutionary transition of bioenergetic chains from the low-potential menaquinones to higher-potential UQ in the proteobacterial phylum, most probably induced by rising levels of dioxygen 2.5 billion years ago. This transition appears to necessarily proceed through bioenergetic ambivalence of the respective organisms, that is, to work both on MK- and on UQ-pools. The establishment of the corresponding low- and high-potential chains was accompanied by duplication and redox optimization of the bc(1) complex or at least of its crucial subunit oxidizing quinols from the pool, the Rieske protein. Evolutionary driving forces rationalizing the empirically observed redox tuning of the chain to the quinone pool are discussed.

Was nitric oxide the first deep electron sink?

Evolutionary histories of enzymes involved in chemiosmotic energy conversion indicate that a strongly oxidizing substrate was available to the last universal common ancestor before the divergence of Bacteria and Archaea. According to palaeogeochemical evidence, O(2) was not present beyond trace amounts on the early Earth. Based on recent phylogenetic, enzymatic and geochemical results, we propose that, in the earliest Archaean, nitric oxide (NO) and its derivatives nitrate and nitrite served as strongly oxidizing substrates driving the evolution of a bioenergetic pathway related to modern dissimilatory denitrification. Aerobic respiration emerged later from within this ancestral pathway via adaptation of the enzyme NO reductase to its new substrate, dioxygen.

Risks from GMOs due to horizontal gene transfer

Horizontal gene transfer (HGT) is the stable transfer of genetic material from one organism to another without reproduction or human intervention. Transfer occurs by the passage of donor genetic material across cellular boundaries, followed by heritable incorporation to the genome of the recipient organism. In addition to conjugation, transformation and transduction, other diverse mechanisms of DNA and RNA uptake occur in nature. The genome of almost every organism reveals the footprint of many ancient HGT events. Most commonly, HGT involves the transmission of genes on viruses or mobile genetic elements. HGT first became an issue of public concern in the 1970s through the natural spread of antibiotic resistance genes amongst pathogenic bacteria, and more recently with commercial production of genetically modified (GM) crops. However, the frequency of HGT from plants to other eukaryotes or prokaryotes is extremely low. The frequency of HGT to viruses is potentially greater, but is restricted by stringent selection pressures. In most cases the occurrence of HGT from GM crops to other organisms is expected to be lower than background rates. Therefore, HGT from GM plants poses negligible risks to human health or the environment.

Enzyme phylogenies as markers for the oxidation state of the environment: the case of respiratory arsenate reductase and related enzymes

Background

Phylogenies of certain bioenergetic enzymes have proved to be useful tools for deducing evolutionary ancestry of bioenergetic pathways and their relationship to geochemical parameters of the environment. Our previous phylogenetic analysis of arsenite oxidase, the molybdopterin enzyme responsible for the biological oxidation of arsenite to arsenate, indicated its probable emergence prior to the Archaea/Bacteria split more than 3 billion years ago, in line with the geochemical fact that arsenite was present in biological habitats on the early Earth. Respiratory arsenate reductase (Arr), another molybdopterin enzyme involved in microbial arsenic metabolism, serves as terminal oxidase, and is thus situated at the opposite end of bioenergetic electron transfer chains as compared to arsenite oxidase. The evolutionary history of the Arr-enzyme has not been studied in detail so far.

Results

We performed a genomic search of genes related to arrA coding for the molybdopterin subunit. The multiple alignment of the retrieved sequences served to reconstruct a neighbor-joining phylogeny of Arr and closely related enzymes. Our analysis confirmed the previously proposed proximity of Arr to the cluster of polysulfide/thiosulfate reductases but also unravels a hitherto unrecognized clade even more closely related to Arr. The obtained phylogeny strongly suggests that Arr originated after the Bacteria/Archaea divergence in the domain Bacteria, and was subsequently laterally distributed within this domain. It further more indicates that, as a result of accumulation of arsenate in the environment, an enzyme related to polysulfide reductase and not to arsenite oxidase has evolved into Arr.

Conclusion

These findings are paleogeochemically rationalized by the fact that the accumulation of arsenate over arsenite required the increase in oxidation state of the environment brought about by oxygenic photosynthesis.

Variation of average g values and effective exchange coupling constants among [2Fe-2S] clusters: a density functional theory study of the impact of localization (trapping forces) versus delocalization (double-exchange) as competing factors

A phenomenological model aimed at rationalizing variations in both average g-tensor values (gav identical with 1/3Sigmaigi ) and effective exchange coupling constants Jeff (defined as two-thirds of the energy difference between the S = 3/2 and S = 1/2 spin states) has been derived in order to describe the great variety of magnetic properties exhibited by reduced [2Fe-2S] clusters in proteins. The key quantity in the present analysis is the ratio Delta E/B computed from two competing terms. Delta Ecomprises various effects that result in trapping-site asymmetries: vibronic coupling and the chemical nature (S/N/O) and conformations of the ligands on the one hand and solvation terms, the hydrogen bonding network, etc., on the other. All of these additive terms (in a "bottom-up" approach) favor valence localization of the reducing electron onto one of the two iron sites. In contrast, the B term is the double-exchange term, which favors electronic delocalization. Both gav and Jeff can be expressed as functions of Delta E/ B. We have also shown that electronic localization generally favors small gav and large Jeff values (while the opposite is true for electronic delocalization) in a comparative study of the spectroscopic features of plant-type ferredoxins (Fd's) and Rieske centers (and related mutants). Two other types of problems were particularly challenging. The first of these involved deprotonated Rieske centers and the xanthine oxidase clusters II, which are characterized by very small Jeff values (40-45 cm (-1) with a J S A. S B model) correlated with unusually large gav values (in the range 1.97-2.01) as a result of an antisymmetric exchange coupling mechanism. The second concerned the analogous Fd's from Clostridium pasteurianum (Cp) and Aquifex aeolicus (Aa). Detailed Mössbauer studies of the C56S mutant of the Cp system revealed a mixture of clusters with valence-localized S = 1/2 and valence-delocalized S = 9/2 ground states. We relied on crystallographic structures of wild-type and mutant Aa Fd's in order to explain such a distribution of spin states.

Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

We use a multigene data set (the mitochondrial locus and nine nuclear gene regions) to test phylogenetic relationships in the South American "lava lizards" (genus Microlophus) and describe a strategy for aligning noncoding sequences that accounts for differences in tempo and class of mutational events. We focus on seven nuclear introns that vary in size and frequency of multibase length mutations (i.e., indels) and present a manual alignment strategy that incorporates insertions and deletions (indels) for each intron. Our method is based on mechanistic explanations of intron evolution that does not require a guide tree. We also use a progressive alignment algorithm (Probabilistic Alignment Kit; PRANK) and distinguishes insertions from deletions and avoids the "gapcost" conundrum. We describe an approach to selecting a guide tree purged of ambiguously aligned regions and use this to refine PRANK performance. We show that although manual alignment is successful in finding repeat motifs and the most obvious indels, some regions can only be subjectively aligned, and there are limits to the size and complexity of a data matrix for which this approach can be taken. PRANK alignments identified more parsimony-informative indels while simultaneously increasing nucleotide identity in conserved sequence blocks flanking the indel regions. When comparing manual and PRANK with two widely used methods (CLUSTAL, MUSCLE) for the alignment of the most length-variable intron, only PRANK recovered a tree congruent at deeper nodes with the combined data tree inferred from all nuclear gene regions. We take this concordance as an objective function of alignment quality and present a strongly supported phylogenetic hypothesis for Microlophus relationships. From this hypothesis we show that (1) a coded indel data partition derived from the PRANK alignment contributed significantly to nodal support and (2) the indel data set permitted detection of significant conflict between mitochondrial and nuclear data partitions, which we hypothesize arose from secondary contact of distantly related taxa, followed by hybridization and mtDNA introgression.

Iron-sulfur protein folds, iron-sulfur chemistry, and evolution

An inventory of unique local protein folds around Fe-S clusters has been derived from the analysis of protein structure databases. Nearly 50 such folds have been identified, and over 90% of them harbor low-potential [2Fe-2S](2+,+) or [4Fe-4S](2+,+) clusters. In contrast, high-potential Fe-S clusters, notwithstanding their structural diversity, occur in only three different protein folds. These observations suggest that the extant population of Fe-S protein folds has to a large extent been shaped in the reducing iron- and sulfur-rich environment that is believed to have predominated on this planet until approximately two billion years ago. High-potential active sites are then surmised to be rarer because they emerged later, in a more oxidizing biosphere, in conditions where iron and sulfide had become poorly available, Fe-S clusters were less stable, and in addition faced competition from heme iron and copper active sites. Among the low-potential Fe-S active sites, protein folds hosting [4Fe-4S](2+,+) clusters outnumber those with [2Fe-2S](2+,+) ones by a factor of 3 at least. This is in keeping with the higher chemical stability and versatility of the tetranuclear clusters, compared with the binuclear ones. It is therefore suggested that, at least while novel Fe-S sites are evolving within proteins, the intrinsic chemical stability of the inorganic moiety may be more important than the stabilizing effect of the polypeptide chain. The discovery rate of novel Fe-S-containing protein folds underwent a sharp increase around 1995, and has remained stable to this day. The current trend suggests that the mapping of the Fe-S fold space is not near completion, in agreement with predictions made for protein folds in general. Altogether, the data collected and analyzed here suggest that the extant structural landscape of Fe-S proteins has been shaped to a large extent by primeval geochemical conditions on one hand, and iron-sulfur chemistry on the other.

Multiple sequence alignment: in pursuit of homologous DNA positions

DNA sequence alignment is a prerequisite to virtually all comparative genomic analyses, including the identification of conserved sequence motifs, estimation of evolutionary divergence between sequences, and inference of historical relationships among genes and species. While it is mere common sense that inaccuracies in multiple sequence alignments can have detrimental effects on downstream analyses, it is important to know the extent to which the inferences drawn from these alignments are robust to errors and biases inherent in all sequence alignments. A survey of investigations into strengths and weaknesses of sequence alignments reveals, as expected, that alignment quality is generally poor for two distantly related sequences and can often be improved by adding additional sequences as stepping stones between distantly related species. Errors in sequence alignment are also found to have a significant negative effect on subsequent inference of sequence divergence, phylogenetic trees, and conserved motifs. However, our understanding of alignment biases remains rudimentary, and sequence alignment procedures continue to be used somewhat like benign formatting operations to make sequences equal in length. Because of the central role these alignments now play in our endeavors to establish the tree of life and to identify important parts of genomes through evolutionary functional genomics, we see a need for increased community effort to investigate influences of alignment bias on the accuracy of large-scale comparative genomics.

Mutational analysis of cytochrome b at the ubiquinol oxidation site of yeast complex III

The cytochrome bc1 complex is a dimeric enzyme of the inner mitochondrial membrane that links electron transfer from ubiquinol to cytochrome c by a protonmotive Q cycle mechanism in which ubiquinol is oxidized at one center in the enzyme, referred to as center P, and ubiquinone is rereduced at a second center, referred to as center N. To better understand the mechanism of ubiquinol oxidation, we have examined catalytic activities and pre-steady-state reduction kinetics of yeast cytochrome bc1 complexes with mutations in cytochrome b that we expected would affect oxidation of ubiquinol. We mutated two residues thought to be involved in proton conduction linked to ubiquinol oxidation, Tyr132 and Glu272, and two residues proposed to be involved in docking ubiquinol into the center P pocket, Phe129 and Tyr279. Substitution of Phe129 by lysine or arginine yielded a respiration-deficient phenotype and lipid-dependent catalytic activity. Increased bypass reactions were detectable for both variants, with F129K showing the more severe effects. Substitution with lysine leads to a disturbed coordination of a b heme as deduced from changes in the midpoint potential and the EPR signature. Removal of the aromatic side chain in position Tyr279 lowers the catalytic activity accompanied by a low level of bypass reactions. Pre-steady-state kinetics of the enzymes modified at Glu272 and Tyr132 confirmed the importance of their functional groups for electron transfer. Altered center N kinetics and activation of ubiquinol oxidation by binding of cytochrome c in the Y132F and E272D enzymes indicate long range effects of these mutations.