Thunder-DDA-PASEF enables high-coverage immunopeptidomics and is boosted by MS2Rescore with MS2PIP timsTOF fragmentation prediction model.

Human leukocyte antigen (HLA) class I peptide ligands (HLAIps) are key targets for developing vaccines and immunotherapies against infectious pathogens or cancer cells. Identifying HLAIps is challenging due to their high diversity, low abundance, and patient individuality. Here, we develop a highly sensitive method for identifying HLAIps using liquid chromatography-ion mobility-tandem mass spectrometry (LC-IMS-MS/MS). In addition, we train a timsTOF-specific peak intensity MS2PIP model for tryptic and non-tryptic peptides and implement it in MS2Rescore (v3) together with the CCS predictor from ionmob. The optimized method, Thunder-DDA-PASEF, semi-selectively fragments singly and multiply charged HLAIps based on their IMS and m/z. Moreover, the method employs the high sensitivity mode and extended IMS resolution with fewer MS/MS frames (300 ms TIMS ramp, 3 MS/MS frames), doubling the coverage of immunopeptidomics analyses, compared to the proteomics-tailored DDA-PASEF (100 ms TIMS ramp, 10 MS/MS frames). Additionally, rescoring boosts the HLAIps identification by 41.7% to 33%, resulting in 5738 HLAIps from as little as one million JY cell equivalents, and 14,516 HLAIps from 20 million. This enables in-depth profiling of HLAIps from diverse human cell lines and human plasma. Finally, profiling JY and Raji cells transfected to express the SARS-CoV-2 spike protein results in 16 spike HLAIps, thirteen of which have been reported to elicit immune responses in human patients.

Ionmob: a Python package for prediction of peptide collisional cross-section values.

MOTIVATION: Including ion mobility separation (IMS) into mass spectrometry proteomics experiments is useful to improve coverage and throughput. Many IMS devices enable linking experimentally derived mobility of an ion to its collisional cross-section (CCS), a highly reproducible physicochemical property dependent on the ion's mass, charge and conformation in the gas phase. Thus, known peptide ion mobilities can be used to tailor acquisition methods or to refine database search results. The large space of potential peptide sequences, driven also by posttranslational modifications of amino acids, motivates an in silico predictor for peptide CCS. Recent studies explored the general performance of varying machine-learning techniques, however, the workflow engineering part was of secondary importance. For the sake of applicability, such a tool should be generic, data driven, and offer the possibility to be easily adapted to individual workflows for experimental design and data processing. RESULTS: We created ionmob, a Python-based framework for data preparation, training, and prediction of collisional cross-section values of peptides. It is easily customizable and includes a set of pretrained, ready-to-use models and preprocessing routines for training and inference. Using a set of ≈21 000 unique phosphorylated peptides and ≈17 000 MHC ligand sequences and charge state pairs, we expand upon the space of peptides that can be integrated into CCS prediction. Lastly, we investigate the applicability of in silico predicted CCS to increase confidence in identified peptides by applying methods of re-scoring and demonstrate that predicted CCS values complement existing predictors for that task. AVAILABILITY AND IMPLEMENTATION: The Python package is available at github: https://github.com/theGreatHerrLebert/ionmob.

Modulation of cellular transcriptome and proteome composition by azidohomoalanine-implications on click chemistry-based secretome analysis.

The analysis of the secretome provides important information on proteins defining intercellular communication and the recruitment and behavior of cells in specific tissues. Especially in the context of tumors, secretome data can support decisions for diagnosis and therapy. The mass spectrometry-based analysis of cell-conditioned media is widely used for the unbiased characterization of cancer secretomes in vitro. Metabolic labeling using azide-containing amino acid analogs in combination with click chemistry facilitates this type of analysis in the presence of serum, preventing serum starvation-induced effects. The modified amino acid analogs, however, are less efficiently incorporated into newly synthesized proteins and may perturb protein folding. Combining transcriptome and proteome analysis, we elucidate in detail the effects of metabolic labeling with the methionine analog azidohomoalanine (AHA) on gene and protein expression. Our data reveal that 15-39% of the proteins detected in the secretome displayed changes in transcript and protein expression induced by AHA labeling. Gene Ontology (GO) analyses indicate that metabolic labeling using AHA leads to induction of cellular stress and apoptosis-related pathways and provide first insights on how this affects the composition of the secretome on a global scale. KEY MESSAGES: Azide-containing amino acid analogs affect gene expression profiles. Azide-containing amino acid analogs influence cellular proteome. Azidohomoalanine labeling induces cellular stress and apoptotic pathways. Secretome consists of proteins with dysregulated expression profiles.

