Search

Background: Thrombocytopenia is a common side effect of cytotoxic chemotherapies, which is often dose-limiting. Predicting an individual's risk is of high clinical importance, as otherwise, a small subgroup of patients limits dosages for the overall population for safety reasons., Methods: We aim to predict individual platelet dynamics using non-linear auto-regressive networks with exogenous inputs (NARX). We consider different architectures of the NARX networks, namely feed-forward networks (FNN) and gated recurrent units (GRU). To cope with the relative sparsity of individual patient data, we employ transfer learning (TL) approaches based on a semi-mechanistic model of hematotoxicity. We use a large data set of patients with high-grade non-Hodgkin's lymphoma to learn the respective models on an individual scale and to compare prediction performances with that of the semi-mechanistic model., Results: Of the examined network models, the NARX with GRU architecture performs best. In comparison to the semi-mechanistic model, the network model can result in a substantial improvement of prediction accuracy for patients with irregular dynamics, given well-spaced measurements. TL improves individual prediction performances., Conclusion: NARX networks can be utilized to predict an individual's thrombotoxic response to cytotoxic chemotherapy treatment. For reasonable model learning, we recommend at least three well-spaced measurements per cycle: at baseline, during the nadir phase and during the recovery phase. We aim at generalizing our approach to other treatment scenarios and blood lineages in the future., (© 2024. The Author(s).)

Purpose: Hematotoxicity is a common side-effect of cytotoxic gastrointestinal (GI) cancer therapies. An unsolved problem is to predict the individual risk therefore to decide on treatment adaptions. We applied an established biomathematical prediction model and primarily evaluated its predictive value in patients undergoing chemotherapy for GI cancers in curative intent., Methods: In a prospective, observational multicenter study on patients with gastro-esophageal or pancreatic cancer (n = 28) receiving myelosuppressive adjuvant or neoadjuvant chemotherapy (FLO(T) or FOLFIRINOX), individual model parameters were learned based on patients' observed laboratory values during the first chemotherapy cycle and further external data resources. Grades of hematotoxicity of subsequent cycles were predicted by model simulation and compared with observed data., Results: The most common high-grade hematological toxicity was neutropenia [19/28 patients (68%)]. For the FLO(T) regimen, individual grades of thrombocytopenia and leukopenia could be well predicted for cycles 2-4, as well as grades of neutropenia for cycle 2. Prediction accuracy for neutropenia in the third and fourth cycle differed by one toxicity grade on average. For the FOLFIRINOX-regimen, thrombocytopenia predictions showed a maximum deviation of one toxicity grade up to the end of therapy (8 cycles). Deviations of predictions were less than one degree on average up to cycle 4 for neutropenia, and up to cycle 6 for leukopenia., Conclusion: The biomathematical model showed excellent short-term and decent long-term prediction performance for all relevant hematological side effects associated with FLO(T)/FOLFIRINOX. Clinical utility of this precision-medicine approach needs to be further investigated in a larger cohort., (© 2023. The Author(s).)

Cytotoxic cancer therapy often results in dose-limiting haematotoxic side effects. Predicting an individual's risk is a major objective in precision medicine of cancer treatment. In this regard, patient heterogeneity presents a significant challenge. In this paper, we explore the use of hypothesis-free machine learning models based on recurrent nonlinear auto-regressive networks with exogenous inputs (NARX) as an approach to achieve this goal. Also, we propose a knowledge transfer approach to ameliorate the issue of sparse individual data, which typically hampers learning of individual networks. We demonstrate the feasibility of our approach based on a virtual patient population generated using a semi-mechanistic model of haematopoiesis and imposing different cytotoxic therapy scenarios on it. Employing different techniques of model optimisation, we derive robust and parsimonious individual networks with good generalisation performances. Moreover, we analyse in detail possible factors influencing the generalisation performance. Results suggest that our transfer learning approach using NARX networks can provide robust predictions of individual patient's response to treatment. As a practical perspective, we apply our approach to individual time series data of two patients with non-Hodgkin's lymphoma., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: M. Scholz received funding from Pfizer Inc. for a project not related to this research., (© 2023 The Authors.)

