FALCON systematically interrogates free fatty acid biology and identifies a novel mediator of lipotoxicity.

Cellular exposure to free fatty acids (FFAs) is implicated in the pathogenesis of obesity-associated diseases. However, there are no scalable approaches to comprehensively assess the diverse FFAs circulating in human plasma. Furthermore, assessing how FFA-mediated processes interact with genetic risk for disease remains elusive. Here, we report the design and implementation of fatty acid library for comprehensive ontologies (FALCON), an unbiased, scalable, and multimodal interrogation of 61 structurally diverse FFAs. We identified a subset of lipotoxic monounsaturated fatty acids associated with decreased membrane fluidity. Furthermore, we prioritized genes that reflect the combined effects of harmful FFA exposure and genetic risk for type 2 diabetes (T2D). We found that c-MAF-inducing protein (CMIP) protects cells from FFA exposure by modulating Akt signaling. In sum, FALCON empowers the study of fundamental FFA biology and offers an integrative approach to identify much needed targets for diverse diseases associated with disordered FFA metabolism.

FALCON systematically interrogates free fatty acid biology and identifies a novel mediator of lipotoxicity.

Cellular exposure to free fatty acids (FFA) is implicated in the pathogenesis of obesity-associated diseases. However, studies to date have assumed that a few select FFAs are representative of broad structural categories, and there are no scalable approaches to comprehensively assess the biological processes induced by exposure to diverse FFAs circulating in human plasma. Furthermore, assessing how these FFA- mediated processes interact with genetic risk for disease remains elusive. Here we report the design and implementation of FALCON (Fatty Acid Library for Comprehensive ONtologies) as an unbiased, scalable and multimodal interrogation of 61 structurally diverse FFAs. We identified a subset of lipotoxic monounsaturated fatty acids (MUFAs) with a distinct lipidomic profile associated with decreased membrane fluidity. Furthermore, we developed a new approach to prioritize genes that reflect the combined effects of exposure to harmful FFAs and genetic risk for type 2 diabetes (T2D). Importantly, we found that c-MAF inducing protein (CMIP) protects cells from exposure to FFAs by modulating Akt signaling and we validated the role of CMIP in human pancreatic beta cells. In sum, FALCON empowers the study of fundamental FFA biology and offers an integrative approach to identify much needed targets for diverse diseases associated with disordered FFA metabolism. Highlights: FALCON (Fatty Acid Library for Comprehensive ONtologies) enables multimodal profiling of 61 free fatty acids (FFAs) to reveal 5 FFA clusters with distinct biological effectsFALCON is applicable to many and diverse cell typesA subset of monounsaturated FAs (MUFAs) equally or more toxic than canonical lipotoxic saturated FAs (SFAs) leads to decreased membrane fluidityNew approach prioritizes genes that represent the combined effects of environmental (FFA) exposure and genetic risk for diseaseC-Maf inducing protein (CMIP) is identified as a suppressor of FFA-induced lipotoxicity via Akt-mediated signaling.

Targeting a Braf/Mapk pathway rescues podocyte lipid peroxidation in CoQ-deficiency kidney disease.

Mutations affecting mitochondrial coenzyme Q (CoQ) biosynthesis lead to kidney failure due to selective loss of podocytes, essential cells of the kidney filter. Curiously, neighboring tubular epithelial cells are spared early in disease despite higher mitochondrial content. We sought to illuminate noncanonical, cell-specific roles for CoQ, independently of the electron transport chain (ETC). Here, we demonstrate that CoQ depletion caused by Pdss2 enzyme deficiency in podocytes results in perturbations in polyunsaturated fatty acid (PUFA) metabolism and the Braf/Mapk pathway rather than ETC dysfunction. Single-nucleus RNA-Seq from kidneys of Pdss2kd/kd mice with nephrotic syndrome and global CoQ deficiency identified a podocyte-specific perturbation of the Braf/Mapk pathway. Treatment with GDC-0879, a Braf/Mapk-targeting compound, ameliorated kidney disease in Pdss2kd/kd mice. Mechanistic studies in Pdss2-depleted podocytes revealed a previously unknown perturbation in PUFA metabolism that was confirmed in vivo. Gpx4, an enzyme that protects against PUFA-mediated lipid peroxidation, was elevated in disease and restored after GDC-0879 treatment. We demonstrate broader human disease relevance by uncovering patterns of GPX4 and Braf/Mapk pathway gene expression in tissue from patients with kidney diseases. Our studies reveal ETC-independent roles for CoQ in podocytes and point to Braf/Mapk as a candidate pathway for the treatment of kidney diseases.

