Regulation of Cell Proliferation and Nrf2-Mediated Antioxidant Defense: Conservation of Keap1 Cysteines and Nrf2 Binding Site in the Context of the Evolution of KLHL Family.

Keap1 (Kelch-like ECH-associated protein 1) is one of the major negative regulators of the transcription factor Nrf2 (nuclear factor erythroid-2-related factor 2), which induces the expression of numerous proteins defending the cell against different stress conditions. Keap1 is generally negatively regulated by post-translational modification (mostly via its cysteine residues) and interaction with other proteins that compete with Nrf2 for binding. Cysteine residues in Keap1 have different effects on protein regulation, as basic residues (Lys, Arg, and His) in close proximity to them increase cysteine modification potential. In this paper, we present an evolutionary analysis of residues involved in both mechanisms of Keap1 regulation in the broader context of the KLHL protein family in vertebrates. We identified the typical domain structure of the KLHL protein family in several proteins outside of this family (namely in KBTBD proteins 2, 3, 4, 6, 7, 8, 12 and 14). We found several cysteines that are flanked by basic residues (namely, C14, C38, C151, C226, C241, C273, C288, C297, C319, and C613) and, therefore, may be considered more susceptible to regulatory modification. The Nrf2 binding site is completely conserved in Keap1 in vertebrates but is absent or located in nonaligned DA and BC loops of the Kelch domain within the KLHL family. The development of specific substrate binding regions could be an evolutionary factor of diversification in the KLHL protein family.

Common Patterns of Hydrolysis Initiation in P-loop Fold Nucleoside Triphosphatases.

The P-loop fold nucleoside triphosphate (NTP) hydrolases (also known as Walker NTPases) function as ATPases, GTPases, and ATP synthases, are often of medical importance, and represent one of the largest and evolutionarily oldest families of enzymes. There is still no consensus on their catalytic mechanism. To clarify this, we performed the first comparative structural analysis of more than 3100 structures of P-loop NTPases that contain bound substrate Mg-NTPs or their analogues. We proceeded on the assumption that structural features common to these P-loop NTPases may be essential for catalysis. Our results are presented in two articles. Here, in the first, we consider the structural elements that stimulate hydrolysis. Upon interaction of P-loop NTPases with their cognate activating partners (RNA/DNA/protein domains), specific stimulatory moieties, usually Arg or Lys residues, are inserted into the catalytic site and initiate the cleavage of gamma phosphate. By analyzing a plethora of structures, we found that the only shared feature was the mechanistic interaction of stimulators with the oxygen atoms of gamma-phosphate group, capable of causing its rotation. One of the oxygen atoms of gamma phosphate coordinates the cofactor Mg ion. The rotation must pull this oxygen atom away from the Mg ion. This rearrangement should affect the properties of the other Mg ligands and may initiate hydrolysis according to the mechanism elaborated in the second article.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity.

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Identification of branched-chain amino acid aminotransferases active towards (R)-(+)-1-phenylethylamine among PLP fold type IV transaminases.

New class IV transaminases with activity towards L-Leu, which is typical of branched-chain amino acid aminotransferases (BCAT), and with activity towards (R)-(+)-1-phenylethylamine ((R)-PEA), which is typical of (R)-selective (R)-amine:pyruvate transaminases, were identified by bioinformatics analysis, obtained in recombinant form, and analyzed. The values of catalytic activities in the reaction with L-Leu and (R)-PEA are comparable to those measured for characteristic transaminases with the corresponding specificity. Earlier, (R)-selective class IV transaminases were found to be active, apart from (R)-PEA, only with some other (R)-primary amines and D-amino acids. Sequences encoding new transaminases with mixed type of activity were found by searching for changes in the conserved motifs of sequences of BCAT by different bioinformatics tools.

COGcollator: a web server for analysis of distant relationships between homologous protein families.

BACKGROUND: The Clusters of Orthologous Groups (COGs) of proteins systematize evolutionary related proteins into specific groups with similar functions. However, the available databases do not provide means to assess the extent of similarity between the COGs. AIM: We intended to provide a method for identification and visualization of evolutionary relationships between the COGs, as well as a respective web server. RESULTS: Here we introduce the COGcollator, a web tool for identification of evolutionarily related COGs and their further analysis. We demonstrate the utility of this tool by identifying the COGs that contain distant homologs of (i) the catalytic subunit of bacterial rotary membrane ATP synthases and (ii) the DNA/RNA helicases of the superfamily 1. REVIEWERS: This article was reviewed by Drs. Igor N. Berezovsky, Igor Zhulin and Yuri Wolf.

COGNAT: a web server for comparative analysis of genomic neighborhoods.