Background: Short-term forecasts of infectious disease burden can contribute to situational awareness and aid capacity planning. Based on best practice in other fields and recent insights in infectious disease epidemiology, one can maximise the predictive performance of such forecasts if multiple models are combined into an ensemble. Here, we report on the performance of ensembles in predicting COVID-19 cases and deaths across Europe between 08 March 2021 and 07 March 2022., Methods: We used open-source tools to develop a public European COVID-19 Forecast Hub. We invited groups globally to contribute weekly forecasts for COVID-19 cases and deaths reported by a standardised source for 32 countries over the next 1-4 weeks. Teams submitted forecasts from March 2021 using standardised quantiles of the predictive distribution. Each week we created an ensemble forecast, where each predictive quantile was calculated as the equally-weighted average (initially the mean and then from 26th July the median) of all individual models' predictive quantiles. We measured the performance of each model using the relative Weighted Interval Score (WIS), comparing models' forecast accuracy relative to all other models. We retrospectively explored alternative methods for ensemble forecasts, including weighted averages based on models' past predictive performance., Results: Over 52 weeks, we collected forecasts from 48 unique models. We evaluated 29 models' forecast scores in comparison to the ensemble model. We found a weekly ensemble had a consistently strong performance across countries over time. Across all horizons and locations, the ensemble performed better on relative WIS than 83% of participating models' forecasts of incident cases (with a total N=886 predictions from 23 unique models), and 91% of participating models' forecasts of deaths (N=763 predictions from 20 models). Across a 1-4 week time horizon, ensemble performance declined with longer forecast periods when forecasting cases, but remained stable over 4 weeks for incident death forecasts. In every forecast across 32 countries, the ensemble outperformed most contributing models when forecasting either cases or deaths, frequently outperforming all of its individual component models. Among several choices of ensemble methods we found that the most influential and best choice was to use a median average of models instead of using the mean, regardless of methods of weighting component forecast models., Conclusions: Our results support the use of combining forecasts from individual models into an ensemble in order to improve predictive performance across epidemiological targets and populations during infectious disease epidemics. Our findings further suggest that median ensemble methods yield better predictive performance more than ones based on means. Our findings also highlight that forecast consumers should place more weight on incident death forecasts than incident case forecasts at forecast horizons greater than 2 weeks., Funding: AA, BH, BL, LWa, MMa, PP, SV funded by National Institutes of Health (NIH) Grant 1R01GM109718, NSF BIG DATA Grant IIS-1633028, NSF Grant No.: OAC-1916805, NSF Expeditions in Computing Grant CCF-1918656, CCF-1917819, NSF RAPID CNS-2028004, NSF RAPID OAC-2027541, US Centers for Disease Control and Prevention 75D30119C05935, a grant from Google, University of Virginia Strategic Investment Fund award number SIF160, Defense Threat Reduction Agency (DTRA) under Contract No. HDTRA1-19-D-0007, and respectively Virginia Dept of Health Grant VDH-21-501-0141, VDH-21-501-0143, VDH-21-501-0147, VDH-21-501-0145, VDH-21-501-0146, VDH-21-501-0142, VDH-21-501-0148. AF, AMa, GL funded by SMIGE - Modelli statistici inferenziali per governare l'epidemia, FISR 2020-Covid-19 I Fase, FISR2020IP-00156, Codice Progetto: PRJ-0695. AM, BK, FD, FR, JK, JN, JZ, KN, MG, MR, MS, RB funded by Ministry of Science and Higher Education of Poland with grant 28/WFSN/2021 to the University of Warsaw. BRe, CPe, JLAz funded by Ministerio de Sanidad/ISCIII. BT, PG funded by PERISCOPE European H2020 project, contract number 101016233. CP, DL, EA, MC, SA funded by European Commission - Directorate-General for Communications Networks, Content and Technology through the contract LC-01485746, and Ministerio de Ciencia, Innovacion y Universidades and FEDER, with the project PGC2018-095456-B-I00. DE., MGu funded by Spanish Ministry of Health / REACT-UE (FEDER). DO, GF, IMi, LC funded by Laboratory Directed Research and Development program of Los Alamos National Laboratory (LANL) under project number 20200700ER. DS, ELR, GG, NGR, NW, YW funded by National Institutes of General Medical Sciences (R35GM119582; the content is solely the responsibility of the authors and does not necessarily represent the official views of NIGMS or the National Institutes of Health). FB, FP funded by InPresa, Lombardy Region, Italy. HG, KS funded by European Centre for Disease Prevention and Control. IV funded by Agencia de Qualitat i Avaluacio Sanitaries de Catalunya (AQuAS) through contract 2021-021OE. JDe, SMo, VP funded by Netzwerk Universitatsmedizin (NUM) project egePan (01KX2021). JPB, SH, TH funded by Federal Ministry of Education and Research (BMBF; grant 05M18SIA). KH, MSc, YKh funded by Project SaxoCOV, funded by the German Free State of Saxony. Presentation of data, model results and simulations also funded by the NFDI4Health Task Force COVID-19 (https://www.nfdi4health.de/task-force-covid-19-2) within the framework of a DFG-project (LO-342/17-1). LP, VE funded by Mathematical and Statistical modelling project (MUNI/A/1615/2020), Online platform for real-time monitoring, analysis and management of epidemic situations (MUNI/11/02202001/2020); VE also supported by RECETOX research infrastructure (Ministry of Education, Youth and Sports of the Czech Republic: LM2018121), the CETOCOEN EXCELLENCE (CZ.02.1.01/0.0/0.0/17-043/0009632), RECETOX RI project (CZ.02.1.01/0.0/0.0/16-013/0001761). NIB funded by Health Protection Research Unit (grant code NIHR200908). SAb, SF funded by Wellcome Trust (210758/Z/18/Z)., Competing Interests: KS, HG, RG, HJ, RN, BP, FS, JD, DW, SA, AU, GG, ER, NR, DS, YW, NW, LW, JT, GO, TS, DT, LP, EK, JM, MB, NL, JM, JS, JW, JF, BL, IR, JB, RW, IV, VT, MS, MZ, CP, BR, NB, SM, LC, GF, IM, DO, PA, AM, VE, AK, DK, LP, BD, ML, SS, JD, SM, VP, GR, BB, GA, NP, GZ, WB, SH, TH, DS, MG, JA, DM, SA, EA, DL, CP, JB, AR, TZ, AK, JZ, FP, FD, MC, GL, PG, BT, FB, GJ, MM, AF, AS, PM, AA, BH, BL, MM, PP, SV, RB, FD, AG, KG, MG, BK, AM, KN, JN, MR, FR, MS, ES, JZ, JK, BP, KH, YK, MS, BP, MB, MF, RI, TK, TO, JB, SF No competing interests declared, PB, HG, SS, BS Affiliated with IEM, Inc. The author has no financial interests to declare

