Multiscale modeling has emerged as a powerful approach to interpret and capitalize on the biological complexity underlying blood vessel growth. We present a multiscale model of angiogenesis that heralds the start of a large scale initiative to integrate related biological models. The goal of the integrative project is to better understand underlying biological mechanisms from the molecular level up through the organ systems level, and test new therapeutic strategies. Model methodology includes ordinary and partial differential equations, stochastic models, complex logical rules, and agent-based architectures. Current modules represent blood flow, oxygen transport, growth factor distribution and signaling, cell sensing, cell movement and cell proliferation. Challenges of integration lie in connecting modules that are diversely designed, seamlessly coordinating feedback, and representing spatial and time scales from ligand-receptor interactions and intracellular signaling, to cell-level movement and cell-matrix interactions, to vessel branching and capillary network formation, to tissue level characteristics, to organ system response. We briefly introduce the individual modules, discuss our approach to integration, present initial results from the coordination of modules, and propose solutions to some critical issues facing angiogenesis multiscale modeling and integration.
Multiscale modeling provides a powerful and quantitative platform for investigating the complexity inherent in intracellular signaling pathways and rationalizing the effects of molecular perturbations on downstream signaling events and ultimately, on the cell phenotype. Here we describe the application of a multiscale modeling scheme to the HER3/ErbB3 receptor tyrosine kinase (RTK) signaling network, which regulates critical cellular processes including proliferation, migration and differentiation. The HER3 kinase is a topic of current interest and investigation, as it has been implicated in mechanisms of resistance to tyrosine kinase inhibition (TKI) of EGFR and HER2 in the treatment of many human malignancies. Moreover, the commonly regarded status of HER3 as a catalytically inactive ‘pseudokinase’ has recently been challenged by our previous study, which demonstrated robust residual kinase activity for HER3. Through our multiscale model, we investigate the most significant molecular interactions that contribute to potential mechanisms of HER3 activity and the physiological relevance of this activity to mechanisms of drug resistance in an ErbB-driven tumor cell in silico. The results of our molecular-scale simulations support the characterization of HER3 as a weakly active kinase that, in contrast to its fully-active ErbB family members, depends upon a unique hydrophobic interface to coordinate the alignment of specific catalytic residues required for its activity. Translating our molecular simulation results of the uniquely active behavior of the HER3 kinase into a physiologically relevant environment, our HER3 signaling model demonstrates that even a weak level of HER3 activity may be sufficient to induce AKT signaling and TKI resistance in the context of an ErbB signaling-dependent tumor cell, and therefore therapeutic targeting of HER3 may represent a superior treatment strategy for specific ErbB-driven cancers.
Viruses are infectious agents that can cause epidemics and pandemics. The understanding of virus formation, evolution, stability, and interaction with host cells is of great importance to the scientific community and public health. Typically, a virus complex in association with its aquatic environment poses a fabulous challenge to theoretical description and prediction. In this work, we propose a differential geometry-based multiscale paradigm to model complex biomolecule systems. In our approach, the differential geometry theory of surfaces and geometric measure theory are employed as a natural means to couple the macroscopic continuum domain of the fluid mechanical description of the aquatic environment from the microscopic discrete domain of the atomistic description of the biomolecule. A multiscale action functional is constructed as a unified framework to derive the governing equations for the dynamics of different scales. We show that the classical Navier-Stokes equation for the fluid dynamics and Newton's equation for the molecular dynamics can be derived from the least action principle. These equations are coupled through the continuum-discrete interface whose dynamics is governed by potential driven geometric flows.
Large chemical and biological systems such as fuel cells, ion channels, molecular motors, and viruses are of great importance to the scientific community and public health. Typically, these complex systems in conjunction with their aquatic environment pose a fabulous challenge to theoretical description, simulation, and prediction. In this work, we propose a differential geometry based multiscale paradigm to model complex macromolecular systems, and to put macroscopic and microscopic descriptions on an equal footing. In our approach, the differential geometry theory of surfaces and geometric measure theory are employed as a natural means to couple the macroscopic continuum mechanical description of the aquatic environment with the microscopic discrete atom-istic description of the macromolecule. Multiscale free energy functionals, or multiscale action functionals are constructed as a unified framework to derive the governing equations for the dynamics of different scales and different descriptions. Two types of aqueous macromolecular complexes, ones that are near equilibrium and others that are far from equilibrium, are considered in our formulations. We show that generalized Navier–Stokes equations for the fluid dynamics, generalized Poisson equations or generalized Poisson–Boltzmann equations for electrostatic interactions, and Newton's equation for the molecular dynamics can be derived by the least action principle. These equations are coupled through the continuum-discrete interface whose dynamics is governed by potential driven geometric flows. Comparison is given to classical descriptions of the fluid and electrostatic interactions without geometric flow based micro-macro interfaces. The detailed balance of forces is emphasized in the present work. We further extend the proposed multiscale paradigm to micro-macro analysis of electrohydrodynamics, electrophoresis, fuel cells, and ion channels. We derive generalized Poisson–Nernst–Planck equations that are coupled to generalized Navier–Stokes equations for fluid dynamics, Newton's equation for molecular dynamics, and potential and surface driving geometric flows for the micro-macro interface. For excessively large aqueous macromolecular complexes in chemistry and biology, we further develop differential geometry based multiscale fluid-electro-elastic models to replace the expensive molecular dynamics description with an alternative elasticity formulation.
Variational principle; Multiscale; Geometric flows; Solvation analysis; Electrostatic analysis; Implicit solvent models; Molecular dynamics; Elasticity; Navier–Stokes equation; Poisson–Boltzmann equation; Nernst–Planck equation
Glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor. Standard treatment for GBM patients is surgery followed by radiotherapy and/or chemotherapy, but tumors always recur. Traditional therapies seem to fail because they eliminate only the bulk of the tumors and spare a population of stem-like cells termed tumor-initiating cells. The stem-like state and preferential activation of DNA damage response in the GBM tumor-initiating cells contribute to their radio-resistance and recurrence. The molecular mechanisms underlying this efficient activation of damage response and maintenance of stem-like state remain elusive. Here we show that RBM14 controls DNA repair pathways and also prevents cell differentiation in GBM spheres, causing radio-resistance. Knockdown of RBM14 affects GBM sphere maintenance and sensitizes radio-resistant GBM cells at the cellular level. We demonstrate that RBM14 knockdown blocks GBM regrowth after irradiation in vivo. In addition, RBM14 stimulates DNA repair by controlling the DNA-PK-dependent non-homologous end-joining (NHEJ) pathway. These results reveal unexpected functions of the RNA-binding protein RBM14 in control of DNA repair and maintenance of tumor-initiating cells. Targeting the RBM14-dependent pathway may prevent recurrence of tumors and eradicate the deadly disease completely.
GBM; tumor-initiating cell; DNA repair; stem cell maintenance; radio-therapy
In vivo and in vitro studies give a paradoxical picture of the actions of the key regulatory factor TGF-β1 in epidermal wound healing with it stimulating migration of keratinocytes but also inhibiting their proliferation. To try to reconcile these into an easily visualized 3D model of wound healing amenable for experimentation by cell biologists, a multiscale model of the formation of a 3D skin epithelium was established with TGF-β1 literature–derived rule sets and equations embedded within it. At the cellular level, an agent-based bottom-up model that focuses on individual interacting units (keratinocytes) was used. This was based on literature-derived rules governing keratinocyte behavior and keratinocyte/ECM interactions. The selection of these rule sets is described in detail in this paper. The agent-based model was then linked with a subcellular model of TGF-β1 production and its action on keratinocytes simulated with a complex pathway simulator. This multiscale model can be run at a cellular level only or at a combined cellular/subcellular level. It was then initially challenged (by wounding) to investigate the behavior of keratinocytes in wound healing at the cellular level. To investigate the possible actions of TGF-β1, several hypotheses were then explored by deliberately manipulating some of these rule sets at subcellular levels. This exercise readily eliminated some hypotheses and identified a sequence of spatial-temporal actions of TGF-β1 for normal successful wound healing in an easy-to-follow 3D model. We suggest this multiscale model offers a valuable, easy-to-visualize aid to our understanding of the actions of this key regulator in wound healing, and provides a model that can now be used to explore pathologies of wound healing.
Mathematical modeling of angiogenesis has been gaining momentum as a means to shed new light on the biological complexity underlying blood vessel growth. A variety of computational models have been developed, each focusing on different aspects of the angiogenesis process and occurring at different biological scales, ranging from the molecular to the tissue levels. Integration of models at different scales is a challenging and currently unsolved problem.
We present an object-oriented module-based computational integration strategy to build a multiscale model of angiogenesis that links currently available models. As an example case, we use this approach to integrate modules representing microvascular blood flow, oxygen transport, vascular endothelial growth factor transport and endothelial cell behavior (sensing, migration and proliferation). Modeling methodologies in these modules include algebraic equations, partial differential equations and agent-based models with complex logical rules. We apply this integrated model to simulate exercise-induced angiogenesis in skeletal muscle. The simulation results compare capillary growth patterns between different exercise conditions for a single bout of exercise. Results demonstrate how the computational infrastructure can effectively integrate multiple modules by coordinating their connectivity and data exchange. Model parameterization offers simulation flexibility and a platform for performing sensitivity analysis.
This systems biology strategy can be applied to larger scale integration of computational models of angiogenesis in skeletal muscle, or other complex processes in other tissues under physiological and pathological conditions.
Despite rapid advances in the field of stem/progenitor cells through experimental studies, relevant modeling approaches have not progressed with a similar pace. Various models have focused on particular aspects of stem cell physiology including gene regulatory networks, gene expression noise and signaling cascades activated by exogenous factors. However, the self-renewal and differentiation of stem cells is driven by the coordinated regulation of events at the subcellular, intercellular and milieu levels. Such events also span multiple time domains from the fast molecular reactions governing gene expression to the slower cell cycle and division. Thus, the development of multiscale computational frameworks for stem cell populations is highly desirable. Multiscale models are expected to aid the design of efficient differentiation strategies and bioprocesses for the generation of therapeutically useful stem cell progeny. Yet, challenges in making these models tractable and pairing those to sufficient experimental data prevent their wide adoption by the stem cell community. Here, we review modeling approaches reported for stem cell populations and associated hurdles.
Stem cells; self-renewal; differentiation; multiscale model; stochastic model; deterministic model; population balance model
Glioblastoma multiforme (GBM) is a highly invasive and chemoradioresistant brain malignancy. Temozolomide (TMZ), a DNA-alkylating agent, is effective against GBM and has become the standard first-line drug. However, the mechanism by which TMZ regulates the progression of GBM remains elusive. Here, we demonstrate that TMZ targets TAp63, a p53 family member, inducing its expression to suppress the progression of human GBM. High levels of TAp63 expression in GBM tissues after TMZ treatment was an indicator of favourable prognosis. In human GBM cells, TMZ-induced TAp63 directly repressed MYC transcription. Activation of this TAp63-MYC pathway by TMZ inhibited human GBM progression both in vitro and in vivo. Furthermore, downregulation of MYC mRNA levels in recurrent GBMs after TMZ treatment correlated with better patient survival. Therefore, our results suggest that the TAp63-mediated transcriptional repression of MYC is a novel pathway regulating TMZ efficacy in GBM.
The epidermal growth factor receptor (EGFR) is frequently overexpressed in many cancers, including non-small cell lung cancer (NSCLC). In silico modeling is considered to be an increasingly promising tool to add useful insights into the dynamics of the EGFR signal transduction pathway. However, most of the previous modeling work focused on the molecular or the cellular level only, neglecting the crucial feedback between these scales as well as the interaction with the heterogeneous biochemical microenvironment.
We developed a multiscale model for investigating expansion dynamics of NSCLC within a two-dimensional in silico microenvironment. At the molecular level, a specific EGFR-ERK intracellular signal transduction pathway was implemented. Dynamical alterations of these molecules were used to trigger phenotypic changes at the cellular level. Examining the relationship between extrinsic ligand concentrations, intrinsic molecular profiles and microscopic patterns, the results confirmed that increasing the amount of available growth factor leads to a spatially more aggressive cancer system. Moreover, for the cell closest to nutrient abundance, a phase-transition emerges where a minimal increase in extrinsic ligand abolishes the proliferative phenotype altogether.
Our in silico results indicate that in NSCLC, in the presence of a strong extrinsic chemotactic stimulus (and depending on the cell's location) downstream EGFR-ERK signaling may be processed more efficiently, thereby yielding a migration-dominant cell phenotype and overall, an accelerated spatio-temporal expansion rate.
Optimizing anticancer therapeutics needs to account for variable drug responses in heterogeneous cell populations within the tumor as well as in organs of toxicity. To address cell heterogeneity, we propose a multiscale modeling approach—from in vitro to preclinical and clinical studies—to develop cell-type–specific pharmacokinetic–pharmacodynamic (PK-PD) models. A physiologically based mechanistic modeling approach integrating data from aqueous solutions, U87 glioma cells, mice, and cancer patients was utilized to characterize the brain disposition of temozolomide (TMZ), the cornerstone of chemotherapy against glioblastoma multiforme. The final model represented intracellular normal brain and brain tumor compartments in which TMZ pH-dependent conversion to the DNA-alkylating species leads to the formation of DNA adducts that serve as an entry point for a PD model. This multiscale protocol can be extended to account for TMZ PK-PD in different cell populations, thus providing a critical tool to personalize TMZ-based chemotherapy on a cell-type–specific basis.
Glioblastoma multiforme (GBM), an aggressive brain tumor of astrocytic/neural stem cell origin, represents one of the most incurable cancers. GBM tumors are highly heterogeneous. However, most tumors contain a subpopulation of cells that display neural stem cell characteristics in vitro and that can generate a new brain tumor upon transplantation in mice. Hence, previously identified molecular pathways regulating neural stem cell biology were found to represent the cornerstone of GBM stem cell self-renewal mechanism. GBM tumors are also notorious for their resistance to radiation therapy. Notably, GBM “cancer stem cells” were also found to be responsible for this radioresistance. Herein, we will analyze the data supporting or not the cancer stem cell model in GBM, overview the current knowledge regarding GBM stem cell self-renewal and radioresistance molecular mechanisms, and discuss the potential therapeutic application of these findings.
polycomb, BMI1; cancer stem cell; glioma; glioblastoma multiforme; astrocyte; astrocytoma; radioresistance; CD133; prominin
The complexity in intracellular signaling mechanisms relevant for the conquest of many diseases resides at different levels of organization with scales ranging from the subatomic realm relevant to catalytic functions of enzymes to the mesoscopic realm relevant to the cooperative association of molecular assemblies and membrane processes. Consequently, the challenge of representing and quantifying functional or dysfunctional modules within the networks remains due to the current limitations in our understanding of mesoscopic biology, i.e., how the components assemble into functional molecular ensembles. A multiscale approach is necessary to treat a hierarchy of interactions ranging from molecular (nm, ns) to signaling (μm, ms) length and time scales, which necessitates the development and application of specialized modeling tools. Complementary to multiscale experimentation (encompassing structural biology, mechanistic enzymology, cell biology, and single molecule studies) multiscale modeling offers a powerful and quantitative alternative for the study of functional intracellular signaling modules. Here, we describe the application of a multiscale approach to signaling mediated by the ErbB1 receptor which constitutes a network hub for the cell’s proliferative, migratory, and survival programs. Through our multiscale model, we mechanistically describe how point-mutations in the ErbB1 receptor can profoundly alter signaling characteristics leading to the onset of oncogenic transformations. Specifically, we describe how the point mutations induce cascading fragility mechanisms at the molecular scale as well as at the scale of the signaling network to preferentially activate the survival factor Akt. We provide a quantitative explanation for how the hallmark of preferential Akt activation in cell-lines harboring the constitutively active mutant ErbB1 receptors causes these cell-lines to be addicted to ErbB1-mediated generation of survival signals. Consequently, inhibition of ErbB1 activity leads to a remarkable therapeutic response in the addicted cell lines.
Glioblastoma multiforme (GBM) is an aggressive brain malignancy characterized by high heterogeneity and invasiveness. It is increasingly accepted that the refractory feature of GBM to current therapies stems from the existence of few tumorigenic cells that sustain tumor growth and spreading, the so-called glioma-initiating cells (GICs). Previous studies showed that cytokines of the bone morphogenetic protein (BMP) family induce differentiation of the GICs, and thus act as tumor suppressors. Molecular pathways that explain this behavior of BMP cytokines remain largely elusive. Here, we show that BMP signaling induces Smad-dependent expression of the transcriptional regulator Snail in a rapid and sustained manner. Consistent with its already established promigratory function in other cell types, we report that Snail silencing decreases GBM cell migration. Consequently, overexpression of Snail increases GBM invasiveness in a mouse xenograft model. Surprisingly, we found that Snail depletes the GBM capacity to form gliomaspheres in vitro and to grow tumors in vivo, both of which are important features shared by GICs. Thus Snail, acting downstream of BMP signaling, dissociates the invasive capacity of GBM cells from their tumorigenic potential.
BMP; glioblastoma multiforme; glioma-initiating cells; invasiveness; snail
Summary: Morpheus is a modeling environment for the simulation and integration of cell-based models with ordinary differential equations and reaction-diffusion systems. It allows rapid development of multiscale models in biological terms and mathematical expressions rather than programming code. Its graphical user interface supports the entire workflow from model construction and simulation to visualization, archiving and batch processing.
Availability and implementation: Binary packages are available at http://imc.zih.tu-dresden.de/wiki/morpheus for Linux, Mac OSX and MS Windows.
Supplementary data are available at Bioinformatics online.
Human glioblastoma (GBM) cells are notorious for their resistance to apoptosis-inducing therapeutics. We have identified lanatoside C as a sensitizer of GBM cells to tumor necrosis factor–related apoptosis-inducing ligand (TRAIL)–induced cell death partly by upregulation of the death receptor 5. We show that lanatoside C sensitizes GBM cells to TRAIL-induced apoptosis in a GBM xenograft model in vivo. Lanatoside C on its own serves as a therapeutic agent against GBM by activating a caspase-independent cell death pathway. Cells treated with lanatoside C showed necrotic cell morphology with absence of caspase activation, low mitochondrial membrane potential, and early intracellular ATP depletion. In conclusion, lanatoside C sensitizes GBM cells to TRAIL-induced cell death and mitigates apoptosis resistance of glioblastoma cells by inducing an alternative cell death pathway. To our knowledge, this is one of the first examples of use of caspase-independent cell death inducers to trigger tumor regression in vivo. Activation of such mechanism may be a useful strategy to counter resistance of cancer cells to apoptosis.
cardiac glycoside; glioblastoma; lanatoside C; non-apoptotic cell death; TRAIL
Tumorigenic potential of glioblastoma multiforme (GBM) cells is, in part, attributable to their undifferentiated (neural stem cell-like) phenotype. Astrocytic differentiation of GBM cells is associated with transcriptional induction of Glial Fibrillary Acidic Protein (GFAP) and repression of Nestin, whereas the reciprocal transcription program operates in undifferentiated GBM cells. The molecular mechanisms underlying the regulation of these transcription programs remain elusive. Here, we show that the transcriptional co-activator p300 was expressed in GBM tumors and cell lines and acted as an activator of the GFAP gene and a repressor of the Nestin gene. On the other hand, Myc (formerly known as c-Myc overrode these p300 functions by repressing the GFAP gene and inducing the Nestin gene in GBM cells. Moreover, RNAi-mediated inhibition of p300 expression significantly enhanced the invasion potential of GBM cells in vitro. Taken together, these data suggest that dedifferentiated/undifferentiated GBM cells are more invasive than differentiated GBM cells. Because invasion is a major cause of GBM morbidity, differentiation therapy may improve the clinical outcome.
glioblastoma multiforme (GBM); GBM stem cell; differentiation; invasion; p300; Myc
Inhibition of VEGF signaling leads to a pro-invasive phenotype in mouse models of glioblastoma multiforme (GBM) and in a subset of GBM patients treated with bevacizumab. Here we demonstrate that vascular endothelial growth factor (VEGF) directly and negatively regulates tumor cell invasion through enhanced recruitment of the protein tyrosine phosphatase 1B (PTP1B) to a MET/VEGFR2 heterocomplex, thereby suppressing HGF-dependent MET phosphorylation and tumor cell migration. Consequently, VEGF blockade restores and increases MET activity in GBM cells in a hypoxia-independent manner, while inducing a program reminiscent of epithelial-to-mesenchymal transition highlighted by a T-cadherin to N-cadherin switch and enhanced mesenchymal features. Inhibition of MET in GBM mouse models blocks mesenchymal transition and invasion provoked by VEGF ablation, resulting in substantial survival benefit.
Adaptation to hypoxia and consequent pro-inflammatory gene expression of prostate and breast carcinomas have been implicated in the progression toward cancer malignant phenotype. Only partial data are available for the human tumor glioblastoma multiforme (GBM). The aim of our study was to analyze the hypoxic and pro-inflammatory microenvironment in GBMs and to demonstrate that in a stem/progenitor cell line derived from human glioblastoma (GBM-SCs), hypoxia activates a coordinated inflammatory response, evidencing an invasive and migratory phenotype.
From each of 10 human solid glioblastomas, clinically and histopathologically characterized, we obtained three surgical samples taken from the center and the periphery of the tumor, and from adjacent host normal tissue. Molecular and morphological analyses were carried out using quantitative real-time PCR and western blot (WB). GBM stem and differentiated cells were incubated under hypoxic conditions and analyzed for pro-inflammatory gene expression and for invasive/migratory behavior.
A panel of selected representative pro-inflammatory genes (RAGE and P2X7R, COX2, NOS2 and, PTX3) were analyzed, comparing tumor, peritumor and host normal tissues. Tumors containing leukocyte infiltrates (as assessed using CD45 immunohistochemistry) were excluded. Selected genes were overexpressed in the central regions of the tumors (i.e. in the more hypoxic areas), less expressed in peripheral regions, and poorly expressed or absent in adjacent normal host tissues. Western blot analysis confirmed that the corresponding pro-inflammatory proteins were also differently expressed. Hypoxic stem cell lines showed a clear time-dependent activation of the entire panel of pro-inflammatory genes as compared to differentiated tumor cells. Biological assays showed that invasive and migratory behavior was strengthened by hypoxia only in GBM stem cells.
In human solid glioblastoma we have observed a coordinated overexpression of a panel of pro-inflammatory genes as compared to host normal tissue. We have also evidenced a similar pattern of overexpressed genes in GBM-SCs after hypoxic treatment, showing also a gain of invasive and migratory function that was lost when these stem cells differentiated. We suggest that, as has been previously described for prostatic and mammary carcinoma, in human glioblastoma acquisition of a proinflammatory phenotype may be relevant for malignant progression.
Glioblastoma multiforme (GBM) has a dismal prognosis despite aggressive therapy. Initial diagnosis and measurement of response to treatment is usually determined by measurement of gadolinium-enhanced tumor volume with magnetic resonance imaging (MRI). Unfortunately, many GBM treatment modalities can cause changes in tumor gadolinium enhancement patterns that mimic tumor progression. The lack of a definitive imaging modality to distinguish posttreatment radiographic imaging changes (PTRIC), including pseudoprogression and radiation necrosis, from true tumor progression presents a major unmet clinical need in the management of GBM patients.
The authors discuss current modalities available for differentiating PTRIC and tumor progression, describe development of an animal model of PTRIC, and consider potential molecular and cellular pathways involved in the development of PTRIC.
An animal model using glioma cells transfected with a luciferase reporter has been developed, and after conventional GBM therapy, this animal model can be evaluated with posttreatment bioluminescence imaging and various MR tumor imaging modalities.
Posttreatment radiographic changes that mimic tumor progression can influence clinicians to make treatment decisions that are inappropriate for the patient's actual clinical condition. Several imaging modalities have been used to try to distinguish PTRIC and true progression, including conventional MRI, perfusion MRI, MR spectroscopy, and positron emission tomography (PET); however, none of these modalities has consistently and reliably distinguished PTRIC from tumor growth. An animal model using glioma cells transfected with a luciferase reporter may enable mechanistic studies to determine causes and potential treatments for PTRIC.
Brain; chemotherapy; imaging; models; pseudoprogression; radiation injury; rat
The morphogenetic movements, and the embryonic phenotypes they ultimately produce, are the consequence of a series of events that involve signaling pathways, cytoskeletal components, and cell- and tissue-level mechanical interactions. In order to better understand how these events work together in the context of amphibian neurulation, an existing multiscale computational model was augmented. Geometric data for this finite element-based mechanical model were obtained from 3D surface reconstructions of live axolotl embryos and serial sections of fixed specimens. Tissue mechanical properties were modeled using cell-based constitutive equations that include internal force generation and cell rearrangement, and equation parameters were adjusted manually to reflect biochemical changes including alterations in Shroom or the planar-cell-polarity pathway. The model indicates that neural tube defects can arise when convergent extension of the neural plate is reduced by as little as 20%, when it is eliminated on one side of the embryo, when neural ridge elevation is disrupted, when tension in the non-neural ectoderm is increased, or when the ectoderm thickness is increased. Where comparable conditions could be induced in Xenopus embryos, good agreement was found, an important step in model validation. The model reveals the neurulating embryo to be a finely tuned biomechanical system.
Increasing the sensitivity of glioblastoma cells to radiation is a promising approach to improve survival in patients with glioblastoma multiforme (GBM). This study aims to determine if serine/threonine phosphatase (protein phosphatase 6 (PP6)) is a molecular target for GBM radiosensitization treatment. The GBM orthotopic xenograft mice model was used in this study. Our data demonstrated that the protein level of PP6 catalytic subunit (PP6c) was upregulated in the GBM tissue from about 50% patients compared with the surrounding tissue or control tissue. Both the in vitro survival fraction of GBM cells and the patient survival time were highly correlated or inversely correlated with PP6c expression (R2=0.755 and −0.707, respectively). We also found that siRNA knockdown of PP6c reduced DNA-dependent protein kinase (DNA-PK) activity in three different GBM cell lines, increasing their sensitivity to radiation. In the orthotopic mice model, the overexpression of PP6c in GBM U87 cells attenuated the effect of radiation treatment, and reduced the survival time of mice compared with the control mice, while the PP6c knocking-down improved the effect of radiation treatment, and increased the survival time of mice. These findings demonstrate that PP6 regulates the sensitivity of GBM cells to radiation, and suggest small molecules disrupting or inhibiting PP6 association with DNA-PK is a potential radiosensitizer for GBM.
PP6; GBM; radiation resistance; DNA-PK
Proton transport is one of the most important and interesting phenomena in living cells. The present work proposes a multiscale/multiphysics model for the understanding of the molecular mechanism of proton transport in transmembrane proteins. We describe proton dynamics quantum mechanically via a density functional approach while implicitly model other solvent ions as a dielectric continuum to reduce the number of degrees of freedom. The densities of all other ions in the solvent are assumed to obey the Boltzmann distribution. The impact of protein molecular structure and its charge polarization on the proton transport is considered explicitly at the atomic level. We formulate a total free energy functional to put proton kinetic and potential energies as well as electrostatic energy of all ions on an equal footing. The variational principle is employed to derive nonlinear governing equations for the proton transport system. Generalized Poisson-Boltzmann equation and Kohn-Sham equation are obtained from the variational framework. Theoretical formulations for the proton density and proton conductance are constructed based on fundamental principles. The molecular surface of the channel protein is utilized to split the discrete protein domain and the continuum solvent domain, and facilitate the multiscale discrete/continuum/quantum descriptions. A number of mathematical algorithms, including the Dirichlet to Neumann mapping, matched interface and boundary method, Gummel iteration, and Krylov space techniques are utilized to implement the proposed model in a computationally efficient manner. The Gramicidin A (GA) channel is used to demonstrate the performance of the proposed proton transport model and validate the efficiency of proposed mathematical algorithms. The electrostatic characteristics of the GA channel is analyzed with a wide range of model parameters. The proton conductances are studied over a number of applied voltages and reference concentrations. A comparison with experimental data verifies the present model predictions and validates the proposed model.
Proton transport; quantum dynamics in continuum; multiscale model; Poisson-Boltzmann equation; generalized Kohn-Sham equation; variational principle
Recurrent glioblastomas (rGBM) invariably relapse after initial response to anti-VEGF therapy. There are two prevailing hypotheses on how these tumors escape antiangiogenic therapy: switch to VEGF-independent angiogenic pathways and vessel co-option. However, direct evidence in rGBM patients is lacking. Thus, we compared molecular, cellular and vascular parameters in autopsy tissues from five rGBM patients who had been treated with the pan-VEGF receptor tyrosine kinase inhibitor cediranib versus seven patients who received no therapy or chemoradiation but no antiangiogenic agents. After cediranib treatment, endothelial proliferation and glomeruloid vessels were decreased, and vessel diameters and perimeters were reduced to levels comparable to the unaffected contralateral brain hemisphere. In addition, tumor endothelial cells expressed molecular markers specific to the blood-brain barrier, indicative of a lack of revascularization despite the discontinuation of therapy. Surprisingly, in cediranib-treated GBM cellular density in the central area of the tumor was lower than in control cases and gradually decreased towards the infiltrating edge, indicative of a change in growth pattern of rGBMs after cediranib treatment, unlike that after chemo-radiation. Finally, cediranib treated GBMs showed high levels of PDGF-C and c-Met expression and infiltration by myeloid cells, which may potentially contribute to resistance to anti-VEGF therapy. In summary, we show that rGBMs switch their growth pattern after anti-VEGF therapy – characterized by lower tumor cellularity in the central area, decreased pseudopalisading necrosis and blood vessels with normal molecular expression and morphology without a second wave of angiogenesis.
Computational models of signal transduction face challenges of scale below the resolution of a single cell. Here, we organize these challenges around three key interfaces for multiscale models of cell signaling: molecules to pathways, pathways to networks, and networks to outcomes. Each interface requires its own set of computational approaches and systems-level data, and no single approach or dataset can effectively bridge all three interfaces. This suggests that realistic “whole-cell” models of signaling will need to agglomerate different model types that span critical intracellular scales. Future multiscale models will be valuable for understanding the impact of signaling mutations or population variants that lead to cellular diseases such as cancer.
Systems biology; Multiscale modeling; Stochastic; Molecular biology