Development of an integrated metabolomic profiling approach for infectious diseases research

Metabolomic profiling offers direct insights into the chemical environment and metabolic pathway activities at sites of human disease. During infection, this environment may receive important contributions from both host and pathogen. Here we apply an untargeted metabolomics approach to identify compounds associated with an E. coli urinary tract infection population. Correlative and structural data from minimally processed samples were obtained using an optimized LC-MS platform capable of resolving ~2300 molecular features. Principal component analysis readily distinguished patient groups and multiple supervised chemometric analyses resolved robust metabolomic shifts between groups. These analyses revealed nine compounds whose provisional structures suggest candidate infection-associated endocrine, catabolic, and lipid pathways. Several of these metabolite signatures may derive from microbial processing of host metabolites. Overall, this study highlights the ability of metabolomic approaches to directly identify compounds encountered by, and produced from, bacterial pathogens within human hosts.

Simultaneous determination by UPLC-ESI-MS of scoparone, capillarisin, rhein, and emodin in rat urine after oral administration of Yin Chen Hao Tang preparation

A simple, sensitive, and validated method was developed for simultaneous determination of scoparone, capillarisin, rhein, and emodin in rat urine by ultra-performance liquid chromatography/electrospray ionization quadruple time-of-flight mass spectrometry (UPLC-MS). The urinary samples were analyzed on an Acquity UPLC BEH C18 1.7 microm 2.1x50 mm column. Scoparone, capillarisin, rhein, and emodin in rat urine were simultaneously analyzed with good separation. The lower limits of detection were 6.0, 9.0, 7.0, and 3.0 ng/mL, and the lower limits of quantification were 20.0, 33.0, 24.0, and 12.0 ng/mL for scoparone, capillarisin, rhein, and emodin, respectively. The intra- and inter-day precisions (RSD) were less than 9%. The intra- and inter-accuracies were found to be in the range of 94.14-104.54% for scoparone, 101.72-107.34% for capillarisin, 95.24-103.59% for rhein, and 101.32-107.82% for emodin at three concentration levels. The absolute recoveries for scoparone, capillarisin, rhein, and emodin were not less than 77.0%. The developed method has been applied to determine scoparone, capillarisin, rhein, and emodin in rat urine after oral administration of Yin Chen Hao Tang preparation, a traditional Chinese medicine formulation widely used in China for treatment of jaundice and liver disorders.

Analysis of the constituents in the rat plasma after oral administration of Yin Chen Hao Tang by UPLC/Q-TOF-MS/MS

A UPLC/Q-TOF-MS/MS method for analyzing the constituents in rat plasma after oral administration of Yin Chen Hao Tang (YCHT), a traditional Chinese medical formula, has been established. The UPLC/MS fingerprints of the samples were established first in vitro and in vivo, with 45 compounds in YCHT and 21 compounds in rat plasma after oral administration of YCHT were detected. Of the 45 detected compounds in vitro, 30 were identified, and all of the 21 compounds detected in rat plasma were identified either by comparing the retention time and mass spectrometry data with that of reference compounds or by mass spectrometry analysis and retrieving the reference literatures. Of the identified 21 compounds in rat plasma, 19 were the original form of compounds absorbed from the 45 detected compounds in vitro, 2 were the metabolites of the compounds existed in YCHT. It is concluded that a rapid and validated method has been developed based on UPLC-MS/MS, which shows high sensitivity and resolution that is more suitable for identifying the bioactive constituents in plasma after oral administration of Chinese herbal medicines, and provides helpful chemical information for further pharmacology and active mechanism research on the Chinese medical formula.

Metabolite profiling and pathway analysis of acute hepatitis rats by UPLC-ESI MS combined with pattern recognition methods

