Translational Modeling Identifies Synergy between Nanoparticle-Delivered miRNA-22 and Standard-of-Care Drugs in Triple Negative Breast Cancer
=============================================================================================================================================

* Prashant Dogra
* Javier Ruiz Ramírez
* Joseph D. Butner
* Maria J. Peláez
* Caroline Chung
* Anupama Hooda-Nehra
* Renata Pasqualini
* Wadih Arap
* Vittorio Cristini
* George A. Calin
* Bulent Ozpolat
* Zhihui Wang

## Abstract

The downregulation of miRNA-22 in triple negative breast cancer (TNBC) is associated with upregulation of eukaryotic elongation 2 factor kinase (eEF2K) protein, which regulates tumor growth, chemoresistance, and tumor immunosurveillance. Moreover, exogenous administration of miRNA-22, loaded in nanoparticles to prevent degradation and improve tumor delivery (termed miRNA-22 nanotherapy), to suppress eEF2K production has shown potential as an investigational therapeutic agent *in vivo*. To evaluate the translational potential of miRNA-22 nanotherapy, we developed a multiscale mechanistic model, calibrated to published *in vivo* data and extrapolated to the human scale, to describe and quantify the pharmacokinetics and pharmacodynamics of miRNA-22 in virtual patient populations. Our analysis revealed the dose-response relationship, suggested optimal treatment frequency for miRNA-22 nanotherapy, and highlighted key determinants of therapy response, from which combination with immune checkpoint inhibitors was identified as a candidate strategy for improving treatment outcomes. More importantly, drug synergy was identified between miRNA-22 and standard-of-care drugs for TNBC, providing a basis for rational therapeutic combinations for improved response.

Keywords
*   Cancer treatment
*   mathematical modeling
*   microRNA
*   pharmacokinetics and pharmacodynamics
*   precision medicine
*   tumor-immune interaction
*   allometry

## Introduction

Triple negative breast cancer (TNBC) accounts for up to 10-12% of all breast cancer cases, and has a 5-year survival that is 8-16% lower than the hormone-receptor positive (HR+) disease subtype (1). Mechanisms to overcome the aggressiveness, histopathological heterogeneity, and prevalence of TNBC in younger women represent major unmet needs in contemporary cancer medicine (2). The severity of TNBC is further aggravated due to the lack of broadly-applicable targeted therapies, and by a high-rate of early metastases to the central nervous system and lungs (3). PARP inhibitors, such as olaparib and talazoparib, have been approved for TNBC with germline BRCA1 or BRCA2 gene mutations, but these are only reported in 15.4% of cases (4). Despite recent advances in developing therapeutics for treating TNBC such as immune checkpoint inhibitor immuotherapy (5), chemotherapy remains the most common recommended systemic regimens for TNBC, even though rapid development of chemoresistance is common (6).

An emerging body of evidence has supported key functional roles of microRNAs (miRNAs) in sustaining tumor proliferation, resisting growth inhibitors and cell death, inducing tumor invasion and metastasis, and promoting angiogenesis (7), suggesting that miRNAs may function as a valuable oncologic therapy target (8-10). Amongst the multitude of miRNAs, miRNA-22 (a chain of non-coding RNA consisting of 22 nucleotides) has been found to play a critical role in cancer initiation and progression processes (11, 12). Indeed, miRNA-22 has been extensively studied as a regulator of tumor suppressor genes like p53 (13) and as a repressor of the oncogene c-Myc (14) in many different cancers, including TNBC (12), hormone-dependent breast cancer (15), and colon cancer (16), and for its roles in metastasis suppression in breast and ovarian cancer (17), in the sensitization of esophageal carcinoma to γ-ray radiation (18). miRNA-22 has been shown to be downregulated in TNBC, which reduces its inhibitory control over the eukaryotic elongation 2 factor kinase (eEF2K), a tumor growth-promoting and chemoresistance-inducing protein (12, 19). Importantly, eEF2K was also found to enhance the expression of PD-L1, and is thus implicated for its role in blocking tumor immunosurveillance (20). Inhibiting these tumorigenic effects of eEF2K via exogenous administration of miRNA-22 represents a potential therapeutic approach to improve the response of TNBC to chemotherapy and/or immunotherapy with immune checkpoint inhibitors. However, naked miRNA has a short half-life due to its vulnerability to plasma ribonucleases, shows limited tumor penetration and cellular uptake due to its negative charge, and has off-target effects due to non-specific delivery (21). To overcome these shortcomings, nanoparticle (NP)-based drug delivery systems are being studied to improve the delivery of miRNAs to tumor cells (21). While the application of nanomaterials in cancer has been promising in improving tumor imaging and drug delivery (4, 22-27), challenges associated with low tumor deliverability due to off-target accumulation and limited tumor penetration continue to limit their clinical success (28). To this end, mechanistic mathematical modeling can be a valuable in-silico tool to help overcome this challenge, by furthering our understanding of NP-mediated miRNA-22 delivery in TNBC *in vivo*.

Mathematical modeling has been used to investigate the mechanisms relevant to tumor response to miRNA-based treatment. For instance, a system of ordinary differential equations (ODEs) with a delay term has been used to study feedback loops between the oncogenes Myc, EF2, and miRNA-17-92 (29). This model was subsequently expanded by integrating nine different mechanisms to evaluate how miRNAs regulate translation (30), and to study how the inactivation of a transcription factor is involved in cardiac dysfunction and cancer (31). In another notable study (32), an energy availability pathway involving miRNA-451 was analyzed in order to elucidate the difference between invasion and proliferation regimes in cancer cells, which was accomplished by combining a pair of ODEs governing the miRNA and glucose concentration with a system of partial differential equations (PDEs) employing transport mechanisms, such as diffusion, chemotaxis, and haptotaxis. Additionally, a signaling pathway relating miRNA-21, miRNA-155, and miRNA-205 to the proliferation and apoptosis of non-small-cell lung cancer cells has been examined with a series of modeling studies (33, 34). A noteworthy feature of the mathematical approach in (34) is the inclusion of a directional migration term, which takes into account the competition for available space between cells under the assumption of a logistic growth rate for cancer cells.

While previous modeling works focused on the molecular interactions of miRNAs associated with their therapeutic outcome, they lacked the inclusion of a viable drug delivery system and the related pharmacokinetics and drug delivery mechanisms required to assess the feasibility of miRNAs as a systemically-deliverable therapy for cancer. Therefore, to support the development of miRNA-22 as a viable therapy for TNBC, here we present a mechanistic mathematical model, formulated as a system of ODEs, to describe the tumoral delivery of systemically administered miRNA-22-loaded NPs and the pharmacodynamics of miRNA-22 in the context of TNBC growth. Our multiscale model incorporates processes pertinent to systemic NP pharmacokinetics, intratumoral transport of NPs, and the known molecular interactions of miRNA-22 with its associated oncogenes to predict TNBC growth dynamics. Following model calibration with *in vivo* data that was then allometrically scaled to humans, we simulated clinically relevant treatment of TNBC with miRNA-22 to obtain the dose-response relationship at the individual and population scales, thus helping to reveal the optimal dose and frequency of treatment for each individual “virtual” patient. Local and global sensitivity analyses of key model parameters revealed the importance of molecular interactions, tumor vascularization, miRNA-22 potency, NP characteristics, and immune checkpoint effects of anti-PD-L1 in governing the outcome of miRNA-22 therapy, thus highlighting some of the key determinants of treatment outome and suggesting the potential benefit of combination with immune checkpoint inhibitors. Drug synergy was identified to occur between miRNA-22 and standard-of-care therapies (including both chemotherapy and immunotherapy) studied in this work. As a result, our mechanistic model may serve as a useful computational means to help design and optimize a therapeutic framework for future clinical trials of miRNA-22.

## Methods

### Mathematical model development

We present a multiscale mechanistic model to simulate the *in vivo* and translational pharmacokinetics (PK) and pharmacodynamics (PD) of NP-mediated miRNA-22 therapy in TNBC, alone or in combination with chemotherapy or immunotherapy (here collectively referred to as *agents*), and thus investigate the factors governing the delivery of agents and their therapeutic efficacy. The model consists of two main compartments (**Figure 1**), represented by the *plasma* and *tumor*, where the latter is sub-compartmentalized into *vascular, interstitial, cellular membrane*, and *cytosolic space*. After injection into the plasma compartment, agents are cleared through various physiological processes, which (along with the volume of distribution of agents) governs their systemic (i.e., plasma) pharmacokinetics. From the plasma compartment, bi-directional, perfusion-mediated delivery characterizes the transport of agents into the tumor vasculature, from where the extravasation of agents across the permeable tumor vasculature introduces them into the tumor interstitium. Once in the interstitium, given the absence of advection due to high interstitial fluid pressure (35, 36), the agents undergo diffusion to reach the interstitium-cell membrane interface, from where based on the type of agent, specific biophysical processes occur to ensure delivery to the target site to invoke the pharmacodynamic effects of the agents. For instance, NPs undergo endocytosis into cancer cells where they release miRNA-22 in cancer cell cytosol, whereas free drugs simply diffuse into the cytosol and antibodies bind to their corresponding cell surface receptors (e.g. PD-L1 in the current context).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.19.21265154/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/F1)

