An in silico model of enterocytic glutamine to citrulline conversion pathway

Enterocyte is one of the main sites of amino acids metabolism and particularly of the citrulline biosynthesis. Working at the cellular scale and applying ordinary differential equations (ODEs) formalism, we have built a mathematical model of the enterocytic glutamine to citrulline conversion in the fasting state. This model enables us to test different physiopathological scenarios of enzyme activity loss. Results from two different approaches were compared: a standard approach (KA) based on the Michaelis-Menten assumptions and an association-dissociation approach (VH) based on the kinetic mass action law. For both approaches, ODEs system was numerically solved using Mathematica™. In both cases, the model correctly predicts the physiological plasma citrulline steady-state, but the two approaches present clear differences for metabolites of enzymes having a complex mechanism, challenging the validity of the KA approach in such cases. When physiopathological scenarios of enzyme activity loss are simulated, both approaches predict a very sharp transition from the physiological citrulline plasma level to the lack of its production: the concentration profiles of these simulations show a clear threshold of which characteristics vary with the involved enzyme. Moreover, amongst all enzymes included in the model, the ornithine aminotransferase (OAT) shows the highest sensitivity in the system whatever the approach used. This model points out the limits of the Michaelis-Menten approach to model complex enzyme mechanisms. It highlights the key role of OAT in the studied citrulline synthesis pathway and also suggests an order of magnitude about the optimal ratio of enzyme concentrations in this pathway.