Coping with time scales in disease systems analysis: application to bone remodeling

Network Pharmacology of Adaptogens in the Assessment of Their Pleiotropic Therapeutic Activity

The reductionist concept, based on the ligand–receptor interaction, is not a suitable model for adaptogens, and herbal preparations affect multiple physiological functions, revealing polyvalent pharmacological activities, and are traditionally used in many conditions. This review, for the first time, provides a rationale for the pleiotropic therapeutic efficacy of adaptogens based on evidence from recent gene expression studies in target cells and where the network pharmacology and systems biology approaches were applied. The specific molecular targets and adaptive stress response signaling mechanisms involved in nonspecific modes of action of adaptogens are identified.

Modeling osteoporosis to design and optimize pharmacological therapies comprising multiple drug types

For the treatment of postmenopausal osteoporosis, several drug classes with different mechanisms of action are available. Since only a limited set of dosing regimens and drug combinations can be tested in clinical trials, it is currently unclear whether common medication strategies achieve optimal bone mineral density gains or are outperformed by alternative dosing schemes and combination therapies that have not been explored so far. Here, we develop a mathematical framework of drug interventions for postmenopausal osteoporosis that unifies fundamental mechanisms of bone remodeling and the mechanisms of action of four drug classes: bisphosphonates, parathyroid hormone analogs, sclerostin inhibitors, and receptor activator of NF-κB ligand inhibitors. Using data from several clinical trials, we calibrate and validate the model, demonstrating its predictive capacity for complex medication scenarios, including sequential and parallel drug combinations. Via simulations, we reveal that there is a large potential to improve gains in bone mineral density by exploiting synergistic interactions between different drug classes, without increasing the total amount of drug administered.

Establishment of a Disease-Drug Trial Model for Postmenopausal Osteoporosis: A Zoledronic Acid Case Study

Costly and lengthy clinical trials hinder the development of safe and effective treatments for postmenopausal osteoporosis. To reduce the expense associated with these trials, we established a mechanistic disease-drug trial model for postmenopausal osteoporosis that can predict phase 3 trial outcome based on short-term bone turnover marker data. To this end, we applied a previously developed model for tibolone to bisphosphonates using zoledronic acid as paradigm compound by (1) linking the mechanistic bone cell interaction model to bone turnover markers as well as bone mineral density in lumbar spine and total hip, (2) employing a mechanistic disease progression function, and (3) accounting for zoledronic acid's mechanism of action. Once developed, we fitted the model to clinical trial data of 581 postmenopausal women receiving (1) 5-mg zoledronic acid in year 1 and saline in year 2, (2) 5-mg zoledronic acid in year 1 and year 2, or (3) placebo (saline), calcium (500 mg), and vitamin D (400 IU). All biomarker data was fitted reasonably well, with no apparent bias or model misspecification. Age, years since menopause, and body mass index at baseline were identified as significant covariates. In the future, the model can be modified to explore the link between short-term biomarkers and fracture risk.

Translational Quantitative Systems Pharmacology in Drug Development: from Current Landscape to Good Practices

Systems pharmacology approaches have the capability of quantitatively linking the key biological molecules relevant to a drug candidate's mechanism of action (drug-induced signaling pathways) to the clinical biomarkers associated with the proposed target disease, thereby quantitatively facilitating its development and life cycle management. In this review, the model attributes of published quantitative systems pharmacology (QSP) modeling for lowering cholesterol, treating salt-sensitive hypertension, and treating rare diseases as well as describing bone homeostasis and related pharmacological effects are critically reviewed with respect to model quality, calibration, validation, and performance. We further reviewed the common practices in optimizing QSP modeling and prediction. Notably, leveraging genetics and genomic studies for model calibration and validation is common. Statistical and quantitative assessment of QSP prediction and handling of model uncertainty are, however, mostly lacking as are the quantitative and statistical criteria for assessing QSP predictions and the covariance matrix of coefficients between the parameters in a validated virtual population. To accelerate advances and application of QSP with consistent quality, a list of key questions is proposed to be addressed when assessing the quality of a QSP model in hopes of stimulating the scientific community to set common expectations. The common expectations as to what constitutes the best QSP modeling practices, which the scientific community supports, will advance QSP modeling in the realm of informed drug development. In the long run, good practices will extend the life cycles of QSP models beyond the life cycles of individual drugs.