Eculizumab treatment alters the proteometabolome beyond the inhibition of complement.

Therapeutic strategies targeting complement have revolutionized the treatment of myasthenia gravis (MG). However, a deeper understanding of complement modulation in the human system is required to improve treatment responses and identify off-target effects shaping long-term outcomes. For this reason, we studied a cohort of patients with MG treated with either eculizumab or azathioprine as well as treatment-naive patients using a combined proteomics and metabolomics approach. This strategy validated known effects of eculizumab on the terminal complement cascade. Beyond that, eculizumab modulated the serum proteometabolome as distinct pathways were altered in eculizumab-treated patients, including the oxidative stress response, mitogen-activated protein kinase signaling, and lipid metabolism with particular emphasis on arachidonic acid signaling. We detected reduced levels of arachidonate 5-lipoxygenase (ALOX5) and leukotriene A4 in eculizumab-treated patients. Mechanistically, ligation of the C5a receptor (C5aR) is needed for ALOX5 metabolism and generation of downstream leukotrienes. As eculizumab prevents cleavage of C5 into C5a, decreased engagement of C5aR may inhibit ALOX5-mediated synthesis of pro-inflammatory leukotrienes. These findings indicate distinct off-target effects induced by eculizumab, illuminating potential mechanisms of action that may be harnessed to improve treatment outcomes.

AKT activity orchestrates marginal zone B cell development in mice and humans.

The signals controlling marginal zone (MZ) and follicular (FO) B cell development remain incompletely understood. Here, we show that AKT orchestrates MZ B cell formation in mice and humans. Genetic models that increase AKT signaling in B cells or abolish its impact on FoxO transcription factors highlight the AKT-FoxO axis as an on-off switch for MZ B cell formation in mice. In humans, splenic immunoglobulin (Ig) D+CD27+ B cells, proposed as an MZ B cell equivalent, display higher AKT signaling than naive IgD+CD27- and memory IgD-CD27+ B cells and develop in an AKT-dependent manner from their precursors in vitro, underlining the conservation of this developmental pathway. Consistently, CD148 is identified as a receptor indicative of the level of AKT signaling in B cells, expressed at a higher level in MZ B cells than FO B cells in mice as well as humans.

Multiomic analysis of papillary thyroid cancers identifies BAIAP2L1-BRAF fusion and requirement of TRIM25, PDE5A and PKCδ for tumorigenesis.

BACKGROUND: Papillary thyroid carcinoma (PTC) is one of the most common forms of thyroid cancer with a cure rate of over 90% after surgery. However, aggressive forms may still occur, and personalized therapeutic strategies are increasingly required. METHODS: We performed integrated genomic and proteomic analysis of PTC tumor samples from patients who did not harbor BRAF or RAS mutations. We validate the analysis and present in-depth molecular analysis of the identified genetic rearrangement by employing biochemical and cell biological assays. Finally, we employ 3D spheroid models, loss of function studies and chemical inhibitors to target the hitherto upregulated factors. The data are analysed with appropriate statistical tests which are mentioned in the legends section. RESULTS: In a 23-year-old patient with thyroiditis, we identified a novel rearrangement leading to a BAIAP2L1-BRAF fusion that transforms immortalized human thyroid cells in a kinase and CC-domain dependent manner. Moreover, quantitative proteomic analysis of the same patient samples revealed the upregulation of several proteins including the Ubiquitin E3 ligase TRIM25, PDE5A, and PKCδ. Further, in a cohort of PTC patients, we observed higher expression of TRIM25 and PKCδ in the tumor and metastatic lesions, when compared to the matched normal tissue. Inhibition of TRIM25, PDE5A and PKCδ with small molecules or RNA interference affected not only viability and proliferation of BAIAP2L1-BRAF transformed cells, but also the viability, growth and invasion of corresponding 3D spheroid cultures. CONCLUSIONS: Apart from unveiling a novel oncogenic BRAF fusion in PTCs, our data may open a novel avenue of therapeutic targeting in human PTCs.

Plasmodium falciparum S-Adenosylmethionine Synthetase Is Essential for Parasite Survival through a Complex Interaction Network with Cytoplasmic and Nuclear Proteins.