BACKGROUND: In prokaryotic genomes, functionally coupled genes can be organized in conserved gene clusters enabling their coordinated regulation. Such clusters could contain one or several operons, which are groups of co-transcribed genes. Those genes that evolved from a common ancestral gene by speciation (i.e. orthologs) are expected to have similar genomic neighborhoods in different organisms, whereas those copies of the gene that are responsible for dissimilar functions (i.e. paralogs) could be found in dissimilar genomic contexts. Comparative analysis of genomic neighborhoods facilitates the prediction of co-regulated genes and helps to discern different functions in large protein families. AIM: We intended, building on the attribution of gene sequences to the clusters of orthologous groups of proteins (COGs), to provide a method for visualization and comparative analysis of genomic neighborhoods of evolutionary related genes, as well as a respective web server. RESULTS: Here we introduce the COmparative Gene Neighborhoods Analysis Tool (COGNAT), a web server for comparative analysis of genomic neighborhoods. The tool is based on the COG database, as well as the Pfam protein families database. As an example, we show the utility of COGNAT in identifying a new type of membrane protein complex that is formed by paralog(s) of one of the membrane subunits of the NADH:quinone oxidoreductase of type 1 (COG1009) and a cytoplasmic protein of unknown function (COG3002). REVIEWERS: This article was reviewed by Drs. Igor Zhulin, Uri Gophna and Igor Rogozin.

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 b6 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 bc1 -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 b6 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.

Modeling of interaction between cytochrome c and the WD domains of Apaf-1: bifurcated salt bridges underlying apoptosome assembly.

BACKGROUND: Binding of cytochrome c, released from the damaged mitochondria, to the apoptotic protease activating factor 1 (Apaf-1) is a key event in the apoptotic signaling cascade. The binding triggers a major domain rearrangement in Apaf-1, which leads to oligomerization of Apaf-1/cytochrome c complexes into an apoptosome. Despite the availability of crystal structures of cytochrome c and Apaf-1 and cryo-electron microscopy models of the entire apoptosome, the binding mode of cytochrome c to Apaf-1, as well as the nature of the amino acid residues of Apaf-1 involved remain obscure. RESULTS: We investigated the interaction between cytochrome c and Apaf-1 by combining several modeling approaches. We have applied protein-protein docking and energy minimization, evaluated the resulting models of the Apaf-1/cytochrome c complex, and carried out a further analysis by means of molecular dynamics simulations. We ended up with a single model structure where all the lysine residues of cytochrome c that are known as functionally-relevant were involved in forming salt bridges with acidic residues of Apaf-1. This model has revealed three distinctive bifurcated salt bridges, each involving a single lysine residue of cytochrome c and two neighboring acidic resides of Apaf-1. Salt bridge-forming amino acids of Apaf-1 showed a clear evolutionary pattern within Metazoa, with pairs of acidic residues of Apaf-1, involved in bifurcated salt bridges, reaching their highest numbers in the sequences of vertebrates, in which the cytochrome c-mediated mechanism of apoptosome formation seems to be typical. CONCLUSIONS: The reported model of an Apaf-1/cytochrome c complex provides insights in the nature of protein-protein interactions which are hard to observe in crystallographic or electron microscopy studies. Bifurcated salt bridges can be expected to be stronger than simple salt bridges, and their formation might promote the conformational change of Apaf-1, leading to the formation of an apoptosome. Combination of structural and sequence analyses provides hints on the evolution of the cytochrome c-mediated apoptosis.

Phylogenomic reconstruction of archaeal fatty acid metabolism.

While certain archaea appear to synthesize and/or metabolize fatty acids, the respective pathways still remain obscure. By analysing the genomic distribution of the key lipid-related enzymes, we were able to identify the likely components of the archaeal pathway of fatty acid metabolism, namely, a combination of the enzymes of bacterial-type ß-oxidation of fatty acids [acyl-coenzyme A (CoA) dehydrogenase, enoyl-CoA hydratase and 3-hydroxyacyl-CoA dehydrogenase] with paralogs of the archaeal acetyl-CoA C-acetyltransferase, an enzyme of the mevalonate biosynthesis pathway. These three ß-oxidation enzymes working in the reverse direction could potentially catalyse biosynthesis of fatty acids, with paralogs of acetyl-CoA C-acetyltransferase performing addition of C2 fragments. The presence in archaea of the genes for energy-transducing membrane enzyme complexes, such as cytochrome bc complex, cytochrome c oxidase and diverse rhodopsins, was found to correlate with the presence of the proposed system of fatty acid biosynthesis. We speculate that because these membrane complexes functionally depend on fatty acid chains, their genes could have been acquired via lateral gene transfer from bacteria only by those archaea that already possessed a system of fatty acid biosynthesis. The proposed pathway of archaeal fatty acid metabolism operates in extreme conditions and therefore might be of interest in the context of biofuel production and other industrial applications.

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.

Open questions on the origin of life at anoxic geothermal fields.

