A combinatorial approach of Proteomics and Systems Biology in unravelling the mechanisms of acute kidney injury (AKI): involvement of NMDA receptor GRIN1 in murine AKI
BMC Systems Biology volume 7, Article number: 110 (2013)
Acute kidney injury (AKI) is a frequent condition in hospitalised patients undergoing major surgery or the critically ill and is associated with increased mortality. Based on the volume of the published literature addressing this condition, reporting both supporting as well as conflicting molecular evidence, it is apparent that a comprehensive analysis strategy is required to understand and fully delineate molecular events and pathways which can be used to describe disease induction and progression as well as lead to a more targeted approach in intervention therapies.
We used a Systems Biology approach coupled with a de-novo high-resolution proteomic analysis of kidney cortex samples from a mouse model of folic acid-induced AKI (12 animals in total) and show comprehensive mapping of signalling cascades, gene activation events and metabolite interference by mapping high-resolution proteomic datasets onto a de-novo hypothesis-free dataspace. The findings support the involvement of the glutamatergic signalling system in AKI, induced by over-activation of the N-methyl-D-aspartate (NMDA)-receptor leading to apoptosis and necrosis by Ca2+-influx, calpain and caspase activation, and co-occurring reactive oxygen species (ROS) production to DNA fragmentation and NAD-rundown. The specific over-activation of the NMDA receptor may be triggered by the p53-induced protein kinase Dapk1, which is a known non-reversible cell death inducer in a neurological context. The pathway mapping is consistent with the involvement of the Renin-Angiotensin Aldosterone System (RAAS), corticoid and TNFα signalling, leading to ROS production and gene activation through NFκB, PPARγ, SMAD and HIF1α trans-activation, as well as p53 signalling cascade activation. Key elements of the RAAS-glutamatergic axis were assembled as a novel hypothetical pathway and validated by immunohistochemistry.
This study shows to our knowledge for the first time in a molecular signal transduction pathway map how AKI is induced, progresses through specific signalling cascades that may lead to end-effects such as apoptosis and necrosis by uncoupling of the NMDA receptor. Our results can potentially pave the way for a targeted pharmacological intervention in disease progression or induction.
A number of tightly controlled and complex processes are performed by heterogeneous cell populations of the kidney, including blood filtration by microvascular endothelial cells and podocytes and secretion and re-absorption by proximal tubular epithelial cells. Acute kidney injury (AKI) is a frequent clinical event associated with a rapid loss of kidney function, leading to high morbidity and mortality  in up to 22% of hospitalised patients . However, it is estimated that one-fifth of AKI that occurs after hospital admission is predictable and avoidable . In light of the economic burden of 0.4% to 0.6% of the total healthcare costs, corresponding to £400 m - £600 m, annually spent on treatment for acute kidney injury in the UK alone , a clear need for development of means to predict and/or early detect and prevent/treat AKI arises.
It has been by now widely accepted that AKI represents a continuum or spectrum of diseases that could be identified at an early stage, rather than the previous terminology of acute renal failure describing an “all or nothing” condition . The recognition of renal impairment at an early stage would allow for an immediate course of action to alleviate symptoms and disrupt the process of functional decline , however this implies that this condition is comprehensively understood on a molecular level to allow for targeted intervention therapies.
AKI can be caused by many different events such as ischemia, vasoconstrictive drugs, exposure to toxins, hypotension linked to sepsis, and obstruction of the urinary tract, and leads to a number of complications, including metabolic acidosis, high potassium levels, uremia, changes in body fluid balance, and damage to and failure of other organs . Molecular hallmarks of AKI are accumulation of free and esterified cholesterol, inflammation and inflammatory response, altered tubule dynamics leading to increased luminal sodium, hypoxia, cellular ATP depletion, renal cell apoptosis and necrosis, and hyperglycaemia, which is also a contributor to AKI (see Table 1).
The factors causing AKI and their interaction with each other are still incompletely understood, which is surprising in light of the vast amount of scientific reports in this area. However, a systematic and in-depth mining of the literature, as employed by a comprehensive Systems Biology approach, coupled with high-resolution mass spectrometry analysis of diseased tissue could pave the way forward. Systems Biology, or integrative biology, uses descriptive terms to understand complex interactions based on the analysis of network behaviour and dynamic aspects. Here, we demonstrate a comprehensive analysis of AKI induced in a mouse model system to investigate and identify key signalling cascades, inferred metabolic alterations and molecules modulated in this highly complex and clinically relevant disease. Nevertheless it should be noted that animal models of AKI might not reflect this clinically critical condition observed in man, where an unobserved laps of mere hours between onset and catastrophic irreparable kidney failure can mean life or death. An acceptable approximation is the model chosen here, where an single overdose of folic acid causes the same end-effect of kidney damage observed in clinically manifested AKI.
We first established a high-resolution proteomic map of protein expression profiles using 6 sham injected controls and 6 folic acid overdosed mice. Folic acid nephropathy is a well established murine model of AKI using a single intraperitoneal injection overdose of folic acid and harvesting of kidneys after 24 hours . Kidney cortex proteins were extracted, trypsin-digested, and analysed by high-resolution liquid chromatography mass spectrometry (LC-MS/MS). After merging the individual independent mass-spectrometric analyses, 6564 non-redundant proteins were detected (Additional file 1: Table S1, and detailed information of data processing can be found in Additional file 2). 14 of these proteins were detected only in control samples, 1050 only in AKI, and 5500 in both sets of conditions. Statistical testing using Mann–Whitney showed that 2521 proteins were statistically significant (p-values <0.05), and when also taking fold-change calculations into considerations then the dataset indicated 1480 molecules to be significantly changed in AKI (p-values <0.05 and fold changes >2). 420 of these proteins were found in both control and AKI samples, 14 only in controls, and 1046 molecules only in AKI samples. After initial database matching and employing gene ontology clustering to assess data quality and overlapping characteristics, the data were mapped onto existing general pathway maps to identify affected cascades and modulatory events. More than 50 signalling cascades (see Additional file 2), mainly due to overlapping features in the existing maps, were initially identified to be of interest in further downstream analysis based on their potential functional relevance. In parallel, we assembled potential pathway maps derived from the pertinent AKI-focussed literature. Through a manual iterative feedback analysis and deep-mining of inferred and prior knowledge (Figure 1), extending the pathway analysis beyond the reported AKI-signalling cascades, and molecule-by-molecule pathway delineation and manual feature look-up, we assembled plausible signalling cascades which led to a combined pathway map (Additional file 3: Figure S1). The same approach was also used to map the metabolic pathways modulated in AKI, as well as probable gene activation cascades based on reported modulation of transcription factors such as NFκB, which were integrated into the AKI model.
Key pathways which were previously reported to be main events in AKI induction, such as the Renin-Angiotensin Aldosterone System (RAAS) and involvement of the TNFα-signalling cascade, could be confirmed in our analysis to be up-regulated. Upstream activators, such as the kallikrein and cathepsin systems, were also induced, suggesting a specific activation of the aforementioned RAAS axis. Additional up-regulated modulatory events, which were per se not associated with the currently perceived molecular model of AKI are the glutamatergic signalling cascades and associated calcium-flux pathways, which have a major detrimental effect on both apoptosis and necrosis. A substantial amount of information is available about glutamate-dependent pathways and signalling events in a non-renal, specifically neurological context, and it was surprising to encounter a considerable level of glutamatergic pathway elements associated with renal dysfunction. The specific involvement under physiological conditions of ionotropic as well as metabotropic glutamate receptors in kidney is currently unknown, however a dysfunction, such as over-stimulation and –activation is expected to lead to the same effects observed in other systems, e.g. uncontrollable calcium-influx and ultimately cell death. This observation is further acerbated by an apparent simultaneous induction of the calcium-flux machinery, involving the calcium-import and –export channels, such as calcium pumps (SERCA and PMCA) as well as ryanodine receptors and calcium-sensitive modulators. A potential assembly of signalling events originating from the RAAS axis and involving the most prominent glutamate-sensitive calcium-channel NMDA receptor is depicted in Figure 2B. As shown, signalling from the renin-induced angiotensin receptor leads to a cascade of known signalling and induction events involving PLC2β, PKCα, Ras, RalA, p38kinase, MSK and activation of the transcription factor SP1. The latter promotes gene activation of the NMDA receptor GRIN1 which however is also dependent on SP3 inhibition. SP3 inhibition can be induced by oxidative stress, reflected in this scheme in the stimulation of NADPH oxidase NOX, predicted to result from the activation of the up-stream kinase PKCα. Oxidative stress also stimulates the death-associated protein kinase Dapk1, which can target the NMDA receptor and leads to a permanent-open state of the channel, allowing calcification of the intracellular environment. This in turn triggers the mitochondrial- as well as nuclear-based apoptotic and necrotic molecular machinery.
To enhance these hypotheses and predictions, we validated specific key elements in this delineated pathway by immunohistochemistry, as shown in Figure 2A. We also tested whether the potentially anti-apoptotic protein Hsp27 is regulated according to the MS-results. This protein is involved in cytoskeletal remodelling, which might also be an involvement of the NMDA receptor under physiological conditions, however based on our pathway model this would imply that apoptosis is to some extent inhibited, whereas necrosis is not. All of the proteins tested, namely PLC2β, Ras, MSK, Grin1, Dapk1 and Hsp27, could be confirmed to be up-regulated, as measured by the proteomic mass-spectrometry approach, thereby further emphasizing the importance of this downstream pathway in the development of the AKI phenotype. Unsurprisingly we could not detect any indication of the involvement of the insulin pathway, either up- or down-stream of the insulin receptor, but instead observed an up-regulation of known inhibitors of insulin-signalling such as Socs3 and Ptprf.
Additionally, our analysis of the glycolysis and glycogenolysis pathways, which are modulatory targets of insulin-signalling, shows a significant up-regulation of enzymes involved in glucose and fructose release, and no indication of an equally activated down-stream metabolism from Fructose-1,6-bisphosphate to pyruvate or lactate. This suggests an accumulation of Fructose-1,6-bisphosphate, which is a known allosteric inhibitor of Adss, an enzyme in the de novo pathway of AMP production. This inhibition would lead to an accumulation of IMP, and in conjunction with the observed up-regulation of 5′(3′)-deoxyribonucleotidase Nt5c, to inosine and ultimately hypoxanthine and xanthine build-up. Those last two stages are a characteristic hallmark of AKI, thereby explaining the previously described observations of an accumulation of the latter two chemical compounds on a molecular level (see Table 1).
AKI has a major impact on survival in hospitalised patients undergoing extensive surgery , and a large number of studies have been undertaken to understand the molecular causes leading to and being involved in this complex disease [10, 11]. Here, we employed a non-biased systematic approach using a combination of Systems Biology, proteomics and de-novo pathway mapping and immunohistochemical validation to elucidate major pathways and signalling cascades in order to contextualise previous findings and uncover novel potential targets for pharmacological intervention. The analysis focuses mainly on very significant changes observed in AKI covering all the known hallmarks using the folic acid-overdose animal model and can be regarded as a first attempt to describe this condition in a molecular mechanistic way.
An initial finding of this mapping effort was the apparent complexity of AKI and associated pathways, where a plethora of signalling cascades appears to be modulated simultaneously. This could be due to many kidney cell types being in an active state of inflammation signalling, apoptosis induction, stress, and other modulatory events, as also supported by the literature . As many cascades can be activated in more than one way, and induction of AKI can occur via multiple stimulation or entry points , it appears plausible that these pathways are merging to similar down-stream targets leading to the observed deleterious events in kidney injury.
Our analysis revealed that a major pathway involved in AKI is the RAAS axis, which has been reported multiple times before [14, 15], further confirming the validity of the approach. Renin activation can occur in several ways, either by kallikrein, cathepsin or other stimuli as indicated in Additional file 3: Figure S1. Many of those activating proteins were markedly up-regulated in AKI.
A further main crucial initiation step leading towards AKI is based on up-stream activation via TNFα, where inhibition of this molecule prevents apoptotic cell death . Mapping of modulated proteins found in this study clearly indicates the involvement of this cascade in AKI. TNF signals via Jnk to inhibit phosphoinositide 3-kinases (PI3Ks), which in turn leads to the release of cathepsin B-containing lysosomes and to renin activation , as well as activation of the pro-apoptotic protein p53 . TNFα may also potentially contribute to RAAS activation as shown in Figure 3. Mounting evidence suggests that the RAAS plays a major role in kidney injury and inflammatory processes , whereas TNFα and TNF-like cytokines are involved in induction of cellular responses such as inflammation as well as induction and progression of apoptotic cell death .
Along these lines, the data obtained also suggest up-regulation of inflammatory pathways, indicated in our analysis to involve activation by the RAAS axis through the NFκB pathway, leading to downstream signalling by interleukins (see Figure 3). These results are supported by an observed role for TLRs in promoting inflammation and tissue injury in AKI, which synchronise their inflammatory signal activations via NFκB . The activation of the RAAS axis also leads to the activation of a number of transcription factors , and our results indicate gene activation cascades involving among others NFκB, PPARγ, SMAD and HIF1α.
The significant changes in levels of e.g. collagens, Ena/VASP, Hsp27/Hspb1 and proteases such as MMP2 were interpreted by the bioinformatic tools to be correlated with cytoskeletal remodelling, extracellular matrix (ECM) degradation, and associated cell detachment, all linked to kidney injury . The cytoskeletal reorganisation by Hsp27 has been shown to inhibit fibrogenesis in obstructive nephropathy .
The marked up-regulation of the molecular clusters associated with apoptotic pathways, including NADPH oxidases (NOX) and P450 monooxygenases, nitric oxide synthase Nos1, calpains Capn3 and 11, caspase Casp8, poly-[ADP-ribose]-polymerase Parp3, 6, and 12 and Serine-protein kinase ATM, indicates activated cell-death cascades linked to prolonged hypoxia. Vasoconstriction, one of the hallmarks of AKI, is directly responsible for hypoxia, which stimulates several proteins such as NOX and P450 monooxygenases and leads to ROS generation, including superoxide and hydrogen peroxide production by xanthine oxidase as a result of metabolic hypoxanthine accumulation . This is acerbated by the notion that kidney injury can also induce the over-production of nitric oxide (NO), which, with superoxide, forms peroxynitrite, and causes cell damage via DNA fragmentation . Over-production of ROS is linked to the induction of apoptosis and necrosis, and a number of pro-apoptotic and necrotic pathways have been shown to exist, some of which also require the activation of the mitochondrial calpain system through sustained Ca2+-influx . The classical apoptotic pathway involves the release of mitochondrial proteins, leading ultimately to activation of caspases . Regulated necrosis also contributes to cell loss  (Figure 3, Additional file 3: Figure S1 and Table 1).
A major novel finding is the modulation of proteins involved in glutamatergic signalling. This is exemplified by the up-regulation of key components of the glutamate response system, such as the metabotropic glutamate receptor Grm1 (also known as mGluR1), and ionotropic glutamate receptors Grin1 (NMDA receptor NR1), Gria2 (AMPA receptor GluR2), and specific scaffolders of these receptors, as well as known down-stream signalling molecules (e.g. CamK2b, PP2B, and Nos1) (Additional file 3: Figure S1). Activation of glutamatergic signalling components is one possible route to activate Ca2+-flux pathways . Podocytes and tubular cells express glutamate receptors  including functional Grin1 NMDA receptors , raising the possibility that neuron-like signalling contributes to glomerular function . Glutamatergic signalling driven by podocytes contributes to the integrity of the glomerular filtration barrier and may regulate glomerular filtration . Additionally, one of the metabotropic glutamate receptors, Grm1, is expressed in glomerular podocytes . Although the stimulus for increased NMDA receptor expression in AKI has not been characterised, hypoxia is known to induce the transcription of glutamatergic signalling components such as the NMDA receptors in non-renal cells through inhibition of the transcription factor SP3. SP3 inhibits transcription by occupying the genomic translation initiation site, and activation of SP1 which in turn occupies the freed site and promotes the transcription of target genes .
Collectively, based on the presented data, we could hypothesize that excessive activation or sustained activation of NMDA receptors in tubular cells or podocytes promotes oxidative stress, mobilization of cation channels, disproportionate Ca2+ influx and overloading, excessive ROS generation, and apoptotic cell death  similar to what was shown to occur in the neurological environment . Indeed tissue glutamate is high in AKI  and pharmacologically blocking this receptor with D-AP5, a synthetic glutamate analogue which binds to the glutamate-binding site of the NMDA receptor and suppresses pore opening, prevented tissue injury . Shown in Additional file 3: Figure S1 are proteins and signalling events, assembled to depict probable cascades driving both, apoptosis and necrosis in kidney tissue after injury. The glutamate receptor Grin1/NMDA-R1, which is responsible for cellular Ca2+-overload leading to irreversible caspase activation and apoptosis, in conjunction with sustained ROS production, and ultimately energy-depletion and necrosis, may be a major contributor to these events.
Dapk1, a gene activation target of p53, can mediate the pro-apoptotic activity of TNFα by inhibition of NFκB signalling . Additionally, the NMDA receptor is a known target of Dapk1, where Dapk1-mediated channel modulation results in a permanently open NMDA receptor, leading to cell death . Remarkably, we also observed an orchestrated up-regulation of relevant scaffolders (Homer1/3, Dlg1, and Dlgap4) which link the glutamate receptors directly to downstream signalling cascades ranging from Ca2+-signalling to phospholipase and adenylate cyclase cascades and interlinking pathways. This further supports the importance of these receptors in AKI.
Other Ca2+ channels, which we observed to be up-regulated in AKI, and that could potentially contribute to Ca2+-overload are Ca2+-import channels such as voltage-sensitive calcium channels VSCC (isoforms Cacna1g/h and Cacnb1) and transient receptor potential channels TRP (isoforms Trpc4, Trpm1/3), intracellular release channels sarcoplasmic/endoplasmic reticulum calcium ATPases SRCA (isoforms Atp2a1/2) which are dependent on ATP hydrolysis, and others (Slc8a3, Ryr1). Simultaneously we also detected up-regulation of Ca2+-efflux channels plasma membrane calcium-transporting ATPases PMCA (isoforms Atp2b1/2/3), which are dependent on available ATP to function. The depletion of intracellular ATP pools could therefore lead to an asymmetric Ca2+-flux and exacerbate the intracellular calcification even further.
We also extracted supportive information from the literature, where pharmacological or other intervention studies reported modulatory effects on kidney injury. The results are summarised in Additional file 4: Table S2, and also indicated in Figure 3. Inhibition of molecules downstream of TNFα and involved in ROS production and hyperglycaemia, as well as specific activation of NFκB cascades, all had an improving or attenuated outcome on AKI prevention. Of particular interest is the inhibition of TNFα signalling by using TNFα antibodies which resulted in prevention of AKI . Inhibition of glutamate signalling via the NMDA receptor, by using the channel blocker D-AP5, has been shown to significantly reduce ischemia/reperfusion injury (I/RI) -induced glomerular and tubular dysfunction , and NMDA-induced injury, which leads to calcification of the intracellular space, was shown to be reduced by using calcium antagonists nimodipine and dantrolene in an animal model of retinal injury . Nimodipine targets and blocks L-type calcium channels, whereas dantrolene inhibits ryanodine receptors. Nimodipine was shown to alleviate tubular cell necrosis in transplantation-induced AKI . Dantrolene in particular gained interest as it was shown to be beneficial in protecting against I/RI in animal models of heart-, brain- and potentially liver-induced I/RI . However, dantrolene was shown to be ineffective in protecting renal function . The same study investigated two other calcium-channel inhibitors, namely nicardipine, which blocks L-type calcium channels such as CACNB1, and TMB-8, blocking inositol 1,4,5-trisphosphate receptors. They could demonstrate that both compounds applied in pre-conditioning prior to I/RI had a positive effect both in vitro on tubular cells by decreasing the apoptotic effect, and in vivo in a rat renal model in kidney functional integrity. The results of those pharmacological studies listed in Additional file 4: Table S2 suggest that AKI in principle can be prevented or its effect diminished, and a combination of various drugs, targeting specific AKI-induced pathways and molecules, could be a way forward in disease prevention or therapeutic intervention.
A potential synopsis of pathways and signalling events evoked during AKI leading to Ca2+-overload and apoptosis as well as necrosis is shown in Figure 4, and a focus centred around the NMDA receptor is shown in Figure 2B. The relevant endogenous NMDA receptor ligand in AKI has not been characterised, however kidney glutamate levels are increased in AKI . To further substantiate this cascade, we validated the up-regulation of key molecules involved in this axis as by immunohistochemistry shown in Figure 2A. The signalling cascade downstream of angiotensin-receptor 1 activates PLC2β which signals to PKCα. Both these proteins were up-regulated in AKI based on the mass-spectrometric analysis and confirmed for PLC2β by immunohistochemistry. PKCα activation leads to ROS production via NOX , and induces gene activation of the NMDA receptor, which we found to be up-regulated by both techniques used. This gene induction step also requires activation of MSK, which we also found to be up-regulated. Additionally, we could show that Hsp27 is up-regulated, indicative of cytoskeletal reorganisation, and we could demonstrate that Dapk1 is also up-regulated. Hsp27 up-regulation may protect from angiotensin II- and high-glucose-induced apoptosis . All these events are potentially linked to both apoptosis and necrosis, and through this up-regulated pathway could be demonstrating a probable route of tubular epithelial fate in AKI.
The tissue proteomics approach identified here may eventually be optimized for its use in diagnostic/staging clinical pathology and the molecular fingerprint identified may be used to explore different molecular types of clinical AKI that may require different therapeutic approaches.
In summary we assembled a global pathway model for AKI based on proteomic changes, and further support the concept that the approach of in-depth Systems Biology coupled with high-resolution proteomic analysis represents a valuable approach to elucidate complex pathways and systems, especially in multifaceted diseases such as AKI. This combined methodology can be applied to search for druggable targets and identify key nodes where molecular/pharmacological interventions are most likely to have an impact on the system as a whole, as well as predictions for disease biomarkers. One such target could be the glutamatergic system, however the NMDA receptor itself might not be an ideal lead due to its involvement in neurotransmission in other areas like cardiac function, but a renally localised combinatorial pharmacological intervention, blocking both glutamatergic signalling and Ca2+-flux could be of beneficial effect.
Folic acid nephropathy is a classical model of AKI . C57/BL6 mice from IFFA-CREDO (Barcelona, Spain) (12- to 14-week-old) received a single i.p. injection of folic acid (Sigma) 250 mg/kg in 0.3 mol/L sodium bicarbonate (AKI animals) or vehicle alone (controls) and mice were killed 24 h later . Kidneys were cold saline perfused in situ before removal. One kidney from each mouse was fixed in buffered formalin, embedded in paraffin and stained with hematoxylin-eosin or used for immunohistochemistry. The kidney cortex from other kidney was snap-frozen in liquid nitrogen for protein studies. The study was approved by the IIS-FJD animal ethics committee and followed Directive 2010/63/EU on the protection of animals used for scientific purposes.
Tissue samples were weighed out and extracted using the Filter Aided Sample Preparation (FASP) method . Briefly, tissue samples were homogenised in SDS-lysis buffer (1:10 sample to buffer ratio) (0.1 M Tris–HCl pH 7.6 supplemented with 4% SDS and 0.1 M DTT) using an Ultra-Turrax T 25 (IKA, Staufen, Germany), incubated at 95°C for 3 minutes and clarified by centrifugation at 16,000 g for 5 min at room temperature. An aliquot of the supernatant was taken and placed in a Micron YM-30 filter device (Millipore, Watford UK). 8 M Urea buffer (UA) was added to the protein extract and then centrifuged at 14,000 g for 15 minutes. This step was repeated twice after which the protein extract was mixed gently for 1 minute with 0.05 M iodoacetamide buffer (IAA) and incubated for a further 20 minutes prior to centrifugation. UA buffer was again added and centrifuged (twice). Ammonium bicarbonate buffer (50 mM of NH4HCO3, pH 8) (ABC) was then added and centrifuged (twice) before incubating overnight with trypsin. The trypsin homogenate was centrifuged and washed with ABC buffer prior to acidification with 10% formic acid. Sample volumes were adjusted to match final concentration of protein prior to analysis by LC-MS/MS.
LC-MS/MS mass spectrometry analysis
Tissue extracts were separated on a Dionex Ultimate 3000 RSLS nano flow system (Dionex, Camberly UK). A 5 μl sample was loaded in 0.1% formic acid and acetonitrile (98:2) onto a Dionex 100 μm × 2 cm, 5 μm C18 nano trap column at a flowrate of 5 μl/min. Elution was performed on an Acclaim PepMap C18 nano column 75 μm × 50 cm, 2 μm, 100 Å with a linear gradient of solvent A, 0.1% formic acid and acetonitrile (98:2) against solvent B, 0.1% formic acid and acetonitrile (20:80) starting at 1% B for 5 minutes rising to 30% at 400 minutes then to 50%B at 480 minutes. The sample was ionized in positive ion mode using a Proxeon nano spray ESI source (Thermo Fisher Hemel UK) and analyzed in an Orbitrap Velos FTMS (Thermo Finnigan, Bremen, Germany). The MS was operated in a data-dependent mode (top 40) to switch between MS and MS/MS acquisition and parent ions were fragmented by collision-induced dissociation (CID). Data files were searched against the IPI mouse non-redundant database (version 3.8.5) using SEQUEST with enzyme specified as trypsin. A fixed modification of carbamidomethylation was set and oxidation of methionine and proline as variable modifications were selected. Mass error windows of 20 ppm and 0.8 Da were allowed for MS and MS/MS, respectively. In SEQUEST, only peptides that showed mass deviation of less than 10 ppm were passed, the peptide data were extracted using high peptide confidence and top one peptide rank filters. Statistical p-value analysis was performed using the Mann Whitney test (implemented in in-house software). Normalisation of the mass spectrometry data was done using Histone H2B. Similar results were obtained when the normalisation was performed against actin (skeletal muscle isoform).
Systems Biology analysis
Data merging was done by Blast-searching the SwissProt database with all 9930 identified molecules individually, or by batch-transfer using the UniProt online tool, to transfer all IPI accession numbers to the SwissProt identifiers, followed by combining all duplicated entries. Mass spectrometry measurement values were adjusted using the average or, in case where the protein was only found in one instance in a sample, by taking only the measured value. This dataset of 6564 entries was then used in the downstream analysis (Additional file 1: Table S1). The p-values were re-calculated using the Mann Whitney test, and the fold changes calculated by using the averages where missing values were ignored. 1480 molecules were significantly changed in AKI using high-stringency cut-off values of p-values <0.05 and fold changes >2. Of those, 1046 were only found in the AKI kidney samples, 14 only in controls, and 420 in both. Gene ontology analysis was performed using the CytoScape plug-in ClueGO. Initial pathway mapping was done using KEGG (online) as well as WikiPathway maps (using PathVisio software), where 28 partially overlapping metabolic and signalling cascades could be identified. Data for molecular properties were obtained from UniProt and the published literature, as well as other databases (PADB, BRENDA, Reactome, IntAct, GeneCards, NCBI RefSeq and Entrez). De-novo pathway maps were constructed by extended manual literature searches. Gene activation cascades were put together by manual literature mining and use of on-line resources such as GeneCards. AKI-specific metabolic maps were created based on information from KEGG, Reactome, BRENDA and UniProt, as well as the published literature. The delineated pathways were then combined into plausible signalling cascades and initial sub-models were manually assembled using the PathVisio software. A final model was established after several re-iterations and literature cross-checks.
Immunohistochemistry was carried out in 5 μm thick paraffin-embedded tissue sections . The primary antibodies were rabbit polyclonal anti-PLCβ2, goat polyclonal anti-HSP27, rat monoclonal anti-H-Ras, rabbit polyclonal anti-MSK1 (1:100, Santa Cruz biotechnology), rabbit polyclonal anti-NMDAR1 (1:100, Abcam) and rabbit polyclonal anti-DAP Kinase 1 (1:100, Sigma). Sections were counterstained with Carazzi’s hematoxylin. Negative controls included incubation with a non-specific immunoglobulin of the same isotype as the primary antibody. These sections were mounted in 90% glycerol/PBS and examined in a blinded manner.
Acute kidney injury
Liquid chromatography mass spectrometry
Renin-Angiotensin Aldosterone System
Reactive oxygen species.
Sharfuddin AA, Molitoris BA: Pathophysiology of ischemic acute kidney injury. Nat Rev Nephrol. 2011, 7: 189-200.
Aitken E, Carruthers C, Gall L, Kerr L, Geddes C, Kingsmore D: Acute kidney injury: outcomes and quality of care. QJM. 2013, 106: 323-332. 10.1093/qjmed/hcs237.
Selby NM, Kolhe NV, McIntyre CW, Monaghan J, Lawson N, Elliott D, Packington R, Fluck RJ: Defining the cause of death in hospitalised patients with acute kidney injury. PLoS One. 2012, 7: e48580-10.1371/journal.pone.0048580.
Prescott AM, Lewington A, O’Donoghue D: Acute kidney injury: top ten tips. Clin Med. 2012, 12: 328-332. 10.7861/clinmedicine.12-4-328.
Moore EM, Bellomo R, Nichol AD: The meaning of acute kidney injury and its relevance to intensive care and anaesthesia. Anaesth Intensive Care. 2012, 40: 929-948.
Borthwick E, Ferguson A: Perioperative acute kidney injury: risk factors, recognition, management, and outcomes. BMJ. 2010, 341: c3365-10.1136/bmj.c3365.
An JN, Lee JP, Jeon HJ, Kim DH, Oh YK, Kim YS, Lim CS: Severe hyperkalemia requiring hospitalization: predictors of mortality. Crit Care. 2012, 16: R225-10.1186/cc11872.
Wen X, Peng Z, Li Y, Wang H, Bishop JV, Chedwick LR, Singbartl K, Kellum JA: One dose of cyclosporine A is protective at initiation of folic acid-induced acute kidney injury in mice. Nephrol Dial Transplant. 2012, 27: 3100-3109. 10.1093/ndt/gfr766.
Ricci Z, Cruz DN, Ronco C: Classification and staging of acute kidney injury: beyond the RIFLE and AKIN criteria. Nat Rev Nephrol. 2011, 7: 201-208. 10.1038/nrneph.2011.14.
McKay DB: Intracellular pattern recognition receptors and renal ischemia. Crit Rev Immunol. 2011, 31: 297-306. 10.1615/CritRevImmunol.v31.i4.20.
Cantaluppi V, Quercia AD, Dellepiane S, Figliolini F, Medica D, De DLM: [New mechanisms and recent insights in the pathogenesis of acute kidney injury (AKI)]. G Ital Nefrol. 2012, 29: 535-547.
Pretorius M, Murray KT, Yu C, Byrne JG, Billings FT, Petracek MR, Greelish JP, Hoff SJ, Ball SK, Mishra V, Body SC, Brown NJ: Angiotensin-converting enzyme inhibition or mineralocorticoid receptor blockade do not affect prevalence of atrial fibrillation in patients undergoing cardiac surgery. Crit Care Med. 2012, 40: 2805-2812. 10.1097/CCM.0b013e31825b8be2.
Efrati S, Berman S, Abu HR, El NR, Chanimov M, Siman-Tov Y, Weissgarten J: Hyperglycaemia, inflammation, RAS activation: three culprits to blame for acute kidney injury emerging in healthy rats during general anaesthesia. Nephrology (Carlton ). 2012, 17: 591-602. 10.1111/j.1440-1797.2012.01638.x.
Jackson EK, Herzer WA, Vyas SJ, Kost CK: Angiotensin II-induced renal vasoconstriction in genetic hypertension. J Pharmacol Exp Ther. 1999, 291: 329-334.
Nangaku M, Fujita T: Activation of the renin-angiotensin system and chronic hypoxia of the kidney. Hypertens Res. 2008, 31: 175-184. 10.1291/hypres.31.175.
Wan B, Hao L, Qiu Y, Sun Z, Cao Q, Zhang Y, Zhu T, Wang H, Zhang Y: Blocking tumor necrosis factor-alpha inhibits folic acid-induced acute renal failure. Exp Mol Pathol. 2006, 81: 211-216. 10.1016/j.yexmp.2006.02.005.
Li JH, Pober JS: The cathepsin B death pathway contributes to TNF plus IFN-gamma-mediated human endothelial injury. J Immunol. 2005, 175: 1858-1866.
Fuchs SY, Adler V, Pincus MR, Ronai Z: MEKK1/JNK signaling stabilizes and activates p53. Proc Natl Acad Sci USA. 1998, 95: 10541-10546. 10.1073/pnas.95.18.10541.
Ruster C, Wolf G: Angiotensin II as a morphogenic cytokine stimulating renal fibrogenesis. J Am Soc Nephrol. 2011, 22: 1189-1199. 10.1681/ASN.2010040384.
Sanchez-Nino MD, Benito-Martin A, Goncalves S, Sanz AB, Ucero AC, Izquierdo MC, Ramos AM, Berzal S, Selgas R, Ruiz-Ortega M, Egido J, Ortiz A: TNF superfamily: a growing saga of kidney injury modulators. Mediators Inflamm. 2010, 2010:
Allam R, Scherbaum CR, Darisipudi MN, Mulay SR, Hagele H, Lichtnekert J, Hagemann JH, Rupanagudi KV, Ryu M, Schwarzenberger C, Hohenstein B, Hugo C, Uhl B, Reichel CA, Krombach F, Monestier M, Liapis H, Moreth K, Schaefer L, Anders HJ: Histones from dying renal cells aggravate kidney injury via TLR2 and TLR4. J Am Soc Nephrol. 2012, 23: 1375-1388. 10.1681/ASN.2011111077.
Camici M: The Nephrotic Syndrome is an immunoinflammatory disorder. Med Hypotheses. 2007, 68: 900-905. 10.1016/j.mehy.2006.04.072.
Lee PT, Chou KJ, Fang HC: Are tubular cells not only victims but also perpetrators in renal fibrosis?. Kidney Int. 2012, 82: 128-130. 10.1038/ki.2012.120.
Vidyasagar A, Reese SR, Hafez O, Huang LJ, Swain WF, Jacobson LM, Torrealba JR, Chammas PE, Wilson NA, Djamali A: Tubular expression of heat-shock protein 27 inhibits fibrogenesis in obstructive nephropathy. Kidney Int. 2013, 83: 84-92. 10.1038/ki.2012.336.
Devarajan P: Update on mechanisms of ischemic acute kidney injury. J Am Soc Nephrol. 2006, 17: 1503-1520. 10.1681/ASN.2006010017.
Sanz AB, Santamaria B, Ruiz-Ortega M, Egido J, Ortiz A: Mechanisms of renal apoptosis in health and disease. J Am Soc Nephrol. 2008, 19: 1634-1642. 10.1681/ASN.2007121336.
Kar P, Samanta K, Shaikh S, Chowdhury A, Chakraborti T, Chakraborti S: Mitochondrial calpain system: an overview. Arch Biochem Biophys. 2010, 495: 1-7. 10.1016/j.abb.2009.12.020.
Linkermann A, De ZF, Weinberg J, Kunzendorf U, Krautwald S: Programmed necrosis in acute kidney injury. Nephrol Dial Transplant. 2012, 27: 3412-3419. 10.1093/ndt/gfs373.
Wang JH, Kelly P: Calcium-calmodulin signalling pathway up-regulates glutamatergic synaptic function in non-pyramidal, fast spiking rat hippocampal CA1 neurons. J Physiol. 2001, 533: 407-422. 10.1111/j.1469-7793.2001.0407a.x.
Armelloni S, Li M, Messa P, Rastaldi MP: Podocytes: a new player for glutamate signaling. Int J Biochem Cell Biol. 2012, 44: 2272-2277. 10.1016/j.biocel.2012.09.014.
Parisi E, Bozic M, Ibarz M, Panizo S, Valcheva P, Coll B, Fernandez E, Valdivielso JM: Sustained activation of renal N-methyl-D-aspartate receptors decreases vitamin D synthesis: a possible role for glutamate on the onset of secondary HPT. Am J Physiol Endocrinol Metab. 2010, 299: E825-E831. 10.1152/ajpendo.00428.2010.
Bozic M, De RJ, Parisi E, Ortega MR, Fernandez E, Valdivielso JM: Glutamatergic signaling maintains the epithelial phenotype of proximal tubular cells. J Am Soc Nephrol. 2011, 22: 1099-1111. 10.1681/ASN.2010070701.
Kottgen M, Benzing T: Strangers on a train: atypical glutamate receptors in the kidney glomerulus. Focus on “Functional NMDA receptors with atypical properties are expressed in podocytes”. Am J Physiol Cell Physiol. 2011, 300: C9-C10. 10.1152/ajpcell.00453.2010.
Gu L, Liang X, Wang L, Yan Y, Ni Z, Dai H, Gao J, Mou S, Wang Q, Chen X, Wang L, Qian J: Functional metabotropic glutamate receptors 1 and 5 are expressed in murine podocytes. Kidney Int. 2012, 81: 458-468. 10.1038/ki.2011.406.
Hota SK, Hota KB, Prasad D, Ilavazhagan G, Singh SB: Oxidative-stress-induced alterations in Sp factors mediate transcriptional regulation of the NR1 subunit in hippocampus during hypoxia. Free Radic Biol Med. 2010, 49: 178-191. 10.1016/j.freeradbiomed.2010.03.027.
Leung JC, Ragland N, Marphis T, Silverstein DM: NMDA agonists and antagonists induce renal culture cell toxicity. Med Chem. 2008, 4: 565-571. 10.2174/157340608786242034.
Hardingham GE: Coupling of the NMDA receptor to neuroprotective and neurodestructive events. Biochem Soc Trans. 2009, 37: 1147-1160. 10.1042/BST0371147.
Keller AK, Jorgensen TM, Vittrup DM, Kjerkegaard UK, Jespersen B, Krag SR, Bibby BM, Stolle LB: Fast detection of renal ischemia in transplanted kidneys with delayed graft function-an experimental study. Transplantation. 2013, 95: 275-279. 10.1097/TP.0b013e318276a1c8.
Yang CC, Chien CT, Wu MH, Ma MC, Chen CF: NMDA receptor blocker ameliorates ischemia-reperfusion-induced renal dysfunction in rat kidneys. Am J Physiol Renal Physiol. 2008, 294: F1433-F1440. 10.1152/ajprenal.00481.2007.
Yoo HJ, Byun HJ, Kim BR, Lee KH, Park SY, Rho SB: DAPk1 inhibits NF-kappaB activation through TNF-alpha and INF-gamma-induced apoptosis. Cell Signal. 2012, 24: 1471-1477. 10.1016/j.cellsig.2012.03.010.
Tu W, Xu X, Peng L, Zhong X, Zhang W, Soundarapandian MM, Balel C, Wang M, Jia N, Zhang W, Lew F, Chan SL, Chen Y, Lu Y: DAPK1 interaction with NMDA receptor NR2B subunits mediates brain damage in stroke. Cell. 2010, 140: 222-234. 10.1016/j.cell.2009.12.055.
Kaya M, Tunc M, Ozdemir T, Altuntas I: Calcium antagonists in N-methyl d-aspartate-induced retinal injury. Graefes Arch Clin Exp Ophthalmol. 2003, 241: 418-422. 10.1007/s00417-002-0610-8.
Fisher M, Grotta J: New uses for calcium channel blockers. Therapeutic implications. Drugs. 1993, 46: 961-975. 10.2165/00003495-199346060-00003.
Boys JA, Toledo AH, Anaya-Prado R, Lopez-Neblina F, Toledo-Pereyra LH: Effects of dantrolene on ischemia-reperfusion injury in animal models: a review of outcomes in heart, brain, liver, and kidney. J Investig Med. 2010, 58: 875-882.
Wu D, Chen X, Ding R, Qiao X, Shi S, Xie Y, Hong Q, Feng Z: Ischemia/reperfusion induce renal tubule apoptosis by inositol 1,4,5-trisphosphate receptor and L-type Ca2+ channel opening. Am J Nephrol. 2008, 28: 487-499. 10.1159/000113107.
Zhang C, Yi F, Xia M, Boini KM, Zhu Q, Laperle LA, Abais JM, Brimson CA, Li PL: NMDA receptor-mediated activation of NADPH oxidase and glomerulosclerosis in hyperhomocysteinemic rats. Antioxid Redox Signal. 2010, 13: 975-986. 10.1089/ars.2010.3091.
Sanchez-Nino MD, Sanz AB, Sanchez-Lopez E, Ruiz-Ortega M, Benito-Martin A, Saleem MA, Mathieson PW, Mezzano S, Egido J, Ortiz A: HSP27/HSPB1 as an adaptive podocyte antiapoptotic protein activated by high glucose and angiotensin II. Lab Invest. 2012, 92: 32-45. 10.1038/labinvest.2011.138.
Sanz AB, Justo P, Sanchez-Nino MD, Blanco-Colio LM, Winkles JA, Kreztler M, Jakubowski A, Blanco J, Egido J, Ruiz-Ortega M, Ortiz A: The cytokine TWEAK modulates renal tubulointerstitial inflammation. J Am Soc Nephrol. 2008, 19: 695-703. 10.1681/ASN.2007050577.
Wisniewski JR, Zougman A, Nagaraj N, Mann M: Universal sample preparation method for proteome analysis. Nat Methods. 2009, 6: 359-362. 10.1038/nmeth.1322.
Sanchez-Nino MD, Sanz AB, Lorz C, Gnirke A, Rastaldi MP, Nair V, Egido J, Ruiz-Ortega M, Kretzler M, Ortiz A: BASP1 promotes apoptosis in diabetic nephropathy. J Am Soc Nephrol. 2010, 21: 610-621. 10.1681/ASN.2009020227.
Grant support: FEDER funds FIS ISCIII-RETIC REDinREN (RD12/0021 and PS09/00447), Comunidad de Madrid (CIFRA S2010/BMD-2378). Salary support: FIS-Sara Borrell to MDSN, Programa Intensificación Actividad Investigadora (ISCIII/Agencia Laín-Entralgo/CM) to AO.
HM is the owner and CEO of Mosaiques GmbH, Hannover, Germany. No other competing interests are declared.
HH wrote the draft of the manuscript and did the Systems Biology analysis, WM did the mass spectrometry measurements, MDSN did the animal model and the immunohistochemistry work, CD contributed to the clinical relevance of modelling and correctness, AV initiated the project, and both AO and HM led the research. All authors contributed to the writing of the manuscript. All authors read and approved the final manuscript.
Holger Husi, Maria Dolores Sanchez-Niño, Alberto Ortiz and Harald Mischak contributed equally to this work.
Electronic supplementary material
Additional file 1: Table S1: Proteomic and bioinformatic analysis of the AKI dataset. Data entries were matched to various databases by accession numbers such as SwissProt, UniProt, EnsEMBL, TrEMBL, PIR, Unigene, Entrez (GI), IPI and NCBI-GeneIndex. Raw data from the six control and six AKI samples, as measured by LC-MS, and derived values such as fold changes and p-values are included. (XLS 5 MB)
Additional file 2: This file contains in-depth methods and results used to generate the molecular pathway map. It also contains the results of GO, interactome and KEGG pathway mapping analysis, as well as the original assembly of pathway cascades used to establish the global pathway shown in Additional file 3: Figure S1. (PDF 11 MB)
Additional file 3: Figure S1: Pathway map in AKI. Molecules implicated in AKI were mapped onto pathways based on prior knowledge of signalling cascades as well as merged de-novo pathways. Fold changes of individual proteins range from green (down-regulation) to white (unchanged) to red (up-regulation). Grey denotes proteins and metabolites without fold-change information. ROS are highlighted with a green surrounding box. The legend is included in the top-left corner. (PDF 3 MB)
Additional file 4: Table S2: Reported therapeutic agents employed in acute kidney injury models. Abbreviated names are indicated in round brackets, and gene names in square brackets. (DOC 438 KB)
About this article
Cite this article
Husi, H., Sanchez-Niño, M.D., Delles, C. et al. A combinatorial approach of Proteomics and Systems Biology in unravelling the mechanisms of acute kidney injury (AKI): involvement of NMDA receptor GRIN1 in murine AKI. BMC Syst Biol 7, 110 (2013). https://doi.org/10.1186/1752-0509-7-110