S-adenosylmethionine synthetase (SAMS) is a key enzyme for the synthesis of the lone methyl donor S-adenosyl methionine (SAM), which is involved in transmethylation reactions and hence required for cellular processes such as DNA, RNA, and histone methylation, but also polyamine biosynthesis and proteostasis. In the human malaria parasite Plasmodium falciparum, PfSAMS is encoded by a single gene and has been suggested to be crucial for malaria pathogenesis and transmission; however, to date, PfSAMS has not been fully characterized. To gain deeper insight into the function of PfSAMS, we generated a conditional gene knockdown (KD) using the glmS ribozyme system. We show that PfSAMS localizes to the cytoplasm and the nucleus of blood-stage parasites. PfSAMS-KD results in reduced histone methylation and leads to impaired intraerythrocytic growth and gametocyte development. To further determine the interaction network of PfSAMS, we performed a proximity-dependent biotin identification analysis. We identified a complex network of 1114 proteins involved in biological processes such as cell cycle control and DNA replication, or transcription, but also in phosphatidylcholine and polyamine biosynthesis and proteasome regulation. Our findings highlight the diverse roles of PfSAMS during intraerythrocytic growth and sexual stage development and emphasize that PfSAMS is a potential drug target.

Gamma Irradiation Triggers Immune Escape in Glioma-Propagating Cells.

Glioblastoma multiforme is the most common and devastating form of brain tumor for which only palliative radio- and chemotherapy exists. Although some clinical studies on vaccination approaches have shown promising efficacy due to their potential to generate long-term immune surveillance against cancer cells, the evasion mechanisms preventing therapy response are largely uncharacterized. Here, we studied the response of glioblastoma-propagating cells (GPCs) to clinically relevant doses of γ radiation. GPCs were treated with 2.5 Gy of γ radiation in seven consecutive cellular passages to select for GPCs with increased colony-forming properties and intrinsic or radiation-induced resistance (rsGPCs). Quantitative proteomic analysis of the cellular signaling platforms of the detergent-resistant membranes (lipid rafts) in GPCs vs. rsGPCs revealed a downregulation of the MHC class I antigen-processing and -presentation machinery. Importantly, the radio-selected GPCs showed reduced susceptibility towards cytotoxic CD8+ T-cell-mediated killing. While previous studies suggested that high-dose irradiation results in enhanced antigen presentation, we demonstrated that clinically relevant sub-lethal fractionated irradiation results in reduced expression of components of the MHC class I antigen-processing and -presentation pathway leading to immune escape.

The caspase-2 substrate p54nrb exhibits a multifaceted role in tumor cell death susceptibility via gene regulatory functions.

Caspase-2 represents an evolutionary conserved caspase, which plays a role in genotoxic stress-induced apoptosis, ageing-related metabolic changes, and in deleting aneuploid cells in tumors. Genetic deletion of caspase-2 leads to increased tumor susceptibility in vivo. The exact downstream signaling mechanism by which caspase-2 accomplishes its specific tumor suppressor functions is not clear. Caspase-2, uniquely among caspases, resides in the nucleus and other cellular compartments. In this study, we identify a nuclear caspase-2 specific substrate, p54nrb, which is selectively cleaved by caspase-2 at D422, leading to disruption of the C-terminal site, the putative DNA binding region of the protein. P54nrb is an RNA and DNA binding protein, which plays a role in RNA editing, transport, and transcriptional regulation of genes. Overexpression of p54nrb is observed in several human tumor types, such as cervix adenocarcinoma, melanoma, and colon carcinoma. In contrast, the loss of p54nrb in tumor cell lines leads to increased cell death susceptibility and striking decrease in tumorigenic potential. By employing high resolution quantitative proteomics, we demonstrate that the loss/cleavage of p54nrb results in altered expression of oncogenic genes, among which the downregulation of the tumorigenic protease cathepsin-Z and the anti-apoptotic gelsolin can be detected universally across three tumor cell types, including adenocarcinoma, melanoma and colon carcinoma. Finally, we demonstrate that p54nrb interacts with cathepsin-Z and gelsolin DNA, but not RNA. Taken together, this study uncovers a so far not understood mechanism of caspase-2 tumor suppressor function in human tumor cells.

Fluorovinylsulfones and -Sulfonates as Potent Covalent Reversible Inhibitors of the Trypanosomal Cysteine Protease Rhodesain: Structure-Activity Relationship, Inhibition Mechanism, Metabolism, and In Vivo Studies.

