Modeling cholesterol metabolism by gene expression profiling in the hippocampus

An important part of the challenge of building models of biochemical reactions is determining reaction rate constants that transform substrates into products. We present a method to derive enzymatic kinetic values from mRNA expression levels for modeling biological networks without requiring further tuning. The core metabolic reactions of cholesterol in the brain, particularly in the hippocampus, were simulated. To build the model the baseline mRNA expression levels of genes involved in cholesterol metabolism were obtained from the Allen Mouse Brain Atlas. The model is capable of replicating the trends of relative cholesterol levels in Alzheimer's and Huntington's diseases; and reliably simulated SLOS, desmosterolosis, and Dhcr14/Lbr knockout studies. A sensitivity analysis correctly uncovers the Hmgcr, Idi2 and Fdft1 sites that regulate cholesterol homeostasis. Overall, our model and methodology can be used to pinpoint key reactions, which, upon manipulation, may predict altered cholesterol levels and reveal insights into potential drug therapy targets under diseased conditions.