If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Department of Microbiology, Immunology and Molecular Genetics, University of California, Los Angeles, Los Angeles, California, USAUCLA-DOE Institute, University of California, Los Angeles, California, USAUCLA Molecular Biology Institute, University of California, Los Angeles, California, USA
Department of Chemistry and Biochemistry, University of California, Los Angeles, California, USAUCLA-DOE Institute, University of California, Los Angeles, California, USAUCLA Molecular Biology Institute, University of California, Los Angeles, California, USADepartment of Biological Chemistry, David Geffen School of Medicine, University of California, Los Angeles, California, USA
Department of Chemistry and Biochemistry, University of California, Los Angeles, California, USAUCLA-DOE Institute, University of California, Los Angeles, California, USAUCLA Molecular Biology Institute, University of California, Los Angeles, California, USA
Syntrophus aciditrophicus is a model syntrophic bacterium that degrades fatty and aromatic acids into acetate, CO2, formate, and H2 that are utilized by methanogens and other hydrogen-consuming microbes. S. aciditrophicus benzoate degradation proceeds by a multistep pathway with many intermediate reactive acyl-coenzyme A species (RACS) that can potentially Nε-acylate lysine residues. Herein, we describe the identification and characterization of acyl-lysine modifications that correspond to RACS in the benzoate degradation pathway. The amounts of modified peptides are sufficient to analyze the post-translational modifications without antibody enrichment, enabling a range of acylations located, presumably, on the most extensively acylated proteins throughout the proteome to be studied. Seven types of acyl modifications were identified, six of which correspond directly to RACS that are intermediates in the benzoate degradation pathway including 3-hydroxypimeloylation, a modification first identified in this system. Indeed, benzoate-degrading enzymes are heavily represented among the acylated proteins. A total of 125 sites were identified in 60 proteins. Functional deacylase enzymes are present in the proteome, indicating a potential regulatory system/mechanism by which S. aciditrophicus modulates acylation. Uniquely, Nε-acyl-lysine RACS are highly abundant in these syntrophic bacteria, raising the compelling possibility that post-translational modifications modulate benzoate degradation in this and potentially other, syntrophic bacteria. Our results outline candidates for further study of how acylations impact syntrophic consortia.
). They degrade a broad range of aliphatic and aromatic acids to hydrogen, formate, CO2, and acetate by cooperation with hydrogen- and/or formate-consuming microbes such as archaea which generate methane (
). Here, the presence of methanogens and/or other hydrogen- and/or formate-consuming partners is required to maintain low hydrogen and formate concentrations, such that syntrophic substrate degradation can occur spontaneously (
). This obligate partnership between members of the microbial community exists because anaerobic degradation is thermodynamically unfavorable when hydrogen or formate levels are high (
). Importantly, degrading syntrophic substrates requires multiple enzymatically catalyzed reactions to be performed on many acyl-coenzyme A (CoA) intermediates. Even when hydrogen and formate levels are low, the essential acyl-CoA oxidations of syntrophic metabolism are weakly exergonic with free-energy changes close to thermodynamic equilibrium (
). The energetic challenges of living at the edge of thermodynamic feasibility make syntrophic microbes interesting models for exploring energy conservation and metabolic regulation (
Syntrophus aciditrophicus is an anaerobic, gram-negative bacterium that degrades fatty, aromatic, and alicyclic (cyclohexane-1-carboxylate) acids to acetate, CO2, formate, and hydrogen when grown in coculture with hydrogen- and/or formate-consuming partner microorganisms. S. aciditrophicus can grow in pure culture with crotonate or benzoate (
Syntrophus aciditrophicus sp. nov., a new anaerobic bacterium that degrades fatty acids and benzoate in syntrophic association with hydrogen-using microorganisms.
Cyclohexa-1,5-diene-1-carbonyl-coenzyme A (CoA) hydratases of Geobacter metallireducens and Syntrophus aciditrophicus: Evidence for a common benzoyl-CoA degradation pathway in facultative and strict anaerobes.
6-Oxocyclohex-1-ene-1-carbonyl-coenzyme A hydrolases from obligately anaerobic bacteria: Characterization and identification of its gene as a functional marker for aromatic compounds degrading anaerobes.
Due to challenges in isolating, manipulating, and cultivating syntrophic microorganisms, questions still remain regarding the regulatory mechanisms involved in the important anaerobic degradations. Oxidizing benzoate and other substrates generates reduced cofactors (nicotinamide adenine dinucleotide reduced and flavin adenine dinucleotide reduced) that, if unregulated, could cause stress and cellular dysfunction (
). Intermediates of substrate metabolism can also stress cells. Various acid substrates in S. aciditrophicus are activated through their respective CoA derivative, which serves as a scaffold through subsequent conversions, until acetate, CO2, hydrogen, and formate are released (
Metabolism of benzoate, cyclohex-1-ene carboxylate, and cyclohexane carboxylate by Syntrophus aciditrophicus strain SB in syntrophic association with H2-using microorganisms.
). The intermediates of these pathways are a series of reactive acyl-CoA species (RACS), with potential to modify the nucleophilic side chain of lysine (
). Given the multiple stressors that S. aciditrophicus cells experience in carbon metabolism and that its enzymes catalyze in either direction, depending on the substrate and environmental conditions (
), it is suspected that previously undescribed pathway regulation is involved. The reversibility of syntrophic metabolism and its correspondingly low energy yields makes us consider that post-translational modifications (PTMs), specifically acylations associated with RACS intermediates, may play a role modulating degradations and syntheses in the cell.
Identifying and characterizing PTMs is most often done using mass spectrometry. At the systems level, proteomics can capture PTMs across all proteins, but a system presenting a wide range of acylations brings analytical complications, including potential sequence misidentifications due to isomeric and/or isobaric combinations (
). Ambiguities must be addressed to confidently ascribe sequences, modifications, and modified residues in data. Increases in theoretical search space and scoring thresholds needed to satisfy false-discovery rate cutoffs for database searching are also considerations for proteomes possessing a broad variety of acylations (
Here, we take a systems-level approach to identify post-translational acylations associated with syntrophic benzoate metabolism by S. aciditrophicus. Without using PTM-specific enrichment procedures, we identify a wide range of acylations by targeting acyl-CoAs that are known to be involved in fatty and aromatic acid degradation in this and other microbes (
Metabolism of benzoate, cyclohex-1-ene carboxylate, and cyclohexane carboxylate by Syntrophus aciditrophicus strain SB in syntrophic association with H2-using microorganisms.
). Acylations previously undocumented across the bacterial domain are identified. An understanding of how the modifications are regulated is currently lacking, although bacterial sirtuins are known to deacylate promiscuously (
) with 20 mM crotonate or 10 mM crotonate plus 2 mM benzoate. S. aciditrophicus was grown in coculture with Methanospirillum hungatei JF1 (ATCC 27890) in 2-L Schott bottles with 1 L of Tanner’s mineral medium with 14 mM benzoate. Tanner’s mineral medium was used for large culture volumes to avoid chemical precipitants that form in large volumes of the basal medium. Tanner’s mineral medium (
in: Hurst C. Crawford R. Garland J. Lipson D. Mills A. Stetzenbach L. Manual of Environmental Microbiology. 3rd Ed. American Society of Microbiology,
Washington, DC2007: 69-78
) contained (per liter) 10 ml of Tanner’s minerals, 5 ml of Tanner’s metals, 10 ml of Tanner’s vitamins, 1 ml of 0.1% resazurin, 3.75 g of NaHCO3, and 20 ml of cysteine sulfide (2.5%) solution (
New approach to the cultivation of methanogenic bacteria: 2 Mercaptoethanesulfonic acid (HS CoM) dependent growth of methanobacterium ruminantium in a pressurized atmosphere.
) were used to prepare media and solutions and to inoculate and sample cultures. The headspace was pressurized with N2/CO2 (80%:20% v/v) to 70 kPa, and cultures were incubated at 37 °C without shaking.
Cultures were grown to mid-log phase and harvested by centrifugation (8000g for 20 min at 4 °C) under strict anaerobic conditions. The cell pellet was washed twice by resuspending in 50 mM anoxic potassium phosphate buffer (pH 7.5) and centrifuging as described previously. The final cell pellet was resuspended in anoxic potassium phosphate buffer, transferred to cryovials, and stored in −80 °C (
New approach to the cultivation of methanogenic bacteria: 2 Mercaptoethanesulfonic acid (HS CoM) dependent growth of methanobacterium ruminantium in a pressurized atmosphere.
Two-Dimensional Polyacrylamide Gel Electrophoresis of S. aciditrophicus
Two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) was performed, and protein-containing spots identified as described in the Supplemental Materials section of the study by James et al. (
). Frozen cell pellets were resuspended in lysis buffer containing 4.0% (v/v) ammonium lauryl sulfate, 0.1% sodium deoxycholate (w/v), 5 mM tris(2-carboxyethyl)phosphine, and 100 mM ammonium bicarbonate (ABC). Proteins were digested via the enhanced filter-aided sample preparation method, similar to that described by Erde et al. (
). Briefly, lysate was buffer-exchanged into 8 M urea, 0.1% (w/v) sodium deoxycholate, and 0.1% (w/v) n-octyl glucoside using a 10-kDa Microcon ultrafiltration unit (Millipore). Proteins were alkylated with 17 mM iodoacetamide for an hour at room temperature, after which they were exchanged into 0.1% (w/v) sodium deoxycholate/0.1% (w/v) n-octyl glucoside in 100 mM ABC. Trypsin (1:100 w/w) was added, and the solution was incubated overnight at 37 °C. Ethyl acetate extraction was used to remove detergents, as previously described (
The resulting peptides were separated off-line via hydrophilic interaction liquid chromatography. Fifty micrograms of peptides (quantified with Pierce Quantitative Fluorometric Peptide Assay) were deposited onto a BioPureSPN MACRO PolyHYDROXYETHYL A column (The Nest Group) in 90% acetonitrile and 150 mM ammonium formate, pH 3. Peptides were eluted in six fractions using ammonium formate buffers sequentially decreasing in organic content: 80% acetonitrile, 78% acetonitrile, 74% acetonitrile, 71% acetonitrile, 40% acetonitrile, and 35.5% acetonitrile. Fractions one and six were combined. No antibody affinity fractionations were performed to enrich acylated peptides.
Liquid Chromatography and Tandem Mass Spectrometry
Peptides were mass-measured and fragmented using high-performance liquid chromatography tandem mass spectrometry. Using an EASY nLC1000 liquid chromatography system (Thermo Scientific), peptides (200 ng) were loaded onto an Acclaim PepMap100 C18 trap column (Thermo Scientific, Product #16-494-6, 75 μm × 2 cm, 100 Å) and separated on an Acclaim PepMap RSLC C18 analytical column (Thermo Scientific, Product #03-251-873, 75 μm × 25 cm, 100 Å). Buffer A (0.1% formic acid) and buffer B (0.1% formic acid in acetonitrile) were mixed and delivered at 300 nl min−1 in a gradient of 3 to 20% B for 62 min, 20 to 30% B for 31 min, 30 to 50% B for 5 min, and 50 to 80% B for 2 min.
A Q Exactive (Thermo Fisher Scientific) quadrupole-orbitrap mass spectrometer was operated using data-dependent acquisition mode. MS scans (m/z 300–1800) were acquired at 70,000 resolution, with an automatic gain control target set to 1E6 and a maximum fill time of 100 ms. The ten most abundant precursor ions were dissociated sequentially using higher energy collisional dissociation at a normalized collisional energy of 27 (unless otherwise indicated), and MS/MS spectra were acquired at 17,500 resolution with an automatic gain control target of 1E5 at a maximum fill time of 80 ms. The mass spectrometry data from 2D gels were deposited to the ProteomeXchange Consortium via the PRIDE partner repository with dataset identifier PXD025631. Shotgun datasets submitted through MassIVE are identified as PXD025603.
Database Searching
Q Exactive ∗.RAW files were analyzed using ProteomeDiscoverer (version 1.4) employing Matrix Science’s Mascot search algorithm (
). A database concatenating UniProt S. aciditrophicus and M. hungatei protein sequences (as of July 8, 2019) to the sequences of common contaminants was used for searches. The total number of entries in this database was 6479. Parameters used for the search were enzyme name, trypsin; maximum missed cleavage sites, two; precursor mass tolerance, 10 ppm; fragment mass tolerance, 0.02 Da; and variable methionine oxidation and cysteine carbamidomethylation. Searches also considered variable acyl-lysine modifications, all of which are shown in Table 1. These modifications were selected from RACS identified in the aromatic acid–degrading pathways of a variety of bacteria. Spectra matched to acylated peptides with Mascot ion scores ≥20 were subjected to manual examination, the criteria of which prioritized the presence of immonium ions and product ions spanning the modification as well as from regions that do not contain the modification. A score of 20 is associated with an expectation value of ≥0.1, chosen to broadly collect spectra of potentially modified peptides.
The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate relevant pathways associated with the acylated proteins identified. UniProt accessions were converted to KEGG identifiers and mapped to the KEGG pathway database using KEGG Mapper (
). To identify Gene Ontology (GO)/KEGG pathways frequently represented by acylated proteins, functional enrichment was performed by the StringApp (version 1.6.0) in Cytoscape (version 3.8.2), and p-values were corrected for multiple testing within each category using the Benjamini–Hochberg procedure (
) from the Uniprot S. aciditrophicus database. The other parameters were set to default values: p-value threshold: 0.000001, the minimum number of occurrences for residue/position pair: 10.
All acylated proteins identified were searched against the STRING database (version 11.0) to reveal potential protein–protein interactions (PPIs) (
). Only protein interactions found among the acylated proteins were selected, removing any external candidate interactions. Only interactions with at least medium confidence (>0.4) were included. The interaction network was visualized in Cytoscape (version 3.8.2).
Cloning, Expression, and Purification of Syn_00042 and Syn_01020
Two S. aciditrophicus sirtuin candidate genes annotated as sir2 family proteins were identified in the genome by homology search. The genes were obtained via polymerase chain reaction (PCR) of genomic DNA, cloned into plasmid pMAPLE21, and expressed in Escherichia coli strain BL21 as previously described by Arbing et al. (
Heterologous expression of mycobacterial Esx complexes in Escherichia coli for structural studies is facilitated by the use of maltose binding protein fusions.
). The recombinant proteins were purified by Ni-affinity chromatography followed by size exclusion (Superdex 75) and Q-Sepharose anion exchange. Purity was greater than 95% as visualized by SDS-PAGE. The apparent protein sizes were consistent with the predicted gene/protein sizes.
Sirtuin Assay (anhydride)
Acyl-insulins were prepared from the corresponding acyl-anhydrides; for example, glutaric anhydride was used for glutaryl-lysine, using methods adapted from Baeza et al. (
). About 25 μmol of anhydride was added to 100 μl of a 1-mg/ml solution of human insulin (Alfa Aesar, J67626) in 100 mM ABC. After incubating at 4 °C for 20 min, the solution was readjusted to pH ∼8 using an ammonium hydroxide solution. Anhydride addition, incubation, and pH adjustment were repeated twice more. Ostensible O-acylation was reversed by adding 50% w/v of hydroxylamine hydrochloride in H2O (adjusted with ammonium hydroxide to pH 7–8). Following overnight, room-temperature incubation, the modified products were buffer-exchanged into 100 mM ABC using 3-kD molecular weight cut-off Amicon spin filters (Millipore).
For the matrix-assisted laser desorption/ionization (MALDI)–based activity assay, an approximately 27 μM solution of acyl-modified insulin was mixed with 0.24 μM recombinantly expressed SYN_00042 and excess oxidized nicotinamide adenine dinucleotide (NAD+) in 100 mM ABC. The solution was incubated at 37 °C for 2 h. Sample aliquots were mixed 1:1 with saturated sinapinic acid in 50% acetonitrile/0.1% trifluoroacetic acid and spotted onto a MALDI sample stage. MALDI mass spectra were obtained with an Applied BioSystems Voyager-DE STR time-of-flight mass spectrometer.
Time-series measurements of glutaryl-insulin masses used a Synapt G2-Si quadrupole time-of-flight mass spectrometer (Waters) with electrospray ionization. The assay was performed at 37 °C for the period of time indicated and quenched with 50% acetonitrile. The glutaryl-insulin was diluted to a final concentration of 10 μM for mass analysis with 20 μM ubiquitin added as internal standard. The solution was delivered by direct infusion.
Experimental Design and Statistical Rationale
S. aciditrophicus cells were grown in pure culture on crotonate as a sole carbon source or with crotonate/benzoate to analyze via 2D-PAGE. Cells were cocultured with M. hungatei and grown on benzoate for mass spectrometry proteomic analysis. Three samples from each culture condition were grown, and for each biological sample, three technical replicates were analyzed. As the priority was to maximize the identification of potentially low stoichiometry PTMs, offline hydrophilic interaction liquid chromatography was used to increase the depth of the identifiable proteome. Resulting data were searched for the acyl modifications shown in Table 1. To identify PTMs, mass spectra matched to acylated peptides with Mascot ion scores ≥20 were subjected to manual examination. Modifications were considered present if the tandem mass spectra included acyl-lysine–associated immonium ions (
). Frequently, spots were composed of multiple proteins, limiting the value of densitometry for evaluating abundance. In some spots, MS identified seven or more proteins with two or more tryptic peptides (or a single tryptic peptide identified from an MS/MS spectrum matched to that in a different spot, known to contain that protein, displaying 2+ tryptic peptides) (supplemental Table S1). Abundant proteins are listed in Table 2. Within the darkest spots, some MS heterogeneity suggestive of small variations in sequence was also observed, potentially reflecting high error rates in transcription and/or translation. Little is known about this latter effect.
Fig. 1Proteomic patterns of S. aciditrophicus as detected by 2D-PAGE. Two-dimensional gel analysis of the S. aciditrophicus proteome for cells grown on benzoate supplemented with crotonate. The x-axis is separated by isoelectric focusing (IEF), whereas the y-axis is separated by SDS-PAGE. Selected proteins present in the spots are listed in Table 2. 2D-PAGE, two-dimensional polyacrylamide gel electrophoresis.
Polyribonucleotide nucleotidyltransferase, electron transfer flavoprotein beta-subunit
Gel spots numbered in Figure 1 were excised, and their identities determined by mass spectrometry on a QStar-XL quadrupole time-of-flight mass spectrometer.
a Denotes proteins identified that were acetylated.
These analyses revealed many instances in which the same protein was identified from multiple spots at the same molecular size but different isoelectric points. This pattern, often seen with glycoproteins (
), suggests heterogeneous PTMs, but we found no evidence of S. aciditrophicus glycoproteins. Spot trains might strike some as reminiscent of carbamylation, but the artifacts were easily ruled out. Trains do not accompany every spot, and MS/MS analyses found protein carbamylation to be insignificant. Some of these proteins were identified as housekeeping proteins (Table 2), while others have roles in the catabolism of benzoate and other aromatic and fatty acids. Given the high abundance of RACS in these pathways, acyl-lysine modifications were considered possible. Acylation will neutralize the positive charge on lysine, or, in the cases of acidic RACS, such as glutaryl-CoA, will switch lysine to a negatively charged site. Both types of acylation reduce protein isoelectric points, shifting migration toward the anode (the left in standard 2D gel images); hence, a distribution of many and differently charged acylations would give a pattern similar to that seen in Figure 1.
Digests for spot 20 (Fig. 1) established that the acetylated and nonacetylated peptides YGTKacPEDLALIR and YGTKPEDLALIR of the acetyl-CoA acetyltransferase SYN01681 were both present. That both acetylated and nonacetylated lysines were present could suggest that migration-compensating acylations or modifications reside at other sites on some of the SYN01681 protein in this spot. From spot 21, peptides TAVGAFGGSLKacGVR and MoxAKacLAPVFK of the acetyl-CoA acetyltransferase SYN02642 were observed bearing acetyl lysines. Nonacylated TAVGAFGGSLK and LAPVFK were also observed in the spot, but at lower signal levels. Two acetylated sites were also observed from SYN01310 in spot 22, from peptides ATEEFSKacQLGK and LMoxAGSIKacK of 3-hydroxyacyl-CoA dehydrogenase. A corresponding unmodified peptide, GYYTSDETFKATEEFSK, was also observed in this spot, although neither unacetylated LMAGSIK nor LMoxAGSIK was found. From the spot 23 proteins SYN_01788 and SYN_02636, acetylated peptides TAAALEALEKacR and YASLPGIMKacAK were detected, respectively. Unacetylated TAAALEALEK was also detected from the digested spot, but not the unacetylated SYN_02636 peptide. Most interestingly, peptide TATGKacIQR, found in spot 19, could be attributed to SYN02896 and/or SYN02898. These benzoate-CoA ligases are 77% identical and are both present in this spot, as ascertained from their unique peptides. Despite strong signals from the acylated peptide, no unacetylated peptides spanning the TATGK segment were recovered. Supplemental Spectra S1 present annotated tandem mass spectra from spots 20 to 24.
Shotgun Proteomics Identifies Many Acyl Modifications in the S. aciditrophicus Proteome
Three biological replicates of protein lysates from S. aciditrophicus/M. hungatei cocultures utilizing benzoate as the carbon source were analyzed in three technical replicates each by liquid chromatography tandem mass spectrometry data-dependent acquisition (supplemental Tables S2 and S3). Analyses for one biological replicate incorporated stepped collision energies to increase the diagnostic immonium ion abundances from acyl-lysines (
). Speculating that peptides acylated by reactive intermediates from benzoate degradation might be observable, we sought evidence by including acylation mass shifts in database searches. Table 1 illustrates the acylation shifts that were considered known intermediates for S. aciditrophicus, for related organisms, and common acylations. A total of 125 sites were identified in 60 different proteins, and each site was modified by at least one of the acyl groups (supplemental Tables S4 and S5; supplemental Spectra S2). Acylated peptides comprise 3.2% of the unique peptides identified (supplemental Table S3). The acetylated and 3-hydroxypimelylated species predominated in both numbers of proteins and numbers of modified sites. All of the PTMs found in each biological replicate are listed in Table 3, along with information about the supporting immonium ions. Acylations corresponding to seven of the 20 searched intermediates were found (Table 3). Modifications were detected in biological replicates reproducibly (supplemental Fig. S1 and supplemental Table S4). Peptides with missed cleavages enabled elution time comparisons between modified peptides and their unmodified counterparts. Twenty-eight of these peptide pairs were identified; in 26 pairs, the acylated peptide eluted later by an average of 9.37 min (supplemental Table S6).
Table 3Summary of acyl-lysine modifications identified in S. aciditrophicus
Modification
Proteins identified
Sites identified
Cyclized immonium ion m/z
Peptides with cyclized immonium ion (%)
Benzoylation
1
2
188.107
100%
3-Hydroxypimelylation
16
20
224.129, 242.139
75%, 55%
Glutarylation
4
6
198.112
100%
Crotonylation
2
2
152.107
100%
3-Hydroxybutyrylation
4
6
170.118
83%
Acetylation
48
104
126.091
94%
Butyrylation
3
3
154.123
33%
For each type of modification, the number of unique proteins bearing the modification and the total number of sites observed to be modified are presented. For the modified sequences, the percentages of MS/MS spectra that contained the diagnostic ions are stated.
PTMs can be misidentified in complex proteomic datasets due to chimeric precursors, isobaric peptides, spectral complexity, and insufficient sequence-related ions (
). The large search space traversed in seeking multiple, variable modifications from shotgun datasets challenges confident assignments. To meet these challenges, we relied upon the presence of immonium ions specific for each putative acyl modification. Immonium and immonium-like ions are strong indicators of acyl-lysine (
). Proportions of putative acyl-modified spectra displaying diagnostic immonium ions are presented in Table 3. We also examined sequences upstream and downstream of each acylated residue to determine if there were sequence preferences for these modifications, and indeed, we found that nearby glycines favored lysine acetylation (supplemental Fig. S3).
Pathway Analysis of Modified Proteins
A functional annotation analysis classified proteins by KEGG pathways and GO. The GO functional classification grouped acylated proteins into three categories: biological process, molecular function, and cellular component. Most acylated proteins mapped onto KEGG pathways had roles in metabolic processes, as expected. Benzoate degradation, oxidative phosphorylation, and carbon metabolism were highly enriched (Fig. 2A). GO enrichment analyses of all acylated proteins revealed biological processes primarily related to the synthesis of metabolites, as well as energy production (Fig. 2B). Processes related to nucleotide metabolism were also highly enriched, in line with recent reports that nucleotide-binding regions increase acylation of proximate lysines (
). Altogether, the findings suggest that these acylated proteins are generally involved in carbon metabolism, degrading aromatic and fatty acids.
Fig. 2Functional enrichment of acylated proteins.A, Gene Ontology biological processes that are enriched. B, enriched KEGG pathway processes. C, a functional protein association network generated using the STRING database. Proteins involved in oxidative phosphorylation are colored in blue. Proteins involved in degrading benzoate are red. KEGG, Kyoto Encyclopedia of Genes and Genomes.
Acyl-lysine modifications have been shown to impact PPIs and affect the formation of enzyme complexes, altering their role in cellular physiology. To determine connections between the acylated proteins and to elucidate their functional PPI networks, we generated an interaction network based on the STRING database, visualized through Cytoscape (Fig. 2C). Within the network, there are two highly interconnected clusters, benzoate degradation and oxidative phosphorylation (Fig. 2C). These strong physiological interactions among all modified proteins led us to further investigate the role of protein acylation in the degradation of aromatic compounds.
Benzoate Degradation Enzymes Are Heavily Acylated
Proteins across many pathways are modified by different acyl groups at varying levels, summarized by the heat map in Figure 3A. Acetylation is found in the widest range of functional pathways, with the highest number of acetyl-lysines found in proteins related to the biosynthesis of secondary metabolites. Other acyl modifications are enriched in benzoate degradation pathway proteins, corresponding to the reactive acyl-CoA intermediates from the pathway (
Metabolism of benzoate, cyclohex-1-ene carboxylate, and cyclohexane carboxylate by Syntrophus aciditrophicus strain SB in syntrophic association with H2-using microorganisms.
). These findings imply that high local concentrations of intermediates may modify proteins spontaneously.
Fig. 3Mapping sites of protein lysine acylation based on COG functional pathway analysis.A, a heatmap indicating the pathways (KEGG Ontology) involved for the acyl-modified protein identified. The color indicates the number of proteins in each category that have the specified acylation. B, modified sites identified in the “upper” and “lower” benzoate pathway. Bar plots illustrate the range of modifications found on associated enzymes. One square indicates one site of modification, and each color represents a different acylation. RACS modification types are indicated in the upper right corner boxed color code. COG, clusters of orthologous groups;KEGG, Kyoto Encyclopedia of Genes and Genomes; RACS, reactive acyl-coenzyme A species.
The largest number and widest variety of acyl modifications were found in benzoate degradation (Fig. 3B). Ten different proteins in the pathway were observed bearing several different acyl modifications. Acylated proteins included benzoate-CoA ligase, the initial step in the pathway. Several of the other acylated proteins are members of gene clusters involved in benzoic acid metabolism (bam genes) that are found in numerous anaerobic delta-proteobacteria (
). The gene products of bamR (SYN_01653), bamQ (SYN_01655), and bamA (SYN_01654), which successively reduce cyclohexa-1,5-diene-1-carboxyl-CoA to 3-hydroxypimelyl-CoA, display modifications. All three gene products are acetylated, those of bamQ and bamA are glutarylated and 3-hydroxypimelylated, and that from bamA is also 3-hydroxybutyrylated and benzoylated. Interestingly, all five modifications for the bamA product occur at the same site, lysine 261 (supplemental Fig. S3 and supplemental Spectra S2). Peptides containing these lysine residues were also observed unmodified.
S. aciditrophicus uses acetyl-CoA synthetase (acs1, SYN_02635) to make adenosine triphosphate (ATP) from acetyl-CoA. Other bacteria, in contrast, use phosphate acetyltransferase and acetate kinase, enzymes that are absent in the S. aciditrophicus genome. Using acs1 to produce ATP is also remarkable because acetyl-CoA synthetases were previously thought to function only in activating acetate to acetyl-CoA, not in the reverse direction (ATP-forming direction) (
). The rates at which ATP is produced when acetyl-CoA is the limiting substrate differ between purified and recombinantly expressed Asc1 (Vmax of 7.5 and 1.2, respectively and Km of 0.41 and 1.34 mM, respectively) (
), suggesting that there are factors associated with the in vivo protein that are not recapitulated when the enzyme is recombinantly expressed. PTMs may be one such factor. Indeed, six different sites on this protein were seen in acylated forms, all of which included acetylation, except for one site that was only butyrylated.
Proteins involved in energy conservation were also heavily acetylated. ATP synthase, Atp1, is a large complex that contains many components. Acetylation sites were identified from several subunits of Atp1, namely on the alpha (SYN_00546), beta (SYN_00544), and gamma (SYN_00545) chains (supplemental Table S5). Two homologs of ATP synthase B (the beta subunit) were also acetylated (SYN_00548 and SYN_00549).
Benzoate-CoA Ligase Reveals a Highly Modified Residue in a Conserved Region
Five of the seven acyl modifications identified in S. aciditrophicus proteins appeared in the sequence TATGKIQR, present in paralogs Bcl1 (SYN_02898) and Bcl2 (SYN_02896), which differ in their activity toward different fatty and aromatic acids. The modified lysine is critical to benzoate-CoA ligase (BCL) function (
Biochemical and structural characterization of the paralogous benzoate CoA ligases from Burkholderia xenovorans LB400: Defining the entry point into the novel benzoate oxidation (box) pathway.
Reversible Nε-lysine acetylation regulates the activity of acyl-CoA synthetases involved in anaerobic benzoate catabolism in Rhodopseudomonas palustris.
), catalyzing the first step of benzoate degradation. Immonium ions present in the spectra shown in Figure 4A confirm that the unique mass shifts do, indeed, result from acylation and are not induced by misidentifications that may result from additive side reactions that occur during sample processing; for example, formylation from exposure to high concentrations of formic acid or 12-Da mass shifts from exposure to formaldehyde (
). This modification site was further investigated to determine what functional role, if any, the PTMs play. Basic Local Alignment Search Tool (BLAST) sequence alignments against closely related BCLs (
Biochemical and structural characterization of the paralogous benzoate CoA ligases from Burkholderia xenovorans LB400: Defining the entry point into the novel benzoate oxidation (box) pathway.
) indicated that the lysine falls in a highly conserved region (Fig. 4B). Previous Burkholderia xenovorans structural studies of benzoate-bound BCL (Bxe_A1419), a homolog with 46% sequence identity to Bcl1 and 44% to Bcl2, indicated that this conserved lysine is located near the benzoate carboxyl (Fig. 4C) and is hypothesized to coordinate to benzoate during enzymatic activity (
Biochemical and structural characterization of the paralogous benzoate CoA ligases from Burkholderia xenovorans LB400: Defining the entry point into the novel benzoate oxidation (box) pathway.
). Acylation both neutralizes the ε-amino side chain and adds a bulky group disrupting coordination to the carboxylate. The Rhodopseudomonas palustris BCL (BadA), with 43% sequence identity to Bcl1 and 42% to Bcl2, has also been found acetylated at that residue (K512). Studies in vitro established that acetylation abolishes enzymatic activity (
Reversible Nε-lysine acetylation regulates the activity of acyl-CoA synthetases involved in anaerobic benzoate catabolism in Rhodopseudomonas palustris.
), suggesting possible feedback inhibition of the enzyme and, in turn, the pathway.
Fig. 4Modified benzoate-CoA ligase (BCL).A, HCD spectra of acylated TATGKIQR from SYN_02896 and SYN_02898. The benzoylated spectrum (top) was collected with a stepped NCE method of 27 V/40 V. Other spectra were collected at 27 NCE. Immonium ions of modifications are denoted in red. B, sequence alignment of known BCLs across microbial systems. The line above the sequence (TATGKIQR) denotes the peptide identified in (A). The asterisk (∗) denotes the corresponding modified lysine residue identified in S. aciditrophicus. C, crystal structure of a closely related BCL from Burkholderia xenovorans, determined by Bains and Boulanger (
Biochemical and structural characterization of the paralogous benzoate CoA ligases from Burkholderia xenovorans LB400: Defining the entry point into the novel benzoate oxidation (box) pathway.
), indicates that the acyl-lysine is proximal to a benzoate bound to the active site. BCL, benzoate-CoA ligase; HCD, higher energy collisional dissociation; NCE, normalized collision energy.
S. aciditrophicus Sirtuins Deacylate Promiscuously
Given the ubiquity of acyl modifications in S. aciditrophicus, an ability to control and/or regulate acyl modifications would seem to be advantageous and critical. Sirtuin proteins have been shown to possess generalized deacylase activity, often promiscuous with respect to both sequence context and acyl modification (
). To determine if a deacylase capability is present, the S. aciditrophicus genome was searched with BLAST to determine if E. coli cobB homologs were present and SYN_00042 and SYN_01020 were found. A sequence alignment between these proteins displays their similarities (supplemental Fig. S4). SYN_00042 and SYN_01020 were recombinantly expressed to assay activity by mass spectrometry (Experimental Procedures).
Insulin, a dual-chain protein with one lysine and two free NH2-termini, was synthetically acylated using acyl anhydrides. Insulin was selected as it only has a single lysine residue and has been well characterized as a standard protein used in intact mass spectrometry work. Acylated insulin was incubated with recombinantly produced SYN_00042 and its cofactor, NAD+. MALDI mass spectra of glutarylated insulin showed a 114-Da decrease in mass after incubation with SYN_00042 (Fig. 5A and supplemental Table S7). The recombinant protein displayed deacylase activity against butyryl and succinyl acylations, as well (supplemental Fig. S5 and supplemental Table S7, A–G), but did not appear to act on shorter acetyl or propionyl chains. A quantitative assay of SYN_00042 gene product deacylase activity on glutaryl-insulin, performed by electrospray ionization MS, clearly shows time-dependent deglutarylase activity (Fig. 5B and supplemental Table S7, H–M). SYN_01020 was expressed in an attempt to assess its function, but did not display reproducible activity.
Fig. 5S. aciditrophicus sirtuin SYN_00042 shows in vitro deacylase activity.A, MALDI-MS of insulin and modified insulin deacylated by sirtuin (Syn_00042). The peak corresponding to “Glutarylated Nε - Lysine” reflects the mass shift from three glutaryl modifications (+342 Da), two at the N-termini and one at the single lysine. B, ESI-MS of the 4+ charge states 0, 2, 15 and 60 min after sirtuin addition. Purple peak clusters (left) correspond to insulin modified at the two N-termini; red clusters identify insulin modified at all three sites. Ammonium adducts (+17 Da) are also present. ESI, electrospray ionization; MALDI, matrix-assisted laser desorption/ionization; MS, mass spectrometry.
PTMs can modulate protein function in response to cellular or environmental changes and may occur spontaneously due to cellular conditions or simply over time (
). Acylation is one class of modification that can be spontaneously induced in the presence of acyl phosphate or RACS or can be enzyme-mediated by a lysine acyltransferase (KAT) (
). Regardless of its origin, acylation can affect enzymatic activity profoundly. In vitro work has shown that modifying lysine side chains near catalytic regions of bacterial enzymes can directly alter function (
). Work in other systems has shown more subtle but equally significant acylation effects: acetylation can disrupt enzyme complexes, thereby altering activity (
). Most acetylated proteins in our system were involved in the synthesis of secondary metabolites (Fig. 3A), which is sensible given that acetyl-CoA is an intermediate metabolite in many of those pathways. Benzoyl- (
) have shown function in some systems under certain conditions, although the studies have been less extensive than those for acetyl-lysine, and few, if any, bacterial studies have reported their occurrence. A 3-hydroxypimelylation has previously been reported in S. aciditrophicus by our laboratory (
). Detecting these modifications across biological replicates without pre-enrichment suggests that they are uniquely prevalent and abundant in this system. While butyryl-CoA is not expected to be an intermediate in the degradation pathway, S. aciditrophicus has genes for butyrate dehydrogenase activity (
Syntrophus aciditrophicus sp. nov., a new anaerobic bacterium that degrades fatty acids and benzoate in syntrophic association with hydrogen-using microorganisms.
Syntrophus aciditrophicus sp. nov., a new anaerobic bacterium that degrades fatty acids and benzoate in syntrophic association with hydrogen-using microorganisms.
To begin to understand how these modifications affect global cellular function requires a comprehensive and unbiased catalog of acylations. Recording a biological system’s acyl modifications comprehensively presents many challenges. Low stoichiometries of acyl modifications in many model systems have required enrichment to identify and characterize acylated peptides (
). Pan-specific antibodies are valuable for enriching specifically modified peptides but have some limitations. Antibody cross-reactivity may lessen the modification specificity of enriched peptides, a problem for Western blotting, but not mass spectrometry. Sequence-dependent binding efficiencies can bias results and challenge quantification strategies. Our syntrophic system provides a unique opportunity to investigate the array of acylations without enrichment, bypassing technical challenges in other systems. The experimental approach described herein using RACS as a predictor for lysine acylations complements antibody enrichment and utilizes diagnostic marker ions to validate putative modifications proposed from an enlarged search space.
Many of the S. aciditrophicus proteins decorated with acyl modifications function to degrade benzoate, a process that could result in carbon or reductive stress, if unregulated. Modulating enzymatic activity through protein acylation would be a fitting and elegant mechanism for global metabolic regulation, given the limited energy available to these cells and the high concentrations of acyl-CoA intermediates, a feature also evidenced by the unusual ability of this species to use acetyl-CoA to form ATP (
). The benzoate degradation pathway is a likely candidate for this level of regulation as its constituents frequently interact with RACS. While we have identified a wide range of acyl modifications in this system, the biological role that these modifications play has yet to be explored. We can, however, provide some insight into the functional effect of certain acylation sites by drawing on similar modifications in other systems. BCL acylation, specifically acetylation, has been identified as a means of negative feedback inhibition, stopping benzoate-CoA ligase activity in R. palustris by directly inhibiting enzymatic catalysis (
Reversible Nε-lysine acetylation regulates the activity of acyl-CoA synthetases involved in anaerobic benzoate catabolism in Rhodopseudomonas palustris.
). The large number of acyl modifications we report at this site may indicate that when its respective acyl-intermediate builds, the PTM can act as a brake to slow or stop degradation, thereby acting as a negative feedback inhibitor to mitigate carbon and reductive stress. BCL also consumes ATP, creating a direct link between the energetic state of the cell and aromatic degradation. Additionally, the multiple sites of acetylation found on ATP synthase subunits further hint at a link between bioenergetics and acylation as it, too, demonstrates a relationship between the reactive metabolites and oxidative phosphorylation.
Beyond energetics, the cellular reductive state may also be regulated by the acylation of enzymes. Benzoate degradation generates nicotinamide adenine dinucleotide reduced and flavin adenine dinucleotide reduced, but excess reducing capacity stresses cells. A buildup of intermediates that subsequently modifies and slows the rates of enzymatic catalysis would safeguard against damaging reductive stress (
). This link is further enforced by the presence of sirtuins, whose activity is non–energy-conserving, and relies upon the oxidized cofactor NAD+ for activity (
), innately linking the extent of protein acylation to cellular redox state. Interestingly, the sirtuin assayed in this system appears to have a bias for longer chain acylations, particularly glutarylation, a trait that is shared with Sirt4 and Sirt5 in mammalian systems (
). Although S. aciditrophicus sirtuins were identified by homology to CobB, their activities deviate notably from those of the E. coli enzyme, as illustrated by their action on only longer acyl groups, not acetyl. Previous studies indicated that sirtuin specificity across species including E. coli CobB, mammalian Sirt5, and Archaeoglobus fulgidus Sir2Af2 demonstrates that while also having deacetylase activity, their selectivity toward long-chain deacylation is dependent on a YxxR motif present in the binding pocket (
Metabolism of benzoate, cyclohex-1-ene carboxylate, and cyclohexane carboxylate by Syntrophus aciditrophicus strain SB in syntrophic association with H2-using microorganisms.
). It is therefore expected that some metabolites would accumulate under high hydrogen concentrations, providing another link between the metabolic conditions of the cell and the observed acyl modifications. Interestingly, the lower benzoate pathway degrading 3-hydroxypimelyl-CoA (Fig. 2C) is conserved across many anaerobes that degrade benzoate and other aromatic compounds (
). The acyl modifications we identified on these proteins suggest that the modifications are likely present in other bacteria sharing the same or similar metabolic pathways, especially under cellular conditions that induce an accumulation of intermediates. S. aciditrophicus has demonstrated a wide array of acyl modifications within its proteome. Further investigation into the function of these acylations presents an opportunity to understand how these essential environmental microbes regulate their metabolism. Other syntrophs also contain critical pathways with RACS intermediates (
), and a wider investigation of other syntrophs may further reveal the role RACS play in bacterial metabolism.
Data Availability
All data were uploaded to ProteomeXchange (proteomexchange.org) either through PRIDE (2D-gel dataset; PXD025631) or through MassIVE (Shotgun dataset; PXD025603).
We are grateful for the 2D-PAGE and mass spectrometry analyses performed by Dr Yanan Yang. We would like to thank Dr Mark Arbing and Annie Shin of the UCLA-DOE Institute Protein Expression Technology Center (supported by grant 447147-EL-21341 from the US Department of Energy) for expression and purification of Syntrophus aciditrophicus SYN_00042 and SYN_01020.
Funding and additional information
Funding from the Department of Energy Office of Science (BER) contract DE-FC-02-02ER63421 (to J. A. L. and R. P. G; UCLA/DOE Institute for Genomics and Proteomics), NIH Ruth L. Kirschstein National Research Service 18 Award (to J. Y. F.; GM007185), NSF Award 1911781 to R. P. G and M. J. M, and NSF Graduate Research Fellowship (to J. Y. F.; DGE-1650604) are acknowledged. Also acknowledged are the National Institute of Health awards GM104610 and GM085402 to R. O. L. and J. A. L. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Author contributions
M. J. M., R. P. G., J. A. L., and R. R. O. L. conceptualization; J. M. M., H. H. N., M. J. M., R. P. G., J. A. L., and R. R. O. L. methodology; J. M. M., J. Y. F., H. H. N., N. Q. W., H. M., K. L. J., and R. R. O. L. investigation; J. M. M., J. Y. F., H. H. N., and R. R. O. L. data curation; J. M. M., J. Y. F., N. Q. W., H. M., K. L. J., M. J. M., R. P. G., and R. R. O. L. formal analysis; J. M. M., J. Y. F., M. J. M., R. P. G., J. A. L., and R. R. O. L. writing–original draft; J. M. M., J. Y. F., N. Q. W., H. M., K. L. J., M. J. M., R. P. G., J. A. L., and R. R. O. L. writing–review and editing.
Syntrophus aciditrophicus sp. nov., a new anaerobic bacterium that degrades fatty acids and benzoate in syntrophic association with hydrogen-using microorganisms.
Cyclohexa-1,5-diene-1-carbonyl-coenzyme A (CoA) hydratases of Geobacter metallireducens and Syntrophus aciditrophicus: Evidence for a common benzoyl-CoA degradation pathway in facultative and strict anaerobes.
6-Oxocyclohex-1-ene-1-carbonyl-coenzyme A hydrolases from obligately anaerobic bacteria: Characterization and identification of its gene as a functional marker for aromatic compounds degrading anaerobes.
Metabolism of benzoate, cyclohex-1-ene carboxylate, and cyclohexane carboxylate by Syntrophus aciditrophicus strain SB in syntrophic association with H2-using microorganisms.
in: Hurst C. Crawford R. Garland J. Lipson D. Mills A. Stetzenbach L. Manual of Environmental Microbiology. 3rd Ed. American Society of Microbiology,
Washington, DC2007: 69-78
New approach to the cultivation of methanogenic bacteria: 2 Mercaptoethanesulfonic acid (HS CoM) dependent growth of methanobacterium ruminantium in a pressurized atmosphere.
Heterologous expression of mycobacterial Esx complexes in Escherichia coli for structural studies is facilitated by the use of maltose binding protein fusions.
Biochemical and structural characterization of the paralogous benzoate CoA ligases from Burkholderia xenovorans LB400: Defining the entry point into the novel benzoate oxidation (box) pathway.
Reversible Nε-lysine acetylation regulates the activity of acyl-CoA synthetases involved in anaerobic benzoate catabolism in Rhodopseudomonas palustris.
00023-8&id=fx1.jpg){kind=link}
00023-8&id=gr1.jpg){kind=link}
00023-8&id=gr2.jpg){kind=link}
00023-8&id=gr3.jpg){kind=link}
00023-8&id=gr4.jpg){kind=link}
00023-8&id=gr5.jpg){kind=link}