Simulation Strategies for Calcium Microdomains and Calcium Noise.

In this article, we present an overview of simulation strategies in the context of subcellular domains where calcium-dependent signaling plays an important role. The presentation follows the spatial and temporal scales involved and represented by each algorithm. As an exemplary cell type, we will mainly cite work done on striated muscle cells, i.e. skeletal and cardiac muscle. For these cells, a wealth of ultrastructural, biophysical and electrophysiological data is at hand. Moreover, these cells also express ubiquitous signaling pathways as they are found in many other cell types and thus, the generalization of the methods and results presented here is straightforward.The models considered comprise the basic calcium signaling machinery as found in most excitable cell types including Ca2+ ions, diffusible and stationary buffer systems, and calcium regulated calcium release channels. Simulation strategies can be differentiated in stochastic and deterministic algorithms. Historically, deterministic approaches based on the macroscopic reaction rate equations were the first models considered. As experimental methods elucidated highly localized Ca2+ signaling events occurring in femtoliter volumes, stochastic methods were increasingly considered. However, detailed simulations of single molecule trajectories are rarely performed as the computational cost implied is too large. On the mesoscopic level, Gillespie's algorithm is extensively used in the systems biology community and with increasing frequency also in models of microdomain calcium signaling. To increase computational speed, fast approximations were derived from Gillespie's exact algorithm, most notably the chemical Langevin equation and the τ-leap algorithm. Finally, in order to integrate deterministic and stochastic effects in multiscale simulations, hybrid algorithms are increasingly used. These include stochastic models of ion channels combined with deterministic descriptions of the calcium buffering and diffusion system on the one hand, and algorithms that switch between deterministic and stochastic simulation steps in a context-dependent manner on the other. The basic assumptions of the listed methods as well as implementation schemes are given in the text. We conclude with a perspective on possible future developments of the field.

Small Molecule Targets TMED9 and Promotes Lysosomal Degradation to Reverse Proteinopathy.

Intracellular accumulation of misfolded proteins causes toxic proteinopathies, diseases without targeted therapies. Mucin 1 kidney disease (MKD) results from a frameshift mutation in the MUC1 gene (MUC1-fs). Here, we show that MKD is a toxic proteinopathy. Intracellular MUC1-fs accumulation activated the ATF6 unfolded protein response (UPR) branch. We identified BRD4780, a small molecule that clears MUC1-fs from patient cells, from kidneys of knockin mice and from patient kidney organoids. MUC1-fs is trapped in TMED9 cargo receptor-containing vesicles of the early secretory pathway. BRD4780 binds TMED9, releases MUC1-fs, and re-routes it for lysosomal degradation, an effect phenocopied by TMED9 deletion. Our findings reveal BRD4780 as a promising lead for the treatment of MKD and other toxic proteinopathies. Generally, we elucidate a novel mechanism for the entrapment of misfolded proteins by cargo receptors and a strategy for their release and anterograde trafficking to the lysosome.

GDC-0879, a BRAFV600E Inhibitor, Protects Kidney Podocytes from Death.

Progressive kidney diseases affect approximately 500 million people worldwide. Podocytes are terminally differentiated cells of the kidney filter, the loss of which leads to disease progression and kidney failure. To date, there are no therapies to promote podocyte survival. Drug repurposing may therefore help accelerate the development of cures in an area of tremendous unmet need. In a newly developed high-throughput screening assay of podocyte viability, we identified the BRAFV600E inhibitor GDC-0879 and the adenylate cyclase agonist forskolin as podocyte-survival-promoting compounds. GDC-0879 protects podocytes from injury through paradoxical activation of the MEK/ERK pathway. Forskolin promotes podocyte survival by attenuating protein biosynthesis. Importantly, GDC-0879 and forskolin are shown to promote podocyte survival against an array of cellular stressors. This work reveals new therapeutic targets for much needed podocyte-protective therapies and provides insights into the use of GDC-0879-like molecules for the treatment of progressive kidney diseases.

A small-molecule inhibitor of TRPC5 ion channels suppresses progressive kidney disease in animal models.

