A key aim of the cardiac Physiome Project is to develop theoretical models to simulate the functional behaviour of the heart under physiological and pathophysiological conditions. Heart function is critically dependent on the delivery of an adequate blood supply to the myocardium via the coronary vasculature. Key to this critical function of the coronary vasculature is system dynamics that emerge via the interactions of the numerous constituent components at a range of spatial and temporal scales. Here, we focus on several components for which theoretical approaches can be applied, including vascular structure and mechanics, blood flow and mass transport, flow regulation, angiogenesis and vascular remodelling, and vascular cellular mechanics. For each component, we summarise the current state of the art in model development, and discuss areas requiring further research. We highlight the major challenges associated with integrating the component models to develop a computational tool that can ultimately be used to simulate the responses of the coronary vascular system to changing demands and to diseases and therapies.
Vascular structure; Mechanics; Haemodynamics; Mass transport; Regulation; Adaptation; Mathematical and computational model; Multi-scale; Cellular mechanics; Integration
The evolution in our understanding of tumor angiogenesis has been the result of pioneering imaging and computational modeling studies spanning the endothelial cell, microvasculature and tissue levels. Many of these primary data on the tumor vasculature are in the form of images from pre-clinical tumor models that provide a wealth of qualitative and quantitative information in many dimensions and across different spatial scales. However, until recently, the visualization of changes in the tumor vasculature across spatial scales remained a challenge due to a lack of techniques for integrating micro- and macroscopic imaging data. Furthermore, the paucity of three-dimensional (3-D) tumor vascular data in conjunction with the challenges in obtaining such data from patients presents a serious hurdle for the development and validation of predictive, multiscale computational models of tumor angiogenesis. In this review, we discuss the development of multiscale models of tumor angiogenesis, new imaging techniques capable of reproducing the 3-D tumor vascular architecture with high fidelity, and the emergence of “image-based models”of tumor blood flow and molecular transport. Collectively, these developments are helping us gain a fundamental understanding of the cellular and molecular regulation of tumor angiogenesis that will benefit the development of new cancer therapies. Eventually, we expect this exciting integration of multiscale imaging and mathematical modeling to have widespread application beyond the tumor vasculature to other diseases involving a pathological vasculature, such as stroke and spinal cord injury.
Angiogenesis; Tumor; Vasculature; Multiscale; Imaging; Mathematical modeling; Computational modeling; Cancer; Tumor microenvironment; Systems biology
The motion of a suspension of red blood cells (RBCs) flowing in a Y-shaped bifurcating microfluidic channel is investigated using a validated low-dimensional RBC (LD-RBC) model based on dissipative particle dynamics (DPD). Specifically, the RBC is represented as a closed torus-like ring of ten colloidal particles, which leads to efficient simulations of blood flow in microcirculation over a wide range of hematocrits. Adaptive no-slip wall boundary conditions were implemented to model hydrodynamic flow within a specific wall structure of diverging 3D microfluidic channels, paying attention to controlling density fluctuations. Plasma skimming and the all-or-nothing phenomenon of RBCs in a bifurcating microfluidic channel have been investigated in our simulations for healthy and diseased blood, including the size of cell-free layer on the daughter branches. The feed hematocrit level in the parent channel has considerable influence on blood-plasma separation. Compared to the blood-plasma separation efficiencies of healthy RBCs, malaria-infected stiff RBCs (iRBCs) have a tendency to travel into the low flowrate daughter branch because of their different initial distribution in the parent channel. Our simulation results are consistent with previously published experimental results and theoretical predictions.
plasma skimming; Zweifach-Fung effect; cell-free layer; malaria; DPD
Peptides are receiving increased attention as therapeutic agents, due to their high binding specificity and versatility to be modified as targeting or carrier molecules. Particularly, peptides with anti-angiogenic activity are of high interest due to their applicability to a wide range of cancers. In this study we investigate the biological activity of two novel antiangiogenic peptides in pre-clinical glioma models. One peptide SP2000 is derived from collagen IV and the other peptide SP3019 belongs to the CXC family. We previously characterized the capacity of SP2000 and SP3019 to inhibit multiple biological endpoints linked to angiogenesis in human endothelial cells in several assays. Here we report additional studies using endothelial cells and focus on the activity of these peptides against human glioma cell growth, migration and adhesion in vitro and growth as tumor xenografts in vivo. We found that SP2000 completely inhibits migration of the glioma cells at 50 μM and SP3019 produced 50% inhibition at 100 μM. Their relative anti-adhesion activities were similar with SP2000 and SP3019 generating 50% adhesion inhibition at 4.9 ± 0.82 μM and 21.3 ± 5.92 μM respectively. In vivo glioma growth inhibition was 63 % for SP2000 and 76% for SP3019 after 2 weeks of administration at daily doses of 10mg/kg and 20 mg/kg, respectively. The direct activity of these peptides against glioma cells in conjunction with their anti-angiogenic activities warrants their further development as either stand-alone agents or in combination with standard cytotoxic or emerging targeted therapies in malignant brain tumors.
Angiogenesis; cancer therapy; endothelial cell; glioblastoma; proliferation; migration; adhesion
Structure-activity relationship (SAR) studies are essential in the generation of peptides with enhanced activity and efficacy as therapeutic agents. In this study we report a SAR study for a family of mimetic peptides derived from type IV collagen with potent anti-angiogenic properties. The SAR study was conducted using a number of validated in vitro assays including cell proliferation, adhesion, migration and tubule formation. We report a critical sequence (NINNV) within this peptide series which is required for the potent anti-angiogenic activity. Detailed amino acid substitutions resulted in peptides with superior efficacy. Specifically, substitutions with Isoleucine at positions twelve and eighteen along with the substitution of the Methionine at position ten with the non-natural amino acid d-Alanine led to an increase in potency by two orders of magnitude over the parent peptide. Several mimetic peptides in this series exhibit a significant improvement of activity over the parent peptide. This improved in vitro activity is expected to correlate with an increase in in vivo activity leading to effective peptides for anti-angiogenic therapy for different disease applications including cancer and age-related macular degeneration.
angiogenesis; peptidomimetics; endothelial cell; cancer; age-related macular degeneration; lymphangiogenesis
Vascular endothelial growth factor (VEGF) signal transduction is central to angiogenesis in development and in pathological conditions such as cancer, retinopathy and ischemic diseases. However, no detailed mass-action models of VEGF receptor signaling have been developed. We constructed and validated the first computational model of VEGFR2 trafficking and signaling, to study the opposing roles of Gab1 and Gab2 in regulation of Akt phosphorylation in VEGF-stimulated endothelial cells. Trafficking parameters were optimized against 5 previously published in vitro experiments, and the model was validated against six independent published datasets. The model showed agreement at several key nodes, involving scaffolding proteins Gab1, Gab2 and their complexes with Shp2. VEGFR2 recruitment of Gab1 is greater in magnitude, slower, and more sustained than that of Gab2. As Gab2 binds VEGFR2 complexes more transiently than Gab1, VEGFR2 complexes can recycle and continue to participate in other signaling pathways. Correspondingly, the simulation results show a log-linear relationship between a decrease in Akt phosphorylation and Gab1 knockdown while a linear relationship was observed between an increase in Akt phosphorylation and Gab2 knockdown. Global sensitivity analysis demonstrated the importance of initial-concentration ratios of antagonistic molecular species (Gab1/Gab2 and PI3K/Shp2) in determining Akt phosphorylation profiles. It also showed that kinetic parameters responsible for transient Gab2 binding affect the system at specific nodes. This model can be expanded to study multiple signaling contexts and receptor crosstalk and can form a basis for investigation of therapeutic approaches, such as tyrosine kinase inhibitors (TKIs), overexpression of key signaling proteins or knockdown experiments.
Vascular endothelial growth factor (VEGF) is a key mediator of angiogenesis, whose effect on cancer growth and development is well characterized. Alternative splicing of VEGF leads to several different isoforms, which are differentially expressed in various tumor types and have distinct functions in tumor blood vessel formation. Many cancer therapies aim to inhibit angiogenesis by targeting VEGF and preventing intracellular signaling leading to tumor vascularization; however, the effects of targeting specific VEGF isoforms have received little attention in the clinical setting. In this work, we investigate the effects of selectively targeting a single VEGF isoform, as compared with inhibiting all isoforms. We utilize a molecular-detailed whole-body compartment model of VEGF transport and kinetics in the presence of breast tumor. The model includes two major VEGF isoforms, VEGF121 and VEGF165, receptors VEGFR1 and VEGFR2, and co-receptors Neuropilin-1 and Neuropilin-2. We utilize the model to predict the concentrations of free VEGF, the number of VEGF/VEGFR2 complexes (considered to be pro-angiogenic), and the receptor occupancy profiles following inhibition of VEGF using isoform-specific anti-VEGF agents. We predict that targeting VEGF121 leads to a 54% and 84% reduction in free VEGF in tumors that secrete both VEGF isoforms or tumors that overexpress VEGF121, respectively. Additionally, 21% of the VEGFR2 molecules in the blood are ligated following inhibition of VEGF121, compared with 88% when both isoforms are targeted. Targeting VEGF121 reduces tumor free VEGF and is an effective treatment strategy. Our results provide a basis for clinical investigation of isoform-specific anti-VEGF agents.
Electronic supplementary material
The online version of this article (doi:10.1208/s12248-012-9363-4) contains supplementary material, which is available to authorized users.
angiogenesis; cancer drug target; computational model; pharmacokinetic model; systems biology
Angiogenesis, the formation of new blood vessels from preexisting blood vessels, is a process that supports tumor growth and metastatic dissemination. Lymphangiogenesis also facilitates metastasis by increasing dissemination through the lymphatic vessels (LVs). Even after treatment with antiangiogenic agents, breast cancer patients are vulnerable to LV-mediated metastasis. We report that a 14-amino acid peptide derived from transmembrane protein 45A shows multimodal inhibition of lymphangiogenesis and angiogenesis in breast cancer. The peptide blocks lymphangiogenic and angiogenic phenotypes of lymphatic and blood endothelial cells induced by tumor-conditioned media prepared from MDA-MB-231 breast cancer cells. The peptide delays growth of MDA-MB-231 tumor xenografts and normalizes tumor-conditioned lymph nodes (LNs). These studies demonstrate the antilymphangiogenic and antiangiogenic potential of the peptide against primary tumors and premetastatic, tumor-conditioned regional LNs. Mechanistically, the peptide blocks vascular endothelial growth factor receptors 2 and 3 (VEGFR2/3) and downstream proteins by binding to neuropilin 1/2 (NRP1/2) and inhibiting VEGFR2/3 and NRP1/2 complex formation in the presence of VEGFA/C.
Angiogenesis, the formation of new blood vessels from existing vasculature, is important in tumor growth and metastasis. A key regulator of angiogenesis is vascular endothelial growth factor (VEGF), which has been targeted in numerous anti-angiogenic therapies aimed at inhibiting tumor angiogenesis. Systems biology approaches, including computational modeling, are useful for understanding this complex biological process and can aid in the development of novel and effective therapeutics that target the VEGF family of proteins and receptors. We have developed a computational model of VEGF transport and kinetics in the tumor-bearing mouse, which includes three-compartments: normal tissue, blood, and tumor. The model simulates human tumor xenografts and includes human (VEGF121 and VEGF165) and mouse (VEGF120 and VEGF164) isoforms. The model incorporates molecular interactions between these VEGF isoforms and receptors (VEGFR1 and VEGFR2), as well as co-receptors (NRP1 and NRP2). We also include important soluble factors: soluble VEGFR1 (sFlt-1) and α-2-macroglobulin. The model accounts for transport via macromolecular transendothelial permeability, lymphatic flow, and plasma clearance. We have fit the model to available in vivo experimental data on the plasma concentration of free VEGF Trap and VEGF Trap bound to mouse and human VEGF in order to estimate the rates at which parenchymal cells (myocytes and tumor cells) and endothelial cells secrete VEGF. Interestingly, the predicted tumor VEGF secretion rates are significantly lower (0.007–0.023 molecules/cell/s, depending on the tumor microenvironment) than most reported in vitro measurements (0.03–2.65 molecules/cell/s). The optimized model is used to investigate the interstitial and plasma VEGF concentrations and the effect of the VEGF-neutralizing agent, VEGF Trap (aflibercept). This work complements experimental studies performed in mice and provides a framework with which to examine the effects of anti-VEGF agents, aiding in the optimization of such anti-angiogenic therapeutics as well as analysis of clinical data. The model predictions also have implications for biomarker discovery with anti-angiogenic therapies.
systems biology; mathematical model; computational model; angiogenesis; tumor xenograft model; anti-angiogenic therapy; cancer
Using modified oxygen needle microelectrodes and intravital videomicroscopy, measurements were made of tissue oxygen tension (PO2) profiles near cortical arterioles and transmural PO2 gradients in the pial arterioles of the rat. Under control conditions, the transmural PO2 gradient averaged 1.17 ± 0.06 mm Hg/μm (mean ± s.e., n = 40). Local arteriolar dilation resulted in a marked decrease in the transmural PO2 gradient to 0.68 ± 0.04 mm Hg/μm (P < 0.001, n = 38). The major finding of this study is a dependence of the transmural PO2 gradient on the vascular tone of the pial arterioles. Using a model of oxygen transport in an arteriole and experimental PO2 profiles, values of radial perivascular and intravascular O2 fluxes were estimated. Our theoretical estimates show that oxygen flux values at the outer surface of the arteriolar wall are approximately 10−5 mL O2/cm2 per sec, independent of the values of the arteriolar wall O2 consumption within a wide range of consumption values. This also means that PO2 transmural gradients for cerebral arterioles are within the limits of 1 to 2 mm Hg/μm. The data lead to the conclusion that O2 consumption of the arteriolar wall is within the range for the surrounding tissue and O2 consumption of the endothelial layer appears to have no substantial impact on the transmural PO2 gradient.
cortical microvessels; O2 transport model; oxygen microelectrodes; tissue PO2 profiles; transmural PO2 gradient
Blood is modeled as a suspension of red blood cells using the dissipative particle dynamics method. The red blood cell membrane is coarse-grained for efficient simulations of multiple cells, yet accurately describes its viscoelastic properties. Blood flow in microtubes ranging from 10 to 40 μm in diameter is simulated in three dimensions for values of hematocrit in the range of 0.15–0.45 and carefully compared with available experimental data. Velocity profiles for different hematocrit values show an increase in bluntness with an increase in hematocrit. Red blood cell center-of-mass distributions demonstrate cell migration away from the wall to the tube center. This results in the formation of a cell-free layer next to the tube wall corresponding to the experimentally observed Fahraeus and Fahraeus–Lindqvist effects. The predicted cell-free layer widths are in agreement with those found in in vitro experiments; the results are also in qualitative agreement with in vivo experiments. However, additional features have to be taken into account for simulating microvascular flow, e.g., the endothelial glycocalyx. The developed model is able to capture blood flow properties and provides a computational framework at the mesoscopic level for obtaining realistic predictions of blood flow in microcirculation under normal and pathological conditions.
apparent viscosity; red blood cell; blood flow resistance; dissipative particle dynamics
Angiogenesis is thoroughly balanced and regulated in health; however, it is dysregulated in many diseases including cancer, age-related macular degeneration, cardiovascular diseases such as coronary and peripheral artery diseases and stroke, abnormal embryonic development, and abnormal wound healing. In addition to angiogenesis, lymphangiogenesis is pivotal for maintaining the immune system, homeostasis of body fluids and lymphoid organs; dysregulated lymphangiogenesis may cause inflammatory diseases and lymph node mediated tumor metastasis. Anti-angiogenic or anti-lymphangiogenic small peptides may play an important role as therapeutic agents normalizing angiogenesis or lymphangiogenesis in disease conditions. Several novel endogenous peptides derived from proteins containing a conserved somatotropin domain have been previously identified with the help of our bioinformatics-based methodology. These somatotropin peptides were screened for inhibition of angiogenesis and lymphangiogenesis using in vitro proliferation, migration, adhesion and tube formation assays with blood and lymphatic endothelial cells. We found that the peptides have the potential for inhibiting both angiogenesis and lymphangiogenesis. Focusing the study on the inhibition of lymphangiogenesis, we found that a peptide derived from the somatotropin conserved domain of transmembrane protein 45A human was the most potent lymphangiogenesis inhibitor, blocking lymphatic endothelial cell migration, adhesion, and tube formation.
Lymphatic endothelial cell; Blood endothelial cell; Endogenous somatotropin peptides; Transmembrane protein 45A human
Breast cancer is one of the most commonly diagnosed malignancies in women. Despite the remarkable success of mammography screening and use of adjuvant systemic therapy, it is estimated that approximately 200,000 new diagnoses will be made this year and 40,000 deaths will occur due to this disease (American Cancer Society). Angiogenesis, the growth of vessels from pre-existing microvasculature, is an essential component of tumor progression and has emerged as a therapeutic modality for anti-angiogenic therapies in cancer.
Here we report in vitro and in vivo findings with a 20 amino acid peptide belonging to the collagen IV family, modified to facilitate possible translation to clinical applications. The two cysteines in its natural peptide progenitor were replaced by L-α-amino-n-butyric acid, a non-natural amino acid. The modified peptide was tested in vitro using endothelial cells and in vivo using mouse orthotopic breast cancer xenograft model with MDA-MB-231 human breast cancer cells. This modified peptide demonstrated no significant changes in activity from the parent peptide; however, because it lacks cysteines, it is more suitable for clinical translation. We also investigated its efficacy in combination with a commonly used chemotherapeutic agent paclitaxel; the inhibition of tumor growth by the peptide was similar to that of paclitaxel alone, but the combination did not exhibit any additional inhibition. We have performed further characterization of the mechanism of action (MOA) for this peptide to identify its target receptors, enhancing its translation potential as an antiangiogenic, non-vascular endothelial growth factor (VEGF) targeting agent for therapy in breast cancer.
angiogenesis; cancer therapy; drug development; integrins; triple-negative; xenograft
Angiogenesis is the growth of new blood vessels from existing vasculature. Excessive vascularization is associated with a number of diseases including cancer. Anti-angiogenic therapies have the potential to stunt cancer progression. Peptides derived from type IV collagen are potent inhibitors of angiogenesis. We wanted to gain a better understanding of collagen IV structure-activity relationships using a ligand-based approach. We developed novel peptide-specific QSAR models to study the activity of the peptides in endothelial cell proliferation, migration, and adhesion inhibition assays. We found that the models produced quantitatively accurate predictions of activity and provided insight into collagen IV derived peptide structure-activity relationships.
The process of oxygen delivery from capillary to muscle fiber is essential for a tissue with variable oxygen demand, such as skeletal muscle. Oxygen distribution in exercising skeletal muscle is regulated by convective oxygen transport in the blood vessels, oxygen diffusion and consumption in the tissue. Spatial heterogeneities in oxygen supply, such as microvascular architecture and hemodynamic variables, had been observed experimentally and their marked effects on oxygen exchange had been confirmed using mathematical models. In this study, we investigate the effects of heterogeneities in oxygen demand on tissue oxygenation distribution using a multiscale oxygen transport model. Muscles are composed of different ratios of the various fiber types. Each fiber type has characteristic values of several parameters, including fiber size, oxygen consumption, myoglobin concentration, and oxygen diffusivity. Using experimentally measured parameters for different fiber types and applying them to the rat extensor digitorum longus muscle, we evaluated the effects of heterogeneous fiber size and fiber type properties on the oxygen distribution profile. Our simulation results suggest a marked increase in spatial heterogeneity of oxygen due to fiber size distribution in a mixed muscle. Our simulations also suggest that the combined effects of fiber type properties, except size, do not contribute significantly to the tissue oxygen spatial heterogeneity. However, the incorporation of the difference in oxygen consumption rates of different fiber types alone causes higher oxygen heterogeneity compared to control cases with uniform fiber properties. In contrast, incorporating variation in other fiber type-specific properties, such as myoglobin concentration, causes little change in spatial tissue oxygenation profiles.
VEGFR surface localization plays a critical role in converting extracellular VEGF signaling towards angiogenic outcomes, and the quantitative characterization of these parameters is critical for advancing computational models; however the levels of these receptors on blood vessels is currently unknown. Therefore our aim is to quantitatively determine the VEGFR localization on endothelial cells from mouse hindlimb skeletal muscles. We contextualize this VEGFR quantification through comparison to VEGFR-levels on cells in vitro. Using quantitative fluorescence we measure and compare the levels of VEGFR1 and VEGFR2 on endothelial cells isolated from C57BL/6 and BALB/c gastrocnemius and tibialis anterior hindlimb muscles. Fluorescence measurements are calibrated using beads with known numbers of phycoerythrin molecules. The data show a 2-fold higher VEGFR1 surface localization relative to VEGFR2 with 2,000–3,700 VEGFR1/endothelial cell and 1,300–2,000 VEGFR2/endothelial cell. We determine that endothelial cells from the highly glycolytic muscle, tibialis anterior, contain 30% higher number of VEGFR1 surface receptors than gastrocnemius; BALB/c mice display ∼17% higher number of VEGFR1 than C57BL/6. When we compare these results to mouse fibroblasts in vitro, we observe high levels of VEGFR1 (35,800/cell) and very low levels of VEGFR2 (700/cell), while in human endothelial cells in vitro, we observe that the balance of VEGFRs is inverted, with higher levels VEGFR2 (5,800/cell) and lower levels of VEGFR1 (1,800/cell). Our studies also reveal significant cell-to-cell heterogeneity in receptor expression, and the quantification of these dissimilarities ex vivo for the first time provides insight into the balance of anti-angiogenic or modulatory (VEGFR1) and pro-angiogenic (VEGFR2) signaling.
Angiogenesis is the formation of neovasculature from preexisting microvessels. Several endogenous proteins regulate the balance of vessel formation and regression in the body including pigment epithelium-derived factor (PEDF), which has been shown to be antiangiogenic and to suppress tumor growth. Using sequence homology and bioinformatics, we previously identified several peptide sequences homologous to an active region of PEDF existing in multiple proteins in the human proteome. These short 11-mer peptides are found in a DEAH box helicase protein, CKIP-1 and caspase 10, and show similar activity in altering endothelial cell adhesion, migration and inducing apoptosis.We tested the peptide derived from DEAH box helicase protein in a triple-negative MDA-MB-231 breast orthotopic xenograft model in severe combined immunodeficient mice and show significant tumor suppression.
Vascular endothelial growth factor (VEGF) is one of the most potent cytokines targeted in anti-angiogenic therapies. Bevacizumab, a recombinant humanized monoclonal antibody to VEGF, is being used clinically in combination with chemotherapy for colorectal, non-small cell lung and breast cancers, and as a single agent for glioblastoma, and is being tested for other types of cancer in numerous clinical trials. It has been reported that the intravenous injection of bevacizumab leads to an increase of plasma VEGF concentration in cancer patients. The mechanism responsible for this counterintuitive increase has not been elucidated, although several hypotheses have been proposed. We use a multiscale systems biology approach to address this problem. We have constructed a whole-body pharmacokinetic model comprising three compartments: blood, normal tissue and tumor tissue. Molecular interactions between VEGF-A family members, their major receptors, the extracellular matrix, and an anti-VEGF ligand are considered for each compartment. Diffusible molecules extravasate, intravasate, are removed from the healthy tissue through the lymphatics, and are cleared from the blood. Our model reproduces the experimentally-observed increase of plasma VEGF following intravenous administration of bevacizumab, and predicts this increase to be a consequence of inter-compartmental exchange of VEGF, the anti-VEGF agent and the VEGF/anti-VEGF complex. Our results suggest that a fraction of the anti-VEGF drug extravasates, allowing the agent to bind the interstitial VEGF. When the complex intravasates (via a combination of lymphatic drainage and microvascular transport of macromolecules) and dissociates in the blood, VEGF is released and the VEGF concentration increases in the plasma. These results provide a new hypothesis on the kinetics of VEGF and on the VEGF distribution in the body caused by anti-angiogenic therapies, as well as their mechanisms of action and could help in designing anti-angiogenic therapies.
angiogenesis; anti-VEGF; bevacizumab; anti-angiogenic therapy; mathematical model
Vascular endothelial growth factor (VEGF) is a potent regulator of angiogenesis, and its role in cancer biology has been widely studied. Many cancer therapies target angiogenesis, with a focus being on VEGF-mediated signaling such as antibodies to VEGF. However, it is difficult to predict the effects of VEGF-neutralizing agents. We have developed a whole-body model of VEGF kinetics and transport under pathological conditions (in the presence of breast tumor). The model includes two major VEGF isoforms VEGF121 and VEGF165, receptors VEGFR1, VEGFR2 and co-receptors Neuropilin-1 and Neuropilin-2. We have added receptors on parenchymal cells (muscle fibers and tumor cells), and incorporated experimental data for the cell surface density of receptors on the endothelial cells, myocytes, and tumor cells. The model is applied to investigate the action of VEGF-neutralizing agents (called "anti-VEGF") in the treatment of cancer.
Through a sensitivity study, we examine how model parameters influence the level of free VEGF in the tumor, a measure of the response to VEGF-neutralizing drugs. We investigate the effects of systemic properties such as microvascular permeability and lymphatic flow, and of drug characteristics such as the clearance rate and binding affinity. We predict that increasing microvascular permeability in the tumor above 10-5 cm/s elicits the undesired effect of increasing tumor interstitial VEGF concentration beyond even the baseline level. We also examine the impact of the tumor microenvironment, including receptor expression and internalization, as well as VEGF secretion. We find that following anti-VEGF treatment, the concentration of free VEGF in the tumor can vary between 7 and 233 pM, with a dependence on both the density of VEGF receptors and co-receptors and the rate of neuropilin internalization on tumor cells. Finally, we predict that free VEGF in the tumor is reduced following anti-VEGF treatment when VEGF121 comprises at least 25% of the VEGF secreted by tumor cells.
This study explores the optimal drug characteristics required for an anti-VEGF agent to have a therapeutic effect and the tumor-specific properties that influence the response to therapy. Our model provides a framework for investigating the use of VEGF-neutralizing drugs for personalized medicine treatment strategies.
Vascular endothelial growth factor (VEGF) is a key regulator of angiogenesis – the growth of new microvessels from existing microvasculature. Angiogenesis is a complex process involving numerous molecular species, and to better understand it, a systems biology approach is necessary. In vivo preclinical experiments in the area of angiogenesis are typically performed in mouse models; this includes drug development targeting VEGF. Thus, to quantitatively interpret such experimental results, a computational model of VEGF distribution in the mouse can be beneficial. In this paper, we present an in silico model of VEGF distribution in mice, determine model parameters from existing experimental data, conduct sensitivity analysis, and test the validity of the model.
The multiscale model is comprised of two compartments: blood and tissue. The model accounts for interactions between two major VEGF isoforms (VEGF120 and VEGF164) and their endothelial cell receptors VEGFR-1, VEGFR-2, and co-receptor neuropilin-1. Neuropilin-1 is also expressed on the surface of parenchymal cells. The model includes transcapillary macromolecular permeability, lymphatic transport, and macromolecular plasma clearance. Simulations predict that the concentration of unbound VEGF in the tissue is approximately 50-fold greater than in the blood. These concentrations are highly dependent on the VEGF secretion rate. Parameter estimation was performed to fit the simulation results to available experimental data, and permitted the estimation of VEGF secretion rate in healthy tissue, which is difficult to measure experimentally. The model can provide quantitative interpretation of preclinical animal data and may be used in conjunction with experimental studies in the development of pro- and anti-angiogenic agents. The model approximates the normal tissue as skeletal muscle and includes endothelial cells to represent the vasculature. As the VEGF system becomes better characterized in other tissues and cell types, the model can be expanded to include additional compartments and vascular elements.
Vascular endothelial growth factor (VEGF) is a family of cytokines for which the dysregulation of expression is involved in many diseases; for some excess VEGF causes pathological hypervascularization, while for others VEGF-induced vascular remodeling may alleviate ischemia and/or hypoxia. Anti-angiogenic therapies attacking the VEGF pathway have begun to live up to their promise for treatment of certain cancers and of age-related macular degeneration. However, the corollary is not yet true: in coronary artery disease and peripheral artery disease, clinical trials of pro-angiogenic VEGF delivery have not, so far, proven successful. The VEGF and VEGF-receptor system is complex, with at least five ligand genes, some encoding multiple protein isoforms and five receptor genes. A systems biology approach to designing pro-angiogenic therapies, using a combination of quantitative experimental approaches and detailed computational models, is essential to deal with this complexity and to understand the effects of drugs targeting the system. This approach allows us to learn from unsuccessful clinical trials and to design and test novel single therapeutics or combinations of therapeutics. Among the parameters that can be varied in order to determine optimal strategy are dosage, timing of multiple doses, route of administration, and the molecular target.
Drug efficacy; Drug transport; Physiome; Protein interaction networks; Systems medicine
Gene therapy research has expanded from its original concept of replacing absent or defective DNA with functional DNA for transcription. Genetic material may be delivered via multiple vectors, including naked plasmid DNA, viruses and even cells with the goal of increasing gene expression; and the targeting of specific tissues or cell types is aimed at decreasing risks of systemic or side effects. As with the development of any drug, there is an amount of empiricism in the choice of gene target, route of administration, dosing and in particular the scaling-up from pre-clinical models to clinical trials. Systems Biology, whose arsenal includes high-throughput experimental and computational studies that account for the complexities of host-disease-therapy interactions, holds significant promise in aiding the development and optimization of gene therapies, including personalized therapies and the identification of biomarkers for success of these strategies. In this review we describe some of the obstacles and successes in gene therapy, using the specific example of growth factor gene delivery to promote angiogenesis and blood vessel remodeling in ischemic diseases; we also make references to anti-angiogenic gene therapy in cancer. The opportunities for Systems Biology and in silico modeling to improve on current outcomes are highlighted.
Mathematical model; computational model; bioinformatics; angiogenesis; gene delivery; coronary artery disease; peripheral artery disease
Angiogenesis is important for many physiological processes, diseases, and also regenerative medicine. Therapies that inhibit the vascular endothelial growth factor (VEGF) pathway have been used in the clinic for cancer and macular degeneration. In cancer applications, these treatments suffer from a “tumor escape phenomenon” where alternative pathways are upregulated and angiogenesis continues. The redundancy of angiogenesis regulation indicates the need for additional studies and new drug targets. We aimed to (i) identify novel and missing angiogenesis annotations and (ii) verify their significance to angiogenesis. To achieve these goals, we integrated the human interactome with known angiogenesis-annotated proteins to identify a set of 202 angiogenesis-associated proteins. Across endothelial cell lines, we found that a significant fraction of these proteins had highly perturbed gene expression during angiogenesis. After treatment with VEGF-A, we found increasing expression of HIF-1α, APP, HIV-1 tat interactive protein 2, and MEF2C, while endoglin, liprin β1 and HIF-2α had decreasing expression across three endothelial cell lines. The analysis showed differential regulation of HIF-1α and HIF-2α. The data also provided additional evidence for the role of endothelial cells in Alzheimer's disease.
Excessive vascularization is a hallmark of many diseases including cancer, rheumatoid arthritis, diabetic nephropathy, pathologic obesity, age-related macular degeneration, and asthma. Compounds that inhibit angiogenesis represent potential therapeutics for many diseases. Karagiannis and Popel (PNAS, 2008) used a bioinformatics approach to idenify more than 100 peptides with sequence homology to known angiogenesis inhibitors. The peptides could be grouped into families by the conserved domain of the proteins they were derived from. The families included type IV collagen fibrils, CXC chemokine ligands, and type I thrombospondin domain-containing proteins. The relationships between these families have received relatively little attention. To investigate these relationships, we approached the problem by placing the families of proteins in the context of the human interactome including >120,000 physical interactions among proteins, genes, and transcripts. We built on a graph theoretic approach to identify proteins that may represent conduits of crosstalk between protein families. We validated these findings by statistical analysis and analysis of a time series gene expression dataset taken during angiogenesis. We identified six proteins at the center of the angiogenesis-associated network including three syndecans, MMP9, CD44 and versican. These findings shed light on the complex signaling networks that govern angiogenesis phenomena.
Peptides have emerged as important therapeutics that are being rigorously tested in angiogenesis-dependent diseases due to their low toxicity and high specificity. Since the discovery of endogenous proteins and protein fragments that inhibit microvessel formation (thrombospondin, endostatin) several peptides have shown promise in pre-clinical and clinical studies for cancer. Peptides have been derived from thrombospondin, collagens, chemokines, coagulation cascade proteins, growth factors, and other classes of proteins and target different receptors. Here we survey recent developments for anti-angiogenic peptides with length not exceeding 50 amino acid residues that have shown activity in pre-clinical models of cancer or have been tested in clinical trials; some of the peptides have been modified and optimized, e.g., through L-to-D and non-natural amino acid substitutions. We highlight technological advances in peptide discovery and optimization including computational and bioinformatics tools and novel experimental techniques.
Angiogenesis; animal model; computational biology; inhibitor; in vitro model; peptidomimetics; tumor vasculature; tumor