Mitigating temozolomide resistance in glioblastoma via DNA damage-repair inhibition.

Glioblastomas are among the most lethal cancers, with a 5 year survival rate below 25%. Temozolomide is typically used in glioblastoma treatment; however, the enzymes alkylpurine-DNA-N-glycosylase (APNG) and methylguanine-DNA-methyltransferase (MGMT) efficiently mediate the repair of DNA damage caused by temozolomide, reducing treatment efficacy. Consequently, APNG and MGMT inhibition has been proposed as a way of overcoming chemotherapy resistance. Here, we develop a mechanistic mathematical model that explicitly incorporates the effects of chemotherapy on tumour cells, including the processes of DNA damage induction, cell arrest and DNA repair. Our model is carefully parametrized and validated, and then used to virtually recreate the response of heteroclonal glioblastomas to dual treatment with temozolomide and inhibitors of APNG/MGMT. Using our mechanistic model, we identify four combination treatment strategies optimized by tumour cell phenotype, and isolate the strategy most likely to succeed in a pre-clinical and clinical setting. If confirmed in clinical trials, these strategies have the potential to offset chemotherapy resistance in patients with glioblastoma and improve overall survival.

A Biomechanical Model of Tumor-Induced Intracranial Pressure and Edema in Brain Tissue.

Brain tumor growth and tumor-induced edema result in increased intracranial pressure (ICP), which, in turn, is responsible for conditions as benign as headaches and vomiting or as severe as seizures, neurological damage, or even death. Therefore, it has been hypothesized that tracking ICP dynamics may offer improved prognostic potential in terms of early detection of brain cancer and better delimitation of the tumor boundary. However, translating such theory into clinical practice remains a challenge, in part because of an incomplete understanding of how ICP correlates with tumor grade. Here, we propose a multiphase mixture model that describes the biomechanical response of healthy brain tissue-in terms of changes in ICP and edema-to a growing tumor. The model captures ICP dynamics within the diseased brain and accounts for the ability/inability of healthy tissue to compensate for this pressure. We propose parameter regimes that distinguish brain tumors by grade, thereby providing critical insight into how ICP dynamics vary by severity of disease. In particular, we offer an explanation for clinically observed phenomena, such as a lack of symptoms in low-grade glioma patients versus a rapid onset of symptoms in those with malignant tumors. Our model also takes into account the effects tumor-derived proteases may have on ICP levels and the extent of tumor invasion. This work represents an important first step toward understanding the mechanisms that underlie the onset of edema and ICP in cancer-afflicted brains. Continued modeling effort in this direction has the potential to make an impact in the field of brain cancer diagnostics.

Synthesis, Characterization, and Function of an RNA-Based Transfection Reagent.

A synthetic 8-mer, amphipathic, trans-acting poly-2'-O-methyluridylic thiophosphate triester RNA element (2'-OMeUtaPS) can be prepared using solid-phase synthesis protocols. 2'-OMeUtaPS efficiently mediates the delivery of uncharged polyA-tailed phosphorodiamidate morpholino (PMO) sequences in HeLa pLuc 705 cells, as evidenced by flow cytometry measurements. In this cell line, 2'-OMeUtaPS-mediated transfection of an antisense polyA-tailed PMO sequence induces alternative splicing of an aberrant luciferase pre-mRNA splice site, leading to restoration of functional luciferase, as quantitatively measured using a typical luciferase assay. 2'-OMeUtaPS is also potent at delivering an uncharged antisense polyA-tailed PMO sequence in muscle cells of the mdx mouse model of muscular dystrophy; targeting the polyA-tailed PMO sequence against a splice site of the pre-mRNA encoding mutated dystrophin triggers an alternate splicing event that results in excision of the mutated exon (exon 23) from the pre-mRNA and production of functional dystrophin, as demonstrated by agarose gel electrophoresis. © 2018 by John Wiley & Sons, Inc.