BACKGROUND & AIMS Metabolomics is comprehensive analysis of low-molecular-weight endogenous metabolites in a biological sample. It could enable mapping of perturbations of early biochemical changes in diseases and hence provide an opportunity to develop predictive biomarkers that could provide valuable insights into the mechanisms of diseases. The aim of this study was to elucidate the changes in endogenous metabolites and to phenotype the metabolic profiling of d-galactosamine (GalN)-inducing acute hepatitis in rats by UPLC-ESI MS. METHODS The systemic biochemical actions of GalN administration (ip, 400 mg/kg) have been investigated in male wistar rats using conventional clinical chemistry, liver histopathology and metabolomic analysis of UPLC- ESI MS of urine. The urine was collected predose (-24 to 0 h) and 0-24, 24-48, 48-72, 72-96 h post-dose. Mass spectrometry of the urine was analysed visually and via conjunction with multivariate data analysis. RESULTS Results demonstrated that there was a time-dependent biochemical effect of GalN dosed on the levels of a range of low-molecular-weight metabolites in urine, which was correlated with developing phase of the GalN-inducing acute hepatitis. Urinary excretion of beta-hydroxybutanoic acid and citric acid was decreased following GalN dosing, whereas that of glycocholic acid, indole-3-acetic acid, sphinganine, n-acetyl-l-phenylalanine, cholic acid and creatinine excretion was increased, which suggests that several key metabolic pathways such as energy metabolism, lipid metabolism and amino acid metabolism were perturbed by GalN. CONCLUSION This metabolomic investigation demonstrates that this robust non-invasive tool offers insight into the metabolic states of diseases.

Multiple sclerosis susceptibility-associated SNPs do not influence disease severity measures in a cohort of Australian MS patients

Recent association studies in multiple sclerosis (MS) have identified and replicated several single nucleotide polymorphism (SNP) susceptibility loci including CLEC16A, IL2RA, IL7R, RPL5, CD58, CD40 and chromosome 12q13–14 in addition to the well established allele HLA-DR15. There is potential that these genetic susceptibility factors could also modulate MS disease severity, as demonstrated previously for the MS risk allele HLA-DR15. We investigated this hypothesis in a cohort of 1006 well characterised MS patients from South-Eastern Australia. We tested the MS-associated SNPs for association with five measures of disease severity incorporating disability, age of onset, cognition and brain atrophy. We observed trends towards association between the RPL5 risk SNP and time between first demyelinating event and relapse, and between the CD40 risk SNP and symbol digit test score. No associations were significant after correction for multiple testing. We found no evidence for the hypothesis that these new MS disease risk-associated SNPs influence disease severity.

Proteolytic activation of Chlamydia trachomatis HTRA is mediated by PDZ1 domain interactions with protease domain loops L3 and LC and beta strand β5

Chlamydia trachomatis is a bacterial pathogen responsible for one of the most prevalent sexually transmitted infections worldwide. Its unique development cycle has limited our understanding of its pathogenic mechanisms. However, CtHtrA has recently been identified as a potential C. trachomatis virulence factor. CtHtrA is a tightly regulated quality control protein with a monomeric structural unit comprised of a chymotrypsin-like protease domain and two PDZ domains. Activation of proteolytic activity relies on the C-terminus of the substrate allosterically binding to the PDZ1 domain, which triggers subsequent conformational change and oligomerization of the protein into 24-mers enabling proteolysis. This activation is mediated by a cascade of precise structural arrangements, but the specific CtHtrA residues and structural elements required to facilitate activation are unknown. Using in vitro analysis guided by homology modeling, we show that the mutation of residues Arg362 and Arg224, predicted to disrupt the interaction between the CtHtrA PDZ1 domain and loop L3, and between loop L3 and loop LD, respectively, are critical for the activation of proteolytic activity. We also demonstrate that mutation to residues Arg299 and Lys160, predicted to disrupt PDZ1 domain interactions with protease loop LC and strand β5, are also able to influence proteolysis, implying their involvement in the CtHtrA mechanism of activation. This is the first investigation of protease loop LC and strand β5 with respect to their potential interactions with the PDZ1 domain. Given their high level of conservation in bacterial HtrA, these structural elements may be equally significant in the activation mechanism of DegP and other HtrA family members.

No association between MTHFR A1298C and MTRR A66G polymorphisms, and MS in an Australian cohort

