Mechanism of inhibition defines CETP activity: a mathematical model for CETP in vitro.

Because cholesteryl ester transfer protein (CETP) inhibition is a potential HDL-raising therapy, interest has been raised in the mechanisms and consequences of CETP activity. To explore these mechanisms and the dynamics of CETP in vitro, a mechanistic mathematical model was developed based upon the shuttle mechanism for lipid transfer. Model parameters were estimated from eight published experimental datasets, and the resulting model captures observed dynamics of CETP in vitro. Simulations suggest the shuttle mechanism yields behaviors consistent with experimental observations. Three key findings predicted from model simulations are: 1) net CE transfer activity from HDL to VLDL and LDL can be significantly altered by changing the balance of homoexchange versus heteroexchange of neutral lipids via CETP; 2) lipemia-induced increases in CETP activity are more likely caused by increases in lipoprotein particle size than particle number; and 3) the inhibition mechanisms of the CETP inhibitors torcetrapib and JTT-705 are significantly more potent than a classic competitive inhibition mechanism with the irreversible binding mechanism having the most robust response. In summary, the model provides a plausible representation of CETP activity in vitro, corroborates strong evidence for the shuttle hypothesis, and provides new insights into the consequences of CETP activity and inhibition on lipoproteins.

Perspectives on mathematical modeling for receptor-mediated processes.

Mathematical modeling is a potent in silico tool that can help investigate, interpret, and predict the behavior of biological systems. The first step is to develop a working hypothesis of the biology. Then by "translating" the biological phenomena into equations, models can harness the power of mathematical analysis techniques to explore the dynamics and interactions of the biological components. Models can be used together with traditional experimental models to help design new experiments, test hypotheses, identify mechanisms, and predict outcomes. This article reviews the process of building, calibrating, and using mathematical models in the context of the kinetics of receptor and signal transduction biology. An example model related to the androgen receptor-mediated regulation of the prostate is presented to illustrate the steps in the modeling process and to highlight the potential for mathematical modeling in this area.

Response to continuous and pulsatile PTH dosing: a mathematical model for parathyroid hormone receptor kinetics.

In this paper, we propose a mathematical model for parathyroid hormone receptor (PTH1R) kinetics, focusing on the receptor's response to PTH dosing to discern bone formation responses from bone resorption. The PTH1R is a major target for new osteoporosis treatments, as pulsatile PTH dosing has been shown to induce net bone formation in both animals and humans, and PTH(1-34) was recently FDA approved for the treatment of post-menopausal osteoporosis. PTH has also been shown to cause net bone loss when given continuously, so that the net action of PTH on bone is dependent on the dosing pattern. We have developed a simplified two-state receptor kinetics model for the PTH1R, based on the concepts of Segel et al., to distinguish the activity of active and inactive receptor and receptor-ligand complexes. The goal is to develop a plausible model of the minimal essential biological relationships necessary for understanding the responses to PTH dosing. A two-state model is able to effectively discriminate between continuous and pulsatile PTH dosing using the active species as surrogates for the downstream anabolic response. For continuous PTH dosing, the model predicts a desensitized system dominated by the inactive receptor and complex, consistent with downstream net bone loss that has been demonstrated experimentally. Using pulsatile PTH dosing, the model system predicts a highly sensitized state dominated by the active receptor and complex, corresponding to net bone formation. These results are consistent with the hypothesis that the kinetics of the receptor plays a critical role in the downstream effects of PTH dosing. Moreover, these results indicate that within a range of biologically relevant PTH doses, the two-state model is able to capture the differential behavior of the system for both continuous and pulsatile PTH dosing. The development of such a model provides a rational basis for developing more biologically extensive models that may support the design of optimal dosing strategies for PTH-based anti-osteoporosis treatments. Moreover, this model provides a unique starting point from which to design experiments investigating PTH receptor biology.

Modeling the interactions between osteoblast and osteoclast activities in bone remodeling.

We propose a mathematical model explaining the interactions between osteoblasts and osteoclasts, two cell types specialized in the maintenance of the bone integrity. Bone is a dynamic, living tissue whose structure and shape continuously evolves during life. It has the ability to change architecture by removal of old bone and replacement with newly formed bone in a localized process called remodeling. The model described here is based on the idea that the relative proportions of immature and mature osteoblasts control the degree of osteoclastic activity. In addition, osteoclasts control osteoblasts differentially depending on their stage of differentiation. Despite the tremendous complexity of the bone regulatory system and its fragmentary understanding, we obtain surprisingly good correlations between the model simulations and the experimental observations extracted from the literature. The model results corroborate all behaviors of the bone remodeling system that we have simulated, including the tight coupling between osteoblasts and osteoclasts, the catabolic effect induced by continuous administration of PTH, the catabolic action of RANKL, as well as its reversal by soluble antagonist OPG. The model is also able to simulate metabolic bone diseases such as estrogen deficiency, vitamin D deficiency, senescence and glucocorticoid excess. Conversely, possible routes for therapeutic interventions are tested and evaluated. Our model confirms that anti-resorptive therapies are unable to partially restore bone loss, whereas bone formation therapies yield better results. The model enables us to determine and evaluate potential therapies based on their efficacy. In particular, the model predicts that combinations of anti-resorptive and anabolic therapies provide significant benefits compared with monotherapy, especially for certain type of skeletal disease. Finally, the model clearly indicates that increasing the size of the pool of preosteoblasts is an essential ingredient for the therapeutic manipulation of bone formation. This model was conceived as the first step in a bone turnover modeling platform. These initial modeling results are extremely encouraging and lead us to proceed with additional explorations into bone turnover and skeletal remodeling.