Pharmacology of bisphosphonates

The biological effects of the bisphosphonates (BPs) as inhibitors of calcification and bone resorption were first described in the late 1960s. In the 50 years that have elapsed since then, the BPs have become the leading drugs for the treatment of skeletal disorders characterized by increased bone resorption, including Paget's disease of bone, bone metastases, multiple myeloma, osteoporosis and several childhood inherited disorders. The discovery and development of the BPs as a major class of drugs for the treatment of bone diseases is a paradigm for the successful journey from "bench to bedside and back again". Several of the leading BPs achieved "blockbuster" status as branded drugs. However, these BPs have now come to the end of their patent life, making them highly affordable. The opportunity for new clinical applications for BPs also exists in other areas of medicine such as ageing, cardiovascular disease and radiation protection. Their use as inexpensive generic medicines is therefore likely to continue for many years to come. Fifty years of research into the pharmacology of bisphosphonates have led to a fairly good understanding about how these drugs work and how they can be used safely in patients with metabolic bone diseases. However, while we seemingly know much about these drugs, a number of key aspects related to BP distribution and action remain incompletely understood. This review summarizes the existing knowledge of the (pre)clinical and translational pharmacology of BPs, and highlights areas in which understanding is lacking.

Dynamics of Bone Cell Interactions and Differential Responses to PTH and Antibody-Based Therapies

We propose a mathematical model describing the dynamics of osteoblasts and osteoclasts in bone remodeling. The goal of this work is to develop an integrated modeling framework for bone remodeling and bone cell signaling dynamics that could be used to explore qualitatively combination treatments for osteoporosis in humans. The model has been calibrated using 57 checks from the literature. Specific global optimization methods based on qualitative objectives have been developed to perform the model calibration. We also added pharmacokinetics representations of three drugs to the model, which are teriparatide (PTH(1-34)), denosumab (a RANKL antibody) and romosozumab (a sclerostin antibody), achieving excellent goodness-of-fit of human clinical data. The model reproduces the paradoxical effects of PTH on the bone mass, where continuous administration of PTH results in bone loss but intermittent administration of PTH leads to bone gain, thus proposing an explanation of this phenomenon. We used the model to simulate different categories of osteoporosis. The main attributes of each disease are qualitatively well captured by the model, for example changes in bone turnover in the disease states. We explored dosing regimens for each disease based on the combination of denosumab and romosozumab, identifying adequate ratios and doses of both drugs for subpopulations of patients in function of categories of osteoporosis and the degree of severity of the disease.

Regulatory Perspectives in Pharmacometric Models of Osteoporosis

Osteoporosis is a disorder of the bones in which they are weakened to the extent that they become more prone to fracture. There are various forms of osteoporosis: some of them are induced by drugs, and others occur as a chronic progressive disorder as an individual gets older. As the median age of the population rises across the world, the chronic form of the bone disease is drawing attention as an important worldwide health issue. Developing new treatments for osteoporosis and comparing them with existing treatments are complicated processes due to current acceptance by regulatory authorities of bone mineral density (BMD) and fracture risk as clinical end points, which require clinical trials to be large, prolonged, and expensive to determine clinically significant impacts in BMD and fracture risk. Moreover, changes in BMD and fracture risk are not always correlated, with some clinical trials showing BMD improvement without a reduction in fractures. More recently, bone turnover markers specific to bone formation and resorption have been recognized that reflect bone physiology at a cellular level. These bone turnover markers change faster than BMD and fracture risk, and mathematically linking the biomarkers via a computational model to BMD and/or fracture risk may help in predicting BMD and fracture risk changes over time during the progression of a disease or when under treatment. Here, we discuss important concepts of bone physiology, osteoporosis, treatment options, mathematical modeling of osteoporosis, and the use of these models by the pharmaceutical industry and the Food and Drug Administration.

Sobol Sensitivity Analysis: A Tool to Guide the Development and Evaluation of Systems Pharmacology Models

A systems pharmacology model typically integrates pharmacokinetic, biochemical network, and systems biology concepts into a unifying approach. It typically consists of a large number of parameters and reaction species that are interlinked based upon the underlying (patho)physiology and the mechanism of drug action. The more complex these models are, the greater the challenge of reliably identifying and estimating respective model parameters. Global sensitivity analysis provides an innovative tool that can meet this challenge. CPT Pharmacometrics Syst. Pharmacol. (2015) 4, 69-79; doi:10.1002/psp4.6; published online 25 February 2015.