Rhodesain is a major cysteine protease of Trypanosoma brucei rhodesiense, a pathogen causing Human African Trypanosomiasis, and a validated drug target. Recently, we reported the development of α-halovinylsulfones as a new class of covalent reversible cysteine protease inhibitors. Here, α-fluorovinylsulfones/-sulfonates were optimized for rhodesain based on molecular modeling approaches. 2d, the most potent and selective inhibitor in the series, shows a single-digit nanomolar affinity and high selectivity toward mammalian cathepsins B and L. Enzymatic dilution assays and MS experiments indicate that 2d is a slow-tight binder (Ki = 3 nM). Furthermore, the nonfluorinated 2d-(H) shows favorable metabolism and biodistribution by accumulation in mice brain tissue after intraperitoneal and oral administration. The highest antitrypanosomal activity was observed for inhibitors with an N-terminal 2,3-dihydrobenzo[b][1,4]dioxine group and a 4-Me-Phe residue in P2 (2e/4e) with nanomolar EC50 values (0.14/0.80 µM). The different mechanisms of reversible and irreversible inhibitors were explained using QM/MM calculations and MD simulations.

Naphthoquinones as Covalent Reversible Inhibitors of Cysteine Proteases-Studies on Inhibition Mechanism and Kinetics.

The facile synthesis and detailed investigation of a class of highly potent protease inhibitors based on 1,4-naphthoquinones with a dipeptidic recognition motif (HN-l-Phe-l-Leu-OR) in the 2-position and an electron-withdrawing group (EWG) in the 3-position is presented. One of the compound representatives, namely the acid with EWG = CN and with R = H proved to be a highly potent rhodesain inhibitor with nanomolar affinity. The respective benzyl ester (R = Bn) was found to be hydrolyzed by the target enzyme itself yielding the free acid. Detailed kinetic and mass spectrometry studies revealed a reversible covalent binding mode. Theoretical calculations with different density functionals (DFT) as well as wavefunction-based approaches were performed to elucidate the mode of action.

Proteogenomics analysis unveils a TFG-RET gene fusion and druggable targets in papillary thyroid carcinomas.

Papillary thyroid cancer (PTC) is the most common type of endocrine malignancy. By RNA-seq analysis, we identify a RET rearrangement in the tumour material of a patient who does not harbour any known RAS or BRAF mutations. This new gene fusion involves exons 1-4 from the 5' end of the Trk fused Gene (TFG) fused to the 3' end of RET tyrosine kinase leading to a TFG-RET fusion which transforms immortalized human thyroid cells in a kinase-dependent manner. TFG-RET oligomerises in a PB1 domain-dependent manner and oligomerisation of TFG-RET is required for oncogenic transformation. Quantitative proteomic analysis reveals the upregulation of E3 Ubiquitin ligase HUWE1 and DUBs like USP9X and UBP7 in both tumor and metastatic lesions, which is further confirmed in additional patients. Expression of TFG-RET leads to the upregulation of HUWE1 and inhibition of HUWE1 significantly reduces RET-mediated oncogenesis.

New Cysteine Protease Inhibitors: Electrophilic (Het)arenes and Unexpected Prodrug Identification for the Trypanosoma Protease Rhodesain.

Electrophilic (het)arenes can undergo reactions with nucleophiles yielding π- or Meisenheimer (σ-) complexes or the products of the SNAr addition/elimination reactions. Such building blocks have only rarely been employed for the design of enzyme inhibitors. Herein, we demonstrate the combination of a peptidic recognition sequence with such electrophilic (het)arenes to generate highly active inhibitors of disease-relevant proteases. We further elucidate an unexpected mode of action for the trypanosomal protease rhodesain using NMR spectroscopy and mass spectrometry, enzyme kinetics and various types of simulations. After hydrolysis of an ester function in the recognition sequence of a weakly active prodrug inhibitor, the liberated carboxylic acid represents a highly potent inhibitor of rhodesain (Ki = 4.0 nM). The simulations indicate that, after the cleavage of the ester, the carboxylic acid leaves the active site and re-binds to the enzyme in an orientation that allows the formation of a very stable π-complex between the catalytic dyad (Cys-25/His-162) of rhodesain and the electrophilic aromatic moiety. The reversible inhibition mode results because the SNAr reaction, which is found in an alkaline solvent containing a low molecular weight thiol, is hindered within the enzyme due to the presence of the positively charged imidazolium ring of His-162. Comparisons between measured and calculated NMR shifts support this interpretation.