Progressive kidney diseases are often associated with scarring of the kidney's filtration unit, a condition called focal segmental glomerulosclerosis (FSGS). This scarring is due to loss of podocytes, cells critical for glomerular filtration, and leads to proteinuria and kidney failure. Inherited forms of FSGS are caused by Rac1-activating mutations, and Rac1 induces TRPC5 ion channel activity and cytoskeletal remodeling in podocytes. Whether TRPC5 activity mediates FSGS onset and progression is unknown. We identified a small molecule, AC1903, that specifically blocks TRPC5 channel activity in glomeruli of proteinuric rats. Chronic administration of AC1903 suppressed severe proteinuria and prevented podocyte loss in a transgenic rat model of FSGS. AC1903 also provided therapeutic benefit in a rat model of hypertensive proteinuric kidney disease. These data indicate that TRPC5 activity drives disease and that TRPC5 inhibitors may be valuable for the treatment of progressive kidney diseases.

Lysine trimethylation regulates 78-kDa glucose-regulated protein proteostasis during endoplasmic reticulum stress.

The up-regulation of chaperones such as the 78-kDa glucose-regulated protein (GRP78, also referred to as BiP or HSPA5) is part of the adaptive cellular response to endoplasmic reticulum (ER) stress. GRP78 is widely used as a marker of the unfolded protein response, associated with sustained ER stress. Here we report the discovery of a proteostatic mechanism involving GRP78 trimethylation in the context of ER stress. Using mass spectrometry-based proteomics, we identified two GRP78 fractions, one homeostatic and one induced by ER stress. ER stress leads to de novo biosynthesis of non-trimethylated GRP78, whereas homeostatic, METTL21A-dependent lysine 585-trimethylated GRP78 is reduced. This proteostatic mechanism, dependent on the posttranslational modification of GRP78, allows cells to differentially regulate specific protein abundance during cellular stress.

Calcium, TRPC channels, and regulation of the actin cytoskeleton in podocytes: towards a future of targeted therapies.

With more than 6,000 new pediatric patients with treatment-resistant nephrotic syndrome in the US each year alone, the unmet need for novel, podocyte-specific therapies is substantial. Recently, the established therapeutic benefit of angiotensin-converting enzyme (ACE) inhibitors and angiotensin receptor blockers (ARB) was used as a starting point to gain insight into the pathomechanism of primary podocytopathies. A calcium (Ca(2+))-mediated pathway has been identified that connects the angiotensin type 1 receptor (AT1R) to podocyte cytoskeletal dynamics, essential for a functioning glomerular filtration barrier. This discovery provided an important missing piece in our understanding of the pathomechanism of filter barrier damage, revealing Ca(2+) signaling as critical for podocyte health and disease. The identification of the two Ca(2+) permeant channels TRPC5 and TRPC6 as mediators of this pathway not only bolstered the importance of podocyte cytoskeleton dynamics but also revealed promising drug targets for treatment-resistant nephrotic syndrome. This review will focus on this novel signaling pathway in primary podocytopathies and its implications for next-generation therapies for glomerular disease.

Exact stochastic simulation of a calcium microdomain reveals the impact of Ca²âº fluctuations on IP3R gating.

In this study, we numerically analyzed the nonlinear Ca(2+)-dependent gating dynamics of a single, nonconducting inositol 1,4,5-trisphosphate receptor (IP3R) channel, using an exact and fully stochastic simulation algorithm that includes channel gating, Ca(2+) buffering, and Ca(2+) diffusion. The IP3R is a ubiquitous intracellular Ca(2+) release channel that plays an important role in the formation of complex spatiotemporal Ca(2+) signals such as waves and oscillations. Dynamic subfemtoliter Ca(2+) microdomains reveal low copy numbers of Ca(2+) ions, buffer molecules, and IP3Rs, and stochastic fluctuations arising from molecular interactions and diffusion do not average out. In contrast to models treating calcium dynamics deterministically, the stochastic approach accounts for this molecular noise. We varied Ca(2+) diffusion coefficients and buffer reaction rates to tune the autocorrelation properties of Ca(2+) noise and found a distinct relation between the autocorrelation time τac, the mean channel open and close times, and the resulting IP3R open probability PO. We observed an increased PO for shorter noise autocorrelation times, caused by increasing channel open times and decreasing close times. In a pure diffusion model the effects become apparent at elevated calcium concentrations, e.g., at [Ca(2+)] = 25 µM, τac = 0.082 ms, the IP3R open probability increased by ≈20% and mean open times increased by ≈4 ms, compared to a zero noise model. We identified the inactivating Ca(2+) binding site of IP3R subunits as the primarily noise-susceptible element of the De Young and Keizer model. Short Ca(2+) noise autocorrelation times decrease the probability of Ca(2+) association and consequently increase IPvR activity. These results suggest a functional role of local calcium noise properties on calcium-regulated target molecules such as the ubiquitous IP3R. This finding may stimulate novel experimental approaches analyzing the role of calcium noise properties on microdomain behavior.