Figure 1. Multiscale mechanistic model.
Model schematic shows key system interactions and variables. The plasma compartment is connected to the tumor compartment, with the latter sub-compartmentalized into vascular, interstitial, and cytosolic compartments. Key transport processes responsible for drug delivery to the tumor cytosol include perfusion, extravasation across tumor vasculature, diffusion across tumor interstitium, and endocytosis. While target receptors of immune checkpoint inhibitors are on the cell surface, the other agents including miRNA-22 and chemotherapeutics act intracellularly. Key signaling pathways relevant to miRNA-22 included in the model are shown in the cytosolic sub-compartment, including eEF2K induced tumor growth and PD-L1 production, miRNA-22 induced inhibition of eEF2K production, suppression of tumor antigenicity by checkpoint PD-L1, eEF2K induced chemoresistance, and induction of tumor death by chemotherapeutic agents.

Following delivery of agents to the target site, the pharmacodynamic component of the model is then engaged such that miRNA-22 acts by inhibiting the production of eEF2K protein, leading to inhibition of tumor growth. Also, because there is growing evidence in the literature that eEF2K induces the production of PD-L1 (20), this pathway was incorporated in the model to explore the engagement of immune checkpoints by eEF2K for tumor survival. The reference chemotherapeutic (i.e., doxorubicin) acts by inducing apoptotic cell death, whereas anti-PD-L1 antibodies act by inhibiting the tumor protective effects of PD-L1. Note that the mechanism of action of anti-PD-L1 antibodies modeled here includes the degradation of PD-L1 protein. Simulations and analysis of the model will help to provide insights into the systemic and tumoral pharmacokinetics of therapeutic agents and nanoparticles, along with the effects on tumor progression. We then use the model as an *in silico* tool to simulate virtual clinical trials in order to explore the effects of patient variability and other system parameters on treatment outcomes with mono- or combination therapies.

The various transport and pharmacological processes described above and shown in **Figure 1** have been formulated into a system of ordinary differential equations (ODEs; Eqs. 1–17) to obtain the temporal evolution of model behaviors of interest, including tumor growth. Equations pertaining to various compartments and biological processes are described below:

*Equation for NP mass kinetics in plasma* (*N*P(*t*)): ![Formula][1]</img>  where *V*B,V (= *f**v* ** B*(*t*)) and *V*P are volumes of tumor vasculature and plasma compartments, respectively; *f**v* is the vascular volume fraction of the tumor; *B*(*t*) represents tumor volume; *Q* represents plasma flow rate; *k*cl represents systemic clearance of NPs; *N* is the injected dose of NPs; and *i* represents the injection times (in weeks) post inoculation of tumor in mice (*i* = 2, 3, 4, 5). Note that *Q* (units, mL · mL−1 · wk−1) obeys the following empirical relationship with tumor volume (*B*(*t*)): *Q* = 2843 · *e*−0.65·*B*(*t*), obtained by fitting a monoexponential function to data from literature (37). Note that for human simulations, *Q* was assumed to be 1512 mL · mL−1 wk−1 irrespective of tumor size (38). Given that tissue density is ∼1 g · mL−1, tumor blood flow rate is provided in the units of mL · mL−1 · wk−1 in our work, which is numerically equivalent to mL · g−1 · wk−1 as typically used in the literature. In addition, NP clearance in mice *k*cl (units, wk−1) varies empirically with NP diameter (*ϕ*NP ; units, cm) as: ![Graphic][2]</img>, obtained by fitting the plasma half-life data of quantum dots of varying sizes from the literature (39, 40). For human simulations, the value of *k*Cl was allometrically scaled (see next Section).

*Equation for NP mass kinetics in tumor vasculature* (*N*V(*t*)): ![Formula][3]</img>  where *P*NP indicates NP permeability across tumor vasculature and *S* is the tumor vascular surface area per unit tumor volume (units, cm2/cm3). In vivo, *S* relates to tumor volume *B*(*t*) as *S* = 0.26 · *e* −4.5·*B*(*t*) + 138 · *e*−0.04·*B*(*t*), obtained empirically from the literature (41), and *P*NP (units, cm · wk−1) is a function of tumor vascular porosity and the ratio of NP size (*ϕ*NP) to tumor vascular pore size (*ϕ*pore; units, cm) (42). For human simulations, *S* was fixed at 135 cm2⁄cm3 irrespective of tumor size (41).

*Equation for NP mass kinetics in tumor interstitium* (*N*V(*t*)): ![Formula][4]</img>  where *D*NP is the diffusivity of NPs in tumor interstitium, and *L* is the characteristic interstitial distance between tumor vessels and cancerous cells, referred to herein as the intercapillary length in the tumor.

*Equation for NP mass kinetics in cancer cell membrane* (*N*M(*t*)): ![Formula][5]</img>  where *k*endo (units, wk−1) is the rate of endocytosis of NPs into tumor cells, and may be obtained by equating work done by the membrane motor proteins against surface tension of cell membrane (43).

*Equation for NP mass kinetics in cancer cell cytosol* (*N**C*(*t*)): ![Formula][6]</img>  where *δ*NP is NP degradation rate.

*Equation for miRNA concentration kinetics in cancer cell cytosol* (*C*M(*t*)): ![Formula][7]</img>  where ![Graphic][8]</img> is the intrinsic production rate of miRNA-22 in the tumor cytosol; *ε*B is the efficiency of tumor on inhibiting miRNA-22 production; *M* is the mass of miRNAs loaded in a single NP; and *V*B,C is the cytosolic volume of tumor (= *f**c* ** f*cy ** B*(*t*)), where *f**c* (= 0.4) is the cancer cell volume fraction of a tumor (44) and *f*cy (= 0.4) is the cytoplasmic volume fraction of a cancer cell (45). *k*rel is release rate of miRNAs from the endocytosed NP; and *δ*M is the degradation rate of miRNAs. Note that the initial condition *C*M(0) is estimated at the trivial steady state of the system when *B*(0) = 0 and no exogenous administration of miRNAs has occured, such that ![Graphic][9]</img>.

*Equation for anti-PD-L1 antibody concentration kinetics in plasma* (*C*Ab,P(*t*)): ![Formula][10]</img>  where DoseAb represents the dose of antibodies; *k*abs is the absorption rate constant of antibodies from the peritoneum into plasma after intraperitoneal (IP) injection; *V*pc is the volume of peritoneal fluid in female mice (0.1 mL) (46); ClAb is the systemic clearance of antibodies; and *j* represents the injection times (in weeks) post inoculation of tumor in mice defined in the set *Q**T* (= 1, 1.71, 2.43, 3.14, 3.86, 4.57). Note that in mouse experiments, antibodies were given IP; therefore, the initial plasma concentration of antibodies *C*0,Ab is zero. However, for human simulations, immunotherapy was administered as an intravenous (IV) bolus injection, therefore the initial plasma concentration of antibodies *C*0,Ab is non-zero and is calculated based on the initial dose (DoseAb) and systemic volume of distribution (*V*D,Ab) of the given antibody. As a result, the first term of Eq. 7 that accounts for systemic absorption of the drug from the peritoneum is removed during human simulations.

*Equation for anti-PD-L1 antibody concentration kinetics in tumor vasculature* (*C*Ab,v(*t*)): ![Formula][11]</img>  where *P*Ab indicates antibody permeability across tumor vasculature.

*Equation for anti-PD-L1 antibody concentration kinetics in tumor interstitium* (*C*Ab,I(*t*)): ![Formula][12]</img>  where *D*Ab is the diffusivity of antibodies in tumor interstitum.

*Equation for anti-PD-L1 antibody concentration kinetics in cancer cell membrane* (*C*Ab,M(*t*)): ![Formula][13]</img>  where, *δ*Ab is the degradation rate of antibodies.

*Equation for doxorubicin concentration in plasma* (*C*D,P(*t*)): ![Formula][14]</img>  where Cldox is the plasma clearance of doxorubicin; *C*0,D is the initial concentration of doxorubicin calculated based on the injected dose of 4mg/kg and the given volume of distribution (**Table 2**); *i* represents the injection times (in weeks) post inoculation of tumor in mice (*i* = 1, 2, 3).

*Equation for doxorubicin concentration in tumor vasculature* (*C*D,V(*t*)): ![Formula][15]</img>  where *J* is the diffusive flux of doxorubicin across the tumor vasculature into tumor interstitium, and is given by ![Graphic][16]</img>. Here,*D* *dox* is the diffusivity of doxorubicin in tumor interstitium and Δ*x* is the thickness of blood capillary wall (5 μm) (47).

*Equation for doxorubicin concentration in tumor interstitium* (*C*D,I(*t*)): ![Formula][17]</img>  *Equation for doxorubicin concentration in tumor cytosolic space* (*C*D,C(*t*)): ![Formula][18]</img>  where *δ*D is the degradation rate of doxorubicin in the cytosolic space.

*Equation for eEF2K concentration kinetics in cancer cell cytosol* (*C*E(*t*)): ![Formula][19]</img>  where ![Graphic][20]</img> is the basal production rate of eEF2K protein in tumor cytosol; *A*B,E is the stimulation factor of tumor effects on eEF2K production; *K*B,E is the Michaelis-Menten constant for tumor effects on eEF2K production; *δ*E is the degradation rate of eEF2K protein; *A*M,E is the stimulation factor of miRNA-22 effects on eEF2K protein degradation; ![Graphic][21]</img> is the half-maximal effective concentration of miRNA-22 for its effect on eEF2K protein degradation; and ![Graphic][22]</img> is the intial concentration of eEF2K protein.

