Search

Short and branched chain fatty acid kinases participate in both bacterial anabolic and catabolic processes, including fermentation, through the reversible, ATP-dependent synthesis of acyl phosphates. This study reports biochemical properties of a predicted butyrate kinase from Desulfovibrio vulgaris str. Hildenborough (DvBuk) expressed heterologously and purified from Escherichia coli. Gel filtration chromatography indicates purified DvBuk is active as a dimer. The optimum temperature and pH for DvBuk activity is 44°C and 7.5, respectively. The enzyme displays enhanced thermal stability in the presence of substrates as observed for similar enzymes. Measurement of kcat and KM for various substrates reveals DvBuk exhibits the highest catalytic efficiencies for butyrate, valerate and isobutyrate. In particular, these measurements reveal this enzyme's apparent high affinity for C4 fatty acids relative to other butyrate kinases. These results have implications on structure and function relationships within the ASKHA superfamily of phosphotransferases, particularly regarding the acyl binding pocket, as well as potential physiological roles for this enzyme in Desulfovibrio vulgaris str. Hildenborough., (© FEMS 2020.)

Although as many as half of all proteins are thought to require a metal cofactor, the metalloproteomes of microorganisms remain relatively unexplored. Microorganisms from different environments are likely to vary greatly in the metals that they assimilate, not just among the metals with well-characterized roles but also those lacking any known function. Herein we investigated the metal utilization of two microorganisms that were isolated from very similar environments and are of interest because of potential roles in the immobilization of heavy metals, such as uranium and chromium. The metals assimilated and their concentrations in the cytoplasm of Desulfovibrio vulgaris strain Hildenborough (DvH) and Enterobacter cloacae strain Hanford (EcH) varied dramatically, with a larger number of metals present in Enterobacter. For example, a total of 9 and 19 metals were assimilated into their cytoplasmic fractions, respectively, and DvH did not assimilate significant amounts of zinc or copper whereas EcH assimilated both. However, bioinformatic analysis of their genome sequences revealed a comparable number of predicted metalloproteins, 813 in DvH and 953 in EcH. These allowed some rationalization of the types of metal assimilated in some cases (Fe, Cu, Mo, W, V) but not in others (Zn, Nd, Ce, Pr, Dy, Hf and Th). It was also shown that U binds an unknown soluble protein in EcH but this incorporation was the result of extracellular U binding to cytoplasmic components after cell lysis.

Motivation: Gene expression profiling technologies can generally produce mRNA abundance data for all genes in a genome. A dearth of proteomic data persists because identification range and sensitivity of proteomic measurements lag behind those of transcriptomic measurements. Using partial proteomic data, it is likely that integrative transcriptomic and proteomic analysis may introduce significant bias. Developing methodologies to accurately estimate missing proteomic data will allow better integration of transcriptomic and proteomic datasets and provide deeper insight into metabolic mechanisms underlying complex biological systems., Results: In this study, we present a non-linear data-driven model to predict abundance for undetected proteins using two independent datasets of cognate transcriptomic and proteomic data collected from Desulfovibrio vulgaris. We use stochastic gradient boosted trees (GBT) to uncover possible non-linear relationships between transcriptomic and proteomic data, and to predict protein abundance for the proteins not experimentally detected based on relevant predictors such as mRNA abundance, cellular role, molecular weight, sequence length, protein length, guanine-cytosine (GC) content and triple codon counts. Initially, we constructed a GBT model using all possible variables to assess their relative importance and characterize the behavior of the predictive model. A strong plateau effect in the regions of high mRNA values and sparse data occurred in this model. Hence, we removed genes in those areas based on thresholds estimated from the partial dependency plots where this behavior was captured. At this stage, only the strongest predictors of protein abundance were retained to reduce the complexity of the GBT model. After removing genes in the plateau region, mRNA abundance, main cellular functional categories and few triple codon counts emerged as the top-ranked predictors of protein abundance. We then created a new tuned GBT model using the five most significant predictors. The construction of our non-linear model consists of a set of serial regression trees models with implicit strength in variable selection. The model provides variable relative importance measures using as a criterion mean square error. The results showed that coefficients of determination for our nonlinear models ranged from 0.393 to 0.582 in both datasets, providing better results than linear regression used in the past. We evaluated the validity of this non-linear model using biological information of operons, regulons and pathways, and the results demonstrated that the coefficients of variation of estimated protein abundance values within operons, regulons or pathways are indeed smaller than those for random groups of proteins., Supplementary Information: Supplementary data are available at Bioinformatics online.

