The systems biology simulation core library.

SUMMARY: Studying biological systems generally relies on computational modeling and simulation, e.g., model-driven discovery and hypothesis testing. Progress in standardization efforts led to the development of interrelated file formats to exchange and reuse models in systems biology, such as SBML, the Simulation Experiment Description Markup Language (SED-ML) or the Open Modeling EXchange format. Conducting simulation experiments based on these formats requires efficient and reusable implementations to make them accessible to the broader scientific community and to ensure the reproducibility of the results. The Systems Biology Simulation Core Library (SBSCL) provides interpreters and solvers for these standards as a versatile open-source API in JavaTM. The library simulates even complex bio-models and supports deterministic Ordinary Differential Equations; Stochastic Differential Equations; constraint-based analyses; recent SBML and SED-ML versions; exchange of results, and visualization of in silico experiments; open modeling exchange formats (COMBINE archives); hierarchically structured models; and compatibility with standard testing systems, including the Systems Biology Test Suite and published models from the BioModels and BiGG databases. AVAILABILITY AND IMPLEMENTATION: SBSCL is freely available at https://draeger-lab.github.io/SBSCL/ and via Maven Central. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Proteomic analysis of hepatic effects of phenobarbital in mice with humanized liver.

Activation of the constitutive androstane receptor (CAR) may induce adaptive but also adverse effects in rodent liver, including the induction of drug-metabolizing enzymes, transient hepatocellular proliferation, and promotion of liver tumor growth. Human relevance of CAR-related adverse hepatic effects is controversially debated. Here, we used the chimeric FRG-KO mouse model with livers largely repopulated by human hepatocytes, in order to study human hepatocytes and their response to treatment with the model CAR activator phenobarbital (PB) in vivo. Mice received an intraperitoneal injection with 50 mg/kg body weight PB or saline, and were sacrificed after 72-144 h. Non-repopulated FRG-KO mice were used as additional control. Comprehensive proteomics datasets were generated by merging data obtained by targeted as well as non-targeted proteomics approaches. For the first time, a novel proteomics workflow was established to comparatively analyze the effects of PB on human and murine proteins within one sample. Analysis of merged proteome data sets and bioinformatics data mining revealed comparable responses in murine and human hepatocytes with respect to nuclear receptor activation and induction of xenobiotic metabolism. By contrast, activation of MYC, a key regulator of proliferation, was predicted only for mouse but not human hepatocytes. Analyses of 5-bromo-2'-deoxyuridine incorporation confirmed this finding. In summary, this study for the first time presents a comprehensive proteomic analysis of CAR-dependent effects in human and mouse hepatocytes from humanized FRG-KO mice. The data support the hypothesis that PB does induce adaptive metabolic responses, but not hepatocellular proliferation in human hepatocytes in vivo.

Matrix and Sampling Effects on Quantification of Protein Biomarkers of Drug-Induced Liver Injury.

Macrophage colony stimulating factor 1 receptor (MCSF1R), osteopontin (OPN), high-mobility group protein B1 (HMGB1), glutamate dehydrogenase (GLDH), keratin 18 (K18), and caspase-cleaved keratin 18 (ccK18) are considered promising mechanistic biomarkers for the diagnosis of drug-induced liver injury. Here, we aim to elucidate the impact of the sample matrix and handling on the quantification of these emerging protein biomarkers. We investigated effects such as time from collection to centrifugation during serum (± gel) or EDTA plasma preparation on two assay platforms: immunoaffinity liquid chromatography mass spectrometric assays and sandwich immunoassays. Furthermore, we measured GLDH activity with an enzymatic activity assay. Matrix effects were observed particularly for HMGB1 and MCSF1R. HMGB1 levels were higher in serum than in plasma, whereas higher concentrations of MCSF1R were observed in plasma than in serum. A comparison of sample collection to centrifugation time ranging from 15 to 60 min demonstrated increasing levels of HMGB1 in serum, while MCSF1R, OPN, GLDH, and ccK18 concentrations remained stable. Additionally, there was a poor correlation in HMGB1 and ccK18 levels between serum and plasma. Considering the observed matrix effects, we recommend plasma as a matrix of choice and cross-study comparison studies to be limited to those using the same matrix.

