Feilim Mac Gabhann

A systems biology perspective on sVEGFR1: its biological function, pathogenic role and therapeutic use

•  Introduction ‐  Angiogenesis in physiology and pathology ‐  Angiogenesis in current medicine ‐  VEGF ligand and receptor system: where does sVEGFR1 fit? ‐  Ligands: the human VEGF family ‐  Membrane‐bound signalling receptors: VEGFRs ‐  Non‐signalling co‐receptors and matrix proteins: HSPGs and NRPs ‐  Soluble receptors: sVEGFR1, sVEGFR2, sNRP1 •  Molecular biology of sVEGFR1 •  Physiological and pathophysiological roles of sVEGFR1 •  Molecular mechanism of sVEGFR1’s anti‐angiogenic potential •  sVEGFR1 as a clinical marker for disease •  Plasma VEGF and sVEGFR1: non‐uniform predictors of angiogenic status across all diseases •  Systems biology perspective: unifying interpretation of plasma angiogenic markers ‐  Baseline heterogeneity in clinical measurements of healthy VEGF and sVEGFR1 levels in plasma ‐  Effect of sVEGFR1 on VEGF bioavailability: VEGF‐sVEGFR1 complexes ‐  Compartmental analysis: biotransport and biodistribution ‐  Pathogenic phenomenon versus compensatory response •  Concluding remarks

Systems Biology of Vascular Endothelial Growth Factors

Several cytokine families have roles in the development, maintenance, and remodeling of the microcirculation. Of these, the vascular endothelial growth factor (VEGF) family is one of the best studied and one of the most complex. Five VEGF ligand genes and five cell‐surface receptor genes are known in the human, and each of these may be transcribed as multiple splice isoforms to generate an extensive family of proteins, many of which are subject to further proteolytic processing. Using the VEGF family as an example, we describe the current knowledge of growth‐factor expression, processing, and transport in vivo. Experimental studies and computational simulations are being used to measure and predict the activity of these molecules, and we describe avenues of research that seek to fill the remaining gaps in our understanding of VEGF family behavior.

A compartment model of VEGF distribution in blood, healthy and diseased tissues

BackgroundAngiogenesis is a process by which new capillaries are formed from pre-existing blood vessels in physiological (e.g., exercise, wound healing) or pathological (e.g., ischemic limb as in peripheral arterial disease, cancer) contexts. This neovascular mechanism is mediated by the vascular endothelial growth factor (VEGF) family of cytokines. Although VEGF is often targeted in anti-angiogenic therapies, there is little knowledge about how its concentration may vary between tissues and the vascular system. A compartment model is constructed to study the VEGF distribution in the tissue (including matrix-bound, cell surface receptor-bound and free VEGF isoforms) and in the blood. We analyze the sensitivity of this distribution to the secretion rate, clearance rate and vascular permeability of VEGF.ResultsWe find that, in a physiological context, VEGF concentration varies approximately linearly with the VEGF secretion rate. VEGF concentration in blood but not in tissue is dependent on the vascular permeability of healthy tissue. Model simulations suggest that relative VEGF increases are similar in blood and tissue during exercise and return to baseline within several hours. In a pathological context (tumor), we find that blood VEGF concentration is relatively insensitive to increased vascular permeability in tumors, to the secretion rate of VEGF by tumors and to the clearance. However, it is sensitive to the vascular permeability in the healthy tissue. Finally, the VEGF distribution profile in healthy tissue reveals that about half of the VEGF is complexed with the receptor tyrosine kinase VEGFR2 and the co-receptor Neuropilin-1. In diseased tissues, this binding can be reduced to 15% while VEGF bound to the extracellular matrix and basement membranes increases.ConclusionThe results are of importance for physiological conditions (e.g., exercise) and pathological conditions (e.g., peripheral arterial disease, coronary artery disease, cancer). This mathematical model can serve as a tool for understanding the VEGF distribution in physiological and pathological contexts as well as a foundation to investigate pro- or anti-angiogenic strategies.

Targeting Neuropilin-1 to Inhibit VEGF Signaling in Cancer: Comparison of Therapeutic Approaches

Angiogenesis (neovascularization) plays a crucial role in a variety of physiological and pathological conditions including cancer, cardiovascular disease, and wound healing. Vascular endothelial growth factor (VEGF) is a critical regulator of angiogenesis. Multiple VEGF receptors are expressed on endothelial cells, including signaling receptor tyrosine kinases (VEGFR1 and VEGFR2) and the nonsignaling co-receptor Neuropilin-1. Neuropilin-1 binds only the isoform of VEGF responsible for pathological angiogenesis (VEGF165), and is thus a potential target for inhibiting VEGF signaling. Using the first molecularly detailed computational model of VEGF and its receptors, we have shown previously that the VEGFR–Neuropilin interactions explain the observed differential effects of VEGF isoforms on VEGF signaling in vitro, and demonstrated potent VEGF inhibition by an antibody to Neuropilin-1 that does not block ligand binding but blocks subsequent receptor coupling. In the present study, we extend that computational model to simulation of in vivo VEGF transport and binding, and predict the in vivo efficacy of several Neuropilin-targeted therapies in inhibiting VEGF signaling: (a) blocking Neuropilin-1 expression; (b) blocking VEGF binding to Neuropilin-1; (c) blocking Neuropilin–VEGFR coupling. The model predicts that blockade of Neuropilin–VEGFR coupling is significantly more effective than other approaches in decreasing VEGF–VEGFR2 signaling. In addition, tumor types with different receptor expression levels respond differently to each of these treatments. In designing human therapeutics, the mechanism of attacking the target plays a significant role in the outcome: of the strategies tested here, drugs with similar properties to the Neuropilin-1 antibody are predicted to be most effective. The tumor type and the microenvironment of the target tissue are also significant in determining therapeutic efficacy of each of the treatments studied.

