Dissecting the regulatory roles of ORM proteins in the sphingolipid pathway of plants

Abstract Sphingolipids are a vital component of plant cellular endomembranes and carry out multiple functional and regulatory roles. Different sphingolipid species confer rigidity to the membrane structure, facilitate trafficking of secretory proteins, and initiate programmed cell death. Although the regulation of the sphingolipid pathway is yet to be uncovered, increasing evidence has pointed to orosomucoid proteins (ORMs) playing a major regulatory role and potentially interacting with a number of components in the pathway, including both enzymes and sphingolipids. However, experimental exploration of new regulatory interactions is time consuming and often infeasible. In this work, a computational approach was taken to address this challenge. A metabolic network of the sphingolipid pathway in plants was reconstructed. The steady-state rates of reactions in the network were then determined through measurements of growth and cellular composition of the different sphingolipids in Arabidopsis seedlings. The Ensemble modeling framework was modified to accurately account for activation mechanisms and subsequently used to generate sets of kinetic parameters that converge to the measured steady-state fluxes in a thermodynamically consistent manner. In addition, the framework was appended with an additional module to automate screening the parameters and to output models consistent with previously reported network responses to different perturbations. By analyzing the network’s response in the presence of different combinations of regulatory mechanisms, the model captured the experimentally observed repressive effect of ORMs on SPT. Furthermore, predictions point to a second regulatory role of ORM proteins, namely as an activator of class II (or LOH1 and LOH3) ceramide synthases. This activating role was found to be modulated by the concentration of free ceramides, where an accumulation of these sphingolipid species dampened the activating effect of ORMs on ceramide synthase. The predictions pave the way for future guided experiments and have implications in engineering crops with higher biotic stress tolerance. Author summary Due to their vital functional and regulatory roles in plant cells, increasing interest has gone into obtaining a complete understanding of the regulatory behavior of the sphingolipid pathway. However, the process of identifying new regulatory interactions is time consuming and often infeasible. To address this issue, ensemble modeling was used as an in silico method to test the ability of different regulatory schemes to predict all known pathway responses in a thermodynamically consistent manner. The analysis resulted in a significant reduction in the number of possible regulatory interactions. Mainly, the model predicts regulatory interactions between ceramides, ORMs, and ceramide synthases (especially class II). This framework can pave the way for biochemists to systematically identify plausible regulatory networks in understudied metabolic networks where knowledge on the underlying regulatory mechanisms is often missing. As future experimental works explore these predictions, an iterative cycle can begin wherein model predictions allow for targeted experiments which in turn generate results that can be reincorporated into the model to further increase prediction accuracy. Such a model-driven approach will significantly reduce the solution space traversed by the experimentalist.