Background: During the COVID-19 pandemic there has been a strong interest in forecasts of the short-term development of epidemiological indicators to inform decision makers. In this study we evaluate probabilistic real-time predictions of confirmed cases and deaths from COVID-19 in Germany and Poland for the period from January through April 2021., Methods: We evaluate probabilistic real-time predictions of confirmed cases and deaths from COVID-19 in Germany and Poland. These were issued by 15 different forecasting models, run by independent research teams. Moreover, we study the performance of combined ensemble forecasts. Evaluation of probabilistic forecasts is based on proper scoring rules, along with interval coverage proportions to assess calibration. The presented work is part of a pre-registered evaluation study., Results: We find that many, though not all, models outperform a simple baseline model up to four weeks ahead for the considered targets. Ensemble methods show very good relative performance. The addressed time period is characterized by rather stable non-pharmaceutical interventions in both countries, making short-term predictions more straightforward than in previous periods. However, major trend changes in reported cases, like the rebound in cases due to the rise of the B.1.1.7 (Alpha) variant in March 2021, prove challenging to predict., Conclusions: Multi-model approaches can help to improve the performance of epidemiological forecasts. However, while death numbers can be predicted with some success based on current case and hospitalization data, predictability of case numbers remains low beyond quite short time horizons. Additional data sources including sequencing and mobility data, which were not extensively used in the present study, may help to improve performance., (© 2022. The Author(s).)