Asymmetric Disulfanylbenzamides as Irreversible and Selective Inhibitors of Staphylococcus aureus Sortase A.

Staphylococcus aureus is one of the most frequent causes of nosocomial and community-acquired infections, with drug-resistant strains being responsible for tens of thousands of deaths per year. S. aureus sortase A inhibitors are designed to interfere with virulence determinants. We have identified disulfanylbenzamides as a new class of potent inhibitors against sortase A that act by covalent modification of the active-site cysteine. A broad series of derivatives were synthesized to derive structure-activity relationships (SAR). In vitro and in silico methods allowed the experimentally observed binding affinities and selectivities to be rationalized. The most active compounds were found to have single-digit micromolar Ki values and caused up to a 66 % reduction of S. aureus fibrinogen attachment at an effective inhibitor concentration of 10 µM. This new molecule class exhibited minimal cytotoxicity, low bacterial growth inhibition and impaired sortase-mediated adherence of S. aureus cells.

Structural and mechanistic insights into the interaction of the circadian transcription factor BMAL1 with the KIX domain of the CREB-binding protein.

The mammalian CLOCK:BMAL1 transcription factor complex and its coactivators CREB-binding protein (CBP)/p300 and mixed-lineage leukemia 1 (MLL1) critically regulate circadian transcription and chromatin modification. Circadian oscillations are regulated by interactions of BMAL1's C-terminal transactivation domain (TAD) with the KIX domain of CBP/p300 (activating) and with the clock protein CRY1 (repressing) as well as by the BMAL1 G-region preceding the TAD. Circadian acetylation of Lys537 within the G-region enhances repressive BMAL1-TAD-CRY1 interactions. Here, we characterized the interaction of the CBP-KIX domain with BMAL1 proteins, including the BMAL1-TAD, parts of the G-region, and Lys537 Tethering the small compound 1-10 in the MLL-binding pocket of the CBP-KIX domain weakened BMAL1 binding, and MLL1-bound KIX did not form a ternary complex with BMAL1, indicating that the MLL-binding pocket is important for KIX-BMAL1 interactions. Small-angle X-ray scattering (SAXS) models of BMAL1 and BMAL1:KIX complexes revealed that the N-terminal BMAL1 G-region including Lys537 forms elongated extensions emerging from the bulkier BMAL1-TAD:KIX core complex. Fitting high-resolution KIX domain structures into the SAXS-derived envelopes suggested that the G-region emerges near the MLL-binding pocket, further supporting a role of this pocket in BMAL1 binding. Additionally, mutations in the second CREB-pKID/c-Myb-binding pocket of the KIX domain moderately impacted BMAL1 binding. The BMAL1(K537Q) mutation mimicking Lys537 acetylation, however, did not affect the KIX-binding affinity, in contrast to its enhancing effect on CRY1 binding. Our results significantly advance the mechanistic understanding of the protein interaction networks controlling CLOCK:BMAL1- and CBP-dependent gene regulation in the mammalian circadian clock.

Enhancing Sensitivity of Microflow-Based Bottom-Up Proteomics through Postcolumn Solvent Addition.

The introduction of more sensitive mass spectrometers allows researchers to adapt front-end liquid chromatography (LC) to individual needs for the analysis of complex proteomes. Where absolute sensitivity is not paramount, it is advantageous to switch from a highly sensitive nanoflow-LC setup, the de facto standard platform in mass-spectrometry (MS)-based discovery proteomics, to a more robust, high-throughput-compatible microflow or conventional-flow setup. To enhance the microflow-LC-MS electrospray process of complex proteomic samples, we tested the effects of different solvents, including 2-propanol, methanol, and acetonitrile, pure or as mixture with dimethyl sulfoxide, which were added postcolumn to the eluting sample. Postcolumn addition of organic solvents strongly enhanced the electrospray efficiency in microflow-LC-MS experiments and improved the sensitivity across the entire gradient and for early eluting peptides by up to 10-fold. Postcolumn solvent addition did not negatively affect chromatographic performance and resulted in an overall 28-36% increase in identifications at both the protein and peptide levels. The presented microflow-LC-MS workflow, including postcolumn solvent addition, can be easily adopted on any LC-MS/MS platform.

The role of TCF3 as potential master regulator in blastemal Wilms tumors.

