Decimal Place Slope,A Fast and Precise Method for Quantifying 13C Incorporation Levels for Detecting the Metabolic Activity of Microbial Species

The metabolic incorporation of stable isotopes such as ¹³C or ¹⁵N into proteins has become a powerful tool for qualitative and quantitative proteome studies. We recently introduced a method that monitors heavy isotope incorporation into proteins and presented data revealing the metabolic activity of various species in a microbial consortium using this technique. To further develop our method using an liquid chromatography (LC)-mass spectrometry (MS)-based approach, we present here a novel approach for calculating the incorporation level of ¹³C into peptides by using the information given in the decimal places of peptide masses obtained by modern high-resolution MS. In the present study, the applicability of this approach is demonstrated using Pseudomonas putida ML2 proteins uniformly labeled via the consumption of [¹³C₆]benzene present in the medium at concentrations of 0, 10, 25, 50, and 100 atom %. The incorporation of ¹³C was calculated on the basis of several labeled peptides derived from one band on an SDS-PAGE gel. The accuracy of the calculated incorporation level depended upon the number of peptide masses included in the analysis, and it was observed that at least 100 peptide masses were required to reduce the deviation below 4 atom %. This accuracy was comparable with calculations of incorporation based on the isotope envelope. Furthermore, this method can be extended to the calculation of the labeling efficiency for a wide range of biomolecules, including RNA and DNA. The technique will therefore allow a highly accurate determination of the carbon flux in microbial consortia with a direct approach based solely on LC-MS.The metabolic incorporation of stable isotopes such as ¹³C or ¹⁵N into proteins has become a powerful component of qualitative and quantitative proteome studies (1). Incorporation of heavy isotopes can be used to analyze microbial processes such as turnover rates and also to help to establish structure-function relationships within microbial communities. Stable isotope probing (SIP¹) techniques based on DNA-SIP (2) and RNA-SIP (3) have been used for this purpose previously. With the introduction of protein-SIP (4), the need for an accurate alternative method for calculating label incorporation into biomolecules arose. Protein-SIP has several advantages compared with DNA/RNA-SIP, the most important being its capacity to detect dynamic levels of incorporation, whereas only labeled or unlabeled states can be categorized by means of DNA/RNA-SIP because of the need to separate ¹³C-DNA/RNA by density gradient centrifugation. Quantitative analysis of ¹³C incorporation levels is of the utmost importance, especially when unraveling carbon fluxes through either microbial communities or food webs with different trophic levels.In contrast to the incorporation of isotopically labeled amino acids, which is often used in quantitative proteomics (5), metabolic labeling by growth substrates and nutrients (e.g. salts) is often imperfect and makes the processing of mass spectrometry (MS) data difficult. For example, when the incorporation of ¹³C exceeds ∼2 atom %, common database search algorithms fail to identify peptides and proteins. The problem can only be managed successfully if a stable, known degree of ¹³C incorporation can be achieved during the experiment (6). Using a low labeling efficiency of roughly 5 atom %, Huttlin et al. (6) chose the altered envelope chain for calculating the incorporation and simultaneously used the signal intensity for a quantitative comparison with the sample that had a natural abundance of ¹³C. Database approaches for peptide identification can cope only with the natural abundance of carbon isotopes; they fail if the incorporation of ¹³C significantly exceeds the natural isotope abundance or if incorporation patterns occur in unpredictable ways (7).The simplest method for determining the incorporation level is to compare the unlabeled average mass of the monoisotopic peptide with the mass of the labeled protein, as estimated by matrix-assisted laser desorption/ionization or electrospray ionization MS (8, 9). A more advanced approach for determining the isotopic mass distribution of peptides is based on the isotopic distribution of the peaks of a peptide envelope (10, 11). Here, for a given isotopomer, the incorporation efficiency is defined as the percentage of incorporated ¹³C atoms with relation to the total number of carbon atoms with the natural isotope abundance (approximately 1.01 atom % ¹³C). As a reference, the theoretical isotopic distribution of a peptide is calculated based upon an algorithm described elsewhere (12). The isotope distribution of both unlabeled and labeled peptides can subsequently be used to calculate the incorporation level. For this method, an Excel spreadsheet (ProSIPQuant.xls) was developed (4). A similar approach, also based on the calculation of isotopic distributions, has been used in other studies (7). In these studies, however, the identification of the peptides is limited to those that have unlabeled counterparts; in addition, an exact calculation can be hampered by overlapping signals coming from additional peaks with similar masses.In the present study, we describe a new way of determining the isotope incorporation level. Our method makes use of characteristic patterns in the digits after the decimal point of the peptide masses generated by high-accuracy instruments such as the linear ion trap LTQ-Orbitrap (Thermo Fisher Scientific, Bremen, Germany). For tryptic peptides, typical regularities in the decimal places of the monoisotopic masses have been observed (13, 14). These observations have been explored in detail for theoretical and experimental data of proteins originating from Helicobacter pylori (15). As a result, a rule called the “half decimal place rule” (HDPR) was defined; it states that the decimal place is nearly half of the first digit for tryptic peptides with masses in the range of 500–1,000 Da. In other words, the exact mass of a peptide is equal to its nominal mass times ∼1.005. Because the difference between ¹²C and ¹³C is slightly greater than 1 Da, exactly 1.0033548378, the decimal places of a tryptic peptide''s mass are shifted in a regular manner by the incorporation level and lead to a significantly increased slope for the digits in the third and fourth place after the decimal point. This shift can be used to estimate the incorporation level of heavy isotopes into the protein. Detecting such shifts requires the highly accurate measurement possible with modern mass spectrometers such as the LTQ-Orbitrap, the Fourier transform ion cyclotron resonance, or the quadrupole time of flight. In this communication, we demonstrate the applicability of this approach using Pseudomonas putida ML2 proteins labeled uniformly via the consumption of [¹³C₆]benzene with five different substrate concentrations (0, 10, 25, 50, and 100 atom % of ¹³C). The ¹³C incorporation was calculated based on several labeled peptides derived from different proteins in one SDS-PAGE band. By these means, we have established a method that allows the determination of ¹³C incorporation into proteins and can be used to assess the metabolic activity of a given species within a mixed community.     