QSP Toolbox: Computational Implementation of Integrated Workflow Components for Deploying Multi-Scale Mechanistic Models

Quantitative systems pharmacology (QSP) modeling has become increasingly important in pharmaceutical research and development, and is a powerful tool to gain mechanistic insights into the complex dynamics of biological systems in response to drug treatment. However, even once a suitable mathematical framework to describe the pathophysiology and mechanisms of interest is established, final model calibration and the exploration of variability can be challenging and time consuming. QSP models are often formulated as multi-scale, multi-compartment nonlinear systems of ordinary differential equations. Commonly accepted modeling strategies, workflows, and tools have promise to greatly improve the efficiency of QSP methods and improve productivity. In this paper, we present the QSP Toolbox, a set of functions, structure array conventions, and class definitions that computationally implement critical elements of QSP workflows including data integration, model calibration, and variability exploration. We present the application of the toolbox to an ordinary differential equations-based model for antibody drug conjugates. As opposed to a single stepwise reference model calibration, the toolbox also facilitates simultaneous parameter optimization and variation across multiple in vitro, in vivo, and clinical assays to more comprehensively generate alternate mechanistic hypotheses that are in quantitative agreement with available data. The toolbox also includes scripts for developing and applying virtual populations to mechanistic exploration of biomarkers and efficacy. We anticipate that the QSP Toolbox will be a useful resource that will facilitate implementation, evaluation, and sharing of new methodologies in a common framework that will greatly benefit the community.

Disease Systems Analysis of Bone Mineral Density and Bone Turnover Markers in Response to Alendronate, Placebo, and Washout in Postmenopausal Women

A previously established mechanism-based disease systems model for osteoporosis that is based on a mathematically reduced version of a model describing the interactions between osteoclast (bone removing) and osteoblast (bone forming) cells in bone remodeling has been applied to clinical data from women (n = 1,379) receiving different doses and treatment regimens of alendronate, placebo, and washout. The changes in the biomarkers, plasma bone-specific alkaline phosphatase activity (BSAP), urinary N-telopeptide (NTX), lumbar spine bone mineral density (BMD), and total hip BMD, were linked to the underlying mechanistic core of the model. The final model gave an accurate description of all four biomarkers for the different treatments. Simulations were used to visualize the dynamics of the underlying network and the natural disease progression upon alendronate treatment and discontinuation. These results complement the previous applications of this mechanism-based disease systems model to data from various treatments for osteoporosis.

Understanding the Behavior of Systems Pharmacology Models Using Mathematical Analysis of Differential Equations: Prolactin Modeling as a Case Study

In this tutorial, we introduce basic concepts in dynamical systems analysis, such as phase‐planes, stability, and bifurcation theory, useful for dissecting the behavior of complex and nonlinear models. A precursor‐pool model with positive feedback is used to demonstrate the power of mathematical analysis. This model is nonlinear and exhibits multiple steady states, the stability of which is analyzed. The analysis offers insight into model behavior and suggests useful parameter regions, which simulations alone could not.

Systems pharmacology - Towards the modeling of network interactions

Mechanism-based pharmacokinetic and pharmacodynamics (PKPD) and disease system (DS) models have been introduced in drug discovery and development research, to predict in a quantitative manner the effect of drug treatment in vivo in health and disease. This requires consideration of several fundamental properties of biological systems behavior including: hysteresis, non-linearity, variability, interdependency, convergence, resilience, and multi-stationarity. Classical physiology-based PKPD models consider linear transduction pathways, connecting processes on the causal path between drug administration and effect, as the basis of drug action. Depending on the drug and its biological target, such models may contain expressions to characterize i) the disposition and the target site distribution kinetics of the drug under investigation, ii) the kinetics of target binding and activation and iii) the kinetics of transduction. When connected to physiology-based DS models, PKPD models can characterize the effect on disease progression in a mechanistic manner. These models have been found useful to characterize hysteresis and non-linearity, yet they fail to explain the effects of the other fundamental properties of biological systems behavior. Recently systems pharmacology has been introduced as novel approach to predict in vivo drug effects, in which biological networks rather than single transduction pathways are considered as the basis of drug action and disease progression. These models contain expressions to characterize the functional interactions within a biological network. Such interactions are relevant when drugs act at multiple targets in the network or when homeostatic feedback mechanisms are operative. As a result systems pharmacology models are particularly useful to describe complex patterns of drug action (i.e. synergy, oscillatory behavior) and disease progression (i.e. episodic disorders). In this contribution it is shown how physiology-based PKPD and disease models can be extended to account for internal systems interactions. It is demonstrated how SP models can be used to predict the effects of multi-target interactions and of homeostatic feedback on the pharmacological response. In addition it is shown how DS models may be used to distinguish symptomatic from disease modifying effects and to predict the long term effects on disease progression, from short term biomarker responses. It is concluded that incorporation of expressions to describe the interactions in biological network analysis opens new avenues to the understanding of the effects of drug treatment on the fundamental aspects of biological systems behavior.