Microdomain calcium fluctuations as a colored noise process.

Calcium ions play a key role in subcellular signaling as localized transients of the intracellular calcium concentration modify the activity of ion channels, enzymes and transcription factors, among others. The intracellular calcium concentration is inherently noisy, as diffusion, the transient binding to and dissociation from buffer molecules and stochastically gating calcium channels contribute to the fluctuations of the local copy number of Ca(2+) ions. We study the properties of the fluctuating calcium concentration in sub-femtoliter volumes using an exact stochastic simulation algorithm and approximations to the exact stochastic solution. It is shown that the time course of the local calcium concentration represents a colored noise process whose autocorrelation time is a function of buffer kinetics and diffusion constants. Using the chemical Langevin description and the excess buffer approximation of the process, fast approximative algorithms and theoretical connections to the Ornstein-Uhlenbeck process are obtained. In a generic example, we show how calcium noise can couple to the dynamics of a single variable moving in a double-well potential, leading to a colored noise induced transition. Our work shows how a multitude of intracellular signaling pathways may be influenced by the inherent stochasticity of calcium signals, a key messenger in virtually any cell type, and how the calcium signal can be implemented efficiently in cellular signaling models.

Simulation strategies for calcium microdomains and calcium-regulated calcium channels.

In this article, we present an overview of simulation strategies in the context of subcellular domains where calcium-dependent signaling plays an important role. The presentation follows the spatial and temporal scales involved and represented by each algorithm. As an exemplary cell type, we will mainly cite work done on striated muscle cells, i.e. skeletal and cardiac muscle. For these cells, a wealth of ultrastructural, biophysical and electrophysiological data is at hand. Moreover, these cells also express ubiquitous signaling pathways as they are found in many other cell types and thus, the generalization of the methods and results presented here is straightforward.The models considered comprise the basic calcium signaling machinery as found in most excitable cell types including Ca(2+) ions, diffusible and stationary buffer systems, and calcium regulated calcium release channels. Simulation strategies can be differentiated in stochastic and deterministic algorithms. Historically, deterministic approaches based on the macroscopic reaction rate equations were the first models considered. As experimental methods elucidated highly localized Ca(2+) signaling events occurring in femtoliter volumes, stochastic methods were increasingly considered. However, detailed simulations of single molecule trajectories are rarely performed as the computational cost implied is too large. On the mesoscopic level, Gillespie's algorithm is extensively used in the systems biology community and with increasing frequency also in models of microdomain calcium signaling. To increase computational speed, fast approximations were derived from Gillespie's exact algorithm, most notably the chemical Langevin equation and the τ-leap algorithm. Finally, in order to integrate deterministic and stochastic effects in multiscale simulations, hybrid algorithms are increasingly used. These include stochastic models of ion channels combined with deterministic descriptions of the calcium buffering and diffusion system on the one hand, and algorithms that switch between deterministic and stochastic simulation steps in a context-dependent manner on the other. The basic assumptions of the listed methods as well as implementation schemes are given in the text. We conclude with a perspective on possible future developments of the field.

Exact and approximate stochastic simulation of intracellular calcium dynamics.

In simulations of chemical systems, the main task is to find an exact or approximate solution of the chemical master equation (CME) that satisfies certain constraints with respect to computation time and accuracy. While Brownian motion simulations of single molecules are often too time consuming to represent the mesoscopic level, the classical Gillespie algorithm is a stochastically exact algorithm that provides satisfying results in the representation of calcium microdomains. Gillespie's algorithm can be approximated via the tau-leap method and the chemical Langevin equation (CLE). Both methods lead to a substantial acceleration in computation time and a relatively small decrease in accuracy. Elimination of the noise terms leads to the classical, deterministic reaction rate equations (RRE). For complex multiscale systems, hybrid simulations are increasingly proposed to combine the advantages of stochastic and deterministic algorithms. An often used exemplary cell type in this context are striated muscle cells (e.g., cardiac and skeletal muscle cells). The properties of these cells are well described and they express many common calcium-dependent signaling pathways. The purpose of the present paper is to provide an overview of the aforementioned simulation approaches and their mutual relationships in the spectrum ranging from stochastic to deterministic algorithms.