Multiple sclerosis (MS) is a complex neurological disease that affects the central nervous system (CNS) resulting in debilitating neuropathology. Pathogenesis is primarily defined by CNS inflammation and demyelination of nerve axons. Methionine synthase reductase (MTRR) is an enzyme that catalyzes the remethylation of homocysteine (Hcy) to methionine via cobalamin and folate dependant reactions. Cobalamin acts as an intermediate methyl carrier between methylenetetrahydrofolate reductase (MTHFR) and Hcy. MTRR plays a critical role in maintaining cobalamin in an active form and is consequently an important determinant of total plasma Hcy (pHcy) concentrations. Elevated intracellular pHcy levels have been suggested to play a role in CNS dysfunction, neurodegenerative, and cerebrovascular diseases. Our investigation entailed the genotyping of a cohort of 140 cases and matched controls for MTRR and MTHFR, by restriction length polymorphism (RFLP) techniques. Two polymorphisms: MTRR A66G and MTHFR A1298C were investigated in an Australian age and gender matched case-control study. No significant allelic frequency difference was observed between cases and controls at the α = 0.05 level (MTRR χ2 = 0.005, P = 0.95, MTHFR χ2 = 1.15, P = 0.28). Our preliminary findings suggest no association between the MTRR A66G and MTHFR A1298C polymorphisms and MS

An examination of MS candidate genes identified as differentially regulated in multiple sclerosis plaque tissue, using absolute and comparative real-time Q-PCR analysis

In our laboratory, we have developed methods in real-time detection and quantitative-polymerase chain reaction (Q-PCR) to analyse the relative levels of gene expression in post mortem brain tissues. We have then applied this method to examine differences in gene activity between normal white matter (NWM) and plaque tissue from multiple sclerosis (MS) patients. Genes were selected based on their association with pathology and through identification by previously conducted global gene expression analysis. Plaque tissue was obtained from secondary progressive (SP) patients displaying chronic active, as well as acute pathologies; while NWM from the same location was obtained from age- and sex-matched controls (normal patients). In this study, we used both SYBR Green I supplementation and commercially available mixes to assess both comparative and absolute levels of gene activity. The results of both methods compared favourably for four of the five genes examined (P < 0.05, Pearsons), while differences in gene expression between chronic active and acute pathologies were also identified. For example, a >50-fold increase in osteopontin (Spp1) and inositol 1-4-5 phosphate 3 kinase B (Itpkb) levels in acute plaques contrasted with the 5-fold or less increase in chronic active plaques (P < 0.05, unpaired t test). By contrast, there was no significant difference in the levels of the MS marker and calcium-dependent protease (Calpain, Capns1) in MS plaque tissue. In summary, Q-PCR analysis using SYBR Green I has allowed us to economically obtain what may be clinically significant information from small amounts of the CNS, providing an opportunity for further clinical investigations.

MS-based metabolomics facilitates the discovery of in vivo functional small molecules with a diversity of biological contexts

In vivo small molecules as necessary intermediates are involved in numerous critical metabolic pathways and biological processes associated with many essential biological functions and events. There is growing evidence that MS-based metabolomics is emerging as a powerful tool to facilitate the discovery of functional small molecules that can better our understanding of development, infection, nutrition, disease, toxicity, drug therapeutics, gene modifications and host-pathogen interaction from metabolic perspectives. However, further progress must still be made in MS-based metabolomics because of the shortcomings in the current technologies and knowledge. This technique-driven review aims to explore the discovery of in vivo functional small molecules facilitated by MS-based metabolomics and to highlight the analytic capabilities and promising applications of this discovery strategy. Moreover, the biological significance of the discovery of in vivo functional small molecules with different biological contexts is also interrogated at a metabolic perspective.

Bioequivalence of two tablet formulations of capecitabine and exploration of age, gender, body surface area, and creatinine clearance as factors influencing systemic exposure in cancer patients

