Proteomic Candidate Biomarkers of Drug-Induced Nephrotoxicity in the Rat

Improved biomarkers of acute nephrotoxicity are coveted by the drug development industry, regulatory agencies, and clinicians. In an effort to identify such biomarkers, urinary peptide profiles of rats treated with two different nephrotoxins were investigated. 493 marker candidates were defined that showed a significant response to cis-platin comparing a cis-platin treated cohort to controls. Next, urine samples from rats that received three consecutive daily doses of 150 or 300 mg/kg gentamicin were examined. 557 potential biomarkers were initially identified; 108 of these gentamicin-response markers showed a clear temporal response to treatment. 39 of the cisplatin-response markers also displayed a clear response to gentamicin. Of the combined 147 peptides, 101 were similarly regulated by gentamicin or cis-platin and 54 could be identified by tandem mass spectrometry. Most were collagen type I and type III fragments up-regulated in response to gentamicin treatment. Based on these peptides, classification models were generated and validated in a longitudinal study. In agreement with histopathology, the observed changes in classification scores were transient, initiated after the first dose, and generally persistent over a period of 10–20 days before returning to control levels. The data support the hypothesis that gentamicin-induced renal toxicity up-regulates protease activity, resulting in an increase in several specific urinary collagen fragments. Urinary proteomic biomarkers identified here, especially those common to both nephrotoxins, may serve as a valuable tool to investigate potential new drug candidates for the risk of nephrotoxicity.


Introduction
The drug development industry, government regulatory agencies, and healthcare professionals are all major stakeholders in the development of biomarkers of drug-induced injury. The advancement of organ-specific biomarkers of drug-induced injury promises each stakeholder improved efficiency and effectiveness in the drug development process that ultimately results in safe and efficacious products coming to the market. The United States Food and Drug Administration (U. S. FDA) has gone as far as describing an official biomarker qualification process [1] to hasten adoption of candidate biomarkers. Kidney injury biomarkers were the logical test case for the biomarker qualification process. Classical, functional biomarkers of kidney injury (blood urea nitrogen, serum creatinine), while accessible (serum, plasma), are not sensitive or specific to etiology or location of injury, leaving considerable room for biomarker improvement. Individual biomarkers have been evaluated and qualified by the U. S. FDA based on validated immune-based assays and data packages put together by industrial consortia [2]. Ongoing evaluation of these biomarkers suggests that specific insults may illicit different biomarker responses and that building biomarker profiles might be the ultimate tool for identifying injury [3][4][5]. As a consequence, there is a specific need for additional biomarkers that enable the generation of such specific biomarker profiles. Ideally, the biomarkers should be easily accessible in a non-invasive way, and should be applicable in animal models, as well as in man.
''Omic'' technologies (genomics, metabolomics, proteomics, etc.) hold the promise to fulfil this need and enable identification of multiple biomarkers that reflect specific types of injury in the kidney. Several proteomics approaches have been described in this context [6][7][8][9]. CE-MS methodology was validated as an analytical tool for the measurement of peptides in rat urine and subsequently used to profile the urinary low-molecular proteome of the rat [10].
In an earlier publication [11], CE-MS was used as a biomarker discovery tool for nephrotoxicity in rats treated with cis-platin. In the study reported here we aimed to identify common and disparate biomarkers of cis-platin-and gentamicin-induced nephrotoxicity by applying CE-MS proteomics in rat urine. The aim of the study was to detect multiple biomarkers that can be efficiently analysed in a non-invasive approach. Such biomarkers could form the basis for specific multi-marker models for displaying drug-induced kidney injury in pre-clinical and clinical application and may have substantial translational value.