Application of Mass Spectrometry-Based Immunoassays for the Species- and Tissue-Specific Quantification of Banned Processed Animal Proteins in Feeds.

Processed Animal Proteins (PAPs) are considered as a sustainable protein source to improve the nutritional profile of feed for livestock and aquaculture. However, the use of these proteins is strongly regulated since the bovine spongiform encephalopathy (BSE) crisis. The reintroduction of nonruminant PAPs for use in aquaculture in 2013 has driven the need for alternative analytical methods to determine the species origin as well as the tissue source (legal or not). The current official methods, light microscopy and polymerase chain reaction, do not fulfill these requirements. Furthermore, future methods need to be quantitative, because the pending zero-tolerance-concept is planned to be replaced by accurate thresholds. Here, we developed a 7-plex mass spectrometry-based immunoassay that is capable of quantifying 0.1% (w/w) ruminant PAP in feed in a tissue- and species-specific way. The workflow comprises a 2 h tryptic digestion of PAPs in suspension, an immunoaffinity enrichment of peptides, and LC-MS/MS-based quantification. In combination with a previously published assay for species identification, we were able to confirm the species and tissue origin of six ring trial samples obtained in former PCR and microscopy proficiency tests. The sensitive, quantitative, species- and tissue-specific character of the developed assays meets the requirements for new methods for PAP detection and can be used in future feed authentication studies.

Mass Spectrometry-Based Immunoassay for the Quantification of Banned Ruminant Processed Animal Proteins in Vegetal Feeds.

The ban of processed animal proteins (PAPs) in feed for farmed animals introduced in 2001 was one of the main EU measures to control the bovine spongiform encephalopathy (BSE) crisis. Currently, microscopy and polymerase chain reaction (PCR) are the official methods for the detection of illegal PAPs in feed. However, the progressive release of the feed ban, recently with the legalization of nonruminant PAPs for the use in aquaculture, requires the development of alternative methods to determine the species origin and the source (legal or not). Additionally, discussions about the need for quantitative tests came up, particularly if the zero-tolerance-concept is replaced by introducing PAP thresholds. To address this issue, we developed and partially validated a multiplex mass spectrometry-based immunoassay to quantify ruminant specific peptides in vegetal cattle feed. The workflow comprises a new sample preparation procedure based on a tryptic digestion of PAPs in suspension, a subsequent immunoaffinity enrichment of the released peptides, and a LC-MS/MS-based analysis for peptide quantification using isotope labeled standard peptides. For the very first time, a mass spectrometry-based method is capable of detecting and quantifying illegal PAPs in animal feed over a concentration range of 4 orders of magnitude with a detection limit in the range of 0.1% to 1% (w/w).

Direct Quantification of Cytochromes P450 and Drug Transporters-A Rapid, Targeted Mass Spectrometry-Based Immunoassay Panel for Tissues and Cell Culture Lysates.

The quantification of drug metabolizing enzymes and transporters has recently been revolutionized on the basis of targeted proteomic approaches. Isotope-labeled peptides are used as standards for the quantification of the corresponding proteins in enzymatically fragmented samples. However, hurdles in these approaches are low throughput and tedious sample prefractionation steps prior to mass spectrometry (MS) readout. We have developed an assay platform using sensitive and selective immunoprecipitation coupled with mass spectrometric readout allowing the quantification of proteins directly from whole cell lysates using less than 20,000 cells per analysis. Peptide group-specific antibodies (triple X proteomics antibodies) enable the enrichment of proteotypic peptides sharing a common terminus. These antibodies were employed to establish a MS-based immunoassay panel for the quantification of 14 cytochrome P450 (P450) enzymes and nine transporters. We analyzed the P450 enzyme and transporter levels in genotyped liver tissue homogenates and microsomes, and in samples from a time course induction experiment in human hepatocytes addressing different induction pathways. For the analysis of P450 enzymes and transporters only a minute amount of sample is required and no prefractionation is necessary, thus the assay platform bears the potential to bridge cell culture model experiments and results from whole organ tissue studies.

Catch and measure-mass spectrometry-based immunoassays in biomarker research.