Purpose: The objective of the study was to assess the bioequivalence of two tablet formulations of capecitabine and to explore the effect of age, gender, body surface area and creatinine clearance on the systemic exposure to capecitabine and its metabolites. Methods: The study was designed as an open, randomized two-way crossover trial. A single oral dose of 2000 mg capecitabine was administered on two separate days to 25 patients with solid tumors. On one day, the patients received four 500-mg tablets of formulation B (test formulation) and on the other day, four 500-mg tablets of formulation A (reference formulation). The washout period between the two administrations was between 2 and 8 days. After each administration, serial blood and urine samples were collected for up to 12 and 24 h, respectively. Unchanged capecitabine and its metabolites were determined in plasma using LC/MS-MS and in urine by NMRS. Results: Based on the primary pharmacokinetic parameter, AUC(0-∞) of 5'-DFUR, equivalence was concluded for the two formulations, since the 90% confidence interval of the estimate of formulation B relative to formulation A of 97% to 107% was within the acceptance region 80% to 125%. There was no clinically significant difference between the t(max) for the two formulations (median 2.1 versus 2.0 h). The estimate for C(max) was 111% for formulation B compared to formulation A and the 90% confidence interval of 95% to 136% was within the reference region 70% to 143%. Overall, these results suggest no relevant difference between the two formulations regarding the extent to which 5'-DFUR reached the systemic circulation and the rate at which 5'-DFUR appeared in the systemic circulation. The overall urinary excretions were 86.0% and 86.5% of the dose, respectively, and the proportion recovered as each metabolite was similar for the two formulations. The majority of the dose was excreted as FBAL (61.5% and 60.3%), all other chemical species making a minor contribution. Univariate and multivariate regression analysis to explore the influence of age, gender, body surface area and creatinine clearance on the log-transformed pharmacokinetic parameters AUC(0-∞) and C(max) of capecitabine and its metabolites revealed no clinically significant effects. The only statistically significant results were obtained for AUC(0-∞) and C(max) of intact drug and for C(max) of FBAL, which were higher in females than in males. Conclusion: The bioavailability of 5'-DFUR in the systemic circulation was practically identical after administration of the two tablet formulations. Therefore, the two formulations can be regarded as bioequivalent. The variables investigated (age, gender, body surface area, and creatinine clearance) had no clinically significant effect on the pharmacokinetics of capecitabine or its metabolites.

Addressing accuracy and precision issues in iTRAQ quantitation

iTRAQ (isobaric tags for relative or absolute quantitation) is a mass spectrometry technology that allows quantitative comparison of protein abundance by measuring peak intensities of reporter ions released from iTRAQ-tagged peptides by fragmentation during MS/MS. However, current data analysis techniques for iTRAQ struggle to report reliable relative protein abundance estimates and suffer with problems of precision and accuracy. The precision of the data is affected by variance heterogeneity: low signal data have higher relative variability; however, low abundance peptides dominate data sets. Accuracy is compromised as ratios are compressed toward 1, leading to underestimation of the ratio. This study investigated both issues and proposed a methodology that combines the peptide measurements to give a robust protein estimate even when the data for the protein are sparse or at low intensity. Our data indicated that ratio compression arises from contamination during precursor ion selection, which occurs at a consistent proportion within an experiment and thus results in a linear relationship between expected and observed ratios. We proposed that a correction factor can be calculated from spiked proteins at known ratios. Then we demonstrated that variance heterogeneity is present in iTRAQ data sets irrespective of the analytical packages, LC-MS/MS instrumentation, and iTRAQ labeling kit (4-plex or 8-plex) used. We proposed using an additive-multiplicative error model for peak intensities in MS/MS quantitation and demonstrated that a variance-stabilizing normalization is able to address the error structure and stabilize the variance across the entire intensity range. The resulting uniform variance structure simplifies the downstream analysis. Heterogeneity of variance consistent with an additive-multiplicative model has been reported in other MS-based quantitation including fields outside of proteomics; consequently the variance-stabilizing normalization methodology has the potential to increase the capabilities of MS in quantitation across diverse areas of biology and chemistry.

Age-related trends in urinary excretion of bisphenol A in Australian children and adults : evidence from a pooled sample study using samples of convenience