Wilms tumors are the most common type of pediatric kidney tumors. While the overall prognosis for patients is favorable, especially tumors that exhibit a blastemal subtype after preoperative chemotherapy have a poor prognosis. For an improved risk assessment and therapy stratification, it is essential to identify the driving factors that are distinctive for this aggressive subtype. In our study, we compared gene expression profiles of 33 tumor biopsies (17 blastemal and 16 other tumors) after neoadjuvant chemotherapy. The analysis of this dataset using the Regulator Gene Association Enrichment algorithm successfully identified several biomarkers and associated molecular mechanisms that distinguish between blastemal and nonblastemal Wilms tumors. Specifically, regulators involved in embryonic development and epigenetic processes like chromatin remodeling and histone modification play an essential role in blastemal tumors. In this context, we especially identified TCF3 as the central regulatory element. Furthermore, the comparison of ChIP-Seq data of Wilms tumor cell cultures from a blastemal mouse xenograft and a stromal tumor provided further evidence that the chromatin states of blastemal cells share characteristics with embryonic stem cells that are not present in the stromal tumor cell line. These stem-cell like characteristics could potentially add to the increased malignancy and chemoresistance of the blastemal subtype. Along with TCF3, we detected several additional biomarkers that are distinctive for blastemal Wilms tumors after neoadjuvant chemotherapy and that may provide leads for new therapeutic regimens.

NF-κB inducing kinase (NIK) is an essential post-transcriptional regulator of T-cell activation affecting F-actin dynamics and TCR signaling.

NF-κB inducing kinase (NIK) is the key protein of the non-canonical NF-κB pathway and is important for the development of lymph nodes and other secondary immune organs. We elucidated the specific role of NIK in T cells using T-cell specific NIK-deficient (NIKΔT) mice. Despite showing normal development of lymphoid organs, NIKΔT mice were resistant to induction of CNS autoimmunity. T cells from NIKΔT mice were deficient in late priming, failed to up-regulate T-bet and to transmigrate into the CNS. Proteomic analysis of activated NIK-/- T cells showed de-regulated expression of proteins involved in the formation of the immunological synapse: in particular, proteins involved in cytoskeleton dynamics. In line with this we found that NIK-deficient T cells were hampered in phosphorylation of Zap70, LAT, AKT, ERK1/2 and PLCγ upon TCR engagement. Hence, our data disclose a hitherto unknown function of NIK in T-cell priming and differentiation.

Chronic intestinal inflammation in mice expressing viral Flip in epithelial cells.

Viruses are present in the intestinal microflora and are currently discussed as a potential causative mechanism for the development of inflammatory bowel disease. A number of viruses, such as Human Herpesvirus-8, express homologs to cellular FLIPs, which are major contributors for the regulation of epithelial cell death. In this study we analyzed the consequences of constitutive expression of HHV8-viral FLIP in intestinal epithelial cells (IECs) in mice. Surprisingly, expression of vFlip disrupts tissue homeostasis and induces severe intestinal inflammation. Moreover vFlipIEC-tg mice showed reduced Paneth cell numbers, associated with excessive necrotic cell death. On a molecular level vFlip expression altered classical and alternative NFκB activation. Blocking of alternative NFκB signaling by deletion of Ikka in vivo largely protected mice from inflammation and Paneth cell loss induced by vFLIP. Collectively, our data provide functional evidence that expression of a single viral protein in IECs can be sufficient to disrupt epithelial homeostasis and to initiate chronic intestinal inflammation.

REGGAE: a novel approach for the identification of key transcriptional regulators.

Motivation: Transcriptional regulators play a major role in most biological processes. Alterations in their activities are associated with a variety of diseases and in particular with tumor development and progression. Hence, it is important to assess the effects of deregulated regulators on pathological processes. Results: Here, we present REGulator-Gene Association Enrichment (REGGAE), a novel method for the identification of key transcriptional regulators that have a significant effect on the expression of a given set of genes, e.g. genes that are differentially expressed between two sample groups. REGGAE uses a Kolmogorov-Smirnov-like test statistic that implicitly combines associations between regulators and their target genes with an enrichment approach to prioritize the influence of transcriptional regulators. We evaluated our method in two different application scenarios, which demonstrate that REGGAE is well suited for uncovering the influence of transcriptional regulators and is a valuable tool for the elucidation of complex regulatory mechanisms. Availability and implementation: REGGAE is freely available at https://regulatortrail.bioinf.uni-sb.de. Supplementary information: Supplementary data are available at Bioinformatics online.