Mass spectrometry-based (MS) methods are effective tools for discovering protein biomarker candidates that can differentiate between physiological and pathophysiological states. Promising candidates are validated in studies comprising large patient cohorts. Here, targeted protein analytics are used to increase sample throughput. Methods involving antibodies, such as sandwich immunoassays or Western blots, are commonly applied at this stage. Highly-specific and sensitive mass spectrometry-based immunoassays that have been established in recent years offer a suitable alternative to sandwich immunoassays for quantifying proteins. Mass Spectrometric ImmunoAssays (MSIA) and Stable Isotope Standards and Capture by Anti-Peptide Antibodies (SISCAPA/iMALDI) are two prominent types of MS-based immunoassays in which the capture is done either at the protein or the peptide level. We present an overview of these emerging types of immunoassays and discuss their suitability for the discovery and validation of protein biomarkers. This article is part of a Special Issue entitled: Biomarkers: A Proteomic Challenge.

Identification of short terminal motifs enriched by antibodies using peptide mass fingerprinting.

MOTIVATION: Mass spectrometry-based protein profiling has become a key technology in biomedical research and biomarker discovery. Sample preparation strategies that reduce the complexity of tryptic digests by immunoaffinity substantially increase throughput and sensitivity in proteomic mass spectrometry. The scarce availability of peptide-specific capture antibodies limits these approaches. Recently antibodies directed against short terminal motifs were found to enrich subsets of peptides with identical terminal sequences. This approach holds the promise of a significant gain in efficiency. TXP (Triple X Proteomics) and context-independent motif specific/global proteome survey binders are variants of this concept. Principally the binding motifs of such antibodies have to be elucidated after generating these antibodies. This entails a substantial effort in the lab, as it requires synthetic peptide libraries and numerous mass spectrometry experiments. RESULTS: We present an algorithm for predicting the antibody-binding motif in a mass spectrum obtained from a tryptic digest of a common cell line after immunoprecipitation. The epitope prediction, based on peptide mass fingerprinting, reveals the most enriched terminal epitopes. The tool provides a P-value for each potential epitope, estimated by sampling random spectra from a peptide database. The second algorithm combines the predicted sequences to more complex binding motifs. A comparison with library screenings shows that the predictions made by the novel methods are reliable and reproducible indicators of the binding properties of an antibody.

Combining ultracentrifugation and peptide termini group-specific immunoprecipitation for multiplex plasma protein analysis.

Blood plasma is a valuable source of potential biomarkers. However, its complexity and the huge dynamic concentration range of its constituents complicate its analysis. To tackle this problem, an immunoprecipitation strategy was employed using antibodies directed against short terminal epitope tags (triple X proteomics antibodies), which allow the enrichment of groups of signature peptides derived from trypsin-digested plasma. Isolated signature peptides are subsequently detected using MALDI-TOF/TOF mass spectrometry. Sensitivity of the immunoaffinity approach was, however, compromised by the presence of contaminant peaks derived from the peptides of nontargeted high abundant proteins. A closer analysis of the enrichment strategy revealed nonspecific peptide binding to the solid phase affinity matrix as the major source of the contaminating peptides. We therefore implemented a sucrose density gradient ultracentrifugation separation step into the procedure. This yielded a 99% depletion of contaminating peptides from a sucrose fraction containing 70% of the peptide-antibody complexes and enabled the detection of the previously undetected low abundance protein filamin-A. Assessment of this novel approach using 15 different triple X proteomics antibodies demonstrated a more consistent detection of a greater number of targeted peptides and a significant reduction in the intensity of nonspecific peptides. Ultracentrifugation coupled with immunoaffinity MS approaches presents a powerful tool for multiplexed plasma protein analysis without the requirement for demanding liquid chromatography separation techniques.

Immunoaffinity proteomics for kidney injury biomarkers in male beagle dogs.