Hypothetical (HyP) and conserved HyP genes account for >30% of sequenced bacterial genomes. For the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough, 347 of the 3634 genes were annotated as conserved HyP (9.5%) along with 887 HyP genes (24.4%). Given the large fraction of the genome, it is plausible that some of these genes serve critical cellular roles. The study goals were to determine which genes were expressed and provide a more functionally based annotation. To accomplish this, expression profiles of 1234 HyP and conserved genes were used from transcriptomic datasets of 11 environmental stresses, complemented with shotgun LC-MS/MS and AMT tag proteomic data. Genes were divided into putatively polycistronic operons and those predicted to be monocistronic, then classified by basal expression levels and grouped according to changes in expression for one or multiple stresses. One thousand two hundred and twelve of these genes were transcribed with 786 producing detectable proteins. There was no evidence for expression of 17 predicted genes. Except for the latter, monocistronic gene annotation was expanded using the above criteria along with matching Clusters of Orthologous Groups. Polycistronic genes were annotated in the same manner with inferences from their proximity to more confidently annotated genes. Two targeted deletion mutants were used as test cases to determine the relevance of the inferred functional annotations.

The gene encoding a MutM from Desulfovibrio vulgaris (Miyazaki F) was cloned and expressed in Escherichia coli. A 5.9-kb DNA fragment, isolated from D. vulgaris (Miyazaki F) by XhoI and PvuII, contained a MutM gene and other open reading frames. The nucleotide sequence of the MutM gene indicated that the protein was composed of 336 amino acids. The amino-acid sequence deduced from the MutM gene was highly homologous with the MutM of other bacteria; however an additional insert consisted of 64 amino acids. An expression system for the MutM gene under the control of the T7 promoter was constructed in E. coli. From the kinetic analysis results, the purified His-tagged MutM showed 8-oxoguanine-DNA glycosylase activity comparable with that of MutM from E. coli. In this study, the amounts of mRNA and protein for MutM were scant in the D. vulgaris (Miyazaki F). MutM activity may be induced by oxidative stress. However, its induction may not be frequently generated because sulfate-reducing bacteria generally grow in anaerobic conditions. MutM might play a role in the protection against the mutagenicity of oxygen when oxygen stress exceeded the capacity of the defense systems against oxygen toxicity.

Mutants of flavin mononucleotide-binding protein (FMN-bp) were made by site-directed mutagenesis to investigate the role of carboxyl-terminal Leu122 of the pairing subunit in controlling redox potentials, binding the prosthetic group, and forming the tertiary and quaternary structure. We compared the oxidation-reduction potentials, FMN-binding properties, and higher structures of wild-type FMN-bp and four mutant proteins (L122Y, L122E, L122K and L122-deleted). We found that the redox potentials were affected by mutations. Also, the affinities of L122E, L122K and L122 deletion mutant apoproteins for FMN were lower than for the wild-type apoprotein, whereas the affinity of L122Y for FMN was increased. Analytical ultracentrifugation showed that the dissociation constants for dimerization of L122E and L122K were larger than for wild-type FMN-bp, whereas the dissociation constants for L122Y and the deletion mutant were lower than for the wild type. Finally, we determined the higher structures of L122Y, L122E and L122K mutants by X-ray crystallography. Our results show that the mutation of Leu122 in FMN-bp changes midpoint potentials, dissociation constants for FMN, and dimer formation, indicating that this residue is important in the pairing subunit.

The modest correlation between mRNA expression and protein abundance in large-scale data sets is explained in part by experimental challenges, such as technological limitations, and in part by fundamental biological factors in the transcription and translation processes. Among various factors affecting the mRNA-protein correlation, the roles of biological factors related to translation are poorly understood. In this study, using experimental mRNA expression and protein abundance data collected from Desulfovibrio vulgaris by DNA microarray and liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) proteomic analysis, we quantitatively examined the effects of several translational-efficiency-related sequence features on mRNA-protein correlation. Three classes of sequence features were investigated according to different translational stages: (i) initiation, Shine-Dalgarno sequences, start codon identity, and start codon context; (ii) elongation, codon usage and amino acid usage; and (iii) termination, stop codon identity and stop codon context. Surprisingly, although it is widely accepted that translation initiation is the rate-limiting step for translation, our results showed that the mRNA-protein correlation was affected the most by the features at elongation stages, i.e., codon usage and amino acid composition (5.3-15.7% and 5.8-11.9% of the total variation of mRNA-protein correlation, respectively), followed by stop codon context and the Shine-Dalgarno sequence (3.7-5.1% and 1.9-3.8%, respectively). Taken together, all sequence features contributed to 15.2-26.2% of the total variation of mRNA-protein correlation. This study provides the first comprehensive quantitative analysis of the mRNA-protein correlation in bacterial D. vulgaris and adds new insights into the relative importance of various sequence features in prokaryotic protein translation.