We have recently reconstructed the 'hatcheries' of the first cells by combining geochemical analysis with phylogenomic scrutiny of the inorganic ion requirements of universal components of modern cells (Mulkidjanian et al. Proc Natl Acad Sci U S A 109:E821-830, 2012). These ubiquitous, and by inference primordial, proteins and functional systems show affinity to and functional requirement for K⁺, Zn²âº, Mn²âº, and phosphate. Thus, protocells must have evolved in habitats with a high K⁺/Na⁺ ratio and relatively high concentrations of Zn, Mn and phosphorous compounds. Geochemical reconstruction shows that the ionic composition conducive to the origin of cells could not have existed in marine settings but is compatible with emissions of vapor-dominated zones of inland geothermal systems. Under an anoxic, CO2-dominated atmosphere, the ionic composition of pools of cool, condensed vapor at anoxic geothermal fields would resemble the internal milieu of modern cells. Such pools would be lined with porous silicate minerals mixed with metal sulfides and enriched in K⁺ ions and phosphorous compounds. Here we address some questions that have appeared in print after the publication of our anoxic geothermal field scenario. We argue that anoxic geothermal fields, which were identified as likely cradles of life by using a top-down approach and phylogenomics analysis, could provide geochemical conditions similar to those which were suggested as most conducive for the emergence of life by the chemists who pursuit the complementary bottom-up strategy.

The role of energy in the emergence of biology from chemistry.

Any scenario of the transition from chemistry to biology should include an "energy module" because life can exist only when supported by energy flow(s). We addressed the problem of primordial energetics by combining physico-chemical considerations with phylogenomic analysis. We propose that the first replicators could use abiotically formed, exceptionally photostable activated cyclic nucleotides both as building blocks and as the main energy source. Nucleoside triphosphates could replace cyclic nucleotides as the principal energy-rich compounds at the stage of the first cells, presumably because the metal chelates of nucleoside triphosphates penetrated membranes much better than the respective metal complexes of nucleoside monophosphates. The ability to exploit natural energy flows for biogenic production of energy-rich molecules could evolve only gradually, after the emergence of sophisticated enzymes and ion-tight membranes. We argue that, in the course of evolution, sodium-dependent membrane energetics preceded the proton-based energetics which evolved independently in bacteria and archaea.

Origin of first cells at terrestrial, anoxic geothermal fields.

All cells contain much more potassium, phosphate, and transition metals than modern (or reconstructed primeval) oceans, lakes, or rivers. Cells maintain ion gradients by using sophisticated, energy-dependent membrane enzymes (membrane pumps) that are embedded in elaborate ion-tight membranes. The first cells could possess neither ion-tight membranes nor membrane pumps, so the concentrations of small inorganic molecules and ions within protocells and in their environment would equilibrate. Hence, the ion composition of modern cells might reflect the inorganic ion composition of the habitats of protocells. We attempted to reconstruct the "hatcheries" of the first cells by combining geochemical analysis with phylogenomic scrutiny of the inorganic ion requirements of universal components of modern cells. These ubiquitous, and by inference primordial, proteins and functional systems show affinity to and functional requirement for K(+), Zn(2+), Mn(2+), and phosphate. Thus, protocells must have evolved in habitats with a high K(+)/Na(+) ratio and relatively high concentrations of Zn, Mn, and phosphorous compounds. Geochemical reconstruction shows that the ionic composition conducive to the origin of cells could not have existed in marine settings but is compatible with emissions of vapor-dominated zones of inland geothermal systems. Under the anoxic, CO(2)-dominated primordial atmosphere, the chemistry of basins at geothermal fields would resemble the internal milieu of modern cells. The precellular stages of evolution might have transpired in shallow ponds of condensed and cooled geothermal vapor that were lined with porous silicate minerals mixed with metal sulfides and enriched in K(+), Zn(2+), and phosphorous compounds.

Characterization of the N-ATPase, a distinct, laterally transferred Na+-translocating form of the bacterial F-type membrane ATPase.

UNLABELLED: An analysis of the distribution of the Na(+)-translocating ATPases/ATP synthases among microbial genomes identified an atypical form of the F(1)F(o)-type ATPase that is present in the archaea Methanosarcina barkeri and M. acetivorans, in a number of phylogenetically diverse marine and halotolerant bacteria and in pathogens Burkholderia spp. In complete genomes, representatives of this form (referred to here as N-ATPase) are always present as second copies, in addition to the typical proton-translocating ATP synthases. The N-ATPase is encoded by a highly conserved atpDCQRBEFAG operon and its subunits cluster separately from the equivalent subunits of the typical F-type ATPases. N-ATPase c subunits carry a full set of sodium-binding residues, indicating that most of these enzymes are Na(+)-translocating ATPases that likely confer on their hosts the ability to extrude Na(+) ions. Other distinctive properties of the N-ATPase operons include the absence of the delta subunit from its cytoplasmic sector and the presence of two additional membrane subunits, AtpQ (formerly gene 1) and AtpR (formerly gene X). We argue that N-ATPases are an early-diverging branch of membrane ATPases that, similarly to the eukaryotic V-type ATPases, do not synthesize ATP. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.