Drug-induced kidney injury (DIKI) is a cause of drug development failure. Dogs represent a common non-rodent animal model in pre-clinical safety studies; however, biomarker assays for detecting nephrotoxicity in dogs are limited. To identify novel proteins and gain insight into the molecular mechanisms involved in DIKI, we developed an assay to evaluate proteomic changes associated with DIKI in male beagle dogs that received nephrotoxic doses of tobramycin for 10 consecutive days. Label-free quantitative discovery proteomics analysis on representative kidney cortex tissues collected on Day 11 showed that the tobramycin-induced kidney injury led to a significant differential regulation of 94 proteins mostly associated with mechanisms of nephrotoxicity such as oxidative stress and proteasome degradation. For verification of the proteomic results, we developed a multiplex peptide-centric immunoaffinity liquid chromatography tandem mass spectrometry assay (IA LC-MS/MS) to evaluate the association of eight DIKI protein biomarker candidates using kidney cortices collected on Day 11 and urine samples collected on Days -4, 1, 3, 7 and 10. The results showed that most biomarkers evaluated were detected in the kidney cortices and their expression profile in tissue aligned with the label-free data. Cystatin C was the most consistent marker regardless of the magnitude of the renal injury while fatty acid-binding protein-4 (FABP4) and kidney injury molecule-1 (KIM-1) were the most affected biomarkers in response to moderate proximal tubular injury in absence of changes in serum-based concentrations of blood urea nitrogen or creatinine. In the urine, clusterin is considered the most consistent biomarker regardless of the magnitude and time of the renal injury. To our knowledge, this is the most comprehensive multiplex assay for the quantitative analysis of mechanism-based proximal tubular injury biomarkers in dogs.

From spots to beads-PTM-peptide bead arrays for the characterization of anti-histone antibodies.

Antibodies that recognize PTMs of histones play a central role in epigenetic proteomic research. Modification-specific antibodies are employed in chromatin immunoprecipitation, for Western blotting and during the immunoprecipitation steps for MS-based global proteomic analyses. Knowledge about the antibodies' off-target binding is essential for the interpretation of experimental data. To address this challenge we developed a fast and cost efficient system for generating peptide bead arrays. We employed this method to establish a bead-based peptide array containing 384 peptides displaying phosphorylated, acetylated, methylated, and citrullinated N-terminal regions of histones H2A, H2B, H3 and H4 and controls. We profiled the binding of 40 PTM-specific antibodies important for epigenetic proteomic research.

Targeting peptide termini, a novel immunoaffinity approach to reduce complexity in mass spectrometric protein identification.

Mass spectrometry and peptide-centric approaches are powerful techniques for the identification of differentially expressed proteins. Despite enormous improvements in MS technologies, sample preparation and efficient fractionation of target analytes are still major bottlenecks in MS-based protein analysis. The complexity of tryptically digested whole proteomes needs to be considerably reduced before low abundance proteins can be effectively analyzed using MS/MS. Sample preparation strategies that use peptide-specific antibodies are able to reduce the complexity of tryptic digests and lead to a substantial increase in throughput and sensitivity; however, the number of peptide-specific capture reagents is low, and consequently immunoaffinity-based approaches are only capable of detecting small sets of protein-derived peptides. In this proof-of-principle study, special anti-peptide antibodies were used to enrich peptides from a complex mixture. These antibodies recognize short amino acid sequences that are found directly at the termini of the peptides. The recognized epitopes consist of three or four amino acids only and include the terminally charged group of the peptide. Because of its limited length, antibodies recognizing the epitope will enrich not only one peptide but a whole class of peptides that share this terminal epitope. In this study, ß-catenin-derived peptides were used to demonstrate that it is possible (i) to effectively generate antibodies that recognize short C-terminal peptide epitopes and (ii) to enrich and identify peptide classes from a complex mixture using these antibodies in an immunoaffinity MS approach. The expected ß-catenin peptides and a set of 38 epitope-containing peptides were identified from trypsin-digested cell lysates. This might be a first step in the development of proteomics applications that are based on the use of peptide class-specific antibodies.

Pesticide mixture effects on liver protein abundance in HepaRG cells.

Toxicological effects of chemicals are mostly tested individually. However, consumers encounter exposure to complex mixtures, for example multiple pesticide residues, by consuming food such as crops, fruits or vegetables. Currently, more than 450 active substances are approved in the European Union, and there is little data on effects after combined exposure to several pesticides. Toxicological animal studies would increase enormously, if pesticide combinations had to be analyzed in vivo. Therefore, in vitro methods addressing this issue are needed. We have developed 32 immunoaffinity-based mass spectrometry assays to investigate the impact of hepatotoxic active substances on liver proteins in human HepaRG cells. Five compounds were selected based on their (dis)similar capability to modulate protein levels, and on their combined use in commercially available formulations. Four binary mixtures were prepared from these five substances and tested in different concentrations over three time points. We applied a novel statistical method to describe deviations from additivity and to detect antagonistic and synergistic effects. The results regarding the abundance of hepatotoxicity-related proteins showed additive behavior for 1323 out of 1427 endpoints tested, while 104 combinatorial effects deviating from additivity, such as antagonism or synergism were observed.