Motivation: Integrated analysis of global scale transcriptomic and proteomic data can provide important insights into the metabolic mechanisms underlying complex biological systems. However, because the relationship between protein abundance and mRNA expression level is complicated by many cellular and physical processes, sophisticated statistical models need to be developed to capture their relationship., Results: In this study, we describe a novel data-driven statistical model to integrate whole-genome microarray and proteomic data collected from Desulfovibrio vulgaris grown under three different conditions. Based on the Poisson distribution pattern of proteomic data and the fact that a large number of proteins were undetected (excess zeros), zero-inflated Poisson (ZIP)-based models were proposed to define the correlation pattern between mRNA and protein abundance. In addition, by assuming that there is a probability mass at zero representing unexpressed genes and expressed proteins that were undetected owing to technical limitations, a Potential ZIP model was established. Two significant improvements introduced by this approach are (1) the predicted protein abundance level values for experimentally detected proteins are corrected by considering their mRNA levels and (2) protein abundance values can be predicted for undetected proteins (in the case of this study, approximately 83% of the proteins in the D.vulgaris genome) for better biological interpretation. We demonstrated the use of these statistical models by comparatively analyzing proteomic and microarray results from D.vulgaris grown on lactate-based versus formate-based media. These models correctly predicted increased expression of Ech hydrogenase and decreased expression of Coo hydrogenase for D.vulgaris grown on formate.

Peculiar attributes revealed by sequencing the genome of Desulfovibrio vulgaris Hildenborough are analyzed, particularly in relation to the presence of a phosphotransferase system (PTS). The PTS is a typical bacterial carbohydrate transport system functioning via group translocation. Novel avenues for investigations are proposed emphasizing the metabolic diversity of D. vulgaris Hildenborough, especially the likely utilization of mannose-type sugars. Comparative analysis with PTS from other Gram-negative and Gram-positive bacteria indicates regulatory functions for the PTS of D. vulgaris Hildenborough, including catabolite repression and inducer exclusion. Chemotaxis towards PTS substrates is considered. Evidence suggests that this organism may not be a strict anaerobic sulfate reducer typical of the ocean, but a versatile organism capable of bidirectional transmigration and adaptation to both water and terrestrial environments.

The response of Desulfovibrio vulgaris Hildenborough (DvH), a sulphate-reducing bacterium, to nitrate stress was examined using quantitative proteomic analysis. DvH was stressed with 105 mM sodium nitrate (NaNO(3)), a level that caused a 50% inhibition in growth. The protein profile of stressed cells was compared with that of cells grown in the absence of nitrate using the iTRAQ peptide labelling strategy and tandem liquid chromatography separation coupled with mass spectrometry (quadrupole time-of-flight) detection. A total of 737 unique proteins were identified by two or more peptides, representing 22% of the total DvH proteome and spanning every functional category. The results indicate that this was a mild stress, as proteins involved in central metabolism and the sulphate reduction pathway were unperturbed. Proteins involved in the nitrate reduction pathway increased. Increases seen in transport systems for proline, glycine-betaine and glutamate indicate that the NaNO(3) exposure led to both salt stress and nitrate stress. Up-regulation observed in oxidative stress response proteins (Rbr, RbO, etc.) and a large number of ABC transport systems as well as in iron-sulphur-cluster-containing proteins, however, appear to be specific to nitrate exposure. Finally, a number of hypothetical proteins were among the most significant changers, indicating that there may be unknown mechanisms initiated upon nitrate stress in DvH.