Trans-acting oligodeoxythymidine phosphorothioate triester reagents for transient transfection optimized and facilitated by a high-throughput microbioreactor system.

A rapid and cost-effective transient transfection method for mammalian cells is essential for screening biopharmaceuticals in early stages of development. A library of 25 amphipathic trans-acting oligodeoxythymidine phosphorothioate triester (dTtaPS) transfection reagents, carrying positively charged and lipophilic groups, has been constructed for this purpose. High-throughput screening of the library, using an imaging cytometer and an automated microbioreactor system, has led to the identification of dTtaPS10+ as a potent transfection reagent. This reagent efficiently delivers a plasmid encoding enhanced green fluorescent protein in adherent HeLa cells while exhibiting low cytotoxicity. The microbioreactor system has been particularly useful for assessing the ability of dTtaPS10+ to deliver a plasmid encoding immunoglobulin IgG1 in a fed-batch serum-free suspension CHO cell culture; dTtaPS10+ -mediated transfection resulted in the production of IgG1 in yields comparable to or better than those obtained with commercial lipid-based transfection reagents under similar conditions. The ability of dTtaPS10+ to deliver plasmids is essentially unaffected by the presence of a silicone-based antifoaming reagent, which is commonly used in bioreactor cell cultures. The transfection efficiency of lyophilized dTtaPS10+ -plasmid complexes has been significantly restored upon aqueous reconstitution when compared to that achieved while using commercial transfection reagent complexes under similar conditions. The results of all experiments underscore the potential of dTtaPS10+ for transient transfection of plasmids into adherent cells and fed-batch serum-free suspension CHO cells and rapid screening of reagents in a microbioreactor system.

A confidence building exercise in data and identifiability: Modeling cancer chemotherapy as a case study.

Mathematical modeling has a long history in the field of cancer therapeutics, and there is increasing recognition that it can help uncover the mechanisms that underlie tumor response to treatment. However, making quantitative predictions with such models often requires parameter estimation from data, raising questions of parameter identifiability and estimability. Even in the case of structural (theoretical) identifiability, imperfect data and the resulting practical unidentifiability of model parameters can make it difficult to infer the desired information, and in some cases, to yield biologically correct inferences and predictions. Here, we examine parameter identifiability and estimability using a case study of two compartmental, ordinary differential equation models of cancer treatment with drugs that are cell cycle-specific (taxol) as well as non-specific (oxaliplatin). We proceed through model building, structural identifiability analysis, parameter estimation, practical identifiability analysis and its biological implications, as well as alternative data collection protocols and experimental designs that render the model identifiable. We use the differential algebra/input-output relationship approach for structural identifiability, and primarily the profile likelihood approach for practical identifiability. Despite the models being structurally identifiable, we show that without consideration of practical identifiability, incorrect cell cycle distributions can be inferred, that would result in suboptimal therapeutic choices. We illustrate the usefulness of estimating practically identifiable combinations (in addition to the more typically considered structurally identifiable combinations) in generating biologically meaningful insights. We also use simulated data to evaluate how the practical identifiability of the model would change under alternative experimental designs. These results highlight the importance of understanding the underlying mechanisms rather than purely using parsimony or information criteria/goodness-of-fit to decide model selection questions. The overall roadmap for identifiability testing laid out here can be used to help provide mechanistic insight into complex biological phenomena, reduce experimental costs, and optimize model-driven experimentation.

An Amphipathic trans-Acting Phosphorothioate DNA Element Delivers Uncharged PNA and PMO Nucleic Acid Sequences in Mammalian Cells.