A targeted transcriptomics approach for the determination of mixture effects of pesticides.

While real-life exposure occurs to complex chemical mixtures, toxicological risk assessment mostly focuses on individual compounds. There is an increasing demand for in vitro tools and strategies for mixture toxicity analysis. Based on a previously established set of hepatotoxicity marker genes, we analyzed mixture effects of non-cytotoxic concentrations of different pesticides in exposure-relevant binary mixtures in human HepaRG hepatocarcinoma cells using targeted transcriptomics. An approach for mixture analysis at the level of a complex endpoint such as a transcript pattern is presented, including mixture design based on relative transcriptomic potencies and similarities. From a mechanistic point of view, goal of the study was to evaluate combinations of chemicals with varying degrees of similarity in order to determine whether differences in mechanisms of action lead to different mixtures effects. Using a model deviation ratio-based approach for assessing mixture effects, it was revealed that most data points are consistent with the assumption of dose addition. A tendency for synergistic effects was only observed at high concentrations of some combinations of the test compounds azoxystrobin, cyproconazole, difenoconazole, propiconazole and thiacloprid, which may not be representative of human real-life exposure. In summary, the findings of our study suggest that, for the pesticide mixtures investigated, risk assessment based on the general assumption of dose addition can be considered sufficiently protective for consumers. The way of data analysis presented in this paper can pave the way for a more comprehensive use of multi-gene expression data in experimental studies related to mixture toxicity.

Immune response to SARS-CoV-2 variants of concern in vaccinated individuals.

SARS-CoV-2 is evolving with mutations in the receptor binding domain (RBD) being of particular concern. It is important to know how much cross-protection is offered between strains following vaccination or infection. Here, we obtain serum and saliva samples from groups of vaccinated (Pfizer BNT-162b2), infected and uninfected individuals and characterize the antibody response to RBD mutant strains. Vaccinated individuals have a robust humoral response after the second dose and have high IgG antibody titers in the saliva. Antibody responses however show considerable differences in binding to RBD mutants of emerging variants of concern and substantial reduction in RBD binding and neutralization is observed against a patient-isolated South African variant. Taken together our data reinforce the importance of the second dose of Pfizer BNT-162b2 to acquire high levels of neutralizing antibodies and high antibody titers in saliva suggest that vaccinated individuals may have reduced transmission potential. Substantially reduced neutralization for the South African variant further highlights the importance of surveillance strategies to detect new variants and targeting these in future vaccines.

SBML2L(A)T(E)X: conversion of SBML files into human-readable reports.

