Network analysis of metabolic enzyme evolution in Escherichia coli

BACKGROUND: The two most common models for the evolution of metabolism are the patchwork evolution model, where enzymes are thought to diverge from broad to narrow substrate specificity, and the retrograde evolution model, according to which enzymes evolve in response to substrate depletion. Analysis of the distribution of homologous enzyme pairs in the metabolic network can shed light on the respective importance of the two models. We here investigate the evolution of the metabolism in E. coli viewed as a single network using EcoCyc. RESULTS: Sequence comparison between all enzyme pairs was performed and the minimal path length (MPL) between all enzyme pairs was determined. We find a strong over-representation of homologous enzymes at MPL 1. We show that the functionally similar and functionally undetermined enzyme pairs are responsible for most of the over-representation of homologous enzyme pairs at MPL 1. CONCLUSIONS: The retrograde evolution model predicts that homologous enzymes pairs are at short metabolic distances from each other. In general agreement with previous studies we find that homologous enzymes occur close to each other in the network more often than expected by chance, which lends some support to the retrograde evolution model. However, we show that the homologous enzyme pairs which may have evolved through retrograde evolution, namely the pairs that are functionally dissimilar, show a weaker over-representation at MPL 1 than the functionally similar enzyme pairs. Our study indicates that, while the retrograde evolution model may have played a small part, the patchwork evolution model is the predominant process of metabolic enzyme evolution

GAPDOCK: A genetic algorithm approach to protein docking in CAPRI round 1

As part of the first Critical Assessment of PRotein Interactions, round 1, we predict the structure of two protein-protein complexes, by using a genetic algorithm, GAPDOCK, in combination with surface complementarity, buried surface area, biochemical information, and human intervention. Among the five models submitted for target 1, HPr phosphocarrier protein (B. subtilis) and the hexameric HPr kinase (L. lactis), the best correctly predicts 17 of 52 interprotein contacts, whereas for target 2, bovine rotavirus VP6 protein-monoclonal antibody, the best model predicts 27 of 52 correct contacts. Given the difficult nature of the targets, these predictions are very encouraging and compare well with those obtained by other methods. Nevertheless, it is clear that there is a need for improved methods for distinguishing between correct and plausible but incorrect complexes. Proteins 2003;52:10-14

Preferential attachment in the evolution of metabolic networks

BACKGROUND: Many biological networks show some characteristics of scale-free networks. Scale-free networks can evolve through preferential attachment where new nodes are preferentially attached to well connected nodes. In networks which have evolved through preferential attachment older nodes should have a higher average connectivity than younger nodes. Here we have investigated preferential attachment in the context of metabolic networks. RESULTS: The connectivities of the enzymes in the metabolic network of Escherichia coli were determined and representatives for these enzymes were located in 11 eukaryotes, 17 archaea and 46 bacteria. E. coli enzymes which have representatives in eukaryotes have a higher average connectivity while enzymes which are represented only in the prokaryotes, and especially the enzymes only present in βγ-proteobacteria, have lower connectivities than expected by chance. Interestingly, the enzymes which have been proposed as candidates for horizontal gene transfer have a higher average connectivity than the other enzymes. Furthermore, It was found that new edges are added to the highly connected enzymes at a faster rate than to enzymes with low connectivities which is consistent with preferential attachment. CONCLUSION: Here, we have found indications of preferential attachment in the metabolic network of E. coli. A possible biological explanation for preferential attachment growth of metabolic networks is that novel enzymes created through gene duplication maintain some of the compounds involved in the original reaction, throughout its future evolution. In addition, we found that enzymes which are candidates for horizontal gene transfer have a higher average connectivity than other enzymes. This indicates that while new enzymes are attached preferentially to highly connected enzymes, these highly connected enzymes have sometimes been introduced into the E. coli genome by horizontal gene transfer. We speculate that E. coli has adjusted its metabolic network to a changing environment by replacing the relatively central enzymes for better adapted orthologs from other prokaryotic species

Diversity in Free Energy Landscape of Proteins with the Same Native Topology

In order to elucidate the role of the native state topology and the stability of subdomains in protein folding, we investigate free energy landscape of human lysozyme, which is composed of two subdomains, by Monte Carlo simulations. A realistic lattice model with Go-like interaction is used. We take the relative interaction strength (stability, in other word) of two subdomains as a variable parameter and study the folding process. A variety of folding process is observed and we obtained a phase diagram of folding in terms of temperature and the relative stability. Experimentally-observed diversity in folding process of c-type lysozimes is thus understood as a consequence of the difference in the relative stability of subdomains.Comment: 4 pages, 3 figure

The membrane bound N-terminal domain of human adenosine diphosphate ribosylation factor-1 (ARF1)

The small G protein adenosine diphosphate ribosylation factor-1 (ARF1) is activated by cell membrane binding of a self-folding N-terminal domain. We present a model of the human ARF1 N-terminal peptide in planar lipid bilayers, determined from neutron lamellar diffraction and circular dichroism data with molecular modelling. This amphipathic domain lies at a shallow membrane depth, ideal for regulation of the ARF1 bio-timer by rapid, reversible membrane binding. The helical region does not elongate upon membrane binding, leaving the connecting flexible linker region's length unchanged.NRC publication: N

Co-populated Conformational Ensembles of β(2)-Microglobulin Uncovered Quantitatively by Electrospray Ionization Mass Spectrometry