Bisphenol A (BPA or 4,4’-(propane-2,2-diyl)diphenol) is a chemical intermediate in the production of polycarbonate and epoxy resins, and used in a wide range of applications. BPA has attracted significant attention in the past decade due to its frequency of detection in human populations worldwide, demonstrated animal toxicity and potential impact on human health, particularly during critical periods of development. The aim of this study was to perform a preliminary assessment of age-related trends in urinary concentration and to estimate daily excretion of BPA in Australian children (aged (>0 – <5 years) and adults (≥15 – <75 years). This was achieved using 79 samples pooled by age and gender, created from 868 individual samples of convenience collected as part of routine, community-based pathology testing. Total BPA was analyzed using online-SPE-LC-MS/MS and detected in all samples with a range of 0.65 – 265 ng/ml. No significant differences were observed between males and females. A urine flow model was constructed from published values and used to provide an estimate of daily excretion per unit bodyweight for each pooled sample. The daily excretion estimates ranged from 26.2 – 18200 ng/kg-d for children; and 20.1 – 165 ng/kg-d for adults. Urinary concentrations and estimated excretion rates were inversely associated with age, and estimated daily excretion rates in infants and young children were significantly higher than in adults (geometric mean: 107 and 47.0 ng/kg-d, respectively). Higher excretion of BPA in children may be explained by their higher food consumption relative to body weight compared to adults and adolescents, and may also reflect alternative exposure pathways and sources. Keywords: bisphenol A, biomonitoring, children, urine flow, Australia

A gene encoding a new cold-active lipase from an Antarctic isolate of Penicillium expansum

Cold-active lipases are of significant interest as biocatalysts in industrial processes. We have identified a lipase that displayed activity towards long carbon-chain-p-nitrophenyl substrates (C12–C18) at 25 °C from the culture supernatant of an Antarctic Penicillium expansum strain assigned P. expansum SM3. Zymography revealed a protein band of around 30 kDa with activity towards olive oil. DNA fragments of a lipase gene designated as lipPE were isolated from the genomic DNA of P. expansum SM3 by genomic walking PCR. Subsequently, the complete genomic lipPE gene was amplified using gene-specific primers designed from the 5′- and 3′-regions. Reverse transcription PCR was used to amplify the lipPE cDNA. The deduced amino acid sequence consisted of 285 residues that included a predicted signal peptide. Three peptides identified by LC/MS/MS analysis of the proteins in the culture supernatant of P. expansum were also present in the deduced amino acid sequence of the lipPE gene suggesting that this gene encoded the lipase identified by initial zymogram activity analysis. Full analysis of the nucleotide and the deduced amino acid sequences indicated that the lipPE gene encodes a novel P. expansum lipase. The lipPE gene was expressed in E. coli for further characterization of the enzyme with a view of assessing its suitability for industrial applications.

Transcript and metabolite profiling for the evaluation of tobacco tree and poplar as feedstock for the bio-based industry

The global demand for food, feed, energy and water poses extraordinary challenges for future generations. It is evident that robust platforms for the exploration of renewable resources are necessary to overcome these challenges. Within the multinational framework MultiBioPro we are developing biorefinery pipelines to maximize the use of plant biomass. More specifically, we use poplar and tobacco tree (Nicotiana glauca) as target crop species for improving saccharification, isoprenoid, long chain hydrocarbon contents, fiber quality, and suberin and lignin contents. The methods used to obtain these outputs include GC-MS, LC-MS and RNA sequencing platforms. The metabolite pipelines are well established tools to generate these types of data, but also have the limitations in that only well characterized metabolites can be used. The deep sequencing will allow us to include all transcripts present during the developmental stages of the tobacco tree leaf, but has to be mapped back to the sequence of Nicotiana tabacum. With these set-ups, we aim at a basic understanding for underlying processes and at establishing an industrial framework to exploit the outcomes. In a more long term perspective, we believe that data generated here will provide means for a sustainable biorefinery process using poplar and tobacco tree as raw material. To date the basal level of metabolites in the samples have been analyzed and the protocols utilized are provided in this article.