SUMMARY: The XML-based Systems Biology Markup Language (SBML) has emerged as a standard for storage, communication and interchange of models in systems biology. As a machine-readable format XML is difficult for humans to read and understand. Many tools are available that visualize the reaction pathways stored in SBML files, but many components, e.g. unit declarations, complex kinetic equations or links to MIRIAM resources, are often not made visible in these diagrams. For a broader understanding of the models, support in scientific writing and error detection, a human-readable report of the complete model is needed. We present SBML2L(A)T(E)X, a Java-based stand-alone program to fill this gap. A convenient web service allows users to directly convert SBML to various formats, including DVI, L(A)T(E)X and PDF, and provides many settings for customization. AVAILABILITY: Source code, documentation and a web service are freely available at (http://www.ra.cs.uni-tuebingen.de/software/SBML2LaTeX).

Induction and repression effects on CYP and transporter protein abundance by azole mixture uptake in rat liver.

Detection of mixture effects is a major challenge in current experimental and regulatory toxicology. Robust markers are needed that are easy to quantify and responsive to chemical stressors in a broad dose range. Several hepatic enzymes and proteins related to drug metabolism like cytochrome-P-450 (CYP) enzymes and transporters have been shown to be responsive to pesticide active substances in a broad dose range and are therefore good candidates to be used as markers for mixture toxicity. Even though they can be well quantified at the mRNA level, quantification on the protein level is challenging because most of these proteins are membrane bound. Here we report the development of mass spectrometry-based assays using triple-x-proteomics (TXP) antibodies in combination with targeted selected ion monitoring (tSIM) to quantify changes of protein levels due to exposure to mixtures of pesticide active substances. Our results indicate that changes on the protein level of CYP1A1, ABCB2, ABCC3 are in line with observations on the mRNA and enzyme activity level and are indicative of mixture effects. Therefore, the tests are promising to reveal effects by chemical mixture effects in toxicological studies in rats.

RNA-protein correlation of liver toxicity markers in HepaRG cells.

The liver is a main target organ for the toxicity of many different compounds. While in general, in vivo testing is still routinely used for assessing the hepatotoxic potential of test chemicals, the use of in vitro models offers advantages with regard to throughput, consumption of resources, and animal welfare aspects. Using the human hepatoma cell line HepaRG, we performed a comparative evaluation of a panel of hepatotoxicity marker mRNAs and proteins after exposure of the cells to 30 different pesticidal active compounds comprising herbizides, fungicides, insecticides, and others. The panel of hepatotoxicity markers included nuclear receptor target genes, key players of fatty acid and bile acid metabolism-related pathways, as well as recently identified biomarkers of drug-induced liver injury. Moreover, marker genes and proteins were identified, for example, S100P, ANXA10, CYP1A1, and CYP7A1. These markers respond with high sensitivity to stimulation with chemically diverse test compounds already at non-cytotoxic concentrations. The potency of the test compounds, determined as an overall parameter of their ability to deregulate marker expression in vitro, was very similar between the mRNA and protein levels. Thus, this study does not only characterize the response of human liver cells to 30 different pesticides but also demonstrates that hepatotoxicity testing in human HepaRG cells yields well comparable results at the mRNA and protein levels. Furthermore, robust hepatotoxicity marker genes and proteins were identified in HepaRG cells.

Species Differentiation and Quantification of Processed Animal Proteins and Blood Products in Fish Feed Using an 8-Plex Mass Spectrometry-Based Immunoassay.

With the reintroduction of nonruminant processed animal proteins (PAPs) for use in aquaculture in 2013, there is a suitable alternative to replace expensive fish meal in fish feed. Nevertheless, since the bovine spongiform encephalopathy (BSE) crisis, the use of PAPs in feed is strictly regulated. To date, light microscopy and polymerase chain reaction are the official methods for proving the absence of illegal PAPs in feed. Due to their limitations, alternative methods for the quantitative species differentiation are needed. To address this issue, we developed and validated an 8-plex mass spectrometry-based immunoassay. The workflow comprises a tryptic digestion of PAPs and blood products in suspension, a cross-species immunoaffinity enrichment of 8 species-specific alpha-2-macroglobulin peptides using a group-specific antibody, and a subsequent analysis by ultrahigh-performance liquid chromatography coupled to tandem mass spectrometry for species identification and quantification. This workflow can be used to quantitatively determine the species origin in future feed authentication studies.

Performance of a Multiplex Serological Helicobacter pylori Assay on a Novel Microfluidic Assay Platform.

Infection with Helicobacter pylori (H. pylori) occurs in 50% of the world population, and is associated with the development of ulcer and gastric cancer. Serological diagnostic tests indicate an H. pylori infection by detecting antibodies directed against H. pylori proteins. In addition to line blots, multiplex assay platforms provide smart solutions for the simultaneous analysis of antibody responses towards several H. pylori proteins. We used seven H. pylori proteins (FliD, gGT, GroEL, HpaA, CagA, VacA, and HP0231) and an H. pylori lysate for the development of a multiplex serological assay on a novel microfluidic platform. The reaction limited binding regime in the microfluidic channels allows for a short incubation time of 35 min. The developed assay showed very high sensitivity (99%) and specificity (100%). Besides sensitivity and specificity, the technical validation (intra-assay CV = 3.7 ± 1.2% and inter-assay CV = 5.5 ± 1.2%) demonstrates that our assay is also a robust tool for the analysis of the H. pylori-specific antibody response. The integration of the virulence factors CagA and VacA allow for the assessment of the risk for gastric cancer development. The short assay time and the performance of the platform shows the potential for implementation of such assays in a clinical setting.