Numerous prediction models of SARS-CoV-2 pandemic were proposed in the past. Unknown parameters of these models are often estimated based on observational data. However, lag in case-reporting, changing testing policy or incompleteness of data lead to biased estimates. Moreover, parametrization is time-dependent due to changing age-structures, emerging virus variants, non-pharmaceutical interventions, and vaccination programs. To cover these aspects, we propose a principled approach to parametrize a SIR-type epidemiologic model by embedding it as a hidden layer into an input-output non-linear dynamical system (IO-NLDS). Observable data are coupled to hidden states of the model by appropriate data models considering possible biases of the data. This includes data issues such as known delays or biases in reporting. We estimate model parameters including their time-dependence by a Bayesian knowledge synthesis process considering parameter ranges derived from external studies as prior information. We applied this approach on a specific SIR-type model and data of Germany and Saxony demonstrating good prediction performances. Our approach can estimate and compare the relative effectiveness of non-pharmaceutical interventions and provide scenarios of the future course of the epidemic under specified conditions. It can be translated to other data sets, i.e., other countries and other SIR-type models.

Disease modelling has had considerable policy impact during the ongoing COVID-19 pandemic, and it is increasingly acknowledged that combining multiple models can improve the reliability of outputs. Here we report insights from ten weeks of collaborative short-term forecasting of COVID-19 in Germany and Poland (12 October-19 December 2020). The study period covers the onset of the second wave in both countries, with tightening non-pharmaceutical interventions (NPIs) and subsequently a decay (Poland) or plateau and renewed increase (Germany) in reported cases. Thirteen independent teams provided probabilistic real-time forecasts of COVID-19 cases and deaths. These were reported for lead times of one to four weeks, with evaluation focused on one- and two-week horizons, which are less affected by changing NPIs. Heterogeneity between forecasts was considerable both in terms of point predictions and forecast spread. Ensemble forecasts showed good relative performance, in particular in terms of coverage, but did not clearly dominate single-model predictions. The study was preregistered and will be followed up in future phases of the pandemic., (© 2021. The Author(s).)

Aims: Thrombocytopenia is a common major side-effect of cytotoxic cancer therapies. A clinically relevant problem is to predict an individual's thrombotoxicity in the next planned chemotherapy cycle in order to decide on treatment adaptation. To support this task, 2 dynamic mathematical models of thrombopoiesis under chemotherapy were proposed, a simple semimechanistic model and a comprehensive mechanistic model. In this study, we assess the performance of these models with respect to existing thrombocytopenia grading schemes., Methods: We consider close-meshed individual time series data of 135 non-Hodgkin's lymphoma patients treated with 6 cycles of CHOP/CHOEP chemotherapies. Individual parameter estimates were derived on the basis of these data considering a varying number of cycles per patient. Parsimony assumptions were applied to optimize parameter identifiability. Models' predictability are assessed by determining deviations of predicted and observed degrees of thrombocytopenia in the next cycles., Results: The mechanistic model results in better agreement of model prediction and individual time series data. Prediction accuracy of future cycle toxicities by the mechanistic model is higher even if the semimechanistic model is provided with data of more cycles for calibration., Conclusion: We successfully established a quantitative and clinically relevant method for assessing prediction performances of biomathematical models of thrombopoiesis under chemotherapy. We showed that the more comprehensive mechanistic model outperforms the semimechanistic model. We aim at implementing the mechanistic model into clinical practice to assess its utility in real life clinical decision-making., (© 2020 The Authors. British Journal of Clinical Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.)

Background: Individualization and patient-specific optimization of treatment is a major goal of modern health care. One way to achieve this goal is the application of high-resolution diagnostics together with the application of targeted therapies. However, the rising number of different treatment modalities also induces new challenges: Whereas randomized clinical trials focus on proving average treatment effects in specific groups of patients, direct conclusions at the individual patient level are problematic. Thus, the identification of the best patient-specific treatment options remains an open question. Systems medicine, specifically mechanistic mathematical models, can substantially support individual treatment optimization. In addition to providing a better general understanding of disease mechanisms and treatment effects, these models allow for an identification of patient-specific parameterizations and, therefore, provide individualized predictions for the effect of different treatment modalities., Results: In the following we describe a software framework that facilitates the integration of mathematical models and computer simulations into routine clinical processes to support decision-making. This is achieved by combining standard data management and data exploration tools, with the generation and visualization of mathematical model predictions for treatment options at an individual patient level., Conclusions: By integrating model results in an audit trail compatible manner into established clinical workflows, our framework has the potential to foster the use of systems-medical approaches in clinical practice. We illustrate the framework application by two use cases from the field of haematological oncology.