We identified the SOD gene in the obligate anaerobic bacterium Desulfovibrio vulgaris (Miyazaki F) and constructed a high-level expression system in Escherichia coli. A 2.6-kbp DNA fragment isolated from D. vulgaris (Miyazaki F) by double digestion with EcoRI and SmaI contained the SOD gene and part of another open reading frame. The amino acid sequence deduced from the SOD gene, which was composed of 238 amino acid residues, showed high homogeneity with iron-containing SOD (Fe-SOD) and predicted that the amino terminus of this protein would carry an export signal peptide. We produced the precursor form of SOD (PSOD) and the mature form of SOD (MSOD), which lacked the putative signal peptide. In E. coli, PSOD was present in insoluble inclusion bodies, and its putative signal peptide was not cleaved. In contrast, MSOD contained one iron per mononer and formed a dimer, which exhibited an SOD activity of 850 U/mg. Furthermore, D. vulgaris soluble extract showed a band of SOD activity on native polyacrylamide gel that migrated to the same point as MSOD. The intracellular localization of SOD and its role in D. vulgaris are also discussed.

We identified a gene encoding a catalase from the anaerobic bacteria Desulfovibrio vulgaris (Miyazaki F), and the expression of its gene in Escherichia coli. The 3.3-kbp DNA fragment isolated from D. vulgaris (Miyazaki F) by double digestion with EcoRI and SalI was found to produce a protein that binds protoheme IX as a prosthetic group in E. coli. This DNA fragment contained a putative open reading frame (Kat) and one part of another open reading frame (ORF-1). The amino acid sequence of the amino terminus of the protein purified from the transformed cells was consistent with that deduced from the nucleotide sequence of Kat in the cloned fragment of D. vulgaris (Miyazaki F) DNA, which may include promoter and regulatory sequences. The nucleotide sequence of Kat indicates that the protein is composed of 479 amino acids per monomer. The recombinant catalase was found to be active in the decomposition of hydrogen peroxide, as are other catalases from aerobic organisms, but its K(m) value was much greater. The hydrogen peroxide stress against D. vulgaris (Miyazaki F) induced the activity for the decomposition of hydrogen peroxide somewhat, so the catalase gene may not work effectively in vivo.

The gene encoding a flavodoxin of Desulfovibrio vulgaris (Miyazaki F) was cloned, and overexpressed in Escherichia coli. A 1.6-kbp DNA fragment, isolated from D. vulgaris (Miyazaki F) by double digestion with SalI and EcoRI, contained the flavodoxin gene and its regulatory region. An expression system for the flavodoxin gene under control of the T7 promoter was constructed in E. coli. The purified protein was soluble and exhibited a characteristic visible absorption spectrum. HPLC analysis of the recombinant flavodoxin revealed the presence of an identical FMN to that found in the native D. vulgaris flavodoxin, and its dissociation constant with FMN was determined to be 0.38 nM. In vitro H2 reduction analysis indicated that the recombinant flavodoxin is active, and its redox potential was determined to be E1 = -434 and E2 = -151 mV using methyl viologen and 2-hydroxy-1,4-naphthoquinone, respectively. Its redox behavior was also examined with the recombinant flavodoxin adsorbed onto a graphite electrode. The mutant, A16E, was also produced, which revealed the feature of a conserved Glu residue at the surface of the molecule.

A comparative analysis of mRNA decay was carried out in Escherichia coli using the wild-type and an isogenic RNase III deletion strain. We have studied the mRNA degradation from the Escherichia coli gene bolA, the Lactococcus lactis biovar diacetylactis citQRP operon and the Desulfovibrio vulgaris Hildenborough gene cyc. As seen by a dramatic stabilization of the specific mRNAs in the mutant strain, RNase III was crucial for the decay process of these three messages. Since RNase III, unlike RNase E, is not essential for bacterial viability we think that there is potential for using RNase III mutant strains to modulate gene expression.

The effect of Escherichia coli ribonuclease II and polynucleotide phosphorylase was analysed on the degradation of Desulfovibrio vulgaris cytochrome c3 (cyc) mRNA. In the absence of these exoribonucleolytic activities, cyc mRNA was stabilised but the two enzymes had a different role in its decay. Surprisingly, a temperature-sensitive mutation in ribonuclease II gave a degradation pattern similar to what had been observed in the absence of endoribonuclease E activity. In an RNase II deletion mutant this was not observed. We propose and verify a model in which the temperature-sensitive ribonuclease II interferes with the action of ribonuclease E.

The genes for the subunits of the Fe-only hydrogenase from Desulfovibrio vulgaris are transcribed as a 1.9 kb mRNA; the operon contains no other genes besides those encoding the two subunits. The transcriptional start site of the operon was mapped. Determination of hydrogenase activity and hydrogenase mRNA levels indicates a growth-phase dependent regulation of hydrogenase expression at transcriptional level. However, it has not yet been possible to localize the sequences required for regulation and expression of the genes.