*Equation for PD-L1 concentration kinetics in cancer cell membrane* (*C*P(*t*)): ![Formula][23]</img>  where ![Graphic][24]</img> is the basal production rate of PD-L1 protein in tumor cytosol; *A*E,P is the stimulation factor of eEF2K effects on PD-L1 production; *K*E,P is the Michaelis-Menten constant for eEF2K effects on PD-L1 production; *δ*P is the degradation rate of PD-L1 protein; *A*Ab,P is the stimulation factor of anti-PD-L1 antibody effects on PD-L1 degradation; ![Graphic][25]</img> is the half-maximal effective concentration of anti-PD-L1 antibody for its effect on PD-L1 protein degradation; and ![Graphic][26]</img> is the intial concentration of PD-L1 protein.

*Equation for tumor volume kinetics* (*B*(*t*)): ![Formula][27]</img>  where *σ* is the tumor growth rate constant; *A*E,B is the stimulation factor representing eEF2K effects on tumor growth; *K*E,B is the Michaelis-Menten constant for eEF2K effects on tumor growth; *B**** is the tumor carrying capacity; *δ*B,I is the death rate of tumors induced by normal immune system functionality (without drug intervention); *ε*P is the efficiency of PD-L1 protein in inhibit ing immune-induced tumor death; *δ*B,C is the death rate of tumors induced by doxorubicin; EC*D* is the half-maximal effective concentration of doxorubicin for its effect on tumor death; *A*E,D is the stimulation factor for eEF2K effects in inducing chemoresistance; *K*E,D is the Michaelis-Menten constant for eEF2K effects in inducing chemoresistance; and *B* is the size of inoculated tumor, i.e., initial condition (equal to a single cell volume for human simulations).

The model was solved numerically as an initial value problem in MATLAB R2018a using the built - in function *ode15s*, and fit to the *in vivo* data from literature (12, 48-50) using the built-in function *lsqcurvefit*. Correlation analysis was then performed between model fits and experimental data to assess the goodness of fit.

### Allometric scaling

For human simulations, the rate constants *k*Cl, *σ, δ*B,I, and *δ*B,C were allometrically scaled from values determined for mice, based on body weights and the standard allometric exponent for rate constants, i.e., –0.25 (51), such that the value of parameter *i* for humans ![Graphic][28]</img> was: ![Graphic][29]</img>, where ![Graphic][30]</img> is the value of parameter *i* for mice, and BW*h* and BW*m* are the body weights assumed for humans (70 kg) and mice (0.02 kg), respectively. We note that a different scaling exponent was used in the above formula for a subset of parameters, as based on published reports, these were: dose and clearance calculations (exponent = 0.75) and volume of distribution calculations (exponent = 1.0) (52, 53).

### Treatment response evaluation