Microbial communities involved in anaerobic degradation of unsaturated or saturated long-chain fatty acids

Anaerobic long-chain fatty acid (LCFA)-degrading bacteria were identified by combining selective enrichment studies with molecular approaches. Two distinct enrichment cultures growing on unsaturated and saturated LCFAs were obtained by successive transfers in medium containing oleate and palmitate, respectively, as the sole carbon and energy sources. Changes in the microbial composition during enrichment were analyzed by denaturing gradient gel electrophoresis (DGGE) profiling of PCR-amplified 16S rRNA gene fragments. Prominent DGGE bands of the enrichment cultures were identified by 16S rRNA gene sequencing. A significant part of the retrieved 16S rRNA gene sequences was most similar to those of uncultured bacteria. Bacteria corresponding to predominant DGGE bands in oleate and palmitate enrichment cultures clustered with fatty acid-oxidizing bacteria within Syntrophomonadaceae and Syntrophobacteraceae families. A low methane yield, corresponding to 9 to 18% of the theoretical value, was observed in the oleate enrichment, and acetate, produced according to the expected stoichiometry, was not further converted to methane. In the palmitate enrichment culture, the acetate produced was completely mineralized and a methane yield of 48 to 70% was achieved from palmitate degradation. Furthermore, the oleate enrichment culture was able to use palmitate without detectable changes in the DGGE profile. However, the palmitate-specialized consortia degraded oleate only after a lag phase of 3 months, after which the DGGE profile had changed. Two predominant bands appeared, and sequence analysis showed affiliation with the Syntrophomonas genus. These bands were also present in the oleate enrichment culture, suggesting that these bacteria are directly involved in oleate degradation, emphasizing possible differences between the degradation of unsaturated and saturated LCFAs.     

Dual selection pressure by drugs and HLA class I-restricted immune responses on human immunodeficiency virus type 1 protease