Identification of the rat urinary peptides indicative of cisplatin and gentamicin induced nephrotoxicity
Two parallel approaches were employed to identify potential biomarkers for drug-induced nephrotoxicity. The study design is graphically depicted in figure 1. In the first approach, CE-MS data generated in a previous study [11] were used as an initial reference to identify biomarkers for cis-platin-induced nephrotoxicity (table S1 sheet 1). The comparison of the urinary peptides and proteins from 14 controls and 25 treated animals resulted in identification of 493 peptides that showed statistically significant changes in distribution (P,0.05 after correction for multiple testing).
In the second approach, rat urine samples from a gentamicin study [12] were used as a second source for nephrotoxicity biomarker discovery. Samples were analyzed using CE-MS and the data were matched against a previously established rat urinary proteome database [10,11]. The compiled proteomic profiles of all samples in this study are shown according to dose and treatment duration in figure S1. Urine samples collected at days 3, 7 or 10 (the first samples taken after treatment) from rats treated once daily for three consecutive days with 150 mg/kg or 300 mg/kg gentamicin (n = 23) or saline (from all time points n = 40) were analyzed as cases and controls, respectively, to form the discovery set (table S1 sheet 2). A total of 557 potential biomarker candidates could be defined (P,0.05). Of this peptide pool, 88 peptides intersected with the 493 peptides profiled in the cis-platin study (figure 1).
Next, proteomic data from untreated animals and animals treated with 150 or 300 mg/kg gentamicin were evaluated and a mean amplitude was calculated for each peptide at each time point (1,2,3,7,10,15,18,22,29,36 and 44 days) for both doses. The distributions of all peptides (493+557) were visually inspected on a graph with amplitude on the y-axis and treatment-time on the xaxis for controls (0 mg), 150 and 300 mg/kg gentamicin to determine regulation over time. This analysis revealed that 39 of the 493 biomarker candidates responsive to cis-platin also displayed a clear response to gentamicin (figures S2). Of the 557 potential candidate peptides from the gentamicin study, an additional 108 peptides showed a clear gentamicin response over the whole observation period (see figures S2). Combining the markers of both analyses provided a list of 147 marker candidates. Some of these peptides showed opposite regulation in the two studies. Only peptides with a change in the same direction were selected as drug induced nephrotoxicity markers (table S2). The mean signal intensities of the resultant 101 potential biomarkers are shown in figure 2.
To obtain sequence information for the nephrotoxicity marker candidates, tandem mass spectrometry was applied. We were able to identify 54 of the 101 marker candidates. Table 1 gives sequence information for the 54 biomarker candidates. The majority of the identified biomarkers were fragments of collagen alpha-1 (I), (III) and alpha-2 (I). Fragments of apolipoprotein A-IV, fibrinogen beta, fetuin-A, actin, hemoglobin subunit beta, inter alpha-trypsin inhibitor, osteopontin, pro-epidermal growth factor, prothrombin, tropomyosin-1 alpha and contrapsin-like protease inhibitors 3 and 6 were also identified.

Evaluation of the distribution of protein fragments of previously defined biomarkers
The U.S. FDA and European Medicines Agency (EMEA) have qualified beta-2-microglobulin, cystatin C, clusterin, kim-1, trefoil factor-3, albumin, total protein, and rpa-1 as biomarkers of acute drug-induced kidney injury in the rat [13]. Recently, researchers there described the association of clusterin, kim-1, albumin, rpa-1, osteopontin, lipocalin-2 (NGAL), alpha GST and mu GST to gentamicin treatment [12,14]. Our rat urinary proteome database which contains naturally occurring peptides (defined by exact mass, migration time, and, if available, exact sequence), but generally not whole proteins (as these are generally not observed in urine), was searched for peptides derived from these biomarkers. We could identify 5 fragments of clusterin, 16 fragments of osteopontin and 11 fragments of albumin in the rat urine. To examine if any of these fragments are altered after treatment with gentamicin the distribution of these 32 peptides in samples of untreated rats (n = 40) were compared to samples of rats treated with 150 and 300 mg/kg gentamicin collected at days 3, 7 and 10 using correction for multiple testing. One clusterin, two albumin, and five osteopontin fragments were significantly altered (P,0.05, see table 2) and showed a direction of regulation identical to that of the total protein.
Generation of high-dimensional biomarker models for drug induced nephrotoxicity in the rat As drug-induced nephrotoxicity represents a highly complex injury, we aimed to capture this complex pathophysiology in a multidimensional biomarker model. To this end, we selected two different sets of biomarker candidates for modeling purposes. First, we used all selected 101 candidate biomarkers for nephrotoxicity as listed in table S2. Second, we used a biomarker set composed of those 54 candidate biomarkers that were sequenced (table 1).
Linear-and support-vector-machine (SVM)-based models for the 101 and 54 marker combinations were generated using the 39 samples from the cis-platin study as a training set. The classification of the cis-platin (training) data using complete takeone-out cross-validation resulted in an AUC of 1.00 with all SVM models. Classification with the linear model based on 54 and 101 markers yielded AUC's of 0.87 and 0.95, respectively. Subsequently, samples from the gentamicin study were classified with all generated models. Classification factors for all samples from the gentamicin study are listed in table S3 and are shown in figure 3 A for the linear models and figure 3 B for the SVM models. The classification factors represent a composite index of signal intensities of all biomarkers included in the model, calculated using linear-or SVM-based algorithms. All models showed a similar response to gentamicin treatment over time. With the exception of the SVM model with 54 peptides, the models showed higher scores immediately after the first drug dose for treated animals in comparison to untreated controls. The highest scores were observed between days 3 and 10, shortly thereafter the scores returned to the level of controls. Classification with the biomarker models correlated with the administered drug dose; such that higher doses resulted in higher classification scores.
In an effort to estimate the value and validity of the urinary peptide biomarkers, as well as the biomarker models developed here, we examined their correlation with the pathophysiological changes observed and with other kidney injury biomarkers. As evident from the data presented in table 3, significant correlation could be established for most variables.

