The mycotoxin Beauvericin is an uncompetitive inhibitor of Cathepsin B

The mycotoxin Beauvericin is an uncompetitive inhibitor of Cathepsin B | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results The mycotoxin Beauvericin is an uncompetitive inhibitor of Cathepsin B Xiaoli Yang , Pablo Cea-Medina , Mohanraj Gopalswamy , Aparna Vaidya , Sonja Schavier , Shixin Oltzen , Sofie Moßner , Anfei Huang , Jing Qi , Doreen Floss , Markus Uhrberg , Holger Gohlke , Stefanie Scheu doi: https://doi.org/10.1101/2025.01.17.633572 Xiaoli Yang 1 Institute of Medical Microbiology and Hospital Hygiene, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Pablo Cea-Medina 2 Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Mohanraj Gopalswamy 2 Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Aparna Vaidya 1 Institute of Medical Microbiology and Hospital Hygiene, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sonja Schavier 1 Institute of Medical Microbiology and Hospital Hygiene, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Shixin Oltzen 3 Institute of Biochemistry and Molecular Biology II, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sofie Moßner 3 Institute of Biochemistry and Molecular Biology II, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Anfei Huang 5 Würzburg Institute of Systems Immunology, University of Würzburg , D-97078 Würzburg, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jing Qi 4 Institute for Transplantation Diagnostics and Cell Therapeutics, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Doreen Floss 3 Institute of Biochemistry and Molecular Biology II, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Markus Uhrberg 4 Institute for Transplantation Diagnostics and Cell Therapeutics, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Holger Gohlke 2 Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany 6 Institute of Bio- and Geosciences (IBG-4: Bioinformatics) , Forschungszentrum Jülich, D-52425 Jülich, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: stefanie.scheu{at}med.uni-rostock.de gohlke{at}uni-duesseldorf.de Stefanie Scheu 1 Institute of Medical Microbiology and Hospital Hygiene, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf , D-40225 Düsseldorf, Germany 7 Institute of Immunology, Rostock University Medical Center , D-18057 Rostock, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: stefanie.scheu{at}med.uni-rostock.de gohlke{at}uni-duesseldorf.de Abstract Full Text Info/History Metrics Preview PDF Abstract Beauvericin (BEA), a cyclic depsipeptide, is a mycotoxin of the enniatin family and the secondary metabolite of various toxigenic fungi. Multiple biological functions of BEA have been well investigated, such as anti-cancer, anti-inflammatory, anti-microbial, and immune-activating functions. In a recent study, we showed that BEA can target Toll-like receptor 4 (TLR4) to induce dendritic cell (DC) activation. In an in silico screen, we identified Cathepsin B (CTSB) as a potential additional interaction partner for BEA, which has been verified recently in a study showing inhibition of CTSB activity by BEA in cell-free assays. The underlying molecular mechanism of BEA-mediated CTSB inhibition remains unknown as do the cellular entities where this inhibition takes place. In this study, we determine the effects of BEA on CTSB within granulocyte-macrophage colony-stimulating factor (GM-CSF)-cultured bone marrow-derived dendritic cells (BMDCs) and human leukemia monocytic cell line THP-1 induced dendritic cells (iDCs). BEA significantly suppresses CTSB activity in both mouse BMDCs and human iDCs. NMR analyses indicate that BEA directly interacts with CTSB. Enzyme kinetics show that BEA can directly inhibit CTSB activity and acts as an uncompetitive inhibitor. Molecular docking analysis revealed a putative binding site for BEA in human CTSB. Collectively, we describe the molecular mechanisms underlying the biological activity of BEA against human CTSB, suggesting that CTSB may be a candidate target for tumor therapy. Introduction CTSB is a member of the lysosomal cysteine protease family, which also includes cathepsin C, F, H, K, L, O, S, V, X, and W ( 1 , 2 ). Structurally, CTSB has the regular V-shaped active site cleft of papain-like enzymes formed by the folding of two distinct domains. In addition, CTSB has a 20-residue insertion, termed occluding loop, that overlaps with the active site, binds the C-terminus of the protein, and induces carboxypeptidase activity ( 1 , 3 ). Unlike other members of the papain family, CTSB has both exopeptidase and endopeptidase activity. Its exopeptidase activity is dependent on the stability of the occluding loop, which is controlled by acidic pH. By contrast, the endopeptidase activity of CTSB is determined by neutral/alkaline pH, which can disrupt the salt bridges between the occluding loop and lead to the exposure of the active site ( 3 , 4 ). CTSB has been shown to play an important role in intracellular proteolysis and remodeling of the extracellular matrix of cellular proteolysis networks ( 3 , 4 ). Under normal physiological conditions, CTSB activity is well controlled and has been implicated in several oncogenic processes such as bone resorption, antigen processing, and protein turnover ( 5 ). However, aberrant expression of the cysteine protease CTSB is observed in the pathogenesis of many diseases, including neurological ( 6 – 8 ), immune ( 9 – 11 ), and cardiac diseases ( 12 , 13 ) as well as cancer ( 14 , 15 ). Numerous studies have shown that CTSB overexpression is correlated with invasive and metastatic cancers because it supports a pro-malignant phenotype by enabling the maintenance of active proteases through active suppression of their inhibitors ( 16 – 18 ). For example, CTSB enhances the activity of matrix metalloproteinases (MMPs) by destroying their inhibitors (e.g., TIMP1 and TIMP2) in human articular chondrocytes and maintaining high levels of MMPs, thereby promoting extracellular matrix (ECM) degradation and angiogenesis ( 19 ). Therefore, targeting CTSB could have significant clinical relevance in the treatment of various cancers. BEA is a cyclic hexadepsipeptide ( 20 ), composed of three units of phenylalanine (Phe) and three units of 2-hydroxyisovaleric acid (Hiv), and is produced by various fungi, such as Beaveria bassiana and Fusarium spp ( 21 , 22 ). The bioactivity of BEA has been studied in the context of anti-cancer ( 23 , 24 ), anti-virus ( 25 , 26 ), anti-inflammatory ( 27 , 28 ), antimicrobial ( 29 – 32 ), immune-regulating ( 33 ), insecticidal ( 34 , 35 ), and pesticidal activity ( 36 , 37 ). We as well as others have identified CTSB as a molecular target of BEA ( 38 ). However, the mode of action of BEA binding to CTSB as well as if this also takes place in treatment-relevant immune cells remains unclear. In this study, BEA is shown to be an uncompetitive inhibitor of CTSB in mouse and human DCs. In addition, the binding site of BEA in CTSB was revealed using a molecular docking approach. Materials and Methods Mice Cells from the bone marrow of C57BL/6N mice were used for GM-CSF culture. No experiments on live animals were performed. Mice were euthanized by cervical dislocation before bone marrow was harvested. The euthanasia method used is in strict accordance with accepted norms of veterinary best practice. Animals were kept under specific pathogen-free conditions in the animal research facility of the Heinrich Heine University of Düsseldorf strictly according to German animal welfare guidelines. BMDC culture and stimulation conditions 2x10 6 mouse bone marrow cells were cultured in non-treated 94x16 mm dishes (Sarstedt, Cat# 82.1473) in 10 ml VLE DMEM (Biochrom, Cat# P04-04515) containing 10% heat-inactivated FCS (Sigma-Aldrich, Cat# F7524), 0.1% 2-mercaptoethanol (Thermo Fischer Scientific, Cat# 31350-010), and GM-CSF and kept for 9 days. GM-CSF cultures were performed as previously described ( 39 ). 10 ml GM-CSF containing medium was added to the plates on day 3. On day 6, 10 ml medium was carefully removed and centrifuged. The cell pellet was resuspended in 10 ml medium and added to the dish. On day 9, BMDCs were seeded on a 24-well plate and subjected to a specified BEA concentration or to 30 μM of the control inhibitor, CA074 (MedChemExpress, Cat# HY-103350), for 24 hours. The stimulated cells were collected for FACS or CTSB activity assays. iDC culture and stimulation conditions iDCs were generated according to the protocol ( 40 ), namely, 1x10 6 THP1 cells were cultured in a T25 flask with 5 ml RPMI (Thermo Fisher Scientific, Cat# 31870-025) including 10% heat-inactivated FCS (Sigma-Aldrich, Cat# F7524), 0.1% 2-mercaptoethanol (Thermo Fischer Scientific, Cat# 31350-010), 1% PenStrep (Biochrom GmbH, Cat#A2212), 1500 IU/ml rhGM-CSF (ImmunoTools, Cat#11343125), and 1500 IU/ml rhIL-4 (ImmunoTools, Cat#11340045) and incubated at 37°C and 5% CO 2 for 5 days. On day 3, 4 ml of the medium was exchanged for fresh medium including 1500 IU/ml rhGM-CSF (ImmunoTools, Cat#11343125) and 1500 IU/ml rhIL-4 (ImmunoTools, Cat#11340045). On day 5, cells are seeded on a 24-well plate and subjected to the specified concentration of BEA and CA074 (30 μM) (MedChemExpress, Cat# HY-103350) for 24 hours; the stimulated cells were collected for FACS or CTSB activity assays. Cell-based CTSB activity assays The stimulated cells were lysed with chilled cell lysis buffer (Abcam, Cat# ab65300), and the cell lysates were subjected to protein quantification by BCA. Cell CTSB activity was measured according to the manufacturer’s protocol (Abcam, Cat# ab65300). Briefly, equal amounts of protein in 50 µl reaction buffer were added to 96 black clear bottom-well plates (Greiner Bio-One, Cat# 675086), followed by the addition of 50 µl cathepsin B reaction buffer. The reaction was then incubated at 37°C for 1 hour with the addition of substrate and protected from light. The output was read on a fluorescence microplate reader at Ex/Em = 400/505 nm. Alternatively, stimulated cells were incubated with Rhodamine 110-(RR)2 substrate (Abcam, Cat# ab270787) at approximately 1/50 for 30 minutes at 37°C, protected from light. After incubation, the cells were stained with Hoechst 33342 (Abcam, Cat# ab270787) at approximately 0.5% for 10 minutes and then analyzed by FACS. STD-NMR analysis NMR experiments were performed at 298 K on Bruker Avance III HD spectrometers operating at 750 MHz, equipped with 5 mm triple resonance TCI ( 1 H, 13 C, 15 N) cryoprobes and shielded z-gradients. BEA compound was purchased (Item No. 11426) from Cayman Chemical, Ann Arbor, Michigan, USA, and recombinant human CTSB (6His tag at the C-terminus) was purchased (Item No. 153921.10) from Biomol GmbH, Hamburg, Germany. Both CTSB and BEA were used without any further purifications. The purity of the compound (>95%) was checked by 1D 1 H NMR of 250 µM BEA compound measured (256 scans) in 90% (v/v) DMSO d6 . ( Figure 2A ). BEA is sparingly soluble in aqueous buffers, so 0.1% (v/v) Tween 20 was used to increase the compound solubility. To perform STD-NMR experiments, 5 µM of CTSB and 250 µM BEA were prepared in 25 mM oxalic acid/sodium oxalate, 5% (v/v) DMSO d6 , pH 5.0, 10% (v/v) D 2 O, and 0.1% (v/v) Tween 20. The STD-NMR spectra were acquired with 4928 scans, and the subtraction of the reference spectrum was performed internally via phase cycling after every scan to obtain the STD effect. Selective saturation of protein resonances (on-resonance spectrum) was performed by irradiating at -0.3 ppm for a total saturation time of 3 s with a delay time (D1) of 5 s. There was no compound signal detected below 0 ppm. For the reference spectrum (off-resonance), the samples were irradiated at −40 ppm. The STD intensity of the aromatic proton STD effect was set to 100% as a reference, and the relative intensities were determined ( 41 ). The assignment of the protons of the ligand was achieved by the analysis of the 1D 1 H NMR and aided by 1 H chemical shift prediction performed by the ChemNMR package, contained within the software ChemDraw (PerkinElmer). Data was processed and analyzed with TopSpin 3.2 (Bruker BioSpin). Sodium 2,2-dimethyl-2-silapentane-5-sulfonate (DSS) was used for chemical shift referencing. Enzyme kinetics Recombinant human CTSB (R&D Biosystems, Cat# 953-CY) was activated in a solution of 25 mM MES pH 5 and 5 mM DTT at a protein concentration of 25 µg/ml, as indicated by the manufacturer. The activation was carried out by incubating the solution in a thermal bath at 26 °C for 30 minutes and then diluting to a concentration of 0.2 µg/ml in a 25 mM MES buffer pH 5. Reactions were carried out in a 96-Well black microplate (Greiner Bio-One, Cat# 675086) in a final volume of 100 µl, using 20 ng of enzyme with increasing concentrations of the fluorogenic substrate Z-Leu-Arg AMC (R&D Biosystems, Cat# ES008) (from 0 to 200 µM) and BEA (from 0 to 25 µM). Progress curves were measured using a TECAN infinite M plex microplate reader with an excitation wavelength of 380 nm and an emission wavelength of 460 nm. All the components, except the enzyme, were loaded onto the plate and were left to reach a temperature of 26 °C. Protein samples were loaded using a multichannel pipette. Plates were shaken using a pulse of 3.5 s with an amplitude of 4.5 mm. Data was collected for a total of 1400 s with a reading interval of 45 s. Reaction rates were calculated from the progress curves using linear regression. Kinetic parameters were derived by fitting a Michaelis-Menten model to a plot of reaction rates at different substrate concentrations as described in equation (1) where [ S ] is the substrate concentration, v is the observed reaction rate or velocity, V max is the maximal velocity of the enzyme, and K M is the Michaelis-Menten constant. To later derive an inhibition constant, we employed a general inhibitor model ( eq. (2) ), as described by Cornish-Bowden in ref. ( 42 ) where [ I ] is the concentration of the inhibitor, K i is the inhibition constant of the [EI] complex and K i ’ is the inhibition constant for the [ESI] complex. The equation can be reformulated as shown in equation (3) which allows us to derive the K i of competitive, non-competitive, and uncompetitive inhibitors. The interception point of curves at different [ S ] on the x-axis is: and the value on the y-axis for the intercept is: For uncompetitive inhibitors the value of K i → ∞, therefore the y-value of the intercept is: To approximate the observed interception point, we identified the (x, y)-values that minimize the sum of the squared distances perpendicular to each substrate curve i where a i , b i , and c i correspond to the coefficients of the curves at different [ S ], derived by fitting to equation (3) . All data manipulation, including model fitting, was performed using the Pandas, NumPy, and SciPy modules of Python 3.10. Modeling and co-solvent probe molecular dynamics simulations of human CTSB Given the high degree of conformational variability shown by human CTSB, we employed four crystal structures as starting points: PDB ID 1CSB, 3CBK, 1HUC, and 3K9M. To avoid biasing the results due to differences in the sequence length of each solved crystal structure, and to fix missing loops and residues, all the structures were remodeled using Modeller ( 43 ), using as a target the sequence of residues Phe74 to Ile339 of human CTSB (UniProtKB ID: P07858) and the respective crystal structure as a template. 30 models were generated for each system, and the best one was selected based on the DOPE score for further work. Each system was packed in a solvent box with either benzene or isopropanol as co-solvent using PACKMOL-Memgen ( 44 ). A ratio of 99 water molecules per co-solvent molecule (benzene or isopropanol) was utilized. Five different random starting configurations of co-solvent distribution were generated for each system. The most likely protonation state for titratable residues was assigned with Propka ( 45 ). Parameter files were generated using tleap in AmberTools ( 46 ), with ff19SB ( 47 ) as protein force field and the co-solvent parameters as distributed alongside PACKMOL-Memgen, TIP3P as water model ( 48 ) and K + and Cl - counter ions with parameters from Joung & Cheatham ( 49 ). All simulations were performed using pmemd.cuda on Amber22 ( 50 ), using a direct space cut-off of 10 Å for long-range electrostatic interactions. All bonds involving hydrogen atoms were constrained using the SHAKE algorithm ( 51 ). The systems were minimized in three consecutive rounds; first using positional restraint forces of 5 kcal mol -1 Å -2 , followed by restraint forces of 0.1 kcal mol -1 Å -2 , and lastly a restraint-free minimization. Each minimization round was carried out using 500 steps of steepest descent followed by 2000 steps of conjugate gradient. The systems were thermalized from 0 to 100 K in a window of 25 ps, using the Langevin thermostat with a collision coefficient of 1.0 ps -1 under NVT conditions. Afterward, the systems were heated to 300 K over a window of 250 ps under NPT conditions using the Berendsen barostat with isotropic position scaling. Both thermalization and pressure adaptation were carried out using positional restraint forces of 1.0 kcal mol -1 Å -2 . Five subsequent 50 ps-long NPT runs were carried out, reducing the positional restraint forces by 0.2 kcal mol -1 Å -2 at each step until reaching zero. Two production runs of 600 ns length were performed for each system using a time step of 2 fs and saving coordinates every 100 ps. This yielded a cumulative simulation time of 1.2 µs per starting configuration. Since five random initial solvent configurations were generated for each co-solvent, the total cumulative sampling time for each co-solvent probe was 6 µs. All simulations were analyzed using cpptraj ( 52 ). The protein was aligned to the first frame of each trajectory with the rmsd command, using the backbone atoms. The density of each co-solvent was calculated using the grid command with a spacing of 0.5 Å. Docking of BEA against the inactive conformation of human CTSB Molecular docking of BEA was performed against the inactive conformation of human CTSB (PDB ID 3PBH), with the propeptide region removed. The protein structure was prepared using the protein preparation wizard, with protonation states assigned by Propka, while the ligand was prepared with LigPrep (LigPrep, Schrödinger, LLC). A grid with dimensions of 35 Å was centered at the midpoint of residues Phe30 and Val33 of the removed propeptide. Due to the large size of the ligand, a rigid docking using a diverse conformational ensemble of BEA structures was used. The conformational ensemble was generated using Schrödinger’s macrocycle conformational sampling tool, using an RMSD cut-off of 0.75 Å and default options ( 53 ). 255 non-redundant conformations were generated. Rigid dockings were carried out using Glide in extra precision mode (XP) ( 54 ) using OPLS_2005 as a force field. Three poses were generated for each conformation of BEA. The resulting binding modes were analyzed using PyMOL. Statistical analysis The GraphPad Prism 9.0 software was used for data analysis. Data are represented as means ± SEM. For analyzing statistical significance between multiple groups, a one-way ANOVA with Dunnett’s multiple comparisons test was used. For analyzing statistical significance for comparisons of more than two groups with two or more stimulations, two-way ANOVA with Sidak’s multiple comparisons test was used. All p values < 0.05 were considered statistically significant. Results BEA suppresses CTSB activity in DCs Our previous study showed that BEA can target TLR4 to induce DC activation ( 33 ). In searching for alternative molecular targets of BEA in DCs, we performed in silico studies and found CTSB as a potential target of BEA. Recently, Silva et al. indicated that BEA potently inhibits human CTSB activity in cell-free assays ( 38 ). We wondered if BEA could inhibit CTSB in immune cells relevant for anti-cancer immunotherapy approaches such as DCs. To this end, BMDCs were generated by culturing bone marrow cells with GM-CSF followed by stimulation with BEA at various concentrations for 16 hours. First, the CTSB enzyme activity of BEA-stimulated BMDCs was determined in the cell lysates. CTSB activity was significantly suppressed by BEA stimulation already at 2.5 µM showing a dose-dependent inhibition of CTSB activity and suggesting that CTSB is a target of BEA in BMDCs ( Figure 1A ). CTSB activity was also detected in cells using a flow cytometry-based assay. Consistently, BEA can inhibit CTSB activity in BMDCs, and significant inhibition was observed starting at 7.5 µM BEA in this assay ( Figure 1B and 1C ). This led us to ask whether BEA could also inhibit human CTSB. To address this question, human monocyte-derived immature dendritic cells (iDCs) were induced from THP-1 cells in the presence of human GM-CSF and IL-4 and then subjected to various concentrations of BEA stimulation for 16 hours. As shown in Figure 1D , BEA can also suppress human CTSB activity with significant inhibition observed starting at a concentration 5 µM BEA. In addition, flow cytometry data also show that BEA significantly inhibits human CTSB activity ( Figure 1E and 1F ). Taken together, BEA inhibits CTSB activity in mouse BMDCs and human iDCs. Download figure Open in new tab Fig 1. Inhibition of CTSB activity by BEA in mouse BMDCs. (A) 1x10 6 BMDCs from C57BL/6N mice were stimulated with the indicated concentrations of BEA for 16 h, CTSB activity was measured in cell lysates. (B, C) 5x10 5 BMDCs from C57BL/6N mice were stimulated with the indicated concentrations of BEA or the known CTSB inhibitor CA-074 for 16 h, CTSB activity was detected in BMDCs by flow cytometry. (D) 1x10 6 iDCs were stimulated with the indicated concentrations of BEA for 16 h, CTSB activity was measured in cell lysates. (E, F) 5x10 5 iDCs were stimulated with the indicated concentrations of BEA or CA-074 for 16 h, CTSB activity was detected in iDCs by flow cytometry. Data is shown as means ± SEM. Results shown are representative of two to three independent experiments. * p < 0.05, ** p < 0.01, **** p < 0.0001, ns: not significant. STD NMR spectroscopy of BEA binding to CTSB To corroborate that the observed biological effects of BEA are mediated by direct interaction with CTSB, we performed saturation transfer difference (STD) NMR experiments using human CTSB. This technique allows identification of the chemical moieties of a ligand that interact directly with a target protein by measuring the shift in proton resonances that are in close proximity ( 41 ). The results show a strong increase in the signal of aromatic protons belonging to one benzene ring and the alkyl protons in one isopropyl group ( Figure 2 ), indicating that this region of the molecule is directly binding to human CTSB. Download figure Open in new tab Fig 2. STD NMR spectra of BEA in the presence of hu CTSB. Chemical structure of BEA along with proton assignment and NMR spectra. (A) 1D 1 H NMR of 250 µM BEA measured (256 scans) in 100 % (v/v) DMSO d6 . (B) Control experiments of only buffer (128 scans), containing 25 mM oxalic acid/sodium oxalate, 5% (v/v) DMSO d6 , pH 5.0, 10% (v/v) D 2 O, and 0.1% (v/v) Tween 20. (C) Reference 1D NMR spectrum (128 scans) and (D) STD NMR spectrum (4928 scans) of BEA in the presence of human CTSB, respectively. Protons making close contacts with the protein interaction site show significant enhancements in their resonances. The intensity of the STD effect of the aromatic protons was set to 100% as a reference, and relative intensities were determined. In particular, aromatic (f) and methyl (a, a’ - ’’) proton groups show strong signal enhancements. For all spectra, the scale for chemical shifts is shown at the bottom. The intensity of the spectra (scale of the y-axis) is adjusted for clarity. Enzyme kinetics of CTSB in the presence of BEA Since the binding of BEA could occur either via the active site of CTSB or an allosteric site, we performed a kinetic characterization of the enzyme at different concentrations of the inhibitor, with the aim of determining whether BEA acts as a competitive, non-competitive, or uncompetitive inhibitor. Human CTSB shows a classical Michaelis-Menten behavior in the absence of an inhibitor ( Figure 3A ). As the inhibitor concentration increases, the maximal enzyme velocity V max decreases while the K M value remains similar within the experimental uncertainty ( Table 1 ), suggesting that BEA does not act as a competitive inhibitor of the enzyme ( 55 ), and, therefore, does not bind directly to the active site of huCTSB. Download figure Open in new tab Fig 3. Kinetic characterization of BEA inhibition reveals an uncompetitive mechanism. (A) Michaelis-Menten curves for human CTSB in the presence of increasing concentrations of BEA. The lines show the model fit. (B) Dixon plot for the kinetic experiment shown in (A) . The lines show the linear fit. (C) Substrate over enzyme velocity plot for the kinetic experiment shown in (A) . The red dashed line shows the closest neighbor intersection point of the curves. The lines show the linear fit. In all plots, the error bars represent the SEM of a triplicate run. View this table: View inline View popup Download powerpoint Table 1. Kinetic parameters derived from Figure 3A . To clarify whether BEA acts via an uncompetitive or mixed inhibition mechanism, we employed the graphical discrimination method proposed by Cornish-Bowden ( 42 ), which consists of plotting the reciprocal of the velocity (1/v) versus the concentration of the inhibitor ([I]), also known as Dixon plot ( 56 ), and the substrate concentration over the velocity ([S]/v) versus the inhibitor concentration, at multiple values of the substrate concentration. For a purely non-competitive inhibitor, both 1/v vs [I] and [S]/v vs. [I] will give straight lines converging in the negative region of the x-axis with y-values close to zero. For a mixed inhibitor, the lines will also converge in the negative region of the x-axis, but with y ≠ 0. Finally, an uncompetitive inhibitor will show parallel lines for the Dixon plot whereas the curves will intersect in the fourth quadrant of the [S]/v vs. [I] plot. Figure 3B and 3C show the plots 1/v or [S]/v versus [I], respectively. In the first plot ( Figure 3B ), the curves do not intersect each other. In the second plot ( Figure 3C ), there are multiple intersections at x-values 0, which strongly suggests that BEA behaves as an uncompetitive inhibitor. The x-value of the intersection point between the curves corresponds to the negative value of the inhibition constant (- K i ). However, the experimental noise in the assays did not allow for the direct derivation of the K i value. Therefore, we approximated the closest neighbor intersection point from all the curves by finding the value that minimizes the total distance to the multiple intersecting points (shown as the cross between the red dashed lines in Figure 3C ). This results in an estimated K i value of 44.3 µM. From equation (2) , it can be deduced that the y-value of the intersection point of the curves is equal to K M / V max . As an internal check, we compared the experimental value of K M / V max in the absence of inhibitor (2843.73 min µg) to the y-value of the closest neighbor intersection point (3116.68 min µg), revealing a difference of 9.15%, which indicates that our approximation yields a K i of equal order of magnitude. Overall, the kinetic characterization of the inhibition of human CTSB reveals that BEA acts as an uncompetitive inhibitor. Structural modelling of CTSB and BEA interactions An uncompetitive inhibition mode implies that the inhibitor binds to an already established enzyme-substrate complex, ruling out the possibility of BEA binding directly to the active site of human CTSB. Apart from the active site, a regulatory site involved in the heparin-mediated activation has been described ( 57 ), but no inhibitory allosteric sites have been reported. To identify alternative putative binding sites for BEA, we carried out molecular dynamics simulations (MD) using multiple starting structures of human CTSB in the presence of two co-solvent probes that mimic the key chemical moieties of BEA found to be in direct contact with the enzyme by STD-NMR, benzene and isopropanol. Performing MD simulations in the presence of small molecular probes has been shown to be a powerful method to identify previously undescribed ligand binding sites ( 58 , 59 ), as smaller probes tend to accumulate in the protein pockets that have the molecular features necessary for the binding of larger ligands. A large loop adjacent to the active site of CTSB, known as the occluding loop, displays a high degree of conformational variability ( 60 ). To avoid biasing our results by the starting conformations of the occluding loop, we performed the co-solvent MD simulations using four different starting conformations modeled from different crystal structures. Figure 4 shows the regions where benzene (orange densities) and isopropanol (cyan densities) accumulate preferentially during the simulation for each one of the starting conformations of human CTSB. Besides small patches of interacting spots on the surface of the protein, for all the starting structures, the largest co-localization of probes occurs adjacent to the occluding loop near the active site. Download figure Open in new tab Fig 4. Co-solvent MD simulations reveal a common accumulation site adjacent to the catalytic site. Molecular dynamics simulations were started using four different initial conformations of human CTSB, each labelled according to the PDB structure that was used as template. The accumulation of co-solvent probes is shown as blue and orange densities for isopropanol and benzene, respectively. The volumes are shown in the same scale across panels. The largest accumulations are highlighted with arrow heads. The catalytic cysteine is shown as purple stick and highlighted with a green circle. Next, we investigated if other known inhibitors of cathepsin could target the edge of the substrate binding site cleft adjacent to the occluding loop, as suggested by the MD simulations. Most of the small molecule inhibitors of human CTSB are covalent inhibitors that directly target the catalytic serine in the active site. Nonetheless, human CTSB as well as other proteases are expressed as autoinhibited zymogens that undergo activation under certain pH and redox conditions ( 61 ). The zymogen form includes a large N-terminal loop that covers the entirety of the substrate binding cleft and inhibits enzyme activity. Notably, the zymogen form of human CTSB (PDB ID 3PBH) contains an autoinhibitory peptide with a hot spot of Phe and Val residues interacting with the edge of the substrate binding cleft near the occluding loop, very similar to the interaction site we identified through MD simulations ( Figure 5A ). Previous work has reported that BEA can inhibit other papain-like enzymes such as papain and cathepsin V ( 62 ). Therefore, we assessed whether the autoinhibitory peptides of these enzymes also display a similar interaction motif in the same spot. Figure 5A shows that the autoinhibitory peptides of papain and cathepsin V have a Phe-Val-rich motif interacting directly with the edge of the substrate binding cleft, in the same position as observed for human CTSB. Download figure Open in new tab Fig 5. Identification of a putative BEA binding site using molecular docking. (A) Different proteins inhibited by BEA (Cathepsin B: blue, Cathepsin V: white, Papain: yellow) display Phe-Val-rich regions in their proenzyme autoinhibitory peptides (Cathepsin B: beige, Cathepsin V: purple, Papain: cyan). The residues are shown as red sticks, the region is highlighted with a grey circle. (B) Docking score distribution for the generated binding poses of BEA. (C) Binding mode consistent with the placement of the autoinhibitory peptide. BEA is shown as yellow sticks, nearby residues are shown as blue sticks, the catalytic cysteine is shown in cyan. The center of mass of the side chains of the autoinhibitory peptide are shown as red spheres. (D) Superimposition of the binding mode shown in C (protein shown as blue cartoon, BEA shown as yellow van der Waals spheres), against the structure of rat CTSB bound to the tripeptide-mimetic covalent inhibitor benzyloxycarbonyl-Arg-Ser(O-Bzl) chloromethylketone (PDB ID 1THE) (protein shown as green cartoon, inhibitor shown as purple van der Waals spheres). Finally, we tested whether BEA could interact with human CTSB without overlapping directly with the catalytic site. For this, we performed molecular docking focusing the search of putative binding poses in the region where the Phe-Val hot spot of the autoinhibitory peptide is located. Using the autoinhibited conformation of human CTSB without the N-terminal inhibitory peptide as a target structure, we generated 141 binding poses using Glide SP ( 54 ). Overall, the docking scores of the generated poses are low, suggesting weak binding ( Figure 5B ). Within the top five poses with the highest score is one binding mode in which BEA is placed such that its benzene ring and isopropyl group occupy the same region as the Phe-Val hot spot of the autoinhibitory peptide and are far apart from the catalytic cysteine residue ( Figure 5C ). To verify if BEA could bind to this region without obstructing access to the catalytic center of the enzyme, we superimposed the predicted binding mode with the crystal structure of a covalent peptidomimetic inhibitor ( Figure 5D ). As evidenced by the van der Waals radius of the atoms of each ligand, the volume of both molecules does not overlap, suggesting that BEA could bind to human CTSB in a manner consistent with the uncompetitive inhibition mechanism. Discussion For BEA, various biological functions have been reported including anti-cancer ( 24 , 63 – 66 ), anti-inflammatory ( 27 , 28 ), anti-viral ( 25 , 26 ), and anti-microbial properties ( 67 ). Furthermore, an immune-activating function of BEA was identified in that BEA activates BMDCs via targeting TLR4, inducing inflammatory cytokine production ( 33 ). In silico prediction tools showed that CTSB is a potential target of BEA. Recently, Silva et al. indicated that BEA acts as a potent inhibitor against human CTSB in cell-free enzymatic activity assays. Here, we revealed that CTSB is active in BMDCs and that BEA significantly inhibits mouse CTSB activity at 2.5 μM in BMDCs after 16 h of stimulation. A similar inhibition of human CTSB activity by BEA was observed in human iDCs, and a significant inhibition was also observed at 2.5 μM. Moreover, through STD-NMR and enzyme kinetics characterizations, we established that BEA binds to CTSB and inhibits it, without displacing the substrate from the active site. Finally, computational structural analysis of CTSB revealed a putative binding site for BEA, where instead of acting as a substrate mimic, it mimics a motif found within the inhibitory propeptide region. CTSB is the first and currently best-characterized member of the C1 family of papain-like, lysosomal cysteine peptidases ( 68 ). It is widely expressed in most cell and tissue types, and is mainly localized in endo/lysosomes and can be released into cytosol in inflammatory diseases ( 69 ). It was shown that elavated expressions of CTSB are observed in a variety of human diseases, such as cancer, neurodegeneration and immunological and developmental disorders. Increasing evidence of the pathophysiological roles and substrates of CTSB was revealed by the establishment of a CTSB knockout mouse ( 70 ). Gonzalez-Leal et al. showed that CTSB deficiency can lead to increased MHC-II molecule expression and can induce IL-12 production by LPS-treated BMDCs ( 71 ). This suggested that BEA might induce increased IL-12 production after inhibition of CTSB activity or deletion of CTSB in BMDCs. Further studies need to be done to verify whether BEA can induce more IL-12 production by BEA-treated BMDCs after inhibition of CTSB activity or deletion of CTSB. We used STD-NMR to map the BEA epitope binding to human CTSB. The results showed that human CTSB interacts through the exposed side chains of phenylalanine and valine residues. Given that these residues are identical to amino acids of the peptide substrate known to mediate the interaction with the active site of human CTSB ( 75 ), one might speculate that BEA acts as a competitive inhibitor. However, the kinetic characterization clearly showed a decrease of V max , hence, discarding the possibility of BEA being a pure competitive inhibitor. Furthermore, graphical analysis via the Dixon plot and the [S]/v vs. [I] plot revealed that BEA acts as an uncompetitive inhibitor under the experimental conditions, which prompted the search for an alternative binding site. Although this inhibition mechanism is rare, it is sometimes explained by supposing that the inhibitor can bind the ES complex but not the free enzyme ( 76 ). MD simulations in the presence of co-solvent probes did not reveal a novel marked allosteric site. Instead, the probes accumulated at the edge of the substrate binding cleft, which does not not contradict the proposed mechanism. Since CTSB has a large substrate binding cleft delimited by flexible loops that can accommodate a wide variety of substrates and inhibitors ( 60 , 75 ), it is conceivable that human CTSB can bind a BEA molecule within the substrate binding cleft, while there is still enough space in the catalytic center to allow for the peptidase activity to take place. To further support the idea that the edge region of the binding cleft of human CTSB effectively interacts with BEA, we analyzed the contacts that the autoinhibitory peptide of the proenzymatic form of human CTSB forms with the rest of the structure. Since the autoinhibitory peptide contains pharmacophoric moieties required to abort human CTSB function, we hypothesized that there is a common binding subregion for both. Structural studies of the proenzymatic form of human CTSB showed that there is indeed a Phe-Val-rich region in the autoinhibitory peptide. Furthermore, by comparing different structures, we showed that the presence of this Phe-Val rich hot spot is common to the enzymes inhibited by BEA. To test whether the identified site could accommodate BEA, we performed molecular docking focusing pose generation onto this region. Our results yielded a single pose, among the top-ranked ones, that is consistent with the space occupied by the autoinhibitory peptide. Furthermore, the observed STD-NMR signal shows that the interaction between BEA and CTSB occurs through the Phe and Val side chains of the peptide, which is consistent with the proposed binding mode. Considering that BEA acts as an uncompetitive inhibitor, conformational rearrangements induced by the presence of the substrate in the active site are required for the binding of BEA ( 75 ). Moreover, it has been shown that different inhibitors can trigger different conformational changes in CTSB due to its high flexibility ( 60 ). In our docking study, we used the autoinhibited conformation of human CTSB as target structure, as it might have the occluding loop in a position to allow the accommodation of BEA. However, molecular docking can be sensitive to structural variations of the protein target ( 77 ). Thus, it is likely that the overall low scores of the poses resulted from using a conformational state that is not similar enough to the conformation induced by the combined presence of the substrate and BEA. Further experimental structural studies are required to unravel the conformational dynamics of BEA binding to human CTSB, and how this hinders the enzyme function. Taken together, our results show that BEA can suppress mouse and human CTSB activity in DCs. Moreover, we define the inhibition mode by showing that BEA can directly interact with human CTSB as an uncompetitive inhibitor. Molecular simulations and docking suggest a binding site at the edge of the active site, which is concordant with known structural data of other CTSB inhibitors. These first insights into the molecular mechanism of how BEA inhibits human CTSB should help further evaluating BEA as a candidate for cancer therapy. Author contributions St.S. conceived and supervised the study together with H.G., and they reviewed and edited the manuscript. X.Y. designed and performed the experiments, analyzed the results, and wrote the manuscript. P.M. designed and performed enzymatic assays, molecular simulations, and docking analyses and wrote the manuscript. M.H. performed NMR experiments and wrote the manuscript. A.V. performed CTSB assays and So.S. generated BMDCs. S.O. and S.M. performed the measurements of CTSB assays. A.H. and J.Q. wrote the manuscript. D.F. and M.U. gave suggestions on the project. Declaration of interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Acknowledgement This work was funded by the Deutsche Forschungsgemeinschaft (German Research Foundation) DFG-270650915/GRK2158. We are grateful for computational support and infrastructure provided by the “Zentrum für Informations- und Medientechnologie” (ZIM) at the Heinrich Heine University Düsseldorf and the computing time provided by the John von Neumann Institute for Computing (NIC) to H.G. on the supercomputer JUWELS at Jülich Supercomputing Centre (JSC) (user ID: VSK33). Footnotes ↵ † These authors share senior authorship References 1. ↵ Cavallo-Medved D , Moin K , Sloane B . Cathepsin B: Basis Sequence: Mouse . AFCS Nat Mol Pages . 2011 ; 2011 . 2. ↵ Kirschke H , Barrett AJ , Rawlings ND . Proteinases 1: lysosomal cysteine proteinases . Protein Profile . 1995 ; 2 ( 14 ): 1581 – 643 . OpenUrl PubMed 3. ↵ Turk V , Stoka V , Vasiljeva O , Renko M , Sun T , Turk B , et al. Cysteine cathepsins: from structure, function and regulation to new frontiers . Biochim Biophys Acta . 2012 ; 1824 ( 1 ): 68 – 88 . OpenUrl CrossRef PubMed Web of Science 4. ↵ Costa MG , Batista PR , Shida CS , Robert CH , Bisch PM , Pascutti PG . How does heparin prevent the pH inactivation of cathepsin B? Allosteric mechanism elucidated by docking and molecular dynamics . Bmc Genomics . 2010 ; 11 Suppl 5 ( Suppl 5 ): S5 . OpenUrl 5. ↵ Chapman HA , Riese RJ , Shi GP . Emerging roles for cysteine proteases in human biology . Annu Rev Physiol . 1997 ; 59 : 63 – 88 . OpenUrl CrossRef PubMed Web of Science 6. ↵ Embury CM , Dyavarshetty B , Lu Y , Wiederin JL , Ciborowski P , Gendelman HE , et al. Cathepsin B Improves ss-Amyloidosis and Learning and Memory in Models of Alzheimer’s Disease . J Neuroimmune Pharmacol . 2017 ; 12 ( 2 ): 340 – 52 . OpenUrl CrossRef PubMed 7. Hook VY , Kindy M , Reinheckel T , Peters C , Hook G . Genetic cathepsin B deficiency reduces beta-amyloid in transgenic mice expressing human wild-type amyloid precursor protein . Biochem Biophys Res Commun . 2009 ; 386 ( 2 ): 284 – 8 . OpenUrl CrossRef PubMed 8. ↵ Kindy MS , Yu J , Zhu H , El-Amouri SS , Hook V , Hook GR . Deletion of the cathepsin B gene improves memory deficits in a transgenic ALZHeimer’s disease mouse model expressing AbetaPP containing the wild-type beta-secretase site sequence . J Alzheimers Dis . 2012 ; 29 ( 4 ): 827 – 40 . OpenUrl CrossRef PubMed 9. ↵ Hashimoto Y , Kakegawa H , Narita Y , Hachiya Y , Hayakawa T , Kos J , et al. Significance of cathepsin B accumulation in synovial fluid of rheumatoid arthritis . Biochem Biophys Res Commun . 2001 ; 283 ( 2 ): 334 – 9 . OpenUrl CrossRef PubMed Web of Science 10. Van Acker GJ , Saluja AK , Bhagat L , Singh VP , Song AM , Steer ML . Cathepsin B inhibition prevents trypsinogen activation and reduces pancreatitis severity . Am J Physiol Gastrointest Liver Physiol . 2002 ; 283 ( 3 ): G794 – 800 . OpenUrl CrossRef PubMed Web of Science 11. ↵ Yoshifuji H , Umehara H , Maruyama H , Itoh M , Tanaka M , Kawabata D , et al. Amelioration of experimental arthritis by a calpain-inhibitory compound: regulation of cytokine production by E-64-d in vivo and in vitro . Int Immunol . 2005 ; 17 ( 10 ): 1327 – 36 . OpenUrl CrossRef PubMed Web of Science 12. ↵ Wu QQ , Xu M , Yuan Y , Li FF , Yang Z , Liu Y , et al. Cathepsin B deficiency attenuates cardiac remodeling in response to pressure overload via TNF-alpha/ASK1/JNK pathway . Am J Physiol Heart Circ Physiol . 2015 ; 308 ( 9 ): H1143 – 54 . OpenUrl CrossRef PubMed 13. ↵ Wuopio J , Hilden J , Bring C , Kastrup J , Sajadieh A , Jensen GB , et al. Cathepsin B and S as markers for cardiovascular risk and all-cause mortality in patients with stable coronary heart disease during 10 years: a CLARICOR trial sub-study . Atherosclerosis . 2018 ; 278 : 97 – 102 . OpenUrl CrossRef PubMed 14. ↵ Palermo C , Joyce JA . Cysteine cathepsin proteases as pharmacological targets in cancer . Trends Pharmacol Sci . 2008 ; 29 ( 1 ): 22 – 8 . OpenUrl CrossRef PubMed 15. ↵ Yan S , Sloane BF . Molecular regulation of human cathepsin B: implication in pathologies . Biol Chem . 2003 ; 384 ( 6 ): 845 – 54 . OpenUrl CrossRef PubMed Web of Science 16. ↵ Beckham TH , Lu P , Cheng JC , Zhao D , Turner LS , Zhang X , et al. Acid ceramidase-mediated production of sphingosine 1-phosphate promotes prostate cancer invasion through upregulation of cathepsin B . Int J Cancer . 2012 ; 131 ( 9 ): 2034 – 43 . OpenUrl CrossRef PubMed 17. Girotti MR , Fernandez M , Lopez JA , Camafeita E , Fernandez EA , Albar JP , et al. SPARC promotes cathepsin B-mediated melanoma invasiveness through a collagen I/alpha2beta1 integrin axis . J Invest Dermatol . 2011 ; 131 ( 12 ): 2438 – 47 . OpenUrl CrossRef PubMed Web of Science 18. ↵ Gopinathan A , Denicola GM , Frese KK , Cook N , Karreth FA , Mayerle J , et al. Cathepsin B promotes the progression of pancreatic ductal adenocarcinoma in mice . Gut . 2012 ; 61 ( 6 ): 877 – 84 . OpenUrl Abstract / FREE Full Text 19. ↵ Kostoulas G , Lang A , Nagase H , Baici A . Stimulation of angiogenesis through cathepsin B inactivation of the tissue inhibitors of matrix metalloproteinases . FEBS Lett . 1999 ; 455 ( 3 ): 286 – 90 . OpenUrl CrossRef PubMed 20. ↵ Peeters H , Zocher R , Kleinkauf H . Synthesis of beauvericin by a multifunctional enzyme . J Antibiot (Tokyo) . 1988 ; 41 ( 3 ): 352 – 9 . OpenUrl CrossRef PubMed 21. ↵ Logrieco A , Moretti A , Castella G , Kostecki M , Golinski P , Ritieni A , et al. Beauvericin production by Fusarium species . Appl Environ Microbiol . 1998 ; 64 ( 8 ): 3084 – 8 . OpenUrl Abstract / FREE Full Text 22. ↵ Peczynska-Czoch W , Urbanczyk MJ . Production of beauvericin by Beauveria bassiana on L-methionine enriched medium . Arch Immunol Ther Exp (Warsz) . 1991 ; 39 ( 1-2 ): 181 – 4 . OpenUrl PubMed 23. ↵ Lin HI , Lee YJ , Chen BF , Tsai MC , Lu JL , Chou CJ , et al. Involvement of Bcl-2 family, cytochrome c and caspase 3 in induction of apoptosis by beauvericin in human non-small cell lung cancer cells . Cancer Lett . 2005 ; 230 ( 2 ): 248 – 59 . OpenUrl CrossRef PubMed 24. ↵ Zhan J , Burns AM , Liu MX , Faeth SH , Gunatilaka AA . Search for cell motility and angiogenesis inhibitors with potential anticancer activity: beauvericin and other constituents of two endophytic strains of Fusarium oxysporum . J Nat Prod . 2007 ; 70 ( 2 ): 227 – 32 . OpenUrl CrossRef PubMed 25. ↵ Al Khoury C , Bashir Z , Tokajian S , Nemer N , Merhi G , Nemer G . In silico evidence of beauvericin antiviral activity against SARS-CoV-2 . Comput Biol Med . 2022 ; 141 : 105171 . OpenUrl CrossRef PubMed 26. ↵ Shin CG , An DG , Song HH , Lee C. Beauvericin and enniatins H, I and MK1688 are new potent inhibitors of human immunodeficiency virus type-1 integrase . J Antibiot (Tokyo) . 2009 ; 62 ( 12 ): 687 – 90 . OpenUrl CrossRef PubMed 27. ↵ Wu XF , Xu R , Ouyang ZJ , Qian C , Shen Y , Wu XD , et al. Beauvericin ameliorates experimental colitis by inhibiting activated T cells via downregulation of the PI3K/Akt signaling pathway . PLoS One . 2013 ; 8 ( 12 ): e83013 . OpenUrl CrossRef PubMed 28. ↵ Yoo S , Kim MY , Cho JY . Beauvericin, a cyclic peptide, inhibits inflammatory responses in macrophages by inhibiting the NF-kappaB pathway . Korean J Physiol Pharmacol . 2017 ; 21 ( 4 ): 449 – 56 . OpenUrl CrossRef PubMed 29. ↵ Meca G , Sospedra I , Soriano JM , Ritieni A , Moretti A , Manes J . Antibacterial effect of the bioactive compound beauvericin produced by Fusarium proliferatum on solid medium of wheat . Toxicon . 2010 ; 56 ( 3 ): 349 – 54 . OpenUrl CrossRef PubMed 30. Nilanonta C , Isaka M , Kittakoop P , Palittapongarnpim P , Kamchonwongpaisan S , Pittayakhajonwut D , et al. Antimycobacterial and antiplasmodial cyclodepsipeptides from the insect pathogenic fungus Paecilomyces tenuipes BCC 1614 . Planta Med . 2000 ; 66 ( 8 ): 756 – 8 . OpenUrl CrossRef PubMed Web of Science 31. Wang Q , Xu L . Beauvericin, a bioactive compound produced by fungi: a short review . Molecules . 2012 ; 17 ( 3 ): 2367 – 77 . OpenUrl CrossRef PubMed 32. ↵ Xu L , Wang J , Zhao J , Li P , Shan T , Wang J , et al. Beauvericin from the endophytic fungus, Fusarium redolens, isolated from Dioscorea zingiberensis and its antibacterial activity . Nat Prod Commun . 2010 ; 5 ( 5 ): 811 – 4 . OpenUrl PubMed 33. ↵ Yang X , Ali S , Zhao M , Richter L , Schafer V , Schliehe-Diecks J , et al. The Mycotoxin Beauvericin Exhibits Immunostimulatory Effects on Dendritic Cells via Activating the TLR4 Signaling Pathway . Front Immunol . 2022 ; 13 : 856230 . OpenUrl CrossRef PubMed 34. ↵ Fornelli F , Minervini F , Logrieco A . Cytotoxicity of fungal metabolites to lepidopteran (Spodoptera frugiperda) cell line (SF-9) . J Invertebr Pathol . 2004 ; 85 ( 2 ): 74 – 9 . OpenUrl CrossRef PubMed 35. ↵ Maniania NK , Portilla MM , Amnulla FM , Mfuti DK , Darie A , Dhiman G , et al. Infectivity of Entomopathogenic Fungal Isolates Against Tarnished Plant Bug Lygus lineolaris (Hemiptera: Miridae) . J Insect Sci . 2022 ; 22 ( 4 ). 36. ↵ Al Khoury C , Nemer N , Nemer G . Beauvericin potentiates the activity of pesticides by neutralizing the ATP-binding cassette transporters in arthropods . Sci Rep . 2021 ; 11 ( 1 ): 10865 . OpenUrl CrossRef PubMed 37. ↵ Shimada A , Fujioka S , Koshino H , Kimura Y . Nematicidal activity of beauvericin produced by the fungus Fusarium bulbicola . Z Naturforsch C J Biosci . 2010 ; 65 ( 3-4 ): 207 – 10 . OpenUrl CrossRef PubMed 38. ↵ Silva TL , Toffano L , Fernandes JB , das Gracas Fernandes da Silva MF , de Sousa LRF , Vieira PC . Mycotoxins from Fusarium proliferatum: new inhibitors of papain-like cysteine proteases . Braz J Microbiol . 2020 ; 51 ( 3 ): 1169 – 75 . OpenUrl CrossRef PubMed 39. ↵ Scheu S , Dresing P , Locksley RM . Visualization of IFNbeta production by plasmacytoid versus conventional dendritic cells under specific stimulation conditions in vivo . Proc Natl Acad Sci U S A . 2008 ; 105 ( 51 ): 20416 – 21 . OpenUrl Abstract / FREE Full Text 40. ↵ Holken JM , Teusch N . The Monocytic Cell Line THP-1 as a Validated and Robust Surrogate Model for Human Dendritic Cells . International journal of molecular sciences . 2023 ; 24 ( 2 ). 41. ↵ Mayer M , Meyer B . Group epitope mapping by saturation transfer difference NMR to identify segments of a ligand in direct contact with a protein receptor . J Am Chem Soc . 2001 ; 123 ( 25 ): 6108 – 17 . OpenUrl CrossRef PubMed Web of Science 42. ↵ Cornish-Bowden A . A simple graphical method for determining the inhibition constants of mixed, uncompetitive and non-competitive inhibitors . Biochem J . 1974 ; 137 ( 1 ): 143 – 4 . OpenUrl Abstract / FREE Full Text 43. ↵ Webb B , Sali A . Protein Structure Modeling with MODELLER . Methods Mol Biol . 2021 ; 2199 : 239 – 55 . OpenUrl CrossRef PubMed 44. ↵ Schott-Verdugo S , Gohlke H . PACKMOL-Memgen: A Simple-To-Use, Generalized Workflow for Membrane-Protein–Lipid-Bilayer System Building . Journal of Chemical Information and Modeling . 2019 ; 59 ( 6 ): 2522 – 8 . OpenUrl CrossRef PubMed 45. ↵ Rostkowski M , Olsson MH , Sondergaard CR , Jensen JH . Graphical analysis of pH-dependent properties of proteins predicted using PROPKA . BMC Struct Biol . 2011 ; 11 : 6 . OpenUrl CrossRef PubMed 46. ↵ Case DA , Aktulga HM , Belfon K , Cerutti DS , Cisneros GA , Cruzeiro VWD , et al. AmberTools . Journal of Chemical Information and Modeling . 2023 ; 63 ( 20 ): 6183 – 91 . OpenUrl CrossRef PubMed 47. ↵ Tian C , Kasavajhala K , Belfon KAA , Raguette L , Huang H , Migues AN , et al. ff19SB: Amino-Acid-Specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution . Journal of Chemical Theory and Computation . 2020 ; 16 ( 1 ): 528 – 52 . OpenUrl CrossRef 48. ↵ Jorgensen WL , Chandrasekhar J , Madura JD , Impey RW , Klein ML . Comparison of simple potential functions for simulating liquid water . The Journal of Chemical Physics . 1983 ; 79 ( 2 ): 926 – 35 . OpenUrl CrossRef PubMed Web of Science 49. ↵ Joung IS , Cheatham TE , III . Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations . The Journal of Physical Chemistry B . 2008 ; 112 ( 30 ): 9020 – 41 . OpenUrl CrossRef PubMed 50. ↵ Salomon-Ferrer R , Götz AW , Poole D , Le Grand S , Walker RC . Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald . Journal of Chemical Theory and Computation . 2013 ; 9 ( 9 ): 3878 – 88 . OpenUrl CrossRef 51. ↵ Miyamoto S , Kollman PA . Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models . Journal of Computational Chemistry . 1992 ; 13 ( 8 ): 952 – 62 . OpenUrl CrossRef PubMed Web of Science 52. ↵ Roe DR , Cheatham TE , 3rd . . PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data . J Chem Theory Comput . 2013 ; 9 ( 7 ): 3084 – 95 . OpenUrl CrossRef PubMed 53. ↵ Sindhikara D , Spronk SA , Day T , Borrelli K , Cheney DL , Posy SL . Improving Accuracy, Diversity, and Speed with Prime Macrocycle Conformational Sampling . Journal of Chemical Information and Modeling . 2017 ; 57 ( 8 ): 1881 – 94 . OpenUrl CrossRef PubMed 54. ↵ Repasky MP , Shelley M , Friesner RA . Flexible ligand docking with Glide . Curr Protoc Bioinformatics . 2007 ;Chapter 8:Unit 8 12. 55. ↵ Ochs RS . Understanding Enzyme Inhibition . Journal of Chemical Education . 2000 ; 77 ( 11 ): 1453 . OpenUrl CrossRef 56. ↵ Dixon M . The Determination of Enzyme Inhibitor Constants . Biochem J . 1953 ; 55 ( 1 ): 170 – 1 . OpenUrl FREE Full Text 57. ↵ Costa MGS , Batista PR , Shida CS , Robert CH , Bisch PM , Pascutti PG . How does heparin prevent the pH inactivation of cathepsin B? Allosteric mechanism elucidated by docking and molecular dynamics . Bmc Genomics . 2010 ; 11 . 58. ↵ Ghanakota P , Carlson HA . Driving Structure-Based Drug Discovery through Cosolvent Molecular Dynamics . Journal of Medicinal Chemistry . 2016 ; 59 ( 23 ): 10383 – 99 . OpenUrl CrossRef PubMed 59. ↵ Gohlke H , Hergert U , Meyer T , Mulnaes D , Grieshaber MK , Smits SHJ , et al. Binding Region of Alanopine Dehydrogenase Predicted by Unbiased Molecular Dynamics Simulations of Ligand Diffusion . Journal of Chemical Information and Modeling . 2013 ; 53 ( 10 ): 2493 – 8 . OpenUrl CrossRef PubMed 60. ↵ Renko M , Pozgan U , Majera D , Turk D . Stefin A displaces the occluding loop of cathepsin B only by as much as required to bind to the active site cleft . Febs J . 2010 ; 277 ( 20 ): 4338 – 45 . OpenUrl CrossRef PubMed 61. ↵ Podobnik M , Kuhelj R , Turk V , Turk D . Crystal structure of the wild-type human procathepsin B at 2.5 Å resolution reveals the native active site of a papain-like cysteine protease zymogen11Edited by R. Huber . Journal of Molecular Biology . 1997 ; 271 ( 5 ): 774 – 88 . OpenUrl CrossRef PubMed Web of Science 62. ↵ Silva TL , Toffano L , Fernandes JB , das Graças Fernandes da Silva MF , de Sousa LRF , Vieira PC . Mycotoxins from Fusarium proliferatum: new inhibitors of papain-like cysteine proteases . Braz J Microbiol . 2020 ; 51 ( 3 ): 1169 – 75 . OpenUrl CrossRef PubMed 63. ↵ Calo L , Fornelli F , Ramires R , Nenna S , Tursi A , Caiaffa MF , et al. Cytotoxic effects of the mycotoxin beauvericin to human cell lines of myeloid origin . Pharmacol Res . 2004 ; 49 ( 1 ): 73 – 7 . OpenUrl CrossRef PubMed 64. Heilos D , Rodriguez-Carrasco Y , Englinger B , Timelthaler G , van Schoonhoven S , Sulyok M , et al. The Natural Fungal Metabolite Beauvericin Exerts Anticancer Activity In Vivo: A Pre-Clinical Pilot Study . Toxins (Basel) . 2017 ; 9 ( 9 ). 65. Jow GM , Chou CJ , Chen BF , Tsai JH . Beauvericin induces cytotoxic effects in human acute lymphoblastic leukemia cells through cytochrome c release, caspase 3 activation: the causative role of calcium . Cancer Lett . 2004 ; 216 ( 2 ): 165 – 73 . OpenUrl CrossRef PubMed 66. ↵ Wang G , Qiao Y , Zhao Y , Song Y , Li M , Jin M , et al. Beauvericin exerts an anti-tumor effect on hepatocellular carcinoma by inducing PI3K/AKT-mediated apoptosis . Arch Biochem Biophys . 2023 ; 745 : 109720 . OpenUrl CrossRef 67. ↵ Wu Q , Patocka J , Nepovimova E , Kuca K . A Review on the Synthesis and Bioactivity Aspects of Beauvericin, a Fusarium Mycotoxin . Front Pharmacol . 2018 ; 9 : 1338 . OpenUrl CrossRef PubMed 68. ↵ Xie Z , Zhao M , Yan C , Kong W , Lan F , Narengaowa , et al. Cathepsin B in programmed cell death machinery: mechanisms of execution and regulatory pathways . Cell Death Dis . 2023 ; 14 ( 4 ): 255 . OpenUrl CrossRef PubMed 69. ↵ Talukdar R , Sareen A , Zhu H , Yuan Z , Dixit A , Cheema H , et al. Release of Cathepsin B in Cytosol Causes Cell Death in Acute Pancreatitis . Gastroenterology . 2016 ; 151 ( 4 ): 747 – 58 e5. OpenUrl CrossRef PubMed 70. ↵ Halangk W , Lerch MM , Brandt-Nedelev B , Roth W , Ruthenbuerger M , Reinheckel T , et al. Role of cathepsin B in intracellular trypsinogen activation and the onset of acute pancreatitis . J Clin Invest . 2000 ; 106 ( 6 ): 773 – 81 . OpenUrl CrossRef PubMed Web of Science 71. ↵ Gonzalez-Leal IJ , Roger B , Schwarz A , Schirmeister T , Reinheckel T , Lutz MB , et al. Cathepsin B in antigen-presenting cells controls mediators of the Th1 immune response during Leishmania major infection . PLoS Negl Trop Dis . 2014 ; 8 ( 9 ): e3194 . OpenUrl CrossRef PubMed 72. Akkari L , Gocheva V , Quick ML , Kester JC , Spencer AK , Garfall AL , et al. Combined deletion of cathepsin protease family members reveals compensatory mechanisms in cancer . Genes Dev . 2016 ; 30 ( 2 ): 220 – 32 . OpenUrl Abstract / FREE Full Text 73. Sevenich L , Schurigt U , Sachse K , Gajda M , Werner F , Muller S , et al. Synergistic antitumor effects of combined cathepsin B and cathepsin Z deficiencies on breast cancer progression and metastasis in mice . Proc Natl Acad Sci U S A . 2010 ; 107 ( 6 ): 2497 – 502 . OpenUrl Abstract / FREE Full Text 74. Vidak E , Javorsek U , Vizovisek M , Turk B . Cysteine Cathepsins and their Extracellular Roles: Shaping the Microenvironment . Cells . 2019 ; 8 ( 3 ). 75. ↵ Schmitz J , Gilberg E , Löser R , Bajorath J , Bartz U , Gütschow M . Cathepsin B: Active site mapping with peptidic substrates and inhibitors . Bioorganic & Medicinal Chemistry . 2019 ; 27 ( 1 ): 1 – 15 . OpenUrl CrossRef PubMed 76. ↵ Cornish-Bowden A . Why is uncompetitive inhibition so rare? FEBS Letters . 1986 ; 203 ( 1 ): 3 – 6 . OpenUrl CrossRef PubMed 77. ↵ Erickson JA , Jalaie M , Robertson DH , Lewis RA , Vieth M . Lessons in Molecular Recognition: The Effects of Ligand and Protein Flexibility on Molecular Docking Accuracy . Journal of Medicinal Chemistry . 2004 ; 47 ( 1 ): 45 – 55 . OpenUrl CrossRef PubMed Web of Science View the discussion thread. Back to top Previous Next Posted January 19, 2025. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following The mycotoxin Beauvericin is an uncompetitive inhibitor of Cathepsin B Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share The mycotoxin Beauvericin is an uncompetitive inhibitor of Cathepsin B Xiaoli Yang , Pablo Cea-Medina , Mohanraj Gopalswamy , Aparna Vaidya , Sonja Schavier , Shixin Oltzen , Sofie Moßner , Anfei Huang , Jing Qi , Doreen Floss , Markus Uhrberg , Holger Gohlke , Stefanie Scheu bioRxiv 2025.01.17.633572; doi: https://doi.org/10.1101/2025.01.17.633572 Share This Article: Copy Citation Tools The mycotoxin Beauvericin is an uncompetitive inhibitor of Cathepsin B Xiaoli Yang , Pablo Cea-Medina , Mohanraj Gopalswamy , Aparna Vaidya , Sonja Schavier , Shixin Oltzen , Sofie Moßner , Anfei Huang , Jing Qi , Doreen Floss , Markus Uhrberg , Holger Gohlke , Stefanie Scheu bioRxiv 2025.01.17.633572; doi: https://doi.org/10.1101/2025.01.17.633572 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Biochemistry Subject Areas All Articles Animal Behavior and Cognition (7616) Biochemistry (17625) Bioengineering (13852) Bioinformatics (41825) Biophysics (21397) Cancer Biology (18524) Cell Biology (25417) Clinical Trials (138) Developmental Biology (13350) Ecology (19858) Epidemiology (2067) Evolutionary Biology (24277) Genetics (15581) Genomics (22459) Immunology (17698) Microbiology (40278) Molecular Biology (17134) Neuroscience (88400) Paleontology (666) Pathology (2823) Pharmacology and Toxicology (4812) Physiology (7632) Plant Biology (15106) Scientific Communication and Education (2042) Synthetic Biology (4281) Systems Biology (9807) Zoology (2266)