Systems biology of pro-angiogenic therapies targeting the VEGF system

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 for 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. Copyright

Increase of Plasma VEGF after Intravenous Administration of Bevacizumab Is Predicted by a Pharmacokinetic Model

Vascular endothelial growth factor (VEGF) is one of the most potent cytokines targeted in antiangiogenic 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 among 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.

Module-based multiscale simulation of angiogenesis in skeletal muscle

BackgroundMathematical 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.ResultsWe 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.ConclusionsThis 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.

A Compartment Model of VEGF Distribution in Humans in the Presence of Soluble VEGF Receptor-1 Acting as a Ligand Trap

Vascular endothelial growth factor (VEGF), through its activation of cell surface receptor tyrosine kinases including VEGFR1 and VEGFR2, is a vital regulator of stimulatory and inhibitory processes that keep angiogenesis--new capillary growth from existing microvasculature--at a dynamic balance in normal physiology. Soluble VEGF receptor-1 (sVEGFR1)--a naturally-occurring truncated version of VEGFR1 lacking the transmembrane and intracellular signaling domains--has been postulated to exert inhibitory effects on angiogenic signaling via two mechanisms: direct sequestration of angiogenic ligands such as VEGF; or dominant-negative heterodimerization with surface VEGFRs. In pre-clinical studies, sVEGFR1 gene and protein therapy have demonstrated efficacy in inhibiting tumor angiogenesis; while in clinical studies, sVEGFR1 has shown utility as a diagnostic or prognostic marker in a widening array of angiogenesis-dependent diseases. Here we developed a novel computational multi-tissue model for recapitulating the dynamic systemic distributions of VEGF and sVEGFR1. Model features included: physiologically-based multi-scale compartmentalization of the human body; inter-compartmental macromolecular biotransport processes (vascular permeability, lymphatic drainage); and molecularly-detailed binding interactions between the ligand isoforms VEGF(121) and VEGF(165), signaling receptors VEGFR1 and VEGFR2, non-signaling co-receptor neuropilin-1 (NRP1), as well as sVEGFR1. The model was parameterized to represent a healthy human subject, whereupon we investigated the effects of sVEGFR1 on the distribution and activation of VEGF ligands and receptors. We assessed the healthy baseline stability of circulating VEGF and sVEGFR1 levels in plasma, as well as their reliability in indicating tissue-level angiogenic signaling potential. Unexpectedly, simulated results showed that sVEGFR1 - acting as a diffusible VEGF sink alone, i.e., without sVEGFR1-VEGFR heterodimerization--did not significantly lower interstitial VEGF, nor inhibit signaling potential in tissues. Additionally, the sensitivity of plasma VEGF and sVEGFR1 to physiological fluctuations in transport rates may partially account for the heterogeneity in clinical measurements of these circulating angiogenic markers, potentially hindering their diagnostic reliability for diseases.

Collateral capillary arterialization following arteriolar ligation in murine skeletal muscle.

Microcirculation (2010) 17, 333–347. doi: 10.1111/j.1549‐8719.2010.00034.x

Quantifying the proteolytic release of extracellular matrix-sequestered VEGF with a computational model.

Background VEGF proteolysis by plasmin or matrix metalloproteinases (MMPs) is believed to play an important role in regulating vascular patterning in vivo by releasing VEGF from the extracellular matrix (ECM). However, a quantitative understanding of the kinetics of VEGF cleavage and the efficiency of cell-mediated VEGF release is currently lacking. To address these uncertainties, we develop a molecular-detailed quantitative model of VEGF proteolysis, used here in the context of an endothelial sprout. Methodology and Findings To study a cells ability to cleave VEGF, the model captures MMP secretion, VEGF-ECM binding, VEGF proteolysis from VEGF165 to VEGF114 (the expected MMP cleavage product of VEGF165) and VEGF receptor-mediated recapture. Using experimental data, we estimated the effective bimolecular rate constant of VEGF165 cleavage by plasmin to be 328 M−1s−1 at 25°C, which is relatively slow compared to typical MMP-ECM proteolysis reactions. While previous studies have implicated cellular proteolysis in growth factor processing, we show that single cells do not individually have the capacity to cleave VEGF to any appreciable extent (less than 0.1% conversion). In addition, we find that a tip cells receptor system will not efficiently recapture the cleaved VEGF due to an inability of cleaved VEGF to associate with Neuropilin-1. Conclusions Overall, VEGF165 cleavage in vivo is likely to be mediated by the combined effect of numerous cells, instead of behaving in a single-cell-directed, autocrine manner. We show that heparan sulfate proteoglycans (HSPGs) potentiate VEGF cleavage by increasing the VEGF clearance time in tissues. In addition, we find that the VEGF-HSPG complex is more sensitive to proteases than is soluble VEGF, which may imply its potential relevance in receptor signaling. Finally, according to our calculations, experimentally measured soluble protease levels are approximately two orders of magnitude lower than that needed to reconcile levels of VEGF cleavage seen in pathological situations.