Cluster Gauss–Newton and CellNOpt Parameter Estimation in a Small Protein Signaling Network of Vorinostat and Bortezomib Pharmacodynamics
- 7 October 2021
- journal article
- research article
- Published by Springer Science and Business Media LLC in The AAPS Journal
- Vol. 23 (6), 1-15
- https://doi.org/10.1208/s12248-021-00640-7
Abstract
Ordinary differential equation (ODE)–based models of signal transduction pathways often contain parameters that are unidentifiable or unmeasurable by experimental data, and calibrating such models to data remains challenging. Here, two efficient parameter estimation methods, cluster Gauss–Newton (CGN) and CellNOpt (CNO), were applied to fit a signaling network model of U266 multiple myeloma cells to the activity dynamics of key proteins in response to vorinostat and/or bortezomib. A logic-based network model was constructed and transformed to 17 ODEs with 79 parameters estimated within broad ranges of biologically plausible values. The top 10% best-fit parameters by both methods had high uncertainties with CV > 50% for the majority of parameters. The root mean square and prediction errors were comparable without statistically significant differences between the two methods. Despite uncertain parameter estimation, protein dynamics after the sequential combination of bortezomib and vorinostat was predicted with reasonable accuracy and precision. Global sensitivity analyses of partial rank correlation coefficients and Sobol sensitivity demonstrated that apoptosis induction was most sensitive to parameters governing the activity of the proteasome–JNK–caspase-8 axis. Simulations revealed that the greatest magnitude of pharmacodynamic drug interactions between bortezomib and vorinostat occurred at caspase-9, AKT, and Bcl-2. Two sequential combinations were explored in silico, and the outcome matched qualitatively with an empirical evaluation of the pharmacodynamic interaction based on cell viability. Overall, the CGN and CNO algorithms performed similarly for this ODE-based network model calibration, and the calibrated model provided meaningful insights into cellular signaling mechanisms in response to pharmacological perturbations.Keywords
This publication has 41 references indexed in Scilit:
- Preclinical screening of histone deacetylase inhibitors combined with ABT-737, rhTRAIL/MD5-1 or 5-azacytidine using syngeneic Vk*MYC multiple myelomaCell Death & Disease, 2013
- CellNOptR: a flexible toolkit to train protein signaling networks to data using multiple logic formalismsBMC Systems Biology, 2012
- Model-based global sensitivity analysis as applied to identification of anti-cancer drug targets and biomarkers of drug resistance in the ErbB2/3 networkEuropean Journal of Pharmaceutical Sciences, 2012
- Integrating literature-constrained and data-driven inference of signalling networksBioinformatics, 2012
- In Vitro and In Vivo Selective Antitumor Activity of a Novel Orally Bioavailable Proteasome Inhibitor MLN9708 against Multiple Myeloma CellsClinical Cancer Research, 2011
- Odefy - From discrete to continuous modelsBMC Bioinformatics, 2010
- Transforming Boolean models to continuous models: methodology and application to T-cell receptor signalingBMC Systems Biology, 2009
- Input–output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic dataMolecular Systems Biology, 2009
- Universally Sloppy Parameter Sensitivities in Systems Biology ModelsPLoS Computational Biology, 2007
- Experimental Design in Systems Biology, Based on Parameter Sensitivity Analysis Using a Monte Carlo Method: A Case Study for the TNFα-Mediated NF-κ B Signal Transduction PathwaySIMULATION, 2003