Background: Thrombocytopenia is a major side-effect of cytotoxic cancer therapies. The aim of precision medicine is to develop therapy modifications accounting for the individual's risk., Methodology/principle Findings: To solve this task, we develop an individualized bio-mechanistic model of the dynamics of bone marrow thrombopoiesis, circulating platelets and therapy effects thereon. Comprehensive biological knowledge regarding cell differentiation, amplification, apoptosis rates, transition times and corresponding regulations are translated into ordinary differential equations. A model of osteoblast/osteoclast interactions was incorporated to mechanistically describe bone marrow support of quiescent cell stages. Thrombopoietin (TPO) as a major regulator is explicitly modelled including pharmacokinetics and-dynamics of TPO injections. Effects of cytotoxic drugs are modelled by transient depletions of proliferating cells. To calibrate the model, we used population data from the literature and close-meshed individual data of N = 135 high-grade non-Hodgkin's lymphoma patients treated with CHOP-like chemotherapies. To limit the number of free parameters, several parsimony assumptions were derived from biological data and tested via Likelihood methods. Heterogeneity of patients was explained by a few model parameters. The over-fitting issue of individual parameter estimation was successfully dealt with a virtual participation of each patient in population-based experiments. The model qualitatively and quantitatively explains a number of biological observations such as the role of osteoblasts in explaining long-term toxic effects, megakaryocyte-mediated feedback on stem cells, bi-phasic stimulation of thrombopoiesis by TPO, dynamics of megakaryocyte ploidies and non-exponential platelet degradation. Almost all individual time series could be described with high precision. We demonstrated how the model can be used to provide predictions regarding individual therapy adaptations., Conclusions: We propose a mechanistic thrombopoiesis model of unprecedented comprehensiveness in both, biological mechanisms considered and experimental data sets explained. Our innovative method of parameter estimation allows robust determinations of individual parameter settings facilitating the development of individual treatment adaptations during chemotherapy., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following Competing interests: Markus Scholz declared a commercial relationship with Merck Serono with respect to consulting, data analysis and modelling.

Background: Prostate cancer (PCa) is a leading cause of cancer death of men worldwide. In hormone-sensitive prostate cancer (HSPC), androgen deprivation therapy (ADT) is widely used, but an eventual failure on ADT heralds the passage to the castration-resistant prostate cancer (CRPC) stage. Because predicting time to failure on ADT would allow improved planning of personal treatment strategy, we aimed to develop a predictive personalization algorithm for ADT efficacy in HSPC patients., Methods: A mathematical mechanistic model for HSPC progression and treatment was developed based on the underlying disease dynamics (represented by prostate-specific antigen; PSA) as affected by ADT. Following fine-tuning by a dataset of ADT-treated HSPC patients, the model was embedded in an algorithm, which predicts the patient's time to biochemical failure (BF) based on clinical metrics obtained before or early in-treatment., Results: The mechanistic model, including a tumor growth law with a dynamic power and an elaborate ADT-resistance mechanism, successfully retrieved individual time-courses of PSA (R(2) = 0.783). Using the personal Gleason score (GS) and PSA at diagnosis, as well as PSA dynamics from 6 months after ADT onset, and given the full ADT regimen, the personalization algorithm accurately predicted the individual time to BF of ADT in 90% of patients in the retrospective cohort (R(2) = 0.98)., Conclusions: The algorithm we have developed, predicting biochemical failure based on routine clinical tests, could be especially useful for patients destined for short-lived ADT responses and quick progression to CRPC. Prospective studies must validate the utility of the algorithm for clinical decision-making., (© 2015 Wiley Periodicals, Inc.)

Inflammation underlies many diseases and is an undesired effect of several therapy modalities. Biomathematical modeling can help unravel the complex inflammatory processes and the mechanisms triggering their emergence. We developed a model for induction of C-reactive protein (CRP), a clinically reliable marker of inflammation, by interleukin (IL)-11, an approved cytokine for treatment of chemotherapy-induced thrombocytopenia. Due to paucity of information on the mechanisms underlying inflammation-induced CRP dynamics, our model was developed by systematically evaluating several models for their ability to retrieve variable CRP profiles observed in IL-11-treated breast cancer patients. The preliminary semi-mechanistic models were designed by non-linear mixed-effects modeling, and were evaluated by various performance criteria, which test goodness-of-fit, parsimony and uniqueness. The best-performing model, a robust population model with minimal inter-individual variability, uncovers new aspects of inflammation dynamics. It shows that CRP clearance is a nonlinear self-controlled process, indicating an adaptive anti-inflammatory reaction in humans. The model also reveals a dual IL-11 effect on CRP elevation, whereby the drug has not only a potent immediate influence on CRP incline, but also a long-term influence inducing elevated CRP levels for several months. Consistent with this, model simulations suggest that periodic IL-11 therapy may result in prolonged low-grade (chronic) inflammation post treatment. Future application of the model can therefore help design improved IL-11 regimens with minimized long-term CRP toxicity. Our study illuminates the dynamics of inflammation and its control, and provides a prototype for progressive modeling of complex biological processes in the medical realm and beyond.