To determine the influence of human immunodeficiency virus type 1 (HIV-1)-specific CD8+ T cells on the development of drug resistance mutations in the HIV-1 protease, we analyzed protease sequences from viruses from a human leukocyte antigen class I (HLA class I)-typed cohort of 94 HIV-1-positive individuals. In univariate statistical analyses (Fisher's exact test), minor and major drug resistance mutations as well as drug-associated polymorphisms showed associations with HLA class I alleles. All correlations with P values of 0.05 or less were considered to be relevant without corrections for multiple tests. A subset of these observed correlations was experimentally validated by enzyme-linked immunospot assays, allowing the definition of 10 new epitopes recognized by CD8+ T cells from patients with the appropriate HLA class I type. Several drug resistance-associated mutations in the protease acted as escape mutations; however, cells from many patients were still able to generate CD8+ T cells targeting the escape mutants. This result presumably indicates the usage of different T-cell receptors by CD8+ T cells targeting these epitopes in these patients. Our results support a fundamental role for HLA class I-restricted immune responses in shaping the sequence of the HIV-1 protease in vivo. This role may have important clinical implications both for the understanding of drug resistance pathways and for the design of therapeutic vaccines targeting drug-resistant HIV-1.     

Carbon Conversion Efficiency and Limits of Productive Bacterial Degradation of Methyl tert-Butyl Ether and Related Compounds

The utilization of the fuel oxygenate methyl tert-butyl ether (MTBE) and related compounds by microorganisms was investigated in a mainly theoretical study based on the Y_ATP concept. Experiments were conducted to derive realistic maintenance coefficients and K_s values needed to calculate substrate fluxes available for biomass production. Aerobic substrate conversion and biomass synthesis were calculated for different putative pathways. The results suggest that MTBE is an effective heterotrophic substrate that can sustain growth yields of up to 0.87 g g⁻¹, which contradicts previous calculation results (N. Fortin et al., Environ. Microbiol. 3:407-416, 2001). Sufficient energy equivalents were generated in several of the potential assimilatory routes to incorporate carbon into biomass without the necessity to dissimilate additional substrate, efficient energy transduction provided. However, when a growth-related kinetic model was included, the limits of productive degradation became obvious. Depending on the maintenance coefficient m_s and its associated biomass decay term b, growth-associated carbon conversion became strongly dependent on substrate fluxes. Due to slow degradation kinetics, the calculations predicted relatively high threshold concentrations, S_min, below which growth would not further be supported. S_min strongly depended on the maximum growth rate μ_max, and b and was directly correlated with the half maximum rate-associated substrate concentration K_s, meaning that any effect impacting this parameter would also change S_min. The primary metabolic step, catalyzing the cleavage of the ether bond in MTBE, is likely to control the substrate flux in various strains. In addition, deficits in oxygen as an external factor and in reduction equivalents as a cellular variable in this reaction should further increase K_s and S_min for MTBE.     

Induction of heme oxygenase-1 and heat shock protein 70 in rat hepatocytes: The role of calcium signaling

Stress response genes including heat shock proteins are induced under a variety of conditions to confer cellular protection. This study investigated the role of calcium signaling in the induction of two stress response genes, heme oxygenase-1/hsp32 and hsp70, in isolated rat hepatocytes. Both genes were induced by cellular glutathione depletion. This induction could be inhibited by BAPTA-AM. Culturing in a calcium-free medium prevented the induction of hsp70 gene expression after glutathione depletion without affecting heme oxygenase-1 gene expression. Thapsigargin increased the gene expression of heme oxygenase-1 but not that of hsp70. Thapsigargin-induced heme oxygenase-1 induction was completely inhibited by BAPTA-AM. Incubation with the Ca²⁺-ionophore A23187 augmented heme oxygenase-1 (two-fold) and hsp70 (5.2-fold) mRNA levels. Our data suggests a significant role of Ca²⁺-dependent pathways in the induction of the two stress genes. An increase in the cytoplasmic Ca²⁺ activity seems to play a key role in the cascade of signaling leading to the induction of the two genes. However, the source of Ca²⁺ that fluxes into the cytoplasm seems to be different. Our data provides evidence for a compartmentalization of calcium fluxes, i.e. the Ca²⁺ flux from intracellular stores (e.g. the endoplasmic reticulum) plays a major role in the induction of heme oxygenase-1. By contrast, Ca²⁺ flux from the extracellular medium seems to be a mechanism initiating the cellular signaling cascade leading to hsp70 gene induction.     

Long-distance dispersal of a wolf, <Emphasis Type="Italic">Canis lupus</Emphasis>, in northwestern Europe

Several mammal species have recolonized their historical ranges across Europe during the last decades. In November 2012, a wolf-looking canid was found dead in Thy National Park (56° 56′ N, 8° 25′ E) in Jutland, Denmark. DNA from this individual and nine German wolves were genotyped using a genome-wide panel of 22,163 canine single nucleotide polymorphism (SNP) markers and compared to existing profiles based on the same marker panel obtained from northeastern Polish (n?=?13) wolves, domestic dogs (n?=?13) and known wolf-dog hybrids (n?=?4). The Thy canid was confirmed to be a wolf from the German-western Polish population, approximately 800 km to the southeast. Access to the German reference database on DNA profiles based on 13 autosomal microsatellites of German wolves made it possible to pinpoint the exact pack origin of the Thy wolf in Saxony, Germany. This was the first documented observation of a wolf in Denmark in 200 years and another example of long-distance dispersal of a carnivore.