A Six-Stage Workflow for Robust Application of Systems Pharmacology

Quantitative and systems pharmacology (QSP) is increasingly being applied in pharmaceutical research and development. One factor critical to the ultimate success of QSP is the establishment of commonly accepted language, technical criteria, and workflows. We propose an integrated workflow that bridges conceptual objectives with underlying technical detail to support the execution, communication, and evaluation of QSP projects.

Top-down and bottom-up modeling in system pharmacology to understand clinical efficacy: An example with NRTIs of HIV-1

A major aim of Systems Pharmacology is to understand clinically relevant mechanisms of action (MOA) of drugs and to use this knowledge in order to optimize therapy. To enable this mission it is necessary to obtain knowledge on how in vitro testable insights translate into clinical efficacy. Mathematical modeling and data integration are essential components to achieve this goal. Two modeling philosophies are prevalent, each of which in isolation is not sufficient to achieve the above described: In a 'top-down' approach, a minimal pharmacokinetic-pharmacodynamic (PK-PD) model is derived from- and fitted to available clinical data. This model may lack interpretability in terms of mechanisms and may only be predictive for scenarios already covered by the data used to derive it. A 'bottom-up' approach builds on mechanistic insights derived from in vitro/ex vivo experiments, which can be conducted under controlled conditions, but may not be fully representative for the in vivo/clinical situation. In this work, we employ both approaches side-by-side to predict the clinical potency (IC₅₀ values) of the nucleoside reverse transcriptase inhibitors (NRTIs) lamivudine, emtricitabine and tenofovir. In the 'top-down' approach, this requires to establish the dynamic link between the intracellularly active NRTI-triphosphates (which exert the effect) and plasma prodrug PK and to subsequently link this composite PK model to viral kinetics. The 'bottom-up' approach assesses inhibition of reverse transcriptase-mediated viral DNA polymerization by the intracellular, active NRTI-triphosphates, which has to be brought into the context of target cell infection. By using entirely disparate sets of data to derive and parameterize the respective models, our approach serves as a means to assess the clinical relevance of the 'bottom-up' approach. We obtain very good qualitative and quantitative agreement between 'top-down' vs. 'bottom-up' predicted IC₅₀ values, arguing for the validity of the 'bottom-up' approach. We noted, however, that the 'top-down' approach is strongly dependent on the sparse and noisy intracellular pharmacokinetic data. All in all, our work provides confidence that we can translate in vitro parameters into measures of clinical efficacy using the 'bottom-up' approach. This may allow to infer the potency of various NRTIs in inhibiting e.g. mutant viruses, to distinguish sources of interaction of NRTI combinations and to assess the efficacy of different NRTIs for repurposing, e.g. for pre-exposure prophylaxis.

Kinetics of drug action in disease states: towards physiology-based pharmacodynamic (PBPD) models