The model was used to study the effect of miRNA-22 nanotherapy, alone or in combination with doxorubicin and/or Atezolizumab, in virtual patients. The effect of therapy on tumor shrinkage was quantified by a metric defined as percent tumor growth inhibition (%TGI), such that %TGI = (1 −*B*treated/*B*control) · 100, where *B*treated and *B*control represent treatment and control tumor volumes at the end of 104 weeks post tumor inception with a single cell. Note that, in the treatment scenario, therapy was initiated 80 weeks post tumor inception, such that treatment was given over 24 weeks (∼6 months). To evaluate tumor response at a population level, we employed a scale analogous to RECIST 1.1 (50), such that treatment response was classified as *progressive disease* (TGI ≤ 0%), *stable disease* (0% < TGI ≤ 10%), *intermediate response* (10% < TGI ≤ 30%), *partial response* (30% < TGI ≤ 50%), and *major response* (TGI > 50%).

### Parameter sensitivity analysis

To investigate the importance of various model parameters in causing tumor shrinkage in patients undergoing treatment with a weekly dose of 0.026 mg/kg miRNA-22 (loaded in NPs), we performed local (LSA) and global (GSA) sensitivity analyses (42, 54-57) by perturbing the parameters of interest (highlighted by a dagger in **Tables 1, 2**) over a range of 0.2× to 5× of their corresponding baseline values.

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/T1)

Table 1. List of biological parameters and initial conditions.

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/T2)

Table 2. List of therapy-related parameters and initial conditions.

LSA involved perturbation of one model parameter at a time at 500 levels between the range of 0.2× to 5× of the baseline value while the other parameters were held constant at baseline. Each parameter was perturbed individually and %TGI was calculated to obtain the qualitative relationship between parameter factor change and %TGI. Alternatively, in GSA, all model parameters of interest were simultaneously perturbed and %TGI calculated for each simulation (i.e., for a given combination of parameter values). Note that, to comprehensively investigate the vast multiparameter space (21 parameters), yet to minimize the number of simulations, Latin hypercube sampling (LHS) (42, 54, 55) was used to obtain 10,000 combinations of parameter values, and 10 such replicates were obtained. Multivariate linear regression analysis was then performed on every replicate, and regression coefficients were determined as a quantitative measure of parameter *sensitivity index* (SI). A distribution of regression coefficients (or SI) was obtained for each parameter, and one-way ANOVA with Tukey’s test was used to rank the parameters in terms of their sensitivity, such that a higher SI represents a greater influence on model output (i.e., %TGI).

### Determination of drug synergy

The Chou-Talalay method (58) was used to identify drug synergy between miRNA-22 and its combination with standard-of-care drugs for TNBC (doxorubicin and/or Atezolizumab). Occurence of drug synergy allows the possibility of using a lower dose of the constituent drugs, which can reduce their adverse effects. The method involves determination of *combination index* (CI), such that CI < 1 is an indicator of existence of drug synergism. To calculate CI, the open source software COMPUSYN (available at [https://www.combosyn.com/](https://www.combosyn.com/)) was used to generate the analysis report, which has been provided in the **Supplementary Information**.

## Results and Discussion

### Model development, calibration, and baseline solution

The multiscale mechanistic model developed to study the PK-PD of NP-mediated miRNA-22 therapy in TNBC, along with other clinically approved treatment modalities, was formulated as a system of ODEs (Eqs. 1-17) and solved numerically as an initial value problem. Some model parameters were known a priori (**Tables 1, 2**), while the rest were estimated through non-linear least squares fitting of the model to published *in vivo* datasets. The selected datasets include longitudinal measurements of tumor volume in mice bearing MDA-MB-231 xenografts under control conditions, or under treatment with one of the following therapies: 0.15 mg/kg (equivalent 4 μg/mouse) IV miRNA-22-loaded NPs once a week (12), 4 mg/kg IV doxorubicin once a week (48), and 5 mg/kg IP anti-PD-L1 immunotherapy (Atezolizumab) once every 5 days (50).

We used the model to simulate the treatment protocols shown in **Figure 2**, and the numerical solutions of tumor volume kinetics were then fit simultaneously to the four datasets to estimate the unknown model parameters (given in **Tables 1, 2**). Additionally, the model solution for eEF2K protein kinetics from the miRNA-22 simulation was fit to the available data (**Figure 2a**). Model fits were in good agreement with the experimental data, as indicated by a strong Pearson correlation (**Figure S1;** *R* > 0.96, *P* < 0.0001). While the various experimental studies used above demonstrated the effects of individual therapies on TNBC progression, the model revealed additional insights into drug (and also NP) pharmacokinetics and molecular interaction dynamics leading to tumor response to the three therapies.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.19.21265154/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/F2)

Figure 2. Model calibration.
Numerical solution of model fit to published *in vivo* data for treatment of MDA-MB-231 tumor-bearing mice with **a**) NP-delivered miRNA-22, **b**) doxorubicin, and **c**) Atezolizumab. Markers represent experimental data. Pearson correlation analysis results goodness of fit of the model are reported in Figure S1.

As shown in **Figure 2a**, simulated NP-mediated miRNA-22 therapy involved periodic IV administration of miRNA-22-loaded NPs into the plasma compartment, from where the NPs were cleared in a size-dependent fashion characterized by *k*cl, and also transported to the tumor vascular sub-compartment in a perfusion-dependent manner governed by the plasma flow rate *Q*. Extravasation of NPs across the leaky tumor vasculature into the tumor interstitium, determined by the vascular permeability-surface area product (*P*NP *· S*), was followed by NP size-dependent diffusion through the tumor interstitium. Remaining NPs were ultimitely endocytosed into the tumor cytosolic sub-compartment (i.e., cancer cell cytosol), followed by NP degradation and release of miRNA-22 into the cancer cell Interior. As a result, miRNA-22-induced inhibition of eEF2K production was observed relative to the control case, which reduced the downstream production of cancer cell transmembrane protein PD-L1. The overall effect of miRNA-22 therapy on alterations in protein expression manisfested as tumor growth inhibition mediated by suppressed induction of tumor growth by eEF2K and increased vulnerability to tumor immunogenicity due to depletion of the immune checkpoint PD-L1.

To explore the therapeutic combinations of miRNA-22 with FDA approved chemotherapies and immunotherapies, we calibrated the model with *in vivo* treatment of TNBC with doxorubicin (**Figure 2b**) and Atezolizumab (**Figure 2c**). In response to doxorubicin therapy (**Figure 2b**), the model showed reduction in tumor growth relative to the control case due to drug concentration-dependent increase in tumor death rate *δ*B,C. This is accompanied by reduced expression levels of eEF2K and PD-L1, but increased basal miRNA-22 expression level. Note that tumor growth has an inhibitory effect on miRNA-22 production (12), but stimulates eEF2K production, which tends to stimulate tumor growth in a complimentary feedback process (12, 19). We next modeled the effect of anti-PD-L1 immunotherapy (Atezolizumab) in a simplistic fashion by targeting PD-L1 degradation rate *δ*P, such that Atezolizumab enhances the degradation of PD-L1 in a drug concentration-dependent manner. As a result, as shown in **Figure 2c**, due to depletion of PD-L1, there is inhibition of tumor growth compared to the control. The parameters estimated as a result of the above model calibrations are given in **Tables 1, 2**.

### Model extrapolation to human scale

To study the translational value of miRNA-22 and associated potential combination therapies, we extrapolated the *in vivo* mechanistic model to human scale, either by substituting known physiological parameters with human values, or by allometric scaling of unknown parameters from mice to humans (see *Allometric Scaling* in **Methods**). In **Figure 3a**, a representative simulation of NP-mediated miRNA-22 therapy in a virtual adult patient (body weight 70 kg) is shown following once a week (QW) IV administration of 0.026 mg/kg miRNA-22 (allometrically scaled dose) for six months, starting 80 weeks after the inception of tumor with a single cell. As shown, eEF2K and PD-L1 levels are suppressed throughout the duration of treatment, thereby leading to ∼29% TGI compared to the control case. The prameters used for the representative simulation are the baseline values shown in **Tables 1, 2**.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.19.21265154/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/F3)

Figure 3. Translational PK-PD of miRNA-22.
**a**) Human extrapolation of *in vivo* mechanistic model simulating treatment with once weekly dose of miRNA-22 for six months. TGI indicates percent tumor growth inhibition. **b**) Dose response curves of a virtual patient under scenarios of once weekly (QW) or once every two weeks (Q2W) dose of miRNA-22 for slow and fast growing tumors. **c**) Effects of inter-individual variability on miRNA-22 therapy outcome for different QW doses, presented on a scale analogous to RECIST 1.1. Black arrow on x-axis indicates the dose of 0.026 mg/kg used for further analysis.

### Dose-response relationship and population variability

To investigate the effect of changes in miRNA-22 doseage and treatment frequency on %TGI, dose-response curves (DRCs) was generated for a representative individual by simulating treatment with miRNA-22 nanotherapy at different doses (0-0.2 mg/kg), either once weekly (QW) or once every two weeks (Q2W) (**Figure 3b**). DRCs were also generated for slow, medium, and fast growing tumors at the two treatment frequencies. Note that the ratio of tumor immunogenicity-induced death rate (e.g., tumor death rate due to normal immune effects without drug intervention) to tumor growth rate (*δ*B,I/*σ*) indiciates how fast a tumor grows; e.g., within the scope of this work, the ratios of 0.99, 0.9, and 0.75 indicate slow, medium, and fast growing tumors, respectively (**Figure S3**). As shown in **Figure 3b**, theraputic response tends to saturate around a dosage of 0.05 mg/kg in these six scenarios, and even beyond 0.026 mg/kg (dose obtained through allometric scaling; see *Allometric Scaling* in **Methods**), tumors do not exhibit significant increase in %TGI, hence 0.026 mg/kg was chosen as the reference dose in humans for further investigation. As for the slow growing tumors, they show a much higher response to therapy (∼two-fold) than their rapidly proliferating counterparts. Additionally, irrespective of the rate of tumor growth, the QW protocol causes greater %TGI than Q2W.

Further, by creating a virtual population of 2,000 patients through LHS of patient-specific parameters between ±50% of their baseline values, the effects of inter-individual variability on %TGI for different doses of QW miRNA-22 were investigated and presented in a manner *analogous* to the RECIST 1.1 classification (51, 59). As shown in **Figures 3c** and **S2**, the patient population showed significant improvement in response with increasing dose up to ∼0.02 mg/kg, such that stable disease (0% < TGI ≤ 10%) cases dropped exponentially, and the population of intermediate responders (10% < TGI ≤ 30%) and major responders (> 50% TGI) grew rapidly. Also, a steady increase was observed in the population of partial responders (30% < TGI ≤ 50%). However, beyond ∼0.02 mg/kg dosage, the population of major responders quickly saturated at a value of ∼20%, while the population of partial responders increased with increasing drug dosage up to ∼0.1 mg/kg, eventually settling at ∼25%; the remaining population (∼55%) primarily consisted of intermediate responsers. Thus, these observations support our use of the allometrically calculated dose of 0.026 mg/kg; this was used as the reference value for further analysis. Note that progressive disease (≤ 0% TGI) was only seen in the no treatment scenario; this indicates that treated tumors do not grow beyond the size of the corresponding control tumors, and hence as per our definition of %TGI, ≤ 0% values are not observed under treatment. These simulations provide quantification of the variation in treatment response that can be expected from physiological variability and tumor heterogeneity on a population scale, and can thus support treatment personalization to maximize patient benefit. Note that the parameters used for virtual population generation are marked by a dagger in **Tables 1**,**2**.

### Parameter sensitivity analysis

For a more complete understanding of the effects of both patient-specific and treatment-related parameters on %TGI following 0.026 mg/kg QW dose of miRNA-22 nanotherapy for six months, starting 80 weeks post initiation of tumor, we performed local (LSA) and global (GSA) sensitivity analyses by perturbing parameters over a range of 0.2× to 5× of the baseline values (42). As shown in **Figure 4a**, GSA ranked the 21 model parameters into eight categories based on their sensitivity indices (using one-way ANOVA and Tukey’s test), out of which we discuss the top five ranking parameter brackets below.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.19.21265154/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/F4)

Figure 4. Parameter sensitivity analysis.
**a**) Violin plot showing results of global sensitivity analysis such that parameters are plotted in a descending order of sensitivity from left to right. SI denotes sensitivity index. Parameters are bracketed based on their ranking obtained from Tukey’s test. **b**) Efects of individual parameters on %TGI is shown via local sensitivity analysis. Note that for both analyses, parameters were perturbed over a range of 0.2x to 5x of the baseline value. Red dot in each curve indicates the %TGI value corresponding to the baseline parameter values.

First, as shown in **Figure 4a**, miRNA-22 degradation rate (*δ*M) stands out for its influence on %TGI, indicating that the stability of the miRNA is critical to ensure therapeutic efficacy, and an increase in degradation rate of miRNA-22 causes reduction in %TGI (as revealed by LSA, **Figure 4b**), thereby reinforcing the need for NP-mediated delivery to protect the cargo until delivered to the cytosol. In the second bracker, tumor-specific parameters controlling PD-L1 and eEF2K protein production (![Graphic][31]</img> and ![Graphic][32]</img>, respectively), and the efficiency of PD-L1 (*ε*P) at inhibiting immune cell-induced tumor death rank second in GSA, suggesting that PD-L1-mediated tumor immunosurveillance blockade and eEF2K-induced tumor proliferation and PD-L1 production that affect the intrinsic tumor growth and death, are important determinants of miRNA-22 efficacy. This suggests that delivering anti-PD-L1 therapies to target high PD-L1 activity can improve treatment outcomes when used in combination combination with miRNA-22.

These parameters are followed by NP size (*ϕ*NP) and NP degradability (*δ*NP). NP characteristics strongly influence the systemic pharmacokinetics of NPs (driven by hepatic and renal clearance (42, 60-62)) and also NP transport to and accumulation within the tumor (driven by extravasation across tumor vasculature, diffusion through tumor interstitium, endocytosis into cancer cells, and metabolism-dependent degradation in the cancer cell cytosol (4, 22)). These NP-specific parameters rank highly for their influence on %TGI due to their role in miRNA-22 delivery to the tumor. Of note, as for the individual effects of NP size (*ϕ*NP), we observed an inverse monotonic trend between %TGI and the investigated parameter values, suggesting that an increase in NP size leads to reduced %TGI (**Figure 4b**). This suggests that while smaller NPs have smaller drug loading capacity, this may be compensated by using larger quantities to deliver the same dose of drug, which can then outperform larger NPs primarily due to better pharmacokinetics and greater tumor penetration. Note that the corresponding number of NPs injected to deliver 0.026 mg/kg miRNA-22 via NPs of different sizes in our study ranged from ∼90 billion (for size 350 nm) to ∼1.5 quadrillion (for size 14 nm), which lies well within the range of values used in preclinical studies and clinical trials (63). While we did not investigate renally clearable NPs (<10 nm) due to lack of reported clinical application for drug delivery, we anticipate poorer performance from such particles, primarily due to their short circulation half-life driven by rapid renal clearance (42, 64). Further, within the same ranking bracket is the parameter governing the induction of eEF2K degradation by miRNA-22 (*A*M,E), suggesting the expected significance of miRNA-22 for eEF2K degradation to inhibit tumor growth.

Parameters in the 4th and 5th ranking brackets include those that control the positive feedback between eEF2K protein and tumor growth (*A*E,B, *A*B,E), the half-maximal effective concentration of miRNA-22 in suppressing eEF2K ![Graphic][33]</img>, the potency of eEF2K in inducing tumor growth (*K*E,B), and importantly, the tumor microvascular surface area (*S*). This suggests that the positive feedback role of eEF2K with tumor growth is secondary relative to its role in immune suppression (these effects are found in the second bracket), and the primairy effective mechanism of action of miRNA-22 is immume suppression, and associate reduction of tumor growth represents a beneficial – but secondary – therapeutic mechanism. The importance of tumor microvascular surface area is attributed to its role in determining rate of extravasation of NPs across tumor vasculature for drug delivery to the cells. However, tumor blood flow rate (*Q*) does not appear to have a significant impact on therapy efficacy.

Together, our simulations find that tumor response is most sensitive to the immune checkpoint effects of PD-L1, miRNA-22-eEF2K interaction, miRNA-22 potency and stability, and NP characteristics. This finding may provide opportunities for patient-specific optimization of NP-mediated miRNA-22 therapy. It also warrants the use of combination therapies, particularly immune checkpoint inhibitors in combination with with chemotherapeutics due to the chemoresistive influence of eEF2K, in order to achieve a better treatment outcome. In light of the ranking obtained through GSA, we understand that LSA may not be required to obtain a ranked order of parameters for their influence on %TGI, because unlike GSA, LSA does not incorporate the interactions between parameters that may influence the outcome, and thus only provides a limited assessment into the sensitivity of parameters. However, LSA can still be used to obtain the empirical relationships between individual parameters and model output, as shown in **Figure 4b**.

### Combination therapies, population variability, and synergy

We then sought to test the effects of combining miRNA-22 with standard-of-care drugs for TNBC, i.e., chemotherapeutics (doxorubicin) and immune checkpoint inhibitors (Atezolizumab) for improvement in %TGI outcome. For these numerical experiments, clinically relevant doses of 2.4 mg/kg Q3W (once every three weeks) doxorubicin and 2 mg/kg Q3W Atezolizumab were simulated for the representative virtual patient (shown in **Figure 3a**) in various combinations with 0.026 mg/kg QW miRNA-22 given for six months, starting 80 weeks after tumor initation. As shown in **Figure 5a**, combining miRNA-22 nanotherapy with the immune checkpoint inhibitor improves the outcome from intermediate response to partial response (30% < TGI ≤ 50%), and combined with doxorubicin it leads to major response (TGI > 50%), which almost reaches complete response when the three modalities are given together. Further, to assess the effects of patient variability and tumor heterogeneity on drug combination outcomes, 2,000 virtual patients were generated as before, and as shown in **Figure 5b**, the three drug combination therapy produced major response in ∼60% of patients, which is three times the patients that showed major response with QW miRNA-22 monotherapy (**Figure 3c**). However, the response worsened when miRNA-22 was given Q2W, either alone or in combination (**Figure 5c**).

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.19.21265154/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/F5)

Figure 5. Combination therapies.
**a**) Effects of QW dose of miRNA-22 (M) alone or in combination with doxorubicin (Dox or D), or Atezolizumab (Atezo or A) on %TGI are show. **b**,**c**) Effects of inter-individual variability on %TGI following treatment with combination therapies is shown when miRNA-22 is given **b**) QW or **c**) Q2W. Note that the other three drugs were administered once every three weeks (Q3W) in all cases.

Finally, our observation that monotherapies without miRNA-22 show stable disease under the given treatment protocols (**Figure 5a**), while two or three drug combinations with miRNA-22 produce significant improvements in treatment outcome, warrants testing for the occurrence of drug synergy between miRNA-22 and other drugs. For this, as shown in **Figure 6a**, %TGI was calculated through model simulations for various doses of three monotherapies (miRNA-22 QW, doxorubicin Q3W, Atezolizumab Q3W) and three combination therapies (miRNA-22 QW + doxorubicin Q3W, miRNA-22 QW + Atezolizumab Q3W, miRNA-22 QW + doxorubicin Q3W + Atezolizumab Q3W), and was used as an input for Chou-Talalay method (58) to calculate the combination indices (CI) of drug combinations. As shown in **Figure 6b**, CI values <1 for the three combinations of miRNA-22 indicate drug synergy with doxorubicin, Atezolizumab, and doxorubicin+Atezolizumab.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.19.21265154/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2021/10/26/2021.10.19.21265154/F6)

Figure 6. Drug synergism.
**a**) Dose-response data generated from model simulations for various monotherapies and combination therapies. Note: Dox indicates doxorubicin and Ate represents Atezolizumab. **b**) Combination indices calculated with the Chou-Talalay method identify drug synergy for various combinations of miRNA-22. CI<1 indicates drug synergism.

## Conclusion

Using a multiscale mechanistic modeling approach, we studied the translational PK-PD of NP-loaded miRNA-22 as a therapeutic for TNBC, alone or in combination with other FDA approved therapeutics. For this, the model was first calibrated with published *in vivo* data involving treatment of MDA-MB-231 tumor-bearing mice with miRNA-22, doxorubicinin, or an immunecheckpoint inhibitor. The calibrated model was extrapolated to the human scale by substituting the physiological parameter values of mice with humans, or by scaling of the parameters with standard allometric techniques. Using the extrapolated model, the dose-response curves and effects of inter-individual variability on treatment outcome were assessed and quantified through a scale analogous to RECIST 1.1. Percent tumor growth inhibition (%TGI) saturated at a dose of 0.05 mg/kg, irrespective of the treatment frequency and doubling time of the tumor. For our translational analysis, a dose of 0.026 mg/kg was used, obtained through allometric scaling of dose for mice. By creating a virtual patient population through perturbation of patient-specific parameters, patient response to variable miRNA doses was quantified, and it was observed that at 0.02 mg/kg, the fraction of patient population showing major response (≥ 50% TGI) to therapy saturated at ∼40%. Within the scope of our computational investigation, further increment in the dose only increased the fraction of partial responders (i.e., patients exhibiting ≥ 30% and < 50% TGI), and appeared to saturate at 0.14 mg/kg with ∼35% patients exhibiting < 30% TGI above that dose.

Parameter sensitivity analysis was conducted to identify the key determinants of %TGI in miRNA-22 nanotherapy. This analysis revealed the significance of miRNA-22-eEF2K interaction, eEF2K-tumor growth feedback loop, tumor vascularization, miRNA-22 potency and stability, NP characteristics, and also the immunecheckpoint effects of PD-L1, which highlights the potential of combination with anti-PD-L1 therapy to improve %TGI. This was supported by numerical experiments involving the combination of NP-loaded miRNA-22 with a clinically used immune checkpoint inhibitor (Atezolizumab) and/or doxorubicin. The triple combination of miRNA-22 with doxorubicin and this immune checkpoint inhibitor lead to almost three-fold increase in population fraction exhibiting major response in comparison to miRNA-22 alone. Importantly, the suspected drug synergy between miRNA-22 and doxorubicin and immunecheckpoint inhibitors was confirmed through the Chou-Talalay combination, whereindices were found to be < 1.

Our analysis, based on a well-calibrated mathematical model extrapolated to the human scale, provides valuable pre-translational, quantitative insights into the limitations, challenges, and opportunities associated with the translation of miRNA-22 nanotherapy for TNBC patients. The ability of the model to explore the effects of patient variability and tumor heterogeneity through parameter perturbation and sampling demonstrates the utility of our in-silico tool to conduct *virtual* clinical trials to assess the effects of anticancer therapeutic agents, which can provide immediate feedback to biologists and clinicians regarding potential problems and their solutions to support the preclinical development and clinical translation of novel therapeutics. The model presented here captures the key processes involved in systemic pharmacokinetics of nanoparticles; however, for a more detailed characterization, we will integrate the tumor compartment with a whole-body physiologically-based pharmacokinetic model in future. Also, spatial tumor heterogeneity, genetic variability, and a more complete tumor microenvironment (with emphasis on immune cells) will be introduced in the tumor compartment of the model to further explore the effects of heterogeneity in drug diffusion barriers, drug resistant cell populations, and tumor immunosurveillance.

## Supporting information

SI [[supplements/265154_file02.pdf]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors.

## Author Contributions

ZW conceived the study and conceptualized the idea. PD, ZW designed the model and the model analysis. PD, JRR developed the model. PD performed the model analysis. PD, JRR, JDB, MJP, GAC, BO collected the data. PD, JRR, JDB, MJP, CC, AHN, RP, WA, VC, GAC, BO, ZW interpreted the results. PD, JRR, JDB, MJP, CC, AHN, RP, WA, VC, GAC, BO, ZW wrote or edited the manuscript.

## Competing Interests

GAC is the scientific founder of Ithax Pharmaceuticals.

## Acknowledgements

The research work was supported by the National Science Foundation grant DMS-1930583 (VC, ZW), National Institutes of Health grants 1R01CA253865 (VC, BO, ZW), 1U01CA196403 (VC, ZW), 1U01CA213759 (VC, BO, ZW), 1R01CA226537 (RP, WA, VC, ZW), 1R01CA222007 (VC, GAC, BO, ZW), U54CA210181 (VC), and the Cockrell Foundation (PD). GAC also acknowledgs the Felix L. Haas Endowed Professorship in Basic Science, and the following grant support: NCI (1R01 CA182905-01, 1R01CA222007-01A1), NIGMS (1R01GM122775-01), DoD (Idea Award W81XWH2110030), a Team DOD grant in Gastric Cancer, a Chronic Lymphocytic Leukemia Moonshot Flagship project, a CLL Global Research Foundation 2019 grant, a CLL Global Research Foundation 2020 grant, a Mathers Foundation grant, an Institutional Research and Development Grant associated with the Brain SPORE 2P50CA127001. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

*   Received October 19, 2021.
*   Revision received October 19, 2021.
*   Accepted October 26, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1. F.M. Howardand  O.I. Olopade. Epidemiology of Triple-Negative Breast Cancer: A Review. The Cancer Journal. 27:8–16 (2021).
    
    

2.  2. G. Bianchini,  J.M. Balko,  I.A. Mayer,  M.E. Sanders, and  L. Gianni. Triple-negative breast cancer: challenges and opportunities of a heterogeneous disease. Nature reviews Clinical oncology. 13:674–690 (2016).
    
    

3.  3. J. Collignon,  L. Lousberg,  H. Schroeder, and  G. Jerusalem. Triple-negative breast cancer: treatment challenges and solutions. Breast Cancer (Dove Med Press). 8:93–107 (2016).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27284266&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

4.  4. S. Goel,  C.A. Ferreira,  P. Dogra,  B. Yu,  C.J. Kutyreff,  C.M. Siamof,  J.W. Engle,  T.E. Barnhart,  V. Cristini, and  Z. Wang. Size-Optimized Ultrasmall Porous Silica Nanoparticles Depict Vasculature-Based Differential Targeting in Triple Negative Breast Cancer. Small (2019).
    
    

5.  5. A.L. Heekeand  A.R. Tan. Checkpoint inhibitor therapy for metastatic triple-negative breast cancer. Cancer metastasis reviews. 40:537–547 (2021).
    
    

6.  6. H. Goncalves, Jr..,  M.R. Guerra,  J.R. Duarte Cintra,  V.A. Fayer,  I.V. Brum, and  M.T. Bustamante Teixeira. Survival Study of Triple-Negative and Non-Triple-Negative Breast Cancer in a Brazilian Cohort. Clinical Medicine Insights Oncology. 12:1179554918790563 (2018).
    
    

7.  7. R. Garzon,  G.A. Calin, and  C.M. Croce. MicroRNAs in Cancer. Annual review of medicine. 60:167–179 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.med.59.053006.104707&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19630570&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268071100013&link_type=ISI) 

8.  8. A.G. Bader,  D. Brown, and  M. Winkler. The Promise of MicroRNA Replacement Therapy. Cancer Research. 70:7027–7030 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI3MC8xOC83MDI3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

9.  9. J.J. Rossi. New Hope for a MicroRNA Therapy for Liver Cancer. Cell. 137:990–992 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2009.05.038&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19524500&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266916400009&link_type=ISI) 

10. 10. C.J. Cheng,  R. Bahal,  I.A. Babar,  Z. Pincus,  F. Barrera,  C. Liu,  A. Svoronos,  D.T. Braddock,  P.M. Glazer,  D.M. Engelman,  W.M. Saltzman, and  F.J. Slack. MicroRNA silencing for cancer therapy targeted to the tumour microenvironment. Nature. 518:107–110 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature13905&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25409146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

11. 11. G.A. Calinand C.M. Croce. MicroRNA-Cancer Connection: The Beginning of a New Tale. Cancer Research. 66:7390–7394 (2006).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI2Ni8xNS83MzkwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

12. 12. A. Gorur,  R. Bayraktar,  C. Ivan,  H.A. Mokhlis,  E. Bayraktar,  N. Kahraman,  D. Karakas,  S. Karamil,  N.N. Kabil,  P. Kanlikilicer,  B. Aslan,  L. Tamer,  Z. Wang,  V. Cristini,  G. Lopez-Berestein,  G. Calin, and  B. Ozpolat. ncRNA therapy with miRNA-22-3p suppresses the growth of triple-negative breast cancer. Molecular therapy Nucleic acids. 23:930–943 (2021).
    
    

13. 13. N. Tsuchiya,  M. Izumiya,  H. Ogata-Kawata,  K. Okamoto,  Y. Fujiwara,  M. Nakai,  A. Okabe,  A.J. Schetter,  E.D. Bowman,  Y. Midorikawa,  Y. Sugiyama,  H. Aburatani,  C.C. Harris, and  H. Nakagama. Tumor suppressor miR-22 determines p53-dependent cellular fate through post-transcriptional regulation of p21. Cancer Res. 71:4628–4639 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1158/1538-7445.AM2011-4628&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21565979&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

14. 14. J. Xiong,  Q. Du, and  Z. Liang. Tumor-suppressive microRNA-22 inhibits the transcription of E-box-containing c-Myc target genes by silencing c-Myc binding protein. Oncogene. 29:4980–4988 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/onc.2010.241&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20562918&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281578700012&link_type=ISI) 

15. 15. J. Xiong,  D. Yu,  N. Wei,  H. Fu,  T. Cai,  Y. Huang,  C. Wu,  X. Zheng,  Q. Du,  D. Lin, and  Z. Liang. An estrogen receptor α suppressor, microRNA-22, is downregulated in estrogen receptor α-positive human breast cancer cell lines and clinical samples. The FEBS Journal. 277:1684–1694 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1742-4658.2010.07594.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20180843&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

16. 16. S. Alvarez-Diaz,  N. Valle,  G. Ferrer-Mayorga,  L. Lombardia,  M. Herrera,  O. Dominguez,  M.F. Segura,  F. Bonilla,  E. Hernando, and  A. Munoz. MicroRNA-22 is induced by vitamin D and contributes to its antiproliferative, antimigratory and gene regulatory effects in colon cancer cells. Hum Mol Genet. 21:2157–2165 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/dds031&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22328083&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

17. 17. J.B. Patel,  H.N. Appaiah,  R.M. Burnett,  P. Bhat-Nakshatri,  G. Wang,  R. Mehta,  S. Badve,  M.J. Thomson,  S. Hammond,  P. Steeg,  Y. Liu, and  H. Nakshatri. Control of EVI-1 oncogene expression in metastatic breast cancer cells through microRNA miR-22. Oncogene. 30:1290–1301 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/onc.2010.510&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21057539&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000288492100004&link_type=ISI) 

18. 18. X.-c. Wang,  Z.-B. Zhang,  Y.-Y. Wang,  H.-Y. Wu,  D.-G. Li,  A.-M. Meng, and  F.-Y. Fan. Increased miRNA-22 expression sensitizes esophageal squamous cell carcinoma to irradiation. Journal of Radiation Research. 54:401–408 (2012).
    
    

19. 19. X. Wang,  J. Xie, and  C.G. Proud. Eukaryotic elongation factor 2 kinase (eEF2K) in cancer. Cancers. 9:162 (2017).
    
    

20. 20. Y. Wu,  J. Xie,  X. Jin,  R.V. Lenchine,  X. Wang,  D.M. Fang,  Z.D. Nassar,  L.M. Butler,  J. Li, and  C.G. Proud. eEF2K enhances expression of PD-L1 by promoting the translation of its mRNA. Biochem J. 477:4367–4381 (2020).
    
    

21. 21. S.W.L. Lee,  C. Paoletti,  M. Campisi,  T. Osaki,  G. Adriani,  R.D. Kamm,  C. Mattu, and  V. Chiono. MicroRNA delivery through nanoparticles. Journal of Controlled Release. 313:80–95 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jconrel.2019.10.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31622695&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

22. 22. S. Goel,  G. Zhang,  P. Dogra,  S. Nizzero,  V. Cristini,  Z. Wang,  Z. Hu,  Z. Li,  X. Liu, and  H. Shen. Sequential deconstruction of composite drug transport in metastatic breast cancer. Science Advances. 6:eaba4498 (2020).
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo4OiJhZHZhbmNlcyI7czo1OiJyZXNpZCI7czoxMzoiNi8yNi9lYWJhNDQ5OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzI2LzIwMjEuMTAuMTkuMjEyNjUxNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

23. 23. G. Brachi,  J. Ruiz-Ramírez,  P. Dogra,  Z. Wang,  V. Cristini,  G. Ciardelli,  R.C. Rostomily,  M. Ferrari,  A.M. Mikheev, and  E. Blanco. Intratumoral injection of hydrogel-embedded nanoparticles enhances retention in glioblastoma. Nanoscale (2020).
    
    

24. 24. H. Hosoya,  A.S. Dobroff,  W.H. Driessen,  V. Cristini,  L.M. Brinker,  F.I. Staquicini,  M. Cardó-Vila,  S. D’Angelo,  F. Ferrara, and  B. Proneth. Integrated nanotechnology platform for tumor-targeted multimodal imaging and therapeutic cargo release. Proceedings of the National Academy of Sciences. 113:1877–1882 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTEzLzcvMTg3NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzI2LzIwMjEuMTAuMTkuMjEyNjUxNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

25. 25. T.A. Brocato,  E.N. Coker,  P.N. Durfee,  Y.S. Lin,  J. Townson,  E.F. Wyckoff,  V. Cristini,  C.J. Brinker, and  Z. Wang. Understanding the Connection between Nanoparticle Uptake and Cancer Treatment Efficacy using Mathematical Modeling. Sci Rep. 8:7538 (2018).
    
    

26. 26. Z. Wang,  R. Kerketta,  Y.L. Chuang,  P. Dogra,  J.D. Butner,  T.A. Brocato,  A. Day,  R. Xu,  H. Shen,  E. Simbawa,  A.S. Al-Fhaid,  S.R. Mahmoud,  S.A. Curley,  M. Ferrari,  E.J. Koay, and  V. Cristini. Theory and Experimental Validation of a Spatio-temporal Model of Chemotherapy Transport to Enhance Tumor Cell Kill. PLoS Comput Biol. 12:e1004969 (2016).
    
    

27. 27.1.   M. Saravananand H. Barabadi
    
     A. Noureddine,  J.D. Butner,  W. Zhu,  P. Naydenkov,  M.J. Peláez,  S. Goel,  Z. Wang,  C.J. Brinker,  V. Cristini, and  P. Dogra. Emerging Lipid-Coated Silica Nanoparticles for Cancer Therapy. In  M. Saravananand H. Barabadi (eds.), Cancer Nanotheranostics: Volume 1, Springer International Publishing, Cham, 2021, pp. 335–361.
    
    

28. 28. S. Wilhelm,  A.J. Tavares,  Q. Dai,  S. Ohta,  J. Audet,  H.F. Dvorak, and  W.C. Chan. Analysis of nanoparticle delivery to tumours. Nature reviews materials. 1:1–12 (2016).
    
    

29. 29. B.D. Aguda,  Y. Kim,  M.G. Piper-Hunter,  A. Friedman, and  C.B. Marsh. MicroRNA regulation of a cancer network: Consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proceedings of the National Academy of Sciences. 105:19678–19683 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA1LzUwLzE5Njc4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

30. 30. N. Morozova,  A. Zinovyev,  N. Nonne,  L.-L. Pritchard,  A.N. Gorban, and  A. Harel-Bellan. Kinetic signatures of microRNA modes of action. RNA. 18:1635–1655 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoicm5hIjtzOjU6InJlc2lkIjtzOjk6IjE4LzkvMTYzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzI2LzIwMjEuMTAuMTkuMjEyNjUxNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

31. 31. M.J. Tindalland  A. Clerk. Modelling Negative Feedback Networks for Activating Transcription Factor 3 Predicts a Dominant Role for miRNAs in Immediate Early Gene Regulation. PLOS Computational Biology. 10:e1003597 (2014).
    
    

32. 32. Y. Kim,  S. Roh,  S. Lawler, and  A. Friedman. miR451 and AMPK mutual antagonism in glioma cell migration and proliferation: a mathematical model. PLoS One. 6:e28293–e28293 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0028293&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22205943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

33. 33. H.-W. Kang,  M. Crawford,  M. Fabbri,  G. Nuovo,  M. Garofalo,  S.P. Nana-Sinkam, and  A. Friedman. A mathematical model for microRNA in lung cancer. PLoS One. 8:e53663–e53663 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0053663&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23365639&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

34. 34. X. Laiand  A. Friedman. Exosomal miRs in Lung Cancer: A Mathematical Model. PLoS One. 11:e0167706 (2016).
    
    

35. 35. T. Brocato,  P. Dogra,  E.J. Koay,  A. Day,  Y.-L. Chuang,  Z. Wang, and  V. Cristini. Understanding Drug Resistance in Breast Cancer with Mathematical Oncology. Current Breast Cancer Reports:1-11 (2014).
    
    

36. 36. P. Dogra,  J.R. Ramirez,  M.J. Pelaez,  Z. Wang,  V. Cristini,  G. Parasher, and  M. Rawat. Mathematical Modeling to Address Challenges in Pancreatic Cancer. Curr Top Med Chem (2019).
    
    

37. 37. S.-G. Kimand  J.J. Ackerman. Quantitative determination of tumor blood flow and perfusion via deuterium nuclear magnetic resonance spectroscopy in mice. Cancer Research. 48:3449–3453 (1988).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI0OC8xMi8zNDQ5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

38. 38. F.M. Waterman,  L. Tupchong,  R.E. Nerlinger, and  J. Matthews. Blood flow in human tumors during local hyperthermia. International Journal of Radiation Oncology* Biology* Physics. 20:1255–1262 (1991).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2045300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991FQ58100013&link_type=ISI) 

39. 39. H.S. Choi,  W. Liu,  P. Misra,  E. Tanaka,  J.P. Zimmer,  B.I. Ipe,  M.G. Bawendi, and  J.V. Frangioni. Renal clearance of nanoparticles. Nature biotechnology. 25:1165 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt1340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17891134&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250226600028&link_type=ISI) 

40. 40. N. Hoshyar,  S. Gray,  H. Han, and  G. Bao. The effect of nanoparticle size on in vivo pharmacokinetics and cellular interaction. Nanomedicine. 11:673–692 (2016).
    
    

41. 41. R.K. Jain. Determinants of tumor blood flow: a review. Cancer Research. 48:2641–2658 (1988).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI0OC8xMC8yNjQxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

42. 42. P. Dogra,  J.D. Butner,  J.R. Ramírez,  Y.-l. Chuang,  A. Noureddine,  C.J. Brinker,  V. Cristini, and  Z. Wang. A mathematical model to predict nanomedicine pharmacokinetics and tumor delivery. Computational and Structural Biotechnology Journal. 18:518–531 (2020).
    
    

43. 43. O. Lunov,  V. Zablotskii,  T. Syrovets,  C. Röcker,  K. Tron,  G.U. Nienhaus, and  T. Simmet. Modeling receptor-mediated endocytosis of polymer-functionalized iron oxide nanoparticles by human macrophages. Biomaterials. 32:547–555 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biomaterials.2010.08.111&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20880574&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285401500023&link_type=ISI) 

44. 44. M. Mikubo,  K. Seto,  A. Kitamura,  M. Nakaguro,  Y. Hattori,  N. Maeda,  T. Miyazaki,  K. Watanabe,  H. Murakami, and  T. Tsukamoto. Calculating the tumor nuclei content for comprehensive cancer panel testing. Journal of Thoracic Oncology. 15:130–137 (2020).
    
    

45. 45. W.-W. Sung,  Y.-M. Lin,  P.-R. Wu,  H.-H. Yen,  H.-W. Lai,  T.-C. Su,  R.-H. Huang,  C.-K. Wen,  C.-Y. Chen, and  C.-J. Chen. High nuclear/cytoplasmic ratio of Cdk1 expression predicts poor prognosis in colorectal cancer patients. BMC cancer. 14:1–7 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2407-14-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24383403&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

46. 46. F. Hartveitand  S. Thunold. Peritoneal fluid volume and the oestrus cycle in mice. Nature. 210:1123–1125 (1966).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/2101123a0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6007179&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

47. 47. T. Stylianopoulos,  K. Soteriou,  D. Fukumura, and  R.K. Jain. Cationic nanoparticles have superior transvascular flux into solid tumors: insights from a mathematical model. Annals of biomedical engineering. 41:68–77 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10439-012-0630-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22855118&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

48. 48. A. Bandyopadhyay,  L. Wang,  J. Agyin,  Y. Tang,  S. Lin,  I.-T. Yeh,  K. De, and  L.-Z. Sun. Doxorubicin in combination with a small TGFβ inhibitor: a potential novel therapy for metastatic breast cancer in mouse models. PLoS One. 5:e10365 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0010365&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20442777&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

49. 49. J.D. Butner,  D. Fuentes,  B. Ozpolat,  G.A. Calin,  X. Zhou,  J. Lowengrub,  V. Cristini, and  Z. Wang. A Multiscale Agent-Based Model of Ductal Carcinoma In Situ. IEEE transactions on bio-medical engineering. 67:1450–1461 (2020).
    
    

50. 50. J. Rios-Doria,  C. Stevens,  C. Maddage,  K. Lasky, and  H.K. Koblish. Characterization of human cancer xenografts in humanized mice. Journal for immunotherapy of cancer. 8: (2020).
    
    

51. 51. A. Lindauer,  C. Valiathan,  K. Mehta,  V. Sriram,  R. De Greef,  J. Elassaiss-Schaap, and  D. De Alwis. Translational pharmacokinetic/pharmacodynamic modeling of tumor growth inhibition supports dose-range selection of the anti–PD-1 antibody pembrolizumab. CPT: pharmacometrics & systems pharmacology. 6:11–20 (2017).
    
    

52. 52. G.B. Westand J.H. Brown. The origin of allometric scaling laws in biology from genomes to ecosystems: towards a quantitative unifying theory of biological structure and organization. Journal of experimental biology. 208:1575–1592 (2005).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiamV4YmlvIjtzOjU6InJlc2lkIjtzOjEwOiIyMDgvOS8xNTc1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

53. 53. J. Mordenti,  S.A. Chen,  J.A. Moore,  B.L. Ferraiolo, and  J.D. Green. Interspecies scaling of clearance and volume of distribution data for five therapeutic proteins. Pharmaceutical research. 8:1351–1359 (1991).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1015836720294&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1798669&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GM79800002&link_type=ISI) 

54. 54. Z. Wang,  T.S. Deisboeck, and  V. Cristini. Development of a sampling-based global sensitivity analysis workflow for multiscale computational cancer models. IET systems biology. 8:191–197 (2014).
    
    

55. 55. Z. Wang,  V. Bordas, and  T. Deisboeck. Identification of critical molecular components in a multiscale cancer model based on the integration of Monte Carlo, resampling, and ANOVA. Frontiers in physiology. 2:35 (2011).
    
    

56. 56. Z. Wang,  V. Bordas,  J. Sagotsky, and  T.S. Deisboeck. Identifying therapeutic targets in a combined EGFR-TGFbetaR signalling cascade using a multiscale agent-based cancer model. Mathematical medicine and biology : a journal of the IMA. 29:95–108 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/imammb/dqq023&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21147846&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 

57. 57. P. Dogra,  J. Ruiz-Ramirez,  K. Sinha,  J.D. Butner,  M.J. Pelaez,  M. Rawat,  V.K. Yellepeddi,  R. Pasqualini,  W. Arap,  H.D. Sostman,  V. Cristini, and  Z. Wang. Innate Immunity Plays a Key Role in Controlling Viral Load in COVID-19: Mechanistic Insights from a Whole-Body Infection Dynamics Model. ACS pharmacology & translational science. 4:248–265 (2021).
    
    

58. 58. T.-C. Chouand P. Talalay. Quantitative analysis of dose-effect relationships: the combined effects of multiple drugs or enzyme inhibitors. Advances in enzyme regulation. 22:27–55 (1984).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0065-2571(84)90007-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6382953&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1984SU37500003&link_type=ISI) 

59. 59. E.A. Eisenhauer,  P. Therasse,  J. Bogaerts,  L.H. Schwartz,  D. Sargent,  R. Ford,  J. Dancey,  S. Arbuck,  S. Gwyther, and  M. Mooney. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). European journal of cancer. 45:228–247 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejca.2008.10.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19097774&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262948300002&link_type=ISI) 

60. 60. P. Dogra,  N.L. Adolphi,  Z. Wang,  Y.-S. Lin,  K.S. Butler,  P.N. Durfee,  J.G. Croissant,  A. Noureddine,  E.N. Coker,  E.L. Bearer,  V. Cristini, and  C.J. Brinker. Establishing the effects of mesoporous silica nanoparticle properties on in vivo disposition using imaging-based pharmacokinetics. Nature Communications. 9:4551 (2018).
    
    

61. 61. P. Dogra,  J.D. Butner,  Y.-l. Chuang,  S. Caserta,  S. Goel,  C.J. Brinker,  V. Cristini, and  Z. Wang. Mathematical modeling in cancer nanomedicine: a review. Biomedical Microdevices. 21:40 (2019).
    
    

62. 62. P. Dogra,  J.D. Butner,  S. Nizzero,  J. Ruiz Ramirez,  A. Noureddine,  M.J. Pelaez,  D. Elganainy,  Z. Yang,  A.D. Le,  S. Goel,  H.S. Leong,  E.J. Koay,  C.J. Brinker,  V. Cristini, and  Z. Wang. Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies. Wiley interdisciplinary reviews Nanomedicine and nanobiotechnology. 12:e1628 (2020).
    
    

63. 63. B. Ouyang,  W. Poon,  Y.-N. Zhang,  Z.P. Lin,  B.R. Kingston,  A.J. Tavares,  Y. Zhang,  J. Chen,  M.S. Valic, and  A.M. Syed. The dose threshold for nanoparticle tumour delivery. Nature materials. 19:1362–1371 (2020).
    
    

64. 64. H.S. Choi,  W. Liu,  P. Misra,  E. Tanaka,  J.P. Zimmer,  B.I. Ipe,  M.G. Bawendi, and  J.V.J.N.b. Frangioni. Renal clearance of quantum dots. 25:1165–1170 (2007).
    
    

65. 65. H. Hashizume,  P. Baluk,  S. Morikawa,  J.W. McLean,  G. Thurston,  S. Roberge,  R.K. Jain, and  D.M. McDonald. Openings between defective endothelial cells explain tumor vessel leakiness. The American journal of pathology. 156:1363–1380 (2000).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0002-9440(10)65006-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10751361&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086413900028&link_type=ISI) 

66. 66. M. Konerding,  E. Fait, and  A. Gaumann. 3D microvascular architecture of pre-cancerous lesions and invasive carcinomas of the colon. British journal of cancer. 84:1354–1362 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1054/bjoc.2001.1809&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11355947&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169042800013&link_type=ISI) 

67. 67. C.-T. Chen,  D.S. Malkus, and  R. Vanderby Jr.. A fiber matrix model for interstitial fluid flow and permeability in ligaments and tendons. Biorheology. 35:103–118 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-355X(99)80001-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10193483&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000077054100001&link_type=ISI) 

68. 68. S. Seré,  B. De Roo,  M. Vervaele,  S. Van Gool,  S. Jacobs,  J.W. Seo, and  J. -P. Locquet. Altering the biodegradation of mesoporous silica nanoparticles by means of experimental parameters and surface functionalization. Journal of Nanomaterials. 2018: (2018).
    
    

69. 69. M.J. Marzi,  F. Ghini,  B. Cerruti,  S. De Pretis,  P. Bonetti,  C. Giacomelli,  M.M. Gorski,  T. Kress,  M. Pelizzola, and  H. Muller. Degradation dynamics of microRNAs revealed by a novel pulse-chase approach. Genome research. 26:554–565 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjI2LzQvNTU0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMjYvMjAyMS4xMC4xOS4yMTI2NTE1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

70. 70. D.L. Gustafson,  J.C. Rastatter,  T. Colombo, and  M.E. Long. Doxorubicin pharmacokinetics: Macromolecule binding, metabolism, and excretion in the context of a physiologic model. Journal of pharmaceutical sciences. 91:1488–1501 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jps.10161&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12115848&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.19.21265154.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000175960000012&link_type=ISI) 

71. 71. Y. Fan,  N.-m. Lin,  L.-h. Luo,  L. Fang,  Z.-y. Huang,  H.-f. Yu, and  F.-q. Wu. Pharmacodynamic and pharmacokinetic study of pegylated liposomal doxorubicin combination (CCOP) chemotherapy in patients with peripheral T-cell lymphomas. Acta Pharmacologica Sinica. 32:408–414 (2011).
    
    

72. 72. M. Reth. Matching cellular dimensions with molecular sizes. Nature immunology. 14:765–767 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ni.2621&link_type=DOI) 

73. 73.Center for Drug Evaluation and Researchand FDA. Clinical Pharmacology Review - Atezolizumab. [https://www.accessdata.fda.gov/drugsatfda\_docs/nda/2016/761034orig1s000clinpharmr.pdf](https://www.accessdata.fda.gov/drugsatfda_docs/nda/2016/761034orig1s000clinpharmr.pdf) (accessed 10/6/2021.
    
    

74. 74. A.D.S. Linhares,  C. Battin,  S. Jutz,  J. Leitner,  C. Hafner,  J. Tobias,  U. Wiedermann,  M. Kundi,  G.J. Zlabinger, and  K. Grabmeier-Pfistershammer. Therapeutic PD-L1 antibodies are more effective than PD-1 antibodies in blocking PD-1/PD-L1 signaling. Scientific reports. 9:1–9 (2019).

 [1]: /embed/graphic-2.gif
 [2]: /embed/inline-graphic-1.gif
 [3]: /embed/graphic-3.gif
 [4]: /embed/graphic-4.gif
 [5]: /embed/graphic-5.gif
 [6]: /embed/graphic-6.gif
 [7]: /embed/graphic-7.gif
 [8]: /embed/inline-graphic-2.gif
 [9]: /embed/inline-graphic-3.gif
 [10]: /embed/graphic-8.gif
 [11]: /embed/graphic-9.gif
 [12]: /embed/graphic-10.gif
 [13]: /embed/graphic-11.gif
 [14]: /embed/graphic-12.gif
 [15]: /embed/graphic-13.gif
 [16]: /embed/inline-graphic-4.gif
 [17]: /embed/graphic-14.gif
 [18]: /embed/graphic-15.gif
 [19]: /embed/graphic-16.gif
 [20]: /embed/inline-graphic-5.gif
 [21]: /embed/inline-graphic-6.gif
 [22]: /embed/inline-graphic-7.gif
 [23]: /embed/graphic-17.gif
 [24]: /embed/inline-graphic-8.gif
 [25]: /embed/inline-graphic-9.gif
 [26]: /embed/inline-graphic-10.gif
 [27]: /embed/graphic-18.gif
 [28]: /embed/inline-graphic-11.gif
 [29]: /embed/inline-graphic-12.gif
 [30]: /embed/inline-graphic-13.gif
 [31]: /embed/inline-graphic-14.gif
 [32]: /embed/inline-graphic-15.gif
 [33]: /embed/inline-graphic-16.gif