An innovative approach to the delivery of uncharged peptide nucleic acids (PNAs) and phosphorodiamidate morpholino (PMO) oligomers in mammalian cells is described and consists of extending the sequence of those oligomers with a short PNA-polyA or PMO-polyA tail. Recognition of the polyA-tailed PNA or PMO oligomers by an amphipathic trans-acting polythymidylic thiophosphate triester element (dTtaPS) results in efficient internalization of those oligomers in several cell lines. The authors' findings indicate that cellular uptake of the oligomers occurs through an energy-dependent mechanism and macropinocytosis appears to be the predominant endocytic pathway used for internalization. The functionality of the internalized oligomers is demonstrated by alternate splicing of the pre-mRNA encoding luciferase in HeLa pLuc 705 cells. Amphipathic phosphorothioate DNA elements may represent a unique class of cellular transporters for robust delivery of uncharged nucleic acid sequences in live mammalian cells. © 2016 by John Wiley & Sons, Inc.

Assessment of the cellular internalization of thermolytic phosphorothioate DNA oligonucleotide prodrugs.

The bioactivity of a CpG-containing phosphorothioate DNA oligonucleotide with thermolytic 2-(N-formyl-N-methylamino)ethyl (fma) thiophosphate groups in mice led us to investigate the parameters affecting the internalization of these thermosensitive DNA prodrugs in various cell lines. Flow cytometry and confocal microscopy analyses indicate that 5'-fluoresceinated fma-phosphorothioate DNA sequences are poorly internalized in Vero, HeLa and GC-2 cells. However, when four fma-thiophosphate groups of a 15-nucleotide long oligothymidylate prodrug are replaced with 3-(N,N-dimethylamino)prop-1-yl thiophosphate functions, internalization of the positively charged prodrug, under physiological conditions, increased fourfold in HeLa and 40-fold in Vero or GC-2 cells. No cytotoxic effects are observed in Vero cells even at an extracellular prodrug concentration of 50 µM over a period of 72 h. Confocal microscopy studies show that internalization of the positively charged oligothymidylate prodrug in Vero cells is time-dependent with early trafficking of the DNA sequence through endosomal vesicles and, eventually, to the nucleus of the cells. Thus, the incorporation of four 3-(N,N-dimethylamino)prop-1-yl thiophosphate groups into thermosentive fma-phosphorothioate DNA prodrugs is an attractive strategy for efficient cellular internalization of these nucleic acid-based drugs for potential therapeutic indications.

A hybrid model of the role of VEGF binding in endothelial cell migration and capillary formation.

Vascular endothelial growth factor (VEGF) is the most studied family of soluble, secreted mediators of endothelial cell migration, survival, and proliferation. VEGF exerts its function by binding to specific tyrosine kinase receptors on the cell surface and transducing the effect through downstream signaling. In order to study the influence of VEGF binding on endothelial cell motion, we develop a hybrid model of VEGF-induced angiogenesis, based on the theory of reinforced random walks. The model includes the chemotactic response of endothelial cells to angiogenic factors bound to cell-surface receptors, rather than approximating this as a function of extracellular chemical concentrations. This allows us to capture biologically observed phenomena such as activation and polarization of endothelial cells in response to VEGF gradients across their lengths, as opposed to extracellular gradients throughout the tissue. We also propose a novel and more biologically reasonable functional form for the chemotactic sensitivity of endothelial cells, which is also governed by activated cell-surface receptors. This model is able to predict the threshold level of VEGF required to activate a cell to move in a directed fashion as well as an optimal VEGF concentration for motion. Model validation is achieved by comparison of simulation results directly with experimental data.

Synthesis and study of cyclic pronucleotides of 5-fluoro-2'-deoxyuridine.

A one-step method for the synthesis of cyclic pronucleotide (cProTide) derivatives of 5-fluoro-2'-deoxyuridine (FdUrd), utilizing a novel phosphoramidating reagent, is described. Stereochemistry at phosphorus was established by NMR studies and modeling. Cytotoxicity data of representative cProTide derivatives of FdUrd are presented. The observed cell-to-cell variations in activity suggests that it is feasible to screen for structural variations in the cProTide moiety favoring metabolic activation in cancer cells, which may lead to an increase in the therapeutic effectiveness of FdUrd. The method described is applicable to all anticancer and antiviral nucleoside analogs having both the 5'- and the 3'-OH groups available for modification, forming cProTide derivatives capable of delivering the 5'-monophosphates to cells.

