BioDiscovery2012; 4: 4; DOI: 10.7750/BioDiscovery.2012.4.4
Copyright: © 2012 Jackson. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, provided the original authors and source are credited.
Received: October 25, 2012; Accepted: October 30, 2012; Available online /Published: October 31, 2012
Keywords: apoe-/- mouse; atherosclerosis; computational pharmacology; disease modelling; systems pharmacology
*Corresponding Author: Robert C. Jackson, e-mail: firstname.lastname@example.org
Conflict of Interests: No potential conflict of interest was disclosed.
The Chan et al. model, based upon the Apoe-/- mouse, includes elements of cholesterol and macrophage trafficking, inflammation, oxidative stress, endothelial function, and thrombosis. It has the ability to predict relationships between biomarker data, pharmacodynamic effects and clinical outcomes. The model is the outcome of a collaboration between Entelos, an in silico modelling and simulation company, and Philip Morris, a tobacco company. A primary motive for developing the model appears to have been a desire to explore the relationship between smoking (and smoking cessation) and heart disease. However, the scope of the model is broad enough to enable it to be used to model the effects of other lifestyle factors, including diet, and of drug treatment. As an example of drug effects, the model is used to predict effects of ezetimibe (which blocks cholesterol absorption from the intestine) on atherosclerotic progression. By publishing the model in Biodiscovery, the authors have agreed to make it available, free of charge, to all researchers.
In the past, drug developers have regarded computational models of complex biological systems with great scepticism. There has been an impression that the systems involved are so complex that any attempt to describe them mathematically must involve simplifying assumptions that were likely to undermine the dynamics of the system being modelled. There was pessimism about the ability to validate such large models (the Chan et al. model contains 94 ordinary differential equations, 524 algebraic equations, and 3,508 parameters). Yet these same drug developers have for eighty years relied heavily upon pharmacokinetic (PK) models in making drug development decisions, and PK models make equally sweeping simplifications. Why the difference? Two reasons: the first is that PK is a generic technology. We can use the same analytical methods to measure plasma levels of an antihypertensive as we use for an antidepressant, whereas the biological or pharmacodynamic (PD) endpoints that we use in disease modelling are different for every therapeutic area and for every drug class. Secondly, the mathematics of disease modelling and PD modelling is more complex. Both these factors still present barriers to wider use of disease modelling, but the barriers are yielding to advances in technology. Development of prognostic and pharmacodynamic biomarkers (and to an increasing extent, whole-body imaging techniques) is making it easier to collect the data required to validate complex models. In a recent review  I discussed the use of PD models of biomarker data in oncology. At present, regulatory approval of anticancer drugs requires clinical endpoints, which for slowly-progressing tumours means drug development times often in excess of ten years. Demonstrating and validating the predictive power of PD biomarkers in the context of computational disease models has the potential to revolutionise clinical drug development.
In recent years, computational disease models have been published covering a wide range of therapeutic areas. In HIV disease, we are presented with a complex interactive system where the immune system attacks the virus, and the virus attacks the immune system. Predicting PD effects of antiviral drugs requires a model that captures these complexities. Because of the very high mutation rates of retroviruses, in the early days of anti-retroviral drug development there was pessimism about whether sustained responses could be achieved in the face of acquired drug resistance. A disease model predicted, correctly, that with the use of multi-drug combinations, and with sufficient treatment intensity, disease progression could be arrested for many years . In oncology, models of the cancer cell cycle  and three-dimensional virtual tumours have been used to predict optimal drug combination schedules . Disease models have been described for diabetes, rheumatoid arthritis, hypertension and skin ageing .
These disease modelling approaches, formerly the province of theoreticians, are beginning to interest the wider drug development community. Quintiles, a major clinical research organization, recently published an extensive report on modelling and simulation practice throughout the drug development process . In this report they note the increasing acceptance by the US Food and Drug Administration (FDA) of modelling and simulation studies in support of applications for marketing approval. As an example, they cite the FDA’s Office of Clinical Pharmacology recommending concentration-QT modelling as a means of evaluating drug potential for QT interval disturbance . Models such as that of Chan et al. are advancing the technology of computational disease modelling. In the long term, they should make possible higher success rates in clinical development, design of rational combination therapies, and tailoring of clinical protocols to individual patients. This is an ambitious agenda, but this paper represents a step in the right direction.
- Chan JR, Vuillaume G, Bever C, Lebrun S, Lietz M, Steffen Y, et al. The Apoe-/- Mouse PhysioLab® Platform: A Validated Physiologically-based Mathematical Model of Atherosclerotic Plaque Progression in the Apoe-/- Mouse. Biodiscovery 2012; 3:2.
- Jackson RC. Pharmacodynamic modelling of biomarker data in oncology. ISRN Pharmacology2012; 2012:590626.
- Jackson RC. A pharmacokinetic-pharmacodynamic model of chemotherapy of human immunodeficiency virus infection that relates development of resistance to treatment intensity. J Pharmacokinet Biopharm1997; 25:713-730.
- Chassagnole C, Jackson RC, Hussain N, Bashir L, Derow C, Savin J, et al. Using a mammalian cell cycle simulation to interpret differential kinase inhibition in anti-tumour pharmaceutical development. BioSystems2006; 83:91-97.
- Orrell D, Fernandez E, Cronier D, Gelbert LM, Fell DA, De Alwis D, et al. Predicting the effect of combination schedules on xenograft tumor using the Virtual Tumor. Proc Amer Assoc Cancer Res2011; (abstract 4918).
- O’Kelly M, Berry S, Walp D, Garrett A. Lifecycle Modeling and Simulation. WWW.QUINTILES.COM (2012).
- Garnett CE, Beasley N, Bharraram VA, Jadhav PR, Madabushi R, Stockbridge N, et al. Concentration-QT relationships play a key role in the evaluation of proarhthmic risk during regulatory review. J Clin Pharmacol2008; 48: 13-18.