Ordered assembly of monomeric human β(2)-microglobulin (β(2)m) into amyloid fibrils is associated with the disorder hemodialysis-related amyloidosis. Previously, we have shown that under acidic conditions (pH <5.0 at 37 °C), wild-type β(2)m assembles spontaneously into fibrils with different morphologies. Under these conditions, β(2)m populates a number of different conformational states in vitro. However, this equilibrium mixture of conformationally different species is difficult to resolve using ensemble techniques such as nuclear magnetic resonance or circular dichroism. Here we use electrospray ionization mass spectrometry to resolve different species of β(2)m populated between pH 6.0 and 2.0. We show that by linear deconvolution of the charge state distributions, the extent to which each conformational ensemble is populated throughout the pH range can be determined and quantified. Thus, at pH 3.6, conditions under which short fibrils are produced, the conformational ensemble is dominated by a charge state distribution centered on the 9+ ions. By contrast, under more acidic conditions (pH 2.6), where long straight fibrils are formed, the charge state distribution is dominated by the 10+ and 11+ ions. The data are reinforced by investigations on two variants of β(2)m (V9A and F30A) that have reduced stability to pH denaturation and show changes in the pH dependence of the charge state distribution that correlate with the decrease in stability measured by tryptophan fluorescence. The data highlight the potential of electrospray ionization mass spectrometry to resolve and quantify complex mixtures of different conformational species, one or more of which may be important in the formation of amyloid

The isolation and characterisation of temperature-dependent Ricin A chain molecules in Saccharomyces cerevisiae.

Ricin is a heterodimeric plant protein that is potently toxic to mammalian cells. Toxicity results from the catalytic depurination of eukaryotic ribosomes by ricin A chain (RTA) that follows toxin endocytosis to, and translocation across, the endoplasmic reticulum (ER) membrane. To ultimately identify proteins required for these later steps in the entry process, it will be useful to express the catalytic subunit within the ER of yeast cells in a manner that initially permits cell growth. A subsequent switch in conditions to provoke innate toxin action would permit only those strains containing defects in genes normally essential for toxin retro-translocation, refolding or degradation to survive. As a route to such a screen, several RTA mutants with reduced catalytic activity have previously been isolated. Here we report the use of Saccharomyces cerevisiae to isolate temperaturedependent mutants of endoplasmic reticulum-targeted RTA. Two such toxin mutants with opposing phenotypes were isolated. One mutant RTA (RTAF108L/L151P) allowed the yeast cells that express it to grow at 37°C while the same cells did not grow at 23ºC. Both mutations were required for temperature-dependent growth. The second toxin mutant (RTAE177D) allowed cells to grow at 23°C but not at 37°C. Interestingly, RTAE177D has been previously reported to have reduced catalytic activity, but this is the first demonstration of a temperature-sensitive phenotype. To provide a more detailed characterisation of these mutants we have investigated their N-glycosylation, stability, catalytic activity and, where appropriate, a three dimensional structure. The potential utility of these mutants is discussed

The Interaction of Atovaquone with the P. carinii Cytochrome bc 1 Complex

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/92453/1/j.1550-7408.2001.tb00505.x.pd

Structure of Tagatose-1,6-bisphosphate Aldolase. Insight into chiral discrimination, mechanism, and specificity of class II aldolases

Tagatose-1,6-bisphosphate aldolase (TBPA) is a tetrameric class II aldolase that catalyzes the reversible condensation of dihydroxyacetone phosphate with glyceraldehyde 3-phosphate to produce tagatose 1,6-bisphosphate. The high resolution (1.45 Å) crystal structure of the Escherichia coli enzyme, encoded by the agaY gene, complexed with phosphoglycolohydroxamate (PGH) has been determined. Two subunits comprise the asymmetric unit, and a crystallographic 2-fold axis generates the functional tetramer. A complex network of hydrogen bonds position side chains in the active site that is occupied by two cations. An unusual Na(+) binding site is created using a interaction with Tyr(183) in addition to five oxygen ligands. The catalytic Zn(2+) is five-coordinate using three histidine nitrogens and two PGH oxygens. Comparisons of TBPA with the related fructose-1,6-bisphosphate aldolase (FBPA) identifies common features with implications for the mechanism. Because the major product of the condensation catalyzed by the enzymes differs in the chirality at a single position, models of FBPA and TBPA with their cognate bisphosphate products provide insight into chiral discrimination by these aldolases. The TBPA active site is more open on one side than FBPA, and this contributes to a less specific enzyme. The availability of more space and a wider range of aldehyde partners used by TBPA together with the highly specific nature of FBPA suggest that TBPA might be a preferred enzyme to modify for use in biotransformation chemistry

Suicide inhibition of alpha-oxamine synthases:structures of the covalent adducts of 8-amino-7-oxononanoate synthase with trifluoroalanine

The suicide inhibition of the α-oxamine synthases by the substrate analog, L-trifluoroalanine was investigated. The inhibition resulted in the formation of a complex with loss of all three fluorine atoms. Decarboxylation and loss of fluoride occurred immediately after aldimine formation. The inherent flexibility could allow the difluorinated intermediate complex to adopt a suitable conformation. Decarboxylation in the normal mechanism occurs after formation of the ketoacid intermediate.link_to_subscribed_fulltex