|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 49, 51554-51560, December 3, 2004
Quantitative Proteomic Analysis of Metabolic Regulation by Copper Ions in Methylococcus capsulatus (Bath)*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Methanotrophs are classified into three distinct types. Type I methanotrophs show disk-like intracellular membranes and use the RuMP (ribulose 5-phosphate) pathway for carbon assimilation. In contrast, Type II methanotrophs possess intracellular membranes in the cell periphery and use the serine pathway to degrade the formaldehyde produced. Methylococcus capsulatus (Bath) is classified as a Type X methanotroph (1), since it shows the physiological properties of both Type I and II methanotrophs, but it develops Type I intracytoplasmic membranes.
Two methane monooxygenases (MMO)1 catalyze the methane oxidation process in methanotrophs, converting methane to methanol. All methanotrophs express a membrane-bound copper-containing particulate MMO (pMMO) (1, 3), while Type X and a few Type II methanotrophic bacteria are capable of producing a second, soluble form (soluble methane monooxygenase, sMMO). Copper ions are known to switch the methane oxidation from sMMO to pMMO. At low copper-to-biomass ratios, the cytoplasmic sMMO is the dominant MMO. When the cells are grown at high copper-to-biomass ratios, pMMO is expressed and produced instead in the plasma membrane.
Because of the unique role of copper ions in regulating the switch between sMMO and pMMO, we have applied cICAT (cleavable isotope coded affinity tag) in combination with off-line two-dimensional LC-MS/MS to quantitatively compare the expression levels of the unique enzymes involved in the featured metabolic pathways under different copper environments. The ICAT technology has been developed for quantitative analysis of proteomic changes in response to cellular perturbations with wide dynamic range and quantitation accuracy (4, 5). More recently, the cICAT reagent has also been developed to eliminate chromatographic isotope effects caused by hydrogen and deuterium (6). Two groups (the National High-Throughput Genome Sequencing Center at the National Yang-Ming University in Taiwan and the Institute for Genomic Research (Rockville, MD) (www.tigr.org/tdb/mdb/mdbinprogress.html)) have been independently making substantial progress toward completing the sequencing of the genome of this methanotroph. These data, although incomplete, have been indispensable to the present proteomic study.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Peptide Separation and PurificationThe peptides were separated by cation exchange chromatography using a 4.6 x 200-mm polysulfoethyl A column (5-µm particles, 300-A pore size; Poly LC, Columbia, MD) at a flow rate of 200 µl/min. Peptides were eluted by a gradient of 0-25% Buffer B over 30 min, followed by 25-100% Buffer B over 20 min (Buffer A: 5 mM K2HPO4, 25% CH3CN, pH = 3.0; Buffer B: 5 mM K2HPO4, 25% CH3CN, 350 mM KCl, pH = 3.0). The elution profile of the cation exchange chromatography determined which fractions were to be analyzed further, and each were individually processed over avidin cartridges (Applied Biosystems), and the affinity tags were cleaved according to the manufacturer's protocol (cICAT kit for protein labeling; Applied Biosystems) to isolate the labeled Cys-containing peptides for micro-capillary liquid chromatography tandem mass spectrometry (µLC-MS/MS) analysis.
µLC-MS/MS AnalysisSamples were loaded using an auto-sampler and sequentially analyzed by µLC-MS/MS. Injections were made on a 1.5 cm x 100 µm trapping and 12 cm x 75 µm separation column packed in-house (Magic C18; Micromass BioResources, Auburn, CA). Peptides were eluted with a linear gradient of 5-40% Buffer B over 180 min at an elution rate of
200 nl/min (Buffer A: 0.1% formic acid in H2O; Buffer B: 0.1% formic acid in acetonitrile). A Hewlett Packard 1100 solvent delivery system with flow splitting (Hewlett Packard, Palo Alto, CA) was used. An LCQ-Deca ion trap mass spectrometer (ThermoFinnigan, San Jose, CA) with an in-house built microspray device was used for all analyses. Peptide fragmentation by collision-induced dissociation was carried out in an automated fashion using the dynamic-exclusion option, and the resultant MS/MS spectra were recorded. The interpretation of the MS/MS data were finally submitted to a suit of software tools for automated data base searching and statistical interpretation of the search results.
Data Analysis and BioinformaticsGenomic sequencing of M. capsulatus (Bath) was performed by the National High-Throughput Genome Sequencing Center at the National Yang-Ming University. The entire shotgun sequencing is not finished; however, we have obtained both the 4x and 8x coverage contigs. EMBOSS 2.8.0 (9), NCBI BLAST 2.2.6 (10), and HMMER 2.3.2 (11) were installed on FreeBSD 4.8 RELEASE with one Intel® Pentium® 4 processor operating in the batch mode. Open reading frames (ORFs) were predicted and translated into protein sequences by the GETORF function of the EMBOSS suite, and their homologues protein names were searched against the Swiss-Prot, TrEMBL, and the NCBI non-redundant data bases with the expectation value of 1E-05 as the criterion. 5485 protein homologues were predicted. Since the homologues of the ORFs were not exactly the ones that gave the best similarity (for example, some hypothetical proteins were predicted by genomic sequencing), all the entries were checked manually.
Automated data base searching using SEQUEST (12) software was performed to identify peptide and protein sequence matches for each recorded MS/MS spectrum. Uninterrupted MS/MS spectra were searched against the M. capsulatus (Bath) genomic sequencing ORF data base. SEQUEST search parameters for light-cICAT-labeled cysteine were set to 277.13, with a +9 differential modification for heavy cICAT-labeled cysteine and +16 for oxidized methionine, mass tolerance 3 Da. Tryptic digested peptides were specified to limit the search results. Identified peptide scores were validated by Peptide Prophet, and peptides that displayed Peptide Prophet scores >0.9 were further analyzed by Protein Prophet to validate identified open reading frames (13, 14). The light to heavy (L/H) ratios resulting from the differential exprssion were calculated and normalized by ASAPRatio (15) (Automated Statistical Analysis on Protein Ratio) program.
Additional function annotation was performed with the NCBI COGNITOR program on the web (www.ncbi.nlm.nih.gov/COG/), and each protein entry was placed into the highest scoring COG group (16). Protein families were classified by HMMER with the Pfam HMM library. Protein domains were analyzed by the PATMATMOTIFS function of EMBOSS 2.8.0 with the PROSITE data base (17). Protein locations were predicted by the PSORT WWW server (18), and the number of transmembrane domains was predicted by the TMAP function of EMBOSS 2.8.0 (19).
| RESULTS |
|---|
|
|
|---|
682 Proteins Revealed in Proteomic StudyProteins expressed in the absence and presence of copper ion, 0 and 30 µM, were labeled with isotopically light (12C: L) and heavy (13C: H) cICAT reagents, respectively. In summary, a total of 682 proteins of M. capsulatus (Bath) were identified with at least one peptide of high confidence value: 429 were from the cytosolic fraction and 253 from the membrane fraction (supplemental Table I). Moreover, it was possible to quantify the change in expression levels of 533 of the 682 identified proteins from the L/H mass ratios of the protein fragments.
87.8% of the identified proteins has significant sequence homology to known proteins in the Swiss-Prot and TrEMBL data bases, including 59 entries (8.7% of the identified proteins) sequenced from M. capsulatus (Bath). The remaining 12.2% identified proteins show very low or no detectable homology to other proteins in the data base and are considered poorly characterized. They may represent unique proteins expressed in M. capsulatus (Bath) or gene products that have not been discovered previously.
To obtain information regarding the subcellular localization of the 682 identified proteins, the identified ORFs were further analyzed by different algorithms. To predict outer and inner membrane proteins in M. capsulatus (Bath), the PSORT algorithm was used. Transmembrane domains (TMDs) were also predicted by the TMAP function of EMBOSS. The analysis showed that the identified proteins were mainly localized in the cytoplasm (412 soluble proteins). Of the 236 proteins predicted to be in the inner membrane fraction, 84 proteins (36%) were predicted to contain two or more transmembrane domains. The fraction of cytosolic proteins with two or more transmembrane domains was 8.7% (36 proteins), indicating that there was some contamination of the cytosolic pool with membrane proteins. The remaining 27 cytosolic proteins were periplasmic, and seven were Gram-negative outer membrane proteins. These results are summarized in Fig. 1.
|
,
, and
subunits of ATP synthase (L/H = 1.27, 1.1, and 0.99, respectively), NADH dehydrogenase C (L/H = 1.09) and D, F, and G (L/H = 1.43) chains, the proton-translocating pyrophosphate synthase (15 TMDs), and the sodium bile acid symporter family protein (10 TMDs, L/H = 0.58). The cation-transporting ATPase (10 TMDs, L/H = 0.84), magnesium-transporting ATPase (E1 E2 family) (8 TMDs), CopF Cu-ATPase (8 TMDs, L/H = 0.66), the O-antigen acetylase (10 TMDs), and the hypothetical protein PA0575 (10 TMDs, L/H = 1.84) were also found, with a large number of transmembrane domains predicted. All unique protein sequences identified were analyzed by the functional annotation program NCBI COGnitor (Fig. 2). 170 proteins (25%), including all the sMMO and pMMO subunits (see below), could not be classified in this manner, because COG relies on sequenced genomes as the classification data base, and there has been no systematic research and genomic sequencing of methanotrophs until the present study. Although TIGR is sequencing the genome of M. capsulatus (Bath) independently, the data available from this source is too limited and inadequate for the large scale proteomic study undertaken here.
|
|
The tungsten-containing aldehyde ferredoxin oxidoreductase (L/H = 0.08) was another protein that was strongly up-regulated by copper. Another strongly induced ORF was identified by BLAST searching as aldehyde dehydrogenase (L/H = 0.23), which is involved in substrate assimilation of a wide variety of aldehydes. The diversity of formaldehyde-related redox proteins may reflect the high amounts of formaldehyde produced at 30 µM copper from methane oxidation by pMMO, due to the high levels of pMMO in the membranes as well as the higher solubility of methane in the intracytoplasmic membranes. The tungsten-containing formate dehydrogenase (L/H = 0/23), the last enzyme in the pathway of methane oxidation, was also identified in the membrane fraction. Since formaldehyde is toxic to the cell, presumably, these enzymes are up-regulated to remove the formaldehyde efficiently. However, the metal content of these "tungsten" enzymes remains to be addressed.
From the genome sequencing data available, we have located all the genes involved in the tricarboxylic acid cycle. Seven of the enzymes, except succinyl-CoA dehydrogenase, were identified in the cICAT study. Succinyl-CoA synthetase
chain (L/H = 0.14) and
chain (L/H = 0.19) (another fragment on another contig was also identified, with a different L/H ratio of 0.43), fumarate hydratase (L/H = 0.31), dihydrolipoamide acetyltranferase (L/H = 0.42), and pyruvate dehydrogenase (the E1 component, L/H = 0.46) were stimulated at 30 µM copper. It has been reported that M. capsulatus (Bath) lacks a complete tricarboxylic acid cycle due to the absence of 2-oxoglutarate dehydrogenase (20); however, the latter has been identified (L/H = 0.59). These results provide further evidence that Type X methylotrophs share some metabolic properties of both Type I and Type II bacteria. Type II methanotrophs use the serine pathway as its major C1 assimilation route with a complete tricarboxylic acid cycle, while Type I methanotrophs do not (1, 20).
It is interesting that some cell wall and capsule synthesis related proteins such as UDP-N-acetylglucosamine acyltransferase (L/H = 0.45), UDP-2,3-diacyl-glucosamine hydrolase (L/H = 0.49), 2-dehydro-3-deoxyphospho-octonate aldolase (L/H = 0.19), UDP-glucose/GDP-mannose dehydrogenase (L/H = 0.27), and 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase (L/H = 0.44) are stimulated at the higher copper ion concentration. With the cell volume increasing to accommodate the intracytoplasmic membranes, it is possible that the synthesis of additional cell wall material is necessary. Methanotrophic bacteria with pMMO also have higher growth yields than those with sMMO (1). The high copper environment is presumably more compatible with the survival of the M. capsulatus (Bath).
Four ion transport proteins are also overexpressed. Three enzymes provide the coenzyme source for other proteins: porphobilinogen deaminase (L/H = 0.39), the fourth enzyme in porphyrin biosynthesis; 6,7-dimethyl-8-ribityllumazine synthase (L/H = 0.43), the enzyme involved in the last step of riboflavin biosynthesis; and quinolinate synthetase (L/H = 0.36), which is involved in NAD biosynthesis.
Hemerythrin (L/H = 0.26) was observed at 30 µM copper in both the cytosolic and membrane fractions. Hemerythrin is an iron-containing oxygen carrier typically found in marine organisms, but the bacterial fusion domain homologue is also known (21). A BLAST search showed that the hemerythrin domain is conserved. We have successfully purified this hemerythrin, and its UV-visible spectrum is identical to hemerythrins from other sources, indicating the presence of iron.2 Interestingly, this is the first bacterial hemerythrin isolated. Since the level of expression of the hemerythrin is significantly enhanced at the higher copper concentration, we surmise that hemerythrin functions as oxygen carrier to transport the dioxygen efficiently throughout the cell body for consumption by the pMMO, particularly at the high levels of pMMO produced at 30 µM copper.
Finally, at 30 µM copper, there are three protein sequences with no significant similarity in the current data base and six hypothetical proteins. The patterns matched by similarity search are summarized in supplemental Table IIb. An ATP/GTP-binding site motif A (P-loop) motif was found in MCB-2007_160 (predicted open reading frame name) (L/H = 0.39); and MCB-2065_178 (L/H = 0.23) possesses three amidation site motifs and one RGD (arginine, glycine, and aspartic acid) tripeptide pattern (cell attachment sequence). Although their exact functions are unclear, these proteins must be quite important, because their expressions are stimulated at the higher copper concentration.
Proteins with Higher Expression Levels at 0 µM CopperFor comparison, we highlight here the proteins that were expressed at higher levels in the absence of copper. The 68 proteins with L/H ratio >2 appear to affect different parts of the cell physiology than the proteins with ratio <0.5. They tend to be proteins related to post-translational modification, protein turnover, chaperones and signal transduction mechanisms (Fig. 2 and supplemental Table IIIa). In fact, no enzymes related to carbohydrate metabolism were found except sMMO (MmoX, MmoY, MmoB, MmoC, with L/H = 7.23, 4.36, 42.96, and 14.6, respectively; MmoD was also identified, but its L/H ratio could not be determined) and glucokinase (L/H = 2.62). The former is responsible for the first step of methane oxidation and has already been well studied (22). Since M. capsulatus (Bath) does not rely on hexose as the carbon source, and there are no reports supporting the ability of M. capsulatus (Bath) to grow on glucose, the function of glucokinase might be involved in glycogen biosynthesis or gluconeogenesis.
Two of the four gene translations (MmoG, MmoS, MmoR, and MmoQ) in the sMMO cluster were also found. MmoG (L/H = 18.71) is a 60-kDa chaperonin homologue involved in the assembly of the sMMO complex, while MmoS (L/H = 0.84) is a histidine kinase hypothesized to function as a copper sensor (23). No significant change was observed in the expression level of the histidine kinase.
Some cell signaling and regulatory proteins were highly expressed when the cells were grown without added copper, including the transcriptional regulatory protein zraR (L/H = 3.23), sensory box protein (L/H = 2.44), LacZ expression regulator (L/H = 2.45), and the probable universal stress protein homologue (L/H = 2.42). Two proteins that control cell division were identified, including the cell division inhibitor MinD (L/H = 2.27) and the cell division protease ftsH (L/H = 4.02). Also observed was the peptide methionine sulfoxide reductase (L/H = 15.22) with its SelR domain, which repairs methionine sulfoxide damaged peptides by reducing them to form methionine in conjunction with reductants such as thioredoxin, thus protecting proteins against oxidative stress. Other redox proteins such as glutathione S-transferase (L/H = 2.75) involved in glutathione-dependent reactions and the seven-iron ferredoxin (L/H = 9.81) were also detected. Two molecular chaperone-related proteins responsible for folding-corrections were identified.
Several key enzymes involved in metabolism were identified with enhanced abundances at 0 µM copper. The formylmethanofuran hydrolase (L/H = 2.18) is responsible for the conversion of formyl-H4MPT (tetrahydromethanopterin) to formic acid in the H4MPT pathway (24). GTP cyclohydrolase I (L/H = 4.02) is involved in the first step of the biosynthesis of tetrahydrofolate, which is a C1 transfer intermediate found in bacteria, some archaea and eukaryotes (24). The ATP phosphoribosyltransferase regulatory subunit (L/H = 8.67), imidazoleglycerol phosphate synthase subunit hisF1 (L/H = 2.81) and histidinol dehydrogenase (L/H = 19.95), are enzymes involved in the first, fifth, and the last steps of histidine biosynthesis, respectively. Both denitrification enzymes were observed: nitrate reductase (large subunit, L/H = 2.88; small subunit, L/H = 2.22) and nitrite reductase (L/H = 3.84).
MCB-2060_183 (L/H = 2.15) was BLAST searched as the macrophage migration inhibitory factor found in various eukaryotic species (25), but the role of this homologue in bacteria has not been reported and remains to be elucidated in future. We have also observed the surface-associated protein precursor (MopE) (L/H = 6.3) when the cells were grown at 0 µM copper concentration. The C-terminal of this protein has been reported to be secreted by M. capsulatus (Bath) into the growth media, with the possibility of some catalytic and/or copper binding function (26), but this protein was not found in enriched copper-containing environments. Peptidoglycan-associated lipoprotein (L/H = 3.39) is thought to play a role in Gram-negative bacterial envelope integrity. Whether or not copper also controls the outer membrane composition is unclear, but this presents an interesting issue.
Finally, we found 16 proteins with no BLAST hits as well as seven hypothetical proteins whose expression levels were enhanced without copper added to the growth medium (see supplemental Table IIIb). The most interesting one is MCB-2041_433 (L/H = 2.05), which consists of one amidation site motif, one ATP/GTP-binding site motif A (P-loop), and 13 cytochrome c motifs. MCB-2043_616 (L/H = 4.92) also shows three amidation sites and one cytochrome c motif. Their function(s) should be related to specific redox reactions when the organism is grown at low copper-to-biomass ratios.
| DISCUSSION |
|---|
|
|
|---|
The C1 Metabolic PathwaysM. capsulatus (Bath) assimilates methane as its single carbon source. The methanol oxidation protein (MoxJ) (five TMDs, L/H = 1.25), the gene of which is located between the genes for the two methanol dehydrogenase subunits, might be involved in assembly of the active methanol dehydrogenase and/or its cofactor pyrrolquinoline quinine in the periplasm (28). Formaldehyde is the most important metabolic intermediate in the metabolism of these methanotrophs, because four metabolic pathways, namely, direct oxidation of formaldehyde, the RuMP pathway, the serine pathway, and the H4MPT pathway, all branch from this metabolite (Fig. 4).
|
Both ribulose-1,5-bisphosphate carboxylase oxygenase (Rubisco) large chain (L/H = 0.8) and phosphoribulokinase (L/H = 0.83) are key enzymes of the RuBP pathway, but M. capsulatus (Bath) does not seem to be able to grow autotrophically (29). The other RuBP pathway enzymes could be found in the FBP aldolase variant of the RuMP pathway. This observation is interesting because it might explain why M. capsulatus (Bath) possesses both cleavage phase variants in the RuMP cycle. It is possible that the RuMP pathway is the evolutionary precursor of the RuBP pathway (30). In methanotrophs, the only carbon source is methane rather than glucose. However, the RuMP pathway provides the fundamental framework for carbohydrate metabolism. All metabolites are assimilated to pyruvate, and the tricarboxylic acid cycle is used to complete the process. So, we can assume that glycolysis is immersed in the entire carbohydrate metabolism network.
It is reported that M. capsulatus (Bath) has an incomplete serine pathway, and on this basis, it has been classified as a Type X methanotroph (1, 20). However, we have located all the genes of the serine pathway and identified all the proteins as well in the present cICAT study, except glycerate kinase. There is supporting evidence that the serine pathway is connected to the RuBP pathway (20), since phosphoglycolate could be generated by the oxygenase activity of Rubisco, and further converted by phosphoglycolate phosphatase (L/H = 1.62) to glycolate, which is then oxidized to glyoxylate by glyoxylate reductase (L/H = 0.51). The ORF for the glyoxylate reductase is downstream of the neighboring gene of serine-glyoxylate aminotransferase.
Finally, all enzymes of the H4MPT pathway were identified. This is the latest C1 transfer pathway reported in proteobacteria since its original discovery in methanogenic archaeon (24, 31, 32). According to our genomic sequencing, the formylmethanofuran hydrolase (L/H = 2.18 and 1.57 for subunit A and B, respectively) and formylmethanofuran-H4MPT formyltransferase (L/H = 1.41) gene cluster is located next to the pmoC2 of the second pmo gene cluster (33), while the genes of the formaldehyde-activating enzyme (L/H = 0.83), putative tetrahydromethanopterin biosynthesis protein (BLAST searched as Swiss-Prot number AAR90367 [GenBank] , and methenyl-H4MPT cyclohydrolase (L/H = 0.93) form another independent cluster in the genome. Since the "formose reaction" was discovered a century and half ago (30), formaldehyde-related compounds have been hypothesized to play a role in early chemical evolution (34-36), and there has been considerable interest in organisms that utilize C1 compounds in biological evolution. The present genomic and proteomic study has allowed us to integrate the four metabolic pathways in methanotrophic carbon assimilation.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Tables 1-3. ![]()

To whom correspondence should be addressed: Inst. of Chemistry, Academia Sinica, 128 Academia Rd. Section 2, Nangkang 115, Taipei, Taiwan. Tel./Fax: 886-2-2789-8654; E-mail: chans{at}chem.sinica.edu.tw.
1 The abbreviations used are: MMO, methane monooxygenase; pMMO, particulate methane monooxygenase; sMMO, soluble methane monooxygenase; BLAST, basic local alignment search tool; C1, one carbon; ICAT, isotope-coded affinity tag; cICAT, cleavable isotope-coded affinity tag; COG, clusters of orthologous groups of proteins; FBP, fructose 1,6-bisphosphate; H4MPT, tetrahydromethanopterin; KDPG, 2-keto-3-deoxy-6-phosphogluconate; L/H ratio, light-to-heavy ratio; ORF, open reading frame; RuMP, ribulose 5-phosphate; RuBP, ribulose bisphosphate; TMD, transmembrane domain; µLC-MS/MS, micro-capillary liquid chromatography tandem mass spectrometry; contig, group of overlapping clones. ![]()
2 W.-C. Kao and S. I. Chan, unpublished result. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. M. R. Lacerda and K. F. Reardon Environmental proteomics: applications of proteome profiling in environmental microbiology and biotechnology Brief Funct Genomic Proteomic, March 11, 2009; (2009) elp005v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. R. McDonald, L. Bodrossy, Y. Chen, and J. C. Murrell Molecular Ecology Techniques for the Study of Aerobic Methanotrophs Appl. Envir. Microbiol., March 1, 2008; 74(5): 1305 - 1315. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |