Combination chemotherapy with a defined
ratio and sequence of drug
release is a clinically established and effective route to treat advanced
solid tumors. In this context, a growing body of literature demonstrates
the potential of hydrogels constructed with chemically modified polysaccharides
as depots for controlled release of chemotherapeutics. Identifying
the appropriate modification in terms of physicochemical properties
of the functional group and its degree of substitution (χ) to
achieve the desired release profile for multiple drugs is, however,
a complex multivariate problem. To address this issue, we have developed
a computational toolbox that models the migration of a drug pair through
a hydrated network of polysaccharide chains modified with hydrophobic
moieties. In this study, we chose doxorubicin (DOX) and Gemcitabine
(GEM) as model drugs, as their synergistic effect against breast cancer
has been thoroughly investigated, and chitosan as the model polymer.
Our model describes how the modification of chitosan chains with acetyl,
butanoyl, and heptanoyl moieties at different values χ governs
both the structure of the hydrogel network and drug migration through
it. Our experimental data confirm the in silico predictions for both
single- and dual-drug release and, most notably, the counterintuitive
inversion of release vs χ that occurs when switching from a
single- to a dual-drug system. Consensus between predicted and experimental
data indicates that acetyl modifications (χ = 32–42%)
and butanoyl modifications (χ = 19–24%) provide synergistic
GEM/DOX release molar ratios (i.e., 5–10). Collectively, these
results demonstrate the potential of this model in guiding the design
of chemotherapeutic hydrogels to combat cancer.