Gerhard Levy started his investigations on the "Kinetics of Drug Action in Disease States" in the fall of 1980. The objective of his research was to study inter-individual variation in pharmacodynamics. To this end, theoretical concepts and experimental approaches were introduced, which enabled assessment of the changes in pharmacodynamics per se, while excluding or accounting for the cofounding effects of concomitant changes in pharmacokinetics. These concepts were applied in several studies. The results, which were published in 45 papers in the years 1984-1994, showed considerable variation in pharmacodynamics. These initial studies on kinetics of drug action in disease states triggered further experimental research on the relations between pharmacokinetics and pharmacodynamics. Together with the concepts in Levy's earlier publications "Kinetics of Pharmacologic Effects" (Clin Pharmacol Ther 7(3): 362-372, 1966) and "Kinetics of pharmacologic effects in man: the anticoagulant action of warfarin" (Clin Pharmacol Ther 10(1): 22-35, 1969), they form a significant impulse to the development of physiology-based pharmacodynamic (PBPD) modeling as novel discipline in the pharmaceutical sciences. This paper reviews Levy's research on the "Kinetics of Drug Action in Disease States". Next it addresses the significance of his research for the evolution of PBPD modeling as a scientific discipline. PBPD models contain specific expressions to characterize in a strictly quantitative manner processes on the causal path between exposure (in terms of concentration at the target site) and the drug effect (in terms of the change in biological function). Pertinent processes on the causal path are: (1) target site distribution, (2) target binding and activation and (3) transduction and homeostatic feedback.

Application of a Systems Pharmacology-Based Placebo Population Model to Analyze Long-Term Data of Postmenopausal Osteoporosis

Osteoporosis is a progressive bone disease characterized by decreased bone mass resulting in increased fracture risk. The objective of this investigation was to test whether a recently developed disease systems analysis model for osteoporosis could describe disease progression in a placebo-treated population from the Early Postmenopausal Intervention Cohort (EPIC) study. First, we qualified the model using a subset from the placebo arm of the EPIC study of 222 women who had similar demographic characteristics as the 149 women from the placebo arm of the original population. Second, we applied the model to all 470 women. Bone mineral density (BMD) dynamics were changed to an indirect response model to describe lumbar spine and total hip BMD in this second population. This updated disease systems analysis placebo model describes the dynamics of all biomarkers in the corresponding datasets to a very good approximation; a good description of an individual placebo response will be valuable for evaluating treatments for osteoporosis.

Modeling effects of dexamethasone on disease progression of bone mineral density in collagen-induced arthritic rats

A mechanism-based model was developed to characterize the crosstalk between proinflammatory cytokines, bone remodeling biomarkers, and bone mineral density (BMD) in collagen-induced arthritic (CIA) rats. Male Lewis rats were divided into five groups: healthy control, CIA control, CIA receiving single 0.225 mg kg⁻¹ subcutaneous (SC) dexamethasone (DEX), CIA receiving single 2.25 mg kg⁻¹ SC DEX, and CIA receiving chronic 0.225 mg kg⁻¹ SC DEX. The CIA rats underwent collagen induction at day 0 and DEX was injected at day 21 post-induction. Disease activity was monitored throughout the study and rats were sacrificed at different time points for blood and paw collection. Protein concentrations of interleukin (IL)-1β, IL-6, receptor activator of nuclear factor kappa-B ligand (RANKL), osteoprotegerin (OPG), and tartrate-resistant acid phosphatase 5b (TRACP-5b) in paws were measured by enzyme-linked immunosorbent assays (ELISA). Disease progression and DEX pharmacodynamic profiles of IL-1β, IL-6, RANKL, and OPG were fitted simultaneously and parameters were sequentially applied to fit the TRACP-5b and BMD data. The model was built according to the mechanisms reported in the literature and modeling was performed using ADAPT 5 software with naïve pooling. Time profiles of IL-1β and IL-6 protein concentrations correlated with their mRNAs. The RANKL and OPG profiles matched previous findings in CIA rats. DEX inhibited the expressions of IL-1β, IL-6, and RANKL, but did not alter OPG. TRACP-5b was also inhibited by DEX. Model predictions suggested that anti-IL-1β therapy and anti-RANKL therapy would result in similar efficacy for prevention of bone loss among the cytokine antagonists.

A semi-mechanistic model of bone mineral density and bone turnover based on a circular model of bone remodeling