Despite its great promise, personalized oncology still faces many hurdles, and it is increasingly clear that targeted drugs and molecular biomarkers alone yield only modest clinical benefit. One reason is the complex relationships between biomarkers and the patient's response to drugs, obscuring the true weight of the biomarkers in the overall patient's response. This complexity can be disentangled by computational models that integrate the effects of personal biomarkers into a simulator of drug-patient dynamic interactions, for predicting the clinical outcomes. Several computational tools have been developed for personalized oncology, notably evidence-based tools for simulating pharmacokinetics, Bayesian-estimated tools for predicting survival, etc. We describe representative statistical and mathematical tools, and discuss their merits, shortcomings and preliminary clinical validation attesting to their potential. Yet, the individualization power of mathematical models alone, or statistical models alone, is limited. More accurate and versatile personalization tools can be constructed by a new application of the statistical/mathematical nonlinear mixed effects modeling (NLMEM) approach, which until recently has been used only in drug development. Using these advanced tools, clinical data from patient populations can be integrated with mechanistic models of disease and physiology, for generating personal mathematical models. Upon a more substantial validation in the clinic, this approach will hopefully be applied in personalized clinical trials, P-trials, hence aiding the establishment of personalized medicine within the main stream of clinical oncology., (© 2014 Wiley Periodicals, Inc.)

Interleukin (IL)-21 is an attractive antitumor agent with potent immunomodulatory functions. Yet thus far, the cytokine has yielded only partial responses in solid cancer patients, and conditions for beneficial IL-21 immunotherapy remain elusive. The current work aims to identify clinically-relevant IL-21 regimens with enhanced efficacy, based on mathematical modeling of long-term antitumor responses. For this purpose, pharmacokinetic (PK) and pharmacodynamic (PD) data were acquired from a preclinical study applying systemic IL-21 therapy in murine solid cancers. We developed an integrated disease/PK/PD model for the IL-21 anticancer response, and calibrated it using selected "training" data. The accuracy of the model was verified retrospectively under diverse IL-21 treatment settings, by comparing its predictions to independent "validation" data in melanoma and renal cell carcinoma-challenged mice (R(2)>0.90). Simulations of the verified model surfaced important therapeutic insights: (1) Fractionating the standard daily regimen (50 µg/dose) into a twice daily schedule (25 µg/dose) is advantageous, yielding a significantly lower tumor mass (45% decrease); (2) A low-dose (12 µg/day) regimen exerts a response similar to that obtained under the 50 µg/day treatment, suggestive of an equally efficacious dose with potentially reduced toxicity. Subsequent experiments in melanoma-bearing mice corroborated both of these predictions with high precision (R(2)>0.89), thus validating the model also prospectively in vivo. Thus, the confirmed PK/PD model rationalizes IL-21 therapy, and pinpoints improved clinically-feasible treatment schedules. Our analysis demonstrates the value of employing mathematical modeling and in silico-guided design of solid tumor immunotherapy in the clinic.

We perform critical-point analysis for three-variable systems that represent essential processes of the growth of the angiogenic tumor, namely, tumor growth, vascularization, and generation of angiogenic factor (protein) as a function of effective vessel density. Two models that describe tumor growth depending on vascular mass and regulation of new vessel formation through a key angiogenic factor are explored. The first model is formulated in terms of ODEs, while the second assumes delays in this regulation, thus leading to a system of DDEs. In both models, the only nontrivial critical point is always unstable, while one of the trivial critical points is always stable. The models predict unlimited growth, if the initial condition is close enough to the nontrivial critical point, and this growth may be characterized by oscillations in tumor and vascular mass. We suggest that angiogenesis per se does not suffice for explaining the observed stabilization of vascular tumor size.