The role of phosphate in the action of thymidine phosphorylase inhibitors: Implications for the catalytic mechanism.

The design and synthesis of 5-fluoro-6-[(2-aminoimidazol-1-yl)methyl]uracil (AIFU), a potent inhibitor of thymidine phosphorylase (TP) with K(i)-values of 11nM (ecTP) and 17nM (hTP), are described. Kinetic studies established that the type of inhibition of TP by AIFU is uncompetitive with respect to inorganic phosphate (or arsenate). The results obtained suggest that AIFU and other zwitterionic thymine analog inhibitors of TP act as transition state analogs, mimicking the anionic thymine leaving group, consistent with an S(N)2-type catalytic mechanism, and anchored by their protonated side chains to the enzyme-bound phosphate by electrostatic and H-bonding interactions.

Modeling the VEGF-Bcl-2-CXCL8 pathway in intratumoral agiogenesis.

Recent experiments show that vascular endothelial growth factor (VEGF) is the crucial mediator of downstream events that ultimately lead to enhanced endothelial cell survival and increased vascular density within many tumors. The newly discovered pathway involves up-regulation of the anti-apoptotic protein Bcl-2, which in turn leads to increased production of interleukin-8 (CXCL8). The VEGF-Bcl-2-CXCL8 pathway suggests new targets for the development of anti-angiogenic strategies including short interfering RNA (siRNA) that silence the CXCL8 gene and small molecule inhibitors of Bcl-2. In this paper, we present and validate a mathematical model designed to predict the effect of the therapeutic blockage of VEGF, CXCL8, and Bcl-2 at different stages of tumor progression. In agreement with experimental observations, the model predicts that curtailing the production of CXCL8 early in development can result in a delay in tumor growth and vascular development; however, it has little effect when applied at late stages of tumor progression. Numerical simulations also show that blocking Bcl-2 up-regulation, either at early stages or after the tumor has fully developed, ensures that both microvascular and tumor cell density stabilize at low values representing growth control. These results provide insight into those aspects of the VEGF-Bcl-2-CXCL8 pathway, which independently and in combination, are crucial mediators of tumor growth and vascular development. Continued quantitative modeling in this direction may have profound implications for the development of novel therapies directed against specific proteins and chemokines to alter tumor progression.

Level of endothelial cell apoptosis required for a significant decrease in microvessel density.

Endothelial cell apoptosis plays a critical role in the disruption of blood vessels mediated by natural inhibitors of angiogenesis and by anti-vascular drugs. However, the proportion of endothelial cells required to mediate a significant decrease in microvessel density is unknown. A system based on an inducible caspase (iCaspase-9) offers a unique opportunity to address this question. The dimerizer drug AP20187 induces apoptosis of human dermal microvascular endothelial cells stably transduced with iCaspase-9 (HDMEC-iCaspase-9), but not control cells (HDMEC-LXSN). Here, we generated blood vessels containing several HDMEC-iCaspase-9:HDMEC-LXSN ratios, and developed a mathematical modeling involving a system of differential equations to evaluate experimentally inaccessible ratios. A significant decrease in capillary sprouts was observed when at least 17% of the endothelial cells underwent apoptosis in vitro. Exposure to vascular endothelial growth factor (VEGF(165)) did not prevent apoptosis of HDMEC-iCaspase-9, but increased the apoptotic requirement for sprout disruption. In vivo experiments showed the requirement of at least 22% apoptotic endothelial cells for a significant decrease in microvascular density. The combined use of biological experimentation with mathematical modeling allowed us to conclude that apoptosis of a relatively small proportion of endothelial cells is sufficient to mediate a significant decrease in microvessel density.