Development of novel therapies for bone diseases can benefit from mathematical models that predict drug effect on bone remodeling biomarkers. Therefore, a bone cycle model (BCM) was developed that takes into consideration the concept of the basic multicellular unit and the dynamic equilibrium of bone remodeling. The model is a closed form cyclical model with four compartments representing resorption, formation, primary mineralization, and secondary mineralization. Equations describing the time course of bone turnover biomarkers were developed using the flow rate of bone cycle units (BCU) between the compartments or the amount of BCU in each compartment. A disease progression model representing bone loss in osteoporosis, a vitamin D and calcium supplementation (placebo) model, and a drug model for antiresorptive treatments were added to the model. Initial model parameter values were derived from published bone turnover data. The BCM accurately described biomarker-time profiles in postmenopausal women receiving either placebo or bisphosphonate treatment. The slow continual increase in bone mineral density (BMD) observed after 1 year of treatment was accurately described when changes in bone turnover were combined with increases in mineralization. For this purpose, the secondary mineralization compartment was replaced by three catenary chain compartments representing increasing mineral content. The refined BCM satisfactorily predicted biomarker profiles after long-term (10-year) bisphosphonate treatment. Furthermore, the model successfully described individual bone turnover markers and BMD results following treatment with denosumab in postmenopausal women. Analyses with this model could be used to optimize dosing regimens and to predict effects of novel osteoporotic treatments.

Systems approaches in risk assessment

Adverse drug events (ADEs) remain a universal problem in drug development, regulatory review, and clinical practice with a substantial financial burden on the global health-care system. Recent advances in molecular and "omics" technologies, along with online databases and bioinformatics, have enabled a more integrative approach to understanding drug-target (protein) interactions, both desirable and undesirable, within a biological system. This has led to the development of systems approaches to risk assessment in an attempt to complement and improve on contemporary observational and predictive strategies for assessing risk. Although still in an evolutionary phase, systems approaches have the potential to markedly advance our understanding of ADEs and ability to predict them. Systems approaches will also move personalized medicine forward by enabling better identification of individual and subgroup risk factors.

Integrated pharmacometrics and systems pharmacology model-based analyses to guide GnRH receptor modulator development for management of endometriosis

Endometriosis is a gynecological condition resulting from proliferation of endometrial-like tissue outside the endometrial cavity. Estrogen suppression therapies, mediated through gonadotropin-releasing hormone (GnRH) modulation, decrease endometriotic implants and diminish associated pain albeit at the expense of bone mineral density (BMD) loss. Our goal was to provide model-based guidance for GnRH-modulating clinical programs intended for endometriosis management. This included developing an estrogen suppression target expected to provide symptomatic relief with minimal BMD loss and to evaluate end points and study durations supportive of efficient development decisions. An existing multiscale model of calcium and bone was adapted to include systematic estrogen pharmacologic effects to describe estrogen concentration-related effects on BMD. A logistic regression fit to patient-level data from three clinical GnRH agonist (nafarelin) studies described the relationship of estrogen with endometrial-related pain. Targeting estradiol between 20 and 40 pg/ml was predicted to provide efficacious endometrial pain response while minimizing BMD effects.

Merging systems biology with pharmacodynamics

The emerging discipline of systems pharmacology aims to combine analysis and computational modeling of cellular regulatory networks with quantitative pharmacology approaches to drive the drug discovery processes, predict rare adverse events, and catalyze the practice of personalized precision medicine. Here, we introduce the concept of enhanced pharmacodynamic (ePD) models, which synergistically combine the desirable features of systems biology and current PD models within the framework of ordinary or partial differential equations. ePD models that analyze regulatory networks involved in drug action can account for a drug's multiple targets and for the effects of genomic, epigenomic, and posttranslational changes on the drug efficacy. This new knowledge can drive drug discovery and shape precision medicine.

Pharmacokinetic/pharmacodynamic modeling in inflammation

Inflammation is an array of immune responses to infection and injury. It results from a complex immune cascade and is the basis of many chronic diseases such as arthritis, diabetes, and cancer. Numerous mathematical models have been developed to describe the disease progression and effects of anti-inflammatory drugs. This review illustrates the state of the art in modeling the effects of diverse drugs for treating inflammation, describes relevant biomarkers amenable to modeling, and summarizes major advantages and limitations of the published pharmacokinetic/ pharmacodynamic (PK/PD) models. Simple direct inhibitory models are often used to describe in vitro effects of anti-inflammatory drugs. Indirect response models are more mechanism based and have been widely applied to the turnover of symptoms and biomarkers. These, along with target-mediated and transduction models, have been successfully applied to capture the PK/PD of many anti-inflammatory drugs and describe disease progression of inflammation. Biologics have offered opportunities to address specific mechanisms of action, and evolve small systems models to quantitatively capture the underlying physiological processes. More advanced mechanistic models should allow evaluation of the roles of some key mediators in disease progression, assess drug interactions, and better translate drug properties from in vitro and animal data to patients.