Cancer treatment is increasingly moving toward personalized medicine, yet clinicians often still rely on "average" patient responses to guide individual care1, 2. Whether dealing with the extreme unpredictability of pediatric neuroblastoma, which can range from spontaneous regression to fatal progression3, 4, or the diagnostic uncertainty of prostate cancer monitoring via blood tests 5, 6, 7, the need for patient-specific "Digital Twins" is clear.
Our computational team has developed a suite of computational tools that combine multiscale modeling, biomechanics, and artificial intelligence algorithms. By integrating imaging and biochemical biomarkers, we aim to provide a "what-if" platform for testing treatments
D. Sainz-DeMena, M.A. Pérez and J.M. García-Aznar
Read the paperA Digital Twin requires high-quality data. We use Dynamic Contrast-Enhanced MRI (DCE-MRI) to map tumor vascularization, but traditional pharmacokinetic models often ignore the physical diffusion of contrast agents between neighboring voxels, leading to errors. In Sainz-DeMena et al. (2024)4, we introduced Physics-Informed Neural Networks (PINNs) to solve this inverse problem. Unlike standard AI, PINNs are constrained by the physical laws of mass conservation. By embedding these Partial Differential Equations (PDEs) into the neural network’s loss function, we retrieve vascularization parameters (Ktrans) with significantly higher accuracy, even from noisy or incomplete clinical data.
D. Camacho-Gomez, I. Sorzabal-Bellido, C. Ortiz-de-Solorzano, J.M. Garcia-Aznar, M.J. Gomez-Benito Read the paper
How do thousands of cells self-organize? In Camacho-Gomez et al. (2023)3, we developed a hybrid framework where a Deep Learning algorithm manages the decision-making of individual agents (cells) in a 3D model.
Instead of fixed rules, the neural network evaluates the simulation in real-time and "evaluates" if a cell should proliferate, stay quiescent, or secrete fluid to form a lumen. This model successfully reproduced the complex morphogenetic patterns of tumor organoids, proving that AI can help us unravel the underlying principles of how cancer cells self-organize.
C. Borau, K. Y. Wertheim, S. Hervas-Raluy, D. Sainz-DeMena, D. Walker, R. Chisholm, P. Richmond, V. Varella, M. Viceconti, A. Montero, E. Gregori-Puigjané, J. Mestres, M. Kasztelnik, J. M. García-Aznar Read the paper
In Borau et al. (2023)1, we presented a massive multiscale orchestrated framework. This "orchestrator" integrates four distinct modules: macro-scale transport, a subcellular machine learning model of drug-protein interactions, a GPU-accelerated agent-based model (ABM), and a macro biomechanical model. This infrastructure allows clinicians to run an 80-day treatment simulation in less than 12 hours on High-Performance Computing (HPC) clusters.
S. Hervas-Raluy, D. Sainz-DeMena, M. J. Gomez-Benito and J. M. García-Aznar Read the paper
In Hervas-Raluy et al. (2024)5, we turned MRI scans into biomechanical Finite Element models. By feeding the models image-based biomarkers like cellularity (from ADC maps) and vascularization (Ktrans), our simulations predicted the dramatic 90% volume reduction in a responsive patient versus the minimal 15% reduction in a resistant one. This provides a window into "what-if" scenarios, such as the potential impact of proangiogenic treatments to improve drug delivery.
D. Camacho-Gomez, C. Borau, J. M. Garcia-Aznar, M. J. Gomez-Benito, M. Girolami and Maria Angeles Perez Read the paper
Moving from the mechanical to the biochemical, our most recent work in Camacho-Gomez et al. (2025)2 addresses prostate cancer monitoring. Current methods often fail to detect tumor growth if Prostate-Specific Antigen (PSA) levels don't rise significantly, a phenomenon known as "hidden growth"2,5,6. We developed a Physics-Informed Machine Learning (PIML) digital twin that reconstructs a patient's tumor growth trajectory directly from PSA blood tests2. This model is geometry-agnostic; it uses voxelized representations that allow the same neural network architecture to be applied across different patients without retraining. In clinical tests, the twin reconstructed 2.5 years of tumor evolution with errors as low as 0.8%, revealing growth even in stable-PSA scenarios5.
These seven works represent a shift toward personalized computational oncology. By bridging the gap between clinical data (MRI, PSA), physics-based modeling, and deep learning, we are moving closer to a future where every patient’s treatmentis guided by a digital simulation of their unique disease.
© 2026 - The Mathematical Oncology Blog