Discussion
Previously, urine obtained from Sprague Dawley rats before and after administration of cis-platin was analyzed to identify biomarkers of drug-induced nephrotoxicity. In a blinded sample set, a set of 34 urinary peptides was validated that demonstrated significant differences between treated and untreated animals [11]. In the present study, all data obtained previously (discovery and validation set) were used as a discovery set to define additional cisplatin-induced nephrotoxicity biomarker candidates. In addition, a set of marker candidates was defined in the urine of gentamicin treated rats. While an overlap of injury markers existed with the two nephrotoxins, it was quite modest. Mitochondrial injury leading to apoptotic and necrotic cell death is common to both nephrotoxins, however, initiating events are thought to be very different [15][16][17]. Furthermore, gentamicin is recognized to have broader nephrotoxic effects impacting the collecting ducts and glomeruli [15] while cis-platin injury is more localized to the proximal and distal tubules especially that segment (S3) of the proximal tubule located in the outer medulla and the corticomedullary junction [18]. Data from the present study further support that substantial differences exist on a molecular level between cis-platin-and gentamicin-induced injury and suggest that the toxic mechanisms of these two nephrotoxins vary considerably and/or target different structures within the nephron. Depicted is the size of the sample cohorts used in the biomarker definition step and the number of selected peptides after statistical comparison of the respective treated and control groups. Numbers of peptides identified in both cis-platin and genatmicin treated rats as markers are denominated in the intersection of the diagrams. The significant marker candidates were visually inspected and only those candidates were selected that showed a clear response to gentamicin. In the next step the direction of regulation was controlled for the markers identified in both studies and selected were only those demonstrating similar (up or down) regulation with both cis-platin and gentamicin. In the last step the number of markers for which the amino acid sequence could be resolved is listed. doi:10.1371/journal.pone.0034606.g001 Individually, the candidate biomarkers showed a clear response to the gentamicin treatment. Biomarker models representing a composite index of all marker intensities based on both cis-platin and gentamicin specific biomarkers were generated. These models performed well in both the cis-platin and gentamicin treated cohorts. Longitudinal data from the gentamicin cohort revealed a dose-dependent response immediately after the first administration of gentamicin. The highest intensity response was observed between days 3 and 10. These results demonstrated that the gentamicin-induced changes in urinary peptides are rapid but return to control levels concurrent with injury resolution. No significant differences between the control group and the gentamicin-treated animals were observed 4 weeks after injury. Consistent with these peptide data, all the urinary protein biomarkers elevated in the original study and the associated histopathology changes [12] demonstrated the same temporal response. The correlation with the different variables was not uniform (e.g. highest correlation with the different histopathological datasets was not consistently observed with one specific biomarker), but on average the highest correlation was detected with the biomarker model based on the 101 biomarkers, again supporting the concept that a multi-marker model is better suited to display complex pathophysiological changes. Furthermore, when examining the data in a hypothesis-driven approach and limiting the dataspace to only previously described biomarkers [12,[19][20][21][22][23][24], we could detect fragments of peptides from three of the eight specific acute kidney injury (AKI) -associated biomarkers (albumin, osteopontin, and clusterin) from the original study.
Using the untargeted MS it was not possible to detect fragments of all eight AKI-associated biomarkers that were previously described. This may be because the protein biomarkers are not represented by specific peptides (e.g. may not fragment in the same areas consistently or may degrade to very small fragments) and hence cannot be detected by an approach targeted at naturally occurring peptides. In addition, we were not able to identify the sequence of all potential peptide biomarkers in this study. This is due to the specific challenges associated with sequencing of naturally occurring peptides. These challenges are described in detail elsewhere [25]. The most prominent hurdles are poor fragmentation and/or post-translational modifications altering the mass, and thereby interfering with sequence assignment. It should also be noted that the sensitivity of untargeted MS is generally below that of targeted immunological methods. Therefore, low abundant prototypic peptides that may represent the biomarkers may be present but below the limit of detection with the MS technology applied in this study.
Altered regulation of human urinary peptides in AKI and tubular injury was previously described [26,27]. All 54 identified sequences in this study were compared with sequenced human urinary biomarkers for AKI [26] and for Fanconi syndrome [27]. Regulation of collagen fragments in human AKI was opposite to that seen in the rat models described here. Generally, urinary collagen fragments were substantially reduced in human AKI while a significant up-regulation was observed in the rat. In addition to inherent species differences, disease etiology and severity could account for some response differences. In the  gentamicin rodent study, experimental kidney injury due to a toxic insult was very apparent upon histopathology examination and was detectable with urinary protein biomarkers. However, the classical serum biomarkers, blood urea nitrogen and serum creatinine, were only minimally elevated and animals never demonstrated morbidity or clinical signs of kidney injury. In the human study, kidney injury was generally observed in conjunction with morbidity and/or more severe alteration of classical biomarkers. Perhaps more significantly, the human study encompassed AKI of largely pre-renal or glomerular etiology while the animal models represented drug-induced, primarily kidney tubular injury. This may explain the difference in regulation of urinary collagen fragments. Reduction of glomerular function has been associated with a reduction of specific collagen fragments in the urine [28,29]. Tubular damage resulting in reduced tubular reabsorption and consequently in an increases of the urinary protein and peptides secretion has been investigated in proteomic experiments [27]. A comparison of rat nephrotoxicity markers and human markers for tubular injury (Fanconi syndrome) [27] reveals similarities including the elevation of urinary collagen fragments. Interestingly, two collagen fragments were observed to have identical cleavage sites in the human and in the rat, indicating upregulation of similar proteases. These data give rise to the hypothesis that the observed alterations in specific urinary peptides may be indicative of disease-associated changes in extracellular remodeling, displayed in the urine by increase in the specific urinary peptides reported here.
The majority of identified biomarkers present in this study were fragments of the collagen chains alpha-1 (I), alpha-1 (III) and alpha-2 (I). Several of the collagen fragments had a PGP-motif at the C-terminus, suggesting these fragments may be generated by matrix metalloproteinase (MMP) activity. Up-regulated MMPs, especially MMP-2 and MMP-9, have been detected following acute kidney injury, mostly in animal models of ischemia [30,31], and up-regulation of MMP-9 was reported to protect from apoptosis in AKI [32]. Deregulation of MMPs, up-regulation of MMP-2 and down-regulation of MMP-9 have also been reported in CKD [33]. In addition, MMP-9 deletion was reported to mitigate vascular lesions (hence insult) after ischemia [34,35]. It appears reasonable that increases in urinary collagen fragments may evolve from increased MMP (and possibly other protease) activity as a result of a toxic insult. Significant up-regulation of peptides derived from the N-terminus of Fibrinogen-beta was also recognized. Recently, a peptide of similar origin was ascribed nephroprotective properties [36]. The nephroprotective capability of the Fibrinogen-beta fragment described here needs to be evaluated. For a fibrinogen alpha chain derived peptide identified in this study as a rat nephrotoxicity marker, a peptide homologue in humans exists which was previously described as prognostic for AKI development [26]. This study also identified as a potential biomarker candidate a fragment of apolipoprotein A-IV that has been described as predictive for progression of chronic kidney diseases [37].
In conclusion, we have identified a panel of urinary peptide biomarkers that are significantly associated with drug-induced Given are molecular mass (in Da), normalized migration time (in min) regulation factor (mean signal intensity of urine samples from 150 and 300 mg/kg gentamicin treated rats collected at days 3, 7 and 10 divided by mean signal intensity of control samples), adjusted p-value (Benjamini and Hochberg), amino acid sequence (modified amino acids: p = hydroxyproline; k = hydroxylysine; m = oxidized methionine) and short protein name of the UniProt database with the position of the firs and last amino acid in parenthesis. doi:10.1371/journal.pone.0034606.t001 nephrotoxicity in two different rat models. The biomarkers correlate with pathophysiology and likely reflect collagen degradation and changes in extracellular matrix turnover associated with increased MMP activity. These biomarkers, and especially the high-dimensional biomarker models, appear to be valuable for the monitoring of early nephrotoxicity in drug safety trials. To further establish their value and validity, we will aim at analyzing their performance in additional studies, and also investigate their distribution in rat models for kidney disease like the ZDF model.

Specimen characteristics
Urine samples from rats in a previously described cis-platin study [11] were analyzed. Rats were given a single i.p. injection of cis-platin (3 or 6 mg/kg). Urine samples were collected at 0, 24, 48 and 72 hrs after dosing. Details on this sample cohort are given in table S1 sheet 1. In addition, urine samples were obtained from rats following gentamicin administration [12]. Saline treated controls were included with animals receiving intramuscular injection of 150 or 300 mg/kg gentamicin once a day for up to 3 consecutive days. Urine samples were collected 24 hours following a single dose (day 1), two consecutive daily doses (day 2), or three consecutive daily doses (day 3). Subsequently, urine samples were collected following a recovery time (no additional treatments) at 7,10,15,18,22,29,36 or 44 days following the first of the three consecutive daily gentamicin doses. For each dose group, 3-4 samples were analyzed at each time point. Three control animals from which no samples were obtained were tested as sentinels. One control and one treated animal were removed from the study due to non-treatment related complications. Table S1 sheet 2 includes all samples, doses and collection days for the gentamicin study.

CE-MS analysis
150 ml urine were mixed with 150 ml of 2 M urea, 100 mM NaCl, 10 mM NH 4 OH containing 0.02% SDS. Samples were ultrafiltered using a Centristat 20 kDa cut-off centrifugal filter (Satorius, Göttingen, Germany) to eliminate high molecular weight compounds. The filtrate was desalted using a NAP-5 gel filtration column (GE Healthcare Bio Sciences, Uppsala, Sweden) to remove urea and electrolytes. The sample was lyophilized in a Christ Speed-Vac RVC 2-18/Alpha 1-2 (Christ, Germany) and  stored at 4uC until use. Shortly before CE-MS analysis, the samples were re-suspended in 10 mL HPLC grade H 2 O. CE-MS analysis was performed using a P/ACE MDQ capillary electrophoresis system (Beckman Coulter, Fullerton, USA) on-line coupled to a MicrOTOF MS (Bruker Daltonic, Bremen, Germany). The ESI-sprayer (Agilent Technologies, Palo Alto, CA, USA) was grounded, and the ion spray interface potential was set 24.5 kV. Data acquisition and MS acquisition methods were automatically controlled by the CE via contact-close-relays. Spectra were accumulated every 3 s, over a range of m/z 350 to 3000. The analytical characteristics of the CE-MS system were extensively investigated by Theodorescu et al. [38] and Kolch et al. [39].

Data processing
Mass spectral ion peaks representing identical molecules at different charge states were deconvoluted into single masses using MosaiquesVisu software [40]. We defined ''rat urinary housekeeping polypeptides'' and calibrated the CE-MS data utilizing 545 migration time data points, 108 mass data points by applying local and global linear regression, respectively. References of 37 highly abundant peptides were used as ''internal standard peptides'' for ion signal intensity (amplitude) calibration using global linear regression. The procedure to use ''internal standard'' for amplitude normalization, was shown to be a reliable method to address both analytical and dilution variances [41]. The resulting peak list characterizes each polypeptide by its molecular mass [Da], normalized migration time [min] and signal intensity. All detected peptides were deposited, matched, and annotated in a Microsoft SQL database, allowing further analysis and comparison of multiple samples. Polypeptides within different samples were considered identical if the mass deviation was lower than 650 ppm for masses ,4.000 Da, for masses between 4.000 and 6.000 Da gradually increasing to 6150 ppm, and 150 ppm for features .6 kDa. Acceptable migration time deviation was ,61 minutes for 19 min, gradually increasing to ,62.5 min at 50 min.

Statistical analysis
The reported P-values were calculated using the natural logarithm transformed intensities and the Wilcoxon test. Only peptides that were found at frequencies .30% in either case or control group were examined. The false discovery rate (FDR) adjustments of Banjamini-Hochberg [42] were employed to correct for multiple testing. Receiver operator characteristic (ROC) curves have been constructed and the area under the ROC curve (AUC) has been calculated using MedCalc version 8.1.1.0 (MedCalc Software, Belgium, www.medcalc.be).

Sequencing
Urine samples were analysed on a Dionex Ultimate 3000 RSLS nano flow system (Dionex, Camberly UK). After loading (5 ml) onto a Dionex 0.1620 mm 5 mm C18 nano trap column at a flowrate of 5 ml/min in 0.1% formic acid and 2% acetonitrile, sample was eluted onto an Acclaim PepMap C18 nano column 75 mm615 cm, 2 mm 100 Å at a flow rate of 0.3 ml/min. The trap and nano flow column were maintained at 35uC. The samples were eluted with a gradient of solvent A: 0.1% formic acid verses solvent B: 80% acetonitrile starting at 5% B rising to 50% B over 100 min.
The eluant was ionized using a Proxeon nano spray ESI source (Thermo Fisher Hemel UK) operating in positive ion mode into an Orbitrap Velos FTMS. Ionization voltage was 2.5 kV and the capillary temperature was 200uC. The mass spectrometer was operated in MS/MS mode scanning from 380 to 2000 amu. The top 10 multiply charged ions were selected from each scan for MS/MS analysis using HCD at 35% collision energy. The resolution of ions in MS1 was 60,000 and 7,500 for HCD MS2. Data files were searched against the IPI rat non-redundant database using the Open Mass Spectrometry Search Algorithm (OMSSA, http://pubchem.ncbi.nlm.nih.gov/omssa) and SE-QUEST (by using Thermo Proteome Discoverer), without any enzyme specificity. No fixed modification and oxidation of methionine and proline as variable modifications were selected. Mass error window of 10 ppm and 0.05 Da were allowed for MS and MS/MS, respectively. Peptide data were extracted using high peptide confidence and top one peptide rank filters. The OMSSA results were further optimized using COMPASS [43], 1% FDR was used as a cut-off value. The correlation between peptide charge at the working pH of 2 and CE-migration time was utilized to minimize false-positive identification rates [44]: Calculated CEmigration time based on the number of basic amino acids was compared to the experimental migration time. Accepted were only those peptides which were found with both search algorithms (OMSSA and SEQUEST) and having a mass deviation below 650 ppm and a CE-migration time deviation below 62 min.

Establishment of biomarker-based classifiers
For generation of disease-specific polypeptide patterns two different algorithms were used: Support vector machine (SVM)based MosaCluster software [45] and a linear combination of logtransformed, normalized data, as described [46,47]. MosaCluster (version 1.7.0) was developed for discrimination between different patient groups in the high-dimensional parameter space by using SVM learning. It generates high dimensional models, which rely on features (biomarkers) displaying statistically significant differences between data from patients with a specific disease to controls or other diseases. Each feature allegorizes one dimension in the ndimensional parameter space [48][49][50][51].
For linear combination, normalized signal intensity values below 1 were substituted with a value of 1 to avoid negative values by log-transformation. The average signal intensity for a specific biomarker over all cases was compared to the average intensity for the biomarker over all controls. To avoid artificial weighting of specific biomarkers in the set due to the difference in observed signal intensities for case and control, the relative distance between the two averages (case and control) was set to a value of 2. This relative distance of signal intensities between the disease and control samples was provided using the formula: A i k is the log-transformed signal intensity of the i th biomarker in the k th sample in either set, mean averages is the average of the mean intensity of all possible markers for test set samples, X X case represents the mean observed signal intensity of the possible biomarker from all X X control cases samples and represents the mean signal intensity of the possible biomarker from the combined control samples. Figure S1 Capillary electrophoresis coupled to mass spectrometry profiling of rat urine. The compiled data sets of urine samples from gentamicin-treated rats at different doses and sampling days are shown. Molecular mass of the analyzed polypeptides (kDa) in logarithmic scale is plotted against CE migration time (min). The mean signal intensity is represented in arbitrary units on the z-axis of the 3D plot. (PDF) Figure S2 Time course of mean signal intensities of all defined (n = 147) nephrotoxicity peptide maker candidates. Mean signal intensities of the respective peptide in urine samples from untreated animals and animals treated once daily for three consecutive days with 150 and 300 mg/kg gentamicin over time are shown. The first 39 diagrams depict the selected cis-platin markers and the subsequent figures the additional 108 gentamicin markers. (PDF)