Reach Us
+44-1625-708989

Medical, Pharma, Engineering, Science, Technology and Business

Department of Mechanical and Aerospace Engineering, University of California, Davis, California, USA

- *Corresponding Author:
- Aldredge CR

Department of Mechanical and Aerospace Engineering

University of California, Davis

California, USA

**Tel:**(530) 752-5016

**E-mail:**[email protected]

**Received Date**: October 25, 2016; **Accepted Date:** October 27, 2016; **Published Date**: November 11, 2016

**Citation: **Aldredge CR, Nargis NN (2016) Effects of Matrix Metalloproteinase on Tumour Growth and Morphology via Haptotaxis. J Bioengineer & Biomedical Sci 6: 207. doi:10.4172/2155-9538.1000207

**Copyright:** © 2016 Aldredge CR, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Bioengineering & Biomedical Science

A tumour invasion model has been developed to examine the influence of diffusible matrix metalloproteinases (MMPs) and the mediated proteolysis of non-diffusible, membrane type MMPs (MT-MMPs) on tumour growth and morphology via haptotaxis. Our results are the first to explore the influence of localized degradation of extracellular matrix (ECM) by MT-MMPs on the morphology of tumours growing in nutrient-rich and nutrient-poor microenvironments. Two-dimensional numerical simulation, using a level-set-based tumour host interface-capturing method, reveals that haptotaxis due to ECM degradation by MMP causes greater instability than with ECM degradation by MT-MMP in low-nutrient environments, even at low proliferation rates; whereas the resulting morphologies are similar for high apoptosis rates. Our simulation results show that while haptotaxis leads to completely different tumour growth rates and morphologies depending on proliferation and apoptosis rates in low-nutrient environment, there are no significant variations when we compare the haptotaxis due to ECM degradation by MMP and MT-MMP, except for low proliferation rates. Focusing on the differences between MMP and MT-MMP mediated effects; our study has important implications in MMP-target validation and MMP-inhibitor-drug development for anti-cancer clinical trials.

**Metalloproteinase**; Tumour; **Morphology**; Proteolysis

Tumour invasion is driven by proliferation and migration into the surrounding tissue. Invasive cells have different characteristics than those of normal cells in that they are less adhesive and more mobile, **mitotic**, and metabolically active than normal cells [1-4]. Among these characteristics, cell motility is a crucial aspect of tumour invasion as cells with motile capabilities can access new nutrient sources and infiltrate the surrounding tissue which is essential for metastasis. Cell motion is commonly described by both random motion of the cells and directed motion stimulated by chemical (nutrient or other growth factors) or structural (ECM) gradients. The gradients are sensed by the migrating cells and influence the direction of movement. Chemo taxis is the **phenomenon** by which the movement of cells is directed in response to a soluble extracellular chemical gradient; i.e., of nutrient or oxygen or other growth factors. Cells produce matrix degrading enzymes (MDE); e.g., matrix **metalloproteinase** (MMPs) and urokinase plasminogen activators that degrade the ECM locally to make room for migration as cells cannot move into regions of the tissue which are too dense. The directed movement of cells in response to gradients of fixed or bound (non-diffusible) chemicals, such as the ECM, is known as haptotaxis.

Most MMPs are secreted as soluble enzymes but six of them are membrane-type MMPs (MT-MMPs) that are associated with the **cell membrane** [5]. Soluble MMPs, such as MMP-2, are secreted as latent pro-enzymes (proMMPs) and the activation of proMMPs often involves MT-MMPs and tissue inhibitors of matrix metalloproteinases (TIMPs). In addition to activating other MMPs, MT-MMPs degrade a number of ECM macromolecules [6-8]. MT1-MMP (MMP14) has emerged as an important collagenase that cancer cells use to degrade and invade in a collagen-rich environment [9,10]. While poorly invasive breast adenocarcinoma cells undergo apoptosis when confronted with a collagen-rich environment, the production of MT1-MMP endows these cells with the capacity to escape from collagen-induced apoptosis [11]. The localized activity of MT1- MMP helps tumour cells to overcome higher collagen density found in some peritumoural stroma [12]. Sabeh, et al. [13] have demonstrated that when cancer cells are faced with structural barriers created in reconstituted gels by covalently cross-linked fibrils of type I collagen, or that exist in the stromal environment of the mammary gland, invasion is dependent on MT1- MMP-mediated **proteolysis**.

Experimental evidence suggests that hypoxia (shortage of oxygen/ nutrient) can contribute importantly to instability and increase the invasive behaviour of a tumour. Instability (i.e. non-circular shape) provides a mechanism for tumour invasion that does not require development of a neo-vasculature to supply essential nutrients. Poplawski et al. and Poplawski et al. [14,15] demonstrated that the development of instability depends primarily on the diffusional limitation parameter (ratio of tumour growth rate to diffusion rate of nutrients) while the morphological details depend on cell-cell adhesion. Their works show that the lack of competition for nutrients (i.e. high nutrient environment) promotes spherical, non-invasive tumours and low concentrations of nutrients in the environment (which cause tumour-cell competition), or cells with a very high substrateconsumption rate generate a fingering instability and irregular, invasive tumours. Their results agree with the in vitro and in vivo experiments conducted in and with other tumour model predictions showing invasive fingered morphology if the nutrient supply is too small. Similar behaviour has been reported when haptotaxis is considered along with the impact of nutrient availability on tumour morphology [16-27].

Apoptosis (programmed cell death) serves as a natural barrier to cancer development and therefore, the ability of cancer cells to evade **apoptosis** is an essential hallmark of cancer [2,28,29]. **Tumour** cells evolve a variety of strategies to limit or circumvent apoptosis and often show weaker response to hypoxia-induced apoptosis [30,31]. Research has revealed how apoptosis is attenuated in some tumours that succeed in progressing to high-grade malignancy states and resistance to therapy [28,29]. Anderson et al. [32] demonstrated that a high necrosis threshold (below which cell death occurs) leads to more fingers in low nutrient environment. Previous studies (reviewed in) also suggests that defects along apoptotic pathways play a pivotal role in cancer development and many novel treatment strategies targeting apoptosis may be used [28].

Deakin et al. [12] presented a continuum model of cancer cell invasion to examine the differences between the invasion mediated by the localized degradation by MT1-MMP and the more extensive degradation by the soluble and diffusible MMP-2 or MMP-9. The model included cancer cell proliferation (with competition for space with the ECM), random diffusion, a hypotactic response to ECM gradients and the production (by cancer cells) and the natural decay of the MDEs. The results showed that invasion through soluble MMPs produces a greater invasion depth and greater amount of degraded ECM than those through MT-MMPs. Watanabe et al. [33] studied the dynamic nature of MT1-MMP at invadopodia and demonstrated the importance of its transient peak (and subsequent degradation of the ECM) in the activity of MT1-MMP followed by steady state activity. This transient activity was due to the inhibition by TIMP-2, and the steady state activity of MT1-MMP. When the transient activity was forcibly suppressed in computer simulations, the ECM degradation was heavily suppressed, indicating the essential role of this transient activity of MT1-MMP in the ECM degradation.

Other mathematical models that have examined the activation of MMP-2 by MT1-MMP include [34-36]. The experiments and simulations conducted by Hoshino et al. [36] established the role of the rapid turnover of MT1-MMP in ECM degradation at invadopodia (protrusions in the cell membrane extended into the ECM). The degree of the reduction in ECM degradation depended on the degree of the reduction in the MT1-MMP turnover. Furthermore, their simulations suggested synergetic contributions of proteolysis activity and the MT1- MMP turnover to ECM degradation because there was a nonlinear and significant reduction in ECM degradation if both factors were reduced simultaneously. But none of these works have considered the effect of haptotaxis due to ECM gradients created by MT-MMP on tumour growth and morphology with the variations of nutrient availability in the tumour microenvironment. In this study we investigate and compare the influences of haptotaxis due to MT-MMP and MMPmediated effects for different proliferation and apoptosis rates in low and high nutrient environments.

Primarily, tumour growth depends on the balance between expansive forces (caused by cell proliferation) and cell-cell adhesion forces (which maintain the tumour’s compactness) [37]. Lownutrient tumour growth and morphology are found to be strongly affected by adhesion and haptotaxis. Increasing cell-cell and cell- ECM adhesion of the microenvironment can help to limit the rate of tumour fragmentation and the extent of tissue invasion [24]. There is a combination of the strength of adhesion and haptotaxis in which finger like shapes are observed [38]. Change in the adhesion properties of cells and the disruption of the hypotactic movement result in less extensively disperses invasion fronts, as confirmed in vitro studies [39,40]. Cell proliferation, which depends on the availability of nutrient in the tumour microenvironment, is the predominate factor for tumour growth. Therefore, the effect of MT-MMP on tumour growth and morphology in nutrient-poor and nutrient-rich tissues with varying extents of adhesion must be understood for the design of drugs targeting haptotaxis for controlling malignant tumour growth and invasion.

In the present work, a continuum model of tumour cell invasion is proposed, based on previous continuum models of solid tumour growth to investigate the effect of haptotaxis on tumour growth in nutrient-poor and nutrient-rich tissues [41,42]. Our model accounts for tumour cell migration due to proliferation, cell-cell and cell-ECM adhesion and increasing gradients of nutrient and ECM. We capture the response of adjacent healthy tissue to tumour growth by solving for velocity and pressure both inside and outside the tumour. Our model enables qualitative predictions regarding the invasive ability of tumour cells. Using two-dimensional numerical simulations, we will present the results of a parameter study of quantitative aspects of the tumour progression, such as the size and shape of tumour with a range of biophysical and taxis parameters, and investigate the causal link between these parameters on tumour growth and morphology. Our model focuses solely on the interactions between the tumour cells and the surrounding tissue, so we have not included angiogenesis. The objectives of the present work are to differentiate the effect of haptotaxis for localized and extensive ECM degradation by MT-MMP and MMP, respectively, on tumour growth and morphology in high and lownutrient environments with varying extents of adhesion.

A one millimetre radius spheroidal tumour contains roughly one million cells. In large scale systems where the cancer cell population is on the order of one million or more, continuum models are more suitable than discrete cell models. In a continuum model, the cell species velocity is obtained from the inertia less momentum conservation equation based on Darcy’s law, representing the instantaneous equilibrium among forces associated with pressure, cell adhesion, elastic forces, forces exchanged with the ECM leading to haptotaxis and chemo taxis due to gradients of nutrient and other growth factors, and other mechanical effects. The continuum models incorporate the adhesion among cancer cells by a surface tension at the tumour surface [20-22,24,37,41,43-50]. Surface tension also models cell-ECM adhesion and the presence of a membrane, or capsule, formed of ECM macromolecules that may encapsulate tumours in vitro and in vivo. This representation of adhesion is relatively simple and indirect, involving no explicit modelling of cell-cell or cell-ECM contact. In these approaches cell to cell and cell to ECM adhesion are modelled as tumour shape stabilizing mechanical forces. This choice is supported by the experimentally observed presence of surface tension at tissue boundaries [51]. The tumour and healthy tissues are modelled as constant-density fluids. A sharp interface separates the tumour and healthy tissue. Cell-cell and cell-ECM adhesion is modelled using a surface-tension force at the tumour interface. The cell velocity is determined by Darcy's law, and growth occurs due to pressure gradients induced by mitosis (cell proliferation) and haptotaxis. Here we model non-necrotic tumours while accounting for apoptosis (programmed cell death).

We consider a rectangular region as shown in **Figure 1** Containing the tumour mass Ω_{T} with boundary Σ_{T} and the non-cancerous tissue Ω_{H} consisting of the ECM, healthy cells, and any other material immediately surrounding the tumour. We choose to focus on the key variables involved in tumour cell growth and invasion; ρ nutrient concentration N, MDE concentration ρ_{M}, ECM density ρ_{E}, cellular velocity field u and solid pressure p. We derive a system of coupled partial differential equations to model tumour growth in the surrounding tissue and present numerical solutions in two spatial dimensions describing the macroscopic dynamics of invasion in different tumour microenvironments due to haptotaxis. The surrounding tissue contains a variety of host cells like fibroblasts, macrophages and blood vessels, all of which have been shown to be important factors in tumour growth [52]. However, these factors are neglected in order to focus on the impact of nutrient concentration and the haptotaxis on tumour growth.

**Nutrient transport**

The net effect of nutrients and growth-promoting factors are characterized by a single vital nutrient (oxygen) that is required for cell survival and mitosis. The concentration of nutrient ρ_{N}, which diffuses through the ECM and is consumed by the tumour cells, satisfies the following reaction-diffusion equation.

Here, D_{N} is the nutrient diffusivity, δ_{N} is the rate of nutrient uptake by the tumour cells (assumed constant) and I_{ΩT} is a characteristic function equal to 1(0) inside (outside) of the tumour region Ω_{T} . Because nutrient diffusion and uptake occur much more quickly than tumour growth, the quasi-steady assumption applies. We assume that the nutrient and nutrient flux are continuous across the tumour boundary. The two cases that we consider are (i) a low-nutrient environment (e.g., no pre-existing blood vasculature in the computation domain) and (ii) a high-nutrient environment (e.g., the presence of a pre-existing network of nutrient-supplying blood vessels in the domain). For case (ii), nutrient delivery by the blood vasculature and uptake by noncancerous cells are assumed to be in balance outside of Ω_{H} so that the blood vasculature sufficiently delivers nutrient so that ρ_{N} is a constant (ρ_{N∞}) at the boundary of the computational domain. So for case (i), we solve the Eq. (1) in the entire domain and for case (ii) we solve it only in the tumour region. We assume that tumour cells uptake nutrient at a greater rate than noncancerous cells and that nutrient uptake is negligible in Ω_{H} [53-55]. Furthermore, we assume that there is no nutrient decay in the entire domain.

**The MDE and ECM density**

MMPs are produced (or activated) by the tumour cells at a constant rate λ_{M} and diffuse freely throughout the tissue. Removal of these enzymes takes place due to natural decay and by deactivation of the enzymes at a constant rate λ_{DM}. The equation governing the evolution of MMP concentration is therefore given by,

Because the diffusion coefficient for MMP (DM) is much smaller than that for oxygen, the full time-dependent diffusion equation is used [30]. As MT-MMP is bound to tumour cells, it moves with the tumour cells with velocity u. So following, (but not considering random cell movement and ECM as physical growth restraint as modelled in) the MT-MMP equation reads [12].

MT-MMPs are produced by the tumour cells at a constant rate λMT and decreased due to natural decay and by deactivation of the enzymes at a constant rate λDMT. The reference value ρ_{M0} /ρ_{MT0} is the maximum possible MMP concentration, assuming that tumour cells can sense the amount of MDE produced and the production of MDE is a function of ρ_{M0} /ρ_{MT0} [42]. The ECM is a complex mixture of macromolecules, but for the sake of simplicity we represent it with a single concentration. Similarly, we characterize MMP and MT-MMP with a single concentration. Although MMPs have multiple roles either promoting or inhibiting tumour progression or metastasis, we only focuses on pro-invasive roles in this paper [56]. In our model, we do not consider MMP-activation process and the effects of TIMP. As cell- ECM interaction is very complex, it is beneficial to investigate each of these effects separately.

Degradation of the ECM by MDEs increases the ability of the tumour to push into the surrounding tissue, both by reducing the mechanical rigidity of the surrounding tissue and by creating extra space for the growing tumour [57]. The ECM is considered as non-motile matter and changes in its distribution are due to its local degradation by MDEs upon contact at a rate λ_{DE}. The ECM density ρ_{E} is given by,

For degradation by MT-MMP, ρ_{M} and ρ_{M0} will be replaced by ρ_{MT} and ρ_{MT0}, respectively. In the far-field (at the boundary of the computational domain), we use the Dirichlet boundary conditions ρ_{M}=0, ρ_{MT}=0, and ρ_{E}=ρ_{E∞} (ECM density at far field). The ECM may deform and gets remodelled in response to pressure and insoluble ECM macro-molecules released by the cells, but we have not considered this.

**Cellular velocity field**

Tumour cells are assumed to move with a single cellular velocity, and cellular motion through the ECM is assumed to be similar to fluid flow in a porous medium. The proliferating tumour cells generate an internal (oncotic) mechanical pressure that also exerts force on the surrounding noncancerous tissue in Ω_{H}. We assume that the motion of the cancer cells is governed by the proliferation pressure and response to nutrient and ECM gradients by chemo taxis and haptotaxis, respectively. Cells thus move with a mass-averaged velocity arising from a generalized Darcy-type constitutive law for velocity from excess forces due to chemo taxis and haptotaxis.

Here p is the proliferation pressure and μ, χ_{N} and χ_{E} are the mobility, chemotaxis and haptotaxis coefficients, respectively. The cellular mobility μ characterizes the overall ability of the tumour tissue to respond to pressure gradients as well as the permeability of the tumour cells. The growth rate of the tumour is characterized by its rate of volume change and is accounted for with the continuity equation. Cell birth and death are in balance in ΩH and so there is no change in the volume of that region. The diffusion of tumour cells is an order of magnitude smaller than advection Zheng et al. [41], so we neglect diffusion of tumour cells. Although the ECM functions directly as a growth restraint for tumour cells, we have not included this role in our model in order to focus on the effect of directional migration (via haptotaxis) on tumour growth and morphology. The haptotaxis coefficients are assumed to be non-negativeand homogeneous in Ω_{T}, but they can be functions of time, space and ECM density. Experiment Zaman et al. [58] has shown that cancer cell motility depends not only on ECM gradients (haptotaxis) but also on ECM density (haptokinesis). However, since our aim here is to focus on gradient driven migration, we effectively choose to ignore the dependence of the cancer cell motility on the ECM density in the absence of ECM gradients. We have also neglected cell movement along ECM fibers (collagen), which is known as contact guidance [59,60].

Using a mass-conservation argument, we get the following equation for the tumour cell density ρ_{T}:

Where λ_{p} and λ_{A} are the rates of volume gain and loss due to cellular mitosis and apoptosis respectively. The tumour mass increases through cell proliferation and decreases through apoptosis for a non-necrotic tumour. As cells prefer to stay bound to each other, we assume that the total tumour cell density ρ_{T} is constant, defining what is known as the close-packing density. Due to the sharp interface between the tumour and host tissue, we assume that ρ_{T} = 0 Ω_{H} in and ρ_{T} = 1 in Ω_{T} so that equation (5) reduces to

We also assume that in the proliferating region, cell-mitosis is proportional to the amount of nutrient present and that apoptosis may occur. There is no proliferation in the host microenvironment. We assume that the normal velocity is continuous across the tumour boundary.

**Mechanical pressure**

As the tumour phase, encompassing cells, interstitial fluid, and extracellular matrix, is treated as porous media, no distinction between interstitial fluid hydrostatic pressure and mechanical pressure due to cell-cell and cell ECM interactions is made. When there is no taxis term, continuity of the normal velocity across the tumour boundary dictates that there is no jump in the normal derivative of p across Σ_{T}, i.e., (∂p/∂n) = 0. Assuming a uniform cell-cell and cell-matrix adhesion throughout the tumour, the adhesive pressure can be incorporated as a jump boundary condition at the tumour-host interface Σ_{T} of the Laplace- Young surface-tension type; i.e., (p)= p_{outer} -p_{inner} = κ_{Γ}, where κ is the mean curvature of the interface and Γ is a constant adhesion parameter analogous to the surface-tension coefficient. Thus, the adhesive forces are modelled by a curvature boundary condition on the interface Σ_{T}. The non-cancerous tissue in Ω_{H} is assumed to be close enough to the tumour to be affected by the pressure changes within the computational domain. The pressure is assumed to satisfy the homogeneous boundary condition, p=0. Experiment has found that the activity of enzymes is cantered almost exclusively in a narrow region adjacent to the tumour surface where it is balanced by the opposing anti-proteases secreted by nearby healthy tissue [4,61,62]. Therefore, we take taxis coefficients as piecewise constant terms i.e., in the tumour region they are positive numbers and outside the tumour, they are set to zero. Thus, pressure satisfies the following normal-gradient jump condition.

Defining the net pressure responsible for tumour growth as pnet, such that p_{net}=p-(χ_{N} ρ_{N}+χ_{E}ρ_{E})/μ Eq. (5) becomes,

Combining Eqs. (7) and (8), we obtain

With the following jump conditions on pressure and its normal gradient across the tumour boundary.

The cell velocity u can be determined from Eq. (9) once the pressure is known.

**Non-dimensionalization**

As the nutrient-concentration equation reveals intrinsic diffusionlength and relaxation-time scales, we non-dimensionalize space and time with L and the cell proliferation (mitosis) time scale τ, where

The length scale L corresponds to the diffusion-penetration length of oxygen in tissue and therefore characterizes the maximum invasion distance at the early stage of tumour invasion. The characteristic tumour pressure p_{T} =L^{2}λ_{p,max/} is that which results in maximum cell proliferation with a cell speed of unity. The field variables are made dimensionless as specified in **Table 1**, while non-dimensional parameters are defined in **(Tables 2-4)**. As reference densities, we use the far-field nutrient concentration ρ_{N∞}, the maximum sustainable MMP and MT-MMP concentration ρ_{M0} and ρ_{MT0}, respectively and the far-field ECM density ρ_{E∞}. The dimensionless governing equations in terms of the nondimensional variables and parameters are given below. We use the same symbols for the non-dimensional variables as those used for the respective dimensional variables defined above.

Dimensionless Variable | Symbol | Units |
---|---|---|

Nutrient concentration | ρN |
ρN∞ |

ECM density | ρE |
ρE∞ |

MDE concentration | ρM |
ρMO |

MT-MMP concentration | ρMT |
rMTO |

Cell velocity | U |
L /λp,max |

Pressure | P net |
ptot |

**Table 1:** Dimensionless-variable summary.

Parameter | Mathematical Expression | Significance | Value |
---|---|---|---|

MDE diffusion coefficient D_{M} |
DM/(L^{2}λ_{P,max}) |
MDE diffusion coefficient relative to characteristic diffusivity L2/τ [43] | 1 |

MDE decay rate λ_{DM} |
λ_{DM}/λ_{P,max} |
MDE decay rate relative to Maximum tumour-cell proliferation rate [43] | 10 |

MDE tumour-cell production rate λ_{M} |
λMr_{T}/(λ_{P},maxr_{MO}) |
MDE production rate relative to maximum tumour-cell proliferation rate [43] | 100 |

ECM degradation rate λ_{DE} |
λDE/λP,max | ECM degradation rate (by MDE) relative to maximum tumour-cell Proliferation rate [43 | 0.01 |

Nutrient chemotaxis coefficient χ_{N} |
χN ρN∞/ λp,max L2 | The extent of chemotaxis of nutrient gradients [73] | 0 |

**Table 2:** Parameter values.

Parameter | Mathematical Expression | Significance | Value |
---|---|---|---|

Surface tension coefficient Γ | ΓLp_{T} |
A measure of cell-cell and cell-ECM adhesion forces [20] | 0.05 0.01 (figure 6, 7 and 10 only) |

Rate of mitosis λ_{P} |
λP /λP,max | Mitosis rate relative to maximum tumour-cell proliferation rate | 1 |

Rate of apoptosis, λ_{A} |
λA/λP,max | Apoptosis rate relative to maximum tumour-cell proliferation rate | 0.1 (low nutrient, Sect. 4.1.1) 0.5 (high nutrient, Sect. 4.1.2) |

Haptotaxis coefficient χ_{E} |
χ_{E} ρ_{E∞}/λp,max L2 |
The extent of haptotaxis of intact-ECM gradients [43] | 0.2, 1, 5 |

**Table 3:** Parameter values.

Parameter | Mathematical Expression | Significance | Value |
---|---|---|---|

Surface tension coefficient Γ | Γ Lp_{T} |
A measure of cell-cell and cell-ECM adhesion forces [20] | 0.05 |

Rate of mitosis λ_{P} |
λP /λP,max | Mitosis rate relative to maximum tumour-cell proliferation rate | 0.1, 0.5, 1 |

Rate of apoptosis, λ_{A} |
λA/λP,max | Apoptosis rate relative to maximum tumour-cell proliferation rate | 0, 0.1 |

Haptotaxis coefficient χ_{E} |
χE ρE∞/λp,max L2 | The extent of haptotaxis of intact-ECM gradients [43] | 5, 10 |

**Table 4:** Non-dimensional parameters.

The following non-dimensional initial and boundary conditions are applied.

The set of governing equations, Eqs (13 and 14), are solved throughout the entire computational domain, on a uniform 2D cartesian mesh, except for the Poisson equation for ρ_{N}, which for case (ii) is solved for the tumour region Ω_{T} only. Because we anticipate frequent and complex morphological tumour-surface changes we use the levelset method, based on a reformulation of a model proposed in [63]. Let φ be an auxiliary level-set function whose zero level set denotes the boundary Ω_{T} of an avascular tumour growing into a surrounding, non cancerous tissue, satisfying φ < 0inside Ω_{T}, φ> 0 outside on the tumour boundary Σ_{T} . We update the position of the interface by solving the PDE,

We solve the Eq. (16) with the semi Lagrangian scheme developed by Aldredge [63] using upwind transient interpolation modelling [64,65]. The scheme provides third-order spatial accuracy and shape preservation. We use the immersed interface method [66] to solve for pressure, which is discontinuous across the tumour-host interface, in accordance with the boundary conditions specified in Eq (15). This requires calculation of normal vectors and curvature at the interface using the level-set function. For accurate normal φat and every curvature time step calculation by we reinitialize by solving the following equation [67,68].

Where φ^{0} is the level set function prior to re-initialization and τ is the re-initialization time, reset to and 0 at the beginning of each reinitialization process. We discretize the temporal derivative using a third-order total-variation-diminishing Runge-Kutta method (TVDRK) and approximate sign(φ^{0}) |grad(φ)| with the fifth-order WENO scheme. However when two interfaces are in close contact, level set functions develop singularities that yield inaccurate normal vectors and curvatures when applying traditional discretization methods. In these cases, we construct the level-set function on a local sub grid where we can accurately calculate the normal vectors and curvature using standard discretization methods [69,70]. As the velocity field near the interface is very sensitive to variations in the curvature, we use a gaussian smoothing filter to damp high frequency perturbations in the speed and interface position, in accordance with [71].

We employ a computational Ω of domain dimension [0,20] × [0,20] and grid size 256 × 256. The initial tumour shape is slightly elliptical and cantered at the domain and is represented by the following equation.

To minimize the influences of the computational-domain boundary, we stop the simulation when the tumour boundary is within four distance units of the boundary of the computational domain. We assume that the ECM is initially not degraded (i.e., ECM=1 and MDE=0 in the entire domain). To assess tumour growth, we calculate the tumour area and surface perimeter at regular intervals. Moreover, following [24], we define a shape parameter S = (Perimeter)^{2}/ π*Ar(4 ea) and length scale LS=2*Area/Perimeter to characterize the morphological characteristics of the tumour. The shape parameter S is a measure of the non-circularity the initial tumour shape; when S=1, the initial tumour shape is a circle. The length scale LS is a measure of the smallest dimension of the tumour. For example, for a rectangular tumour with width W and length L, LS = LW/ (L + W) and LS ≈ W if W<L.

**Effect of haptotaxis due to ECM gradients from MT-MMP mediated proteolysis**

In this section we examine the influences of haptotaxis due to localized ECM degradation by MT-MMP with variations in nutrient availability in the tumour microenvironment subject to the variations in surface tension, proliferation and apoptosis rates. We also compare these effects with MMP-mediated effects.

Effect of MT-MMP in a nutrient-poor microenvironment: Here, we investigate the influence MT-MMP mediated proteolysis on tumour growth and morphology in a nutrient-poor environment, which has not yet been reported in the literature. **Figure 2** shows the evolution of the tumour morphology with time for haptotaxis coefficient χE=1 for both MMP and MT-MMP mediated proteolysis. Relevant parameter values are listed in **(Tables 2 and 3)**. During early times the growth is similar for both (t=48) and forms a neck by creating bulbs on each side of the y-axis that grow in time. The influence of MMP-mediated proteolysis on the growth rate becomes pronounced, as corrugation of the tumour surface continues and bulges are formed on the tumour surface symmetric about both axes (at t=80), which grow in time while MTMMP creates no bulge along x-axis. The contours of MDEs and ECM are shown in **Figure 3** where we can see the diffusible and non-diffusible MDEs and the resulting ECM degradations. The ECM gradients near the boundary of the tumour are responsible for these differences in tumour morphologies via haptotaxis, as shown in **(Figure 4)**. We notice that gradients are largest on the opposing sides of the bulge (regions A and C in **Figure 4a** for MMP, indicating further hapotaxis induced growth of the bulge which are smooth for MT-MMP despite having larger gradients than that for MMP **(Figure 4c)**. Similarly the growth of top bulges for both MMP and MT-MMP can be explained by **Figure 4b and 4d**.

Qualitatively, similar results are found for five-fold increase of the haptotaxis coefficient for MT-MMP whereas the tumour develops four large, symmetric bulbs with very narrow necks, which continue to spread in the computational domain and eventually spilt into pieces for MMP. The effect of an increase in the haptotaxis coefficient on tumour growth and morphology is shown in **(Figure 5)**. As seen in the plot of tumour area versus time **(Figure 5a)**, the higher the rate of haptotaxis, the greater variations for MMP and MT-MMP. This trend is also reflected in the plot of tumour perimeters with time, presented in **(Figure 5b)**. Growth is faster χ for MMP and more distinguishable as E increases.

It is not known whether decrease in adhesion can cause similar growth and morphologies due to MMP and MT-MMP mediated proteolysis in nutrient-poor microenvironments. Therefore, we decrease the surface tension by 10-fold for the simulation results shown in **(Figure 6)**. We get more fingers for MMP compared to that for MT-MMP (c.f., **Figure 6**). The difference between the growth rates and morphological details also increase as surface tension decrease as shown in **Figure 7** comparing with **Figure 5**. Therefore, the results of our simulations show that MT-MMP mediated proteolysis has less effect on tumour growth and morphology compared to that of MMP mediated proteolysis in a hypoxic microenvironment. For similar cases, the tumour develops more fingers for MMP than MT-MMP. So MMP mediated proteolysis leads to more tumour invasion.

**Effect of MT-MMP in a nutrient-rich microenvironment:** In this section, we investigate the influence MT-MMP mediated proteolysis on tumour growth and morphology in a nutrient-rich environment, which has not yet been addressed in the literature. Relevant parameter values are given in **(Tables 2 and 3)**. The evolution of the tumour morphology with time for both proteases is displayed in **Figure 8** for χ_{E} =1. The tumour develops an H-shaped structure in this case, also observed in Zheng, et al. [41] for no taxis movement although the tumour grows faster due to the combined effect of cell proliferation and haptotaxis. The tumour morphologies are very similar for both MMP and MTMMP even after an increase in the haptotaxis coefficient by a factor of 5 (not shown here). The only difference being that the buds tend to bend inward more for MT-MMP. Although the buds approach each other, they never reconnect. Rather, they continue to grow toward the boundary of the computational domain.

**Figure 9** shows the effect of haptotaxis on tumour growth and morphology for both proteases in high nutrient environment. We note that the growth is faster for MT-MMP for χ_{E} =0.2 and χ_{E} =1 unlike for the case of tumour growth in a nutrient-poor environment. However, for E =5, both proteases gives comparable rate of growth. This indicates lower haptotactic sensitivity for nutrient-rich tumour growth. Still, qualitatively, the tumour morphologies for both proteases in nutrient-rich environment are very similar to that found for a tumour with no taxis even with increase in haptotaxis coefficients.

Now, we investigate the effect of decrease in cell adhesion on tumour growth dynamics and morphologies in nutrient-rich environment. If we decrease the surface tension by ten-fold, we get no change in morphologies for MMP and MT-MMP. However, we notice that MMP gives faster growth as the coefficient of haptotaxis χ_{E} =5, whereas for χ_{E} =0.2 and 1, their growth rates are very similar (c.f. **Figure 10**). Growth is faster for MT-MMP when χ_{E} is small, but for larger χ_{E}, comparable growth. So for nutrient rich environment, the surface tension and haptotaxis coefficient mainly determine the growth rates for MMP and MT-MMP. The results of our simulations, therefore, clearly show that the effects of MMP and MT-MMP on tumour growth and morphologies are not significantly different for a tumour growing in a nutrient-rich environment.

The effect of MT-MMP at low proliferation rate and high apoptosis rate in hypoxic environment: Here, we investigate the influence MT MMP mediated proteolysis on tumour growth and morphology for low proliferation rate in a nutrient-poor environment, which has not yet been reported in the literature. Relevant parameter values are given in **(Tables 2 and 4)**. We investigate the effect of haptotaxis under the same conditions described in section 4.1.1, except with a haptotaxis coefficient, χ_{E} =10 and proliferation rate, λ_{p} = 0.1 and with no apoptosis. For the case of the larger proliferation rate, λ_{p} = 0.5, haptotaxis results in stable λ tumour growth. However, for p = 0.1, **Figure 11** shows that the tumour grows into four symmetric buds for MMP in constrast to the MT-MMP case where two symmetric buds are formed. The effects of haptotaxis on tumour growth and morphology are presented in **(Figure 12)**. When χ_{E} =5, the tumour area, perimeter and length scale all increase linearly with time; very similar to behaviour found for the case without taxis. However, the shape parameter is constant, indicating stable circular growth. When taxis coefficient is increased twice the tumour grows linearly at early times but much more rapidly at later times. The LS first increases, then decreases with time, indicating narrowing of the fingers. Similarly, the shape parameters are nearly constant at early times but increase rapidly with decreasing length scale, an indication of tumour surface instability. Thus at low proliferation rates, in the absence of apoptosis, MMP gives faster growth and more fingers than MT-MMP. Similar results we reported in section 4.1.1, where apoptosis and high proliferation rates are considered. So in the absence of necrosis, MMP gives faster growth and more fingers in lownutrient environment irrespective of proliferation rate. Similar tumour morphologies for both MDEs are obtained for higher apoptosis rates. Thus, the difference between the effects of MMP and MT-MMP on tumour morphology is not significant for high apoptosis rates.

The results of our simulation clearly show that there is no significant variation in tumour growth dynamics and shapes when we consider the haptotaxis due to diffusible MMP and non-diffusible MT-MMPmediated ECM degradation, except for low proliferation rate in lownutrient environment, whereas there are significant variations due to haptotaxis in nutrient-poor microenvironment depending on cell proliferation and apoptosis rates.

A better understanding of how the tumour environment affects cancer progression should provide new targets for the isolation and destruction of cancer cells via interference with the complex crosstalk established between cancer cells, host cells, and their surrounding extracellular matrix [5]. So we have proposed a two-dimensional continuum model for the growth and invasion of tumour cells into healthy tissue that focuses on four key components implicated in the invasion process: tumour cells, ECM, MDE (MMPs or MT-MMPs) and nutrient which are initially homogeneous in the tumour microenvironment. We consider both diffusible MMPs and non-diffusible MT-MMP mediated proteolysis and their separate effects on tumour growth via haptotaxis. Our model is the first to examine the effect of haptotaxis for localized ECM degradation by MT-MMP on tumour growth and morphology in high and low-nutrient environments and differentiate the influences of haptotaxis for various cell proliferation and apoptosis rates in hypoxic environment for soluble and non-soluble MMPs.

With the inclusion of haptotaxis, tumour growth and morphology are substantially characterized by whether the tumour is growing in a nutrient-rich or nutrient-poor microenvironment. The invasive fingering morphologies are found in nutrient-poor microenvironments, as a result of haptotaxis, consistent with the results of previous hybrid models whereas in nutrient-rich environments there is no significant influence of haptotaxis on fingered tumour morphology. These results show that ECM degradation by MMP leads to greater instability than that for MT-MMP in low-nutrient environments, even at low proliferation rates; whereas for high apoptosis rates, these both lead to similar unstable morphologies. In summary, our study shows that while haptotaxis lead to completely different tumour growth rates and morphologies, depending on proliferation and apoptosis rates, in lownutrient environment as demonstrated earlier, there are no significant variations when we compare haptotaxis due to ECM degradation by diffusible MMP and surface bound MT-MMP except for low proliferation rate [71,72].

The authors would like to acknowledge many fruitful discussions with Professor Robert Guy and the guidance and mentorship with which he provided the first author during the course of this work.

- Koochekpour S, Merzak A, Pilkington G (1995) Vascular endothelial growth-factor production is stimulated in response to growth factors in human glioma cells. Oncology reports 2: 1059-1061.
- Hanahan D, Weinberg RA (2000) The hallmarks of cancer. Cell 100: 57-70.
- Gerisch A, Chaplain MA (2008) Mathematical modelling of cancer cell invasion of tissue local and non-local models and the effect of adhesion. J Theor Biol 250: 684-704.
- Turner S, Sherratt JA (2002) Intercellular adhesion and cancer invasion a discrete simulation using the extended Potts model. Journal of theoretical biology 216: 85-100.
- Sounni NE, Noel A (2013) Targeting the tumour microenvironment for cancer therapy. Clinical chemistry 59: 85-93.
- D'ortho MP, Will H, Atkinson S, Butler G, Messent A, et al. (1997) Membrane?type matrix metalloproteinases 1 and 2 exhibit broad spectrum proteolytic capacities comparable to many matrix metalloproteinases. European journal of biochemistry 250: 751-757
- Pei D, Weiss SJ (1996) Transmembrane deletion mutants of the membrane type matrix metalloproteinase 1 process progelatinase A and express intrinsic matrix-degrading activity. The Journal of biological chemistry 271: 9135-9140.
- Ohuchi E, Imai K, Fujii Y, Sato H, Seiki M, el al. (1997) Membrane type 1 matrix metalloproteinase digests interstitial collagens and other extracellular matrix macromolecules. The Journal of biological chemistry 272: 2446-2451.
- Poincloux R, Lizarraga F, Chavrier P (2009) Matrix invasion by tumour cells a focus on MT1-MMP trafficking to invadopodia. Journal of cell science 122: 3015-3024.
- Sabeh F, Li XY, Saunders TL, Rowe RG, Weiss SJ (2009) Secreted versus membrane-anchored collagenases relative roles in fibroblast dependent collagenolysis and invasion. The Journal of biological chemistry 284: 23001-23011.
- Maquoi E, Assent D, Detilleux J, Pequeux C, Foidart JM (2012) MT1-MMP protects breast carcinoma cells against type I collagen-induced apoptosis. Oncogene 31: 480-493.
- Deakin NE, Chaplain MA (2013) Mathematical modelling of cancer invasion the role of membrane-bound matrix metalloproteinases. Frontiers in oncology 3: 70.
- Sabeh F, Shimizu-Hirota R, Weiss SJ (2009) Protease-dependent versus independent cancer cell invasion programs: three dimensional amoeboid movement revisited. The Journal of cell biology 185: 11-19.
- Poplawski NJ, Agero U, Gens JS, Swat M, Glazier JA, et al. (2009) Front instabilities and invasiveness of simulated avascular tumours. Bulletin of mathematical biology 71: 1189-1227.
- Poplawski NJ, Shirinifard A, Agero U, Gens JS, Swat M, et al. (2010) Front instabilities and invasiveness of simulated 3D avascular tumours. PloS one 5.
- Pennacchietti S, Michieli P, Galluzzo M, Mazzone M, Giordano S, et al. (2003) Hypoxia promotes invasive growth by transcriptional activation of the met protooncogene. Cancer cell 3: 347-361.
- Erler JT, Bennewith KL, Nicolau M, Dornhofer N, Kong C, et al. (2006) Lysyl oxidase is essential for hypoxia-induced metastasis. Nature 440: 1222-1226.
- Casanovas O, Hicklin DJ, Bergers G, Hanahan D (2005) Drug resistance by evasion of anti-angiogenic targeting of VEGF signalling in late stage pancreatic islet tumours. Cancer cell 8: 299-309.
- Friedl P, Wolf K (2003) Tumour cell invasion and migration: diversity and escape mechanisms. Nature reviews cancer 3:362-374.
- Cristini V, Lowengrub J, Nie Q (2003) Non-linear simulation of tumour growth. Journal of mathematical biology 46:191-224.
- Cristini V (2005) Morphologic Instability and Cancer Invasion. Clin Cancer Res 11: 6772-6779.
- Frieboes HB, Zheng X, Sun CH, Tromberg B, Gatenby R, et al. (2006) An integrated computational/experimental model of tumour invasion. Cancer research 66: 1597-1604.
- Alarcon T, Byrne HM, Maini PK (2003) A cellular automaton model for tumour growth in homogeneous environment. J Theor Biol 225: 257-274.
- Macklin P, Lowengrub J (2007) Nonlinear simulation of the effect of microenvironment on tumour growth. Journal of theoretical biology 245: 677-704.
- Cristini V, Li X, Lowengrub JS, Wise SM (2008) Nonlinear simulations of solid tumour growth using a mixture model: invasion and branching. Journal of mathematical biology 58: 723-763.
- Anderson AR, Weaver AM, Cummings PT, Quaranta V (2006) Tumour Morphology and Phenotypic Evolution Driven by Selective Pressure from the Microenvironment. Cell 127: 905-915.
- Anderson AR, Rejniak KA, Gerlee P, Quaranta V (2009) Microenvironment driven invasion: a multiscale multimodel investigation. J Math Biol 58: 579-624.
- Adams JM, Cory S (2007) The Bcl-2 apoptotic switch in cancer development and therapy. Oncogene 26: 1324-1337.
- Lowe SW, Cepero E, Evan G (2004) Intrinsic tumour suppression. Nature 432: 307-15.
- Hanahan D, Weinberg RA (2011) Hallmarks of cancer: the next generation. Cell 144: 646-674.
- Lowe SW, Lin AW (2000) Apoptosis in cancer. Carcinogenesis 21: 485-495.
- Wong RS (2011) Apoptosis in cancer: from pathogenesis to treatment. J Exp Clin Cancer Res 30: 87.
- Watanabe A, Hosino D, Koshikawa N, Seiki M, Suzuki T, et al. (2013) Critical role of transient activity of MT1-MMP for ECM degradation in invadopodia. PLoS computational biology 9: e1003086.
- Karagiannis ED, Popel AS (2004) A theoretical model of type I collagen proteolysis by matrix metalloproteinase (MMP) 2 and membrane type 1 MMP in the presence of tissue inhibitor of metalloproteinase 2. The Journal of biological chemistry 279: 39105-391114.
- Donze A, Fanchon E, Gattepaille LM, Maler O, Tracqui P (2011) Robustness analysis and behavior discrimination in enzymatic reaction networks. PLoS One 6: e24246.
- Hoshino D, Koshikawa N, Suzuki T, Quaranta V, Weaver AM, et al. (2012) Establishment and validation of computational model for MT1-MMP dependent ECM degradation and intervention strategies. PLoS computational biology 8: e1002479.
- Byrne HM, Chaplain MAJ (1996) Modelling the role of cell-cell adhesion in the growth and development of carcinoma. Math Comput Model 24: 1-17.
- Jeon J, Quaranta V, Cummings PT (2010) An off-lattice hybrid discrete-continuum model of tumour growth and invasion, Biophysical journal 98: 37-47.
- Scianna M, Preziosi L (2012) A hybrid model describing different morphologies of tumour invasion fronts. Mathematical modelling of natural phenomena 7: 78-104.
- Giverso C, Scianna M, Preziosi L, Lo Buono N, Funaro A (2010) Individual cell-based model for in-vitro mesothelial invasion of ovarian cancer. Math Model Nat Phenom 5: 203-223.
- Zheng X, Wise S, Cristini V (2005) Nonlinear simulation of tumour necrosis, neo-vascularization and tissue invasion via an adaptive finite-element/level-set method. Bulletin of mathematical biology 67: 211-259.
- Macklin P, McDougall S, Anderson AR, Chaplain MA, Cristini V, et al. (2009) Multiscale modelling and nonlinear simulation of vascular tumour growth. J Math Biol 58: 765-798.
- Frieboes HB, Lowengrub JS, Wise S, Zheng X, Macklin P, et al. (2007) Computer simulation of glioma growth and morphology. NeuroImage 37: S59-S70.
- Pham K, Frieboes HB, Cristini V, Lowengrub J (2011) Predictions of tumour morphological stability and evaluation against experimental observations. J R Soc Interface 8: 16-29.
- Carmen Calzada M, Camacho G, Fernández-Cara E, Marín M (2011) Fictitious domains and level sets for moving boundary problems. Applications to the numerical simulation of tumour growth. J Comput Phys 230: 1335-1358.
- Wu M, Frieboes HB, McDougall SR, Chaplain MAJ, Cristini V, et al. (2013) The effect of interstitial pressure on tumour growth: Coupling with the blood and lymphatic vascular systems. J Theor Biol 320: 131-151.
- Lee JJ, Huang J, England CG, McNally LR, Frieboes HB (2013) Predictive Modeling of In Vivo Response to Gemcitabine in Pancreatic Cancer. PLoS computational biology 9: e1003231.
- Ovadia J, Nie Q (2014) Numerical Methods for Two-Dimensional Stem Cell Tissue Growth. Journal of scientific computing 58: 149-175.
- Greenspan HP (1976) On the growth and stability of cell cultures and solid tumours. J Theor Biol 56: 229-42.
- Macklin P, Mumenthaler S, Lowengrub J (2013) Modeling multiscale necrotic and calcified tissue biomechanics in cancer patients: application to ductal carcinoma in situ (DCIS). Multiscale Computer Modeling in Biomechanics and Biomedical Engineering, Springer 14: 349-380.
- Gonzalez-Rodriguez D, Guevorkian K, Douezan S, Brochard-Wyart F (2012) Soft matter models of developing tissues and tumours. Science 338: 910-917.
- Rubin H (2003) Micro-environmental regulation of the initiated cell. Adv Cancer Res 90: 1-62.
- Vaupel P, Kallinowski F, Okunieff P (1989) Blood flow, oxygen and nutrient supply, and metabolic microenvironment of human tumours: a review. Cancer research 49: 6449-6465.
- Garber K (2004) Energy boost: the Warburg effect returns in a new theory of cancer. Journal of the National Cancer Institute 96: 1805-1806.
- Esteban MA, Maxwell PH (2005) HIF, a missing link between metabolism and cancer. Nature medicine 11: 1047-1048.
- Stetler-Stevenson WG, Gavil NV (2014) Normalization of the tumour microenvironment: evidence for tissue inhibitor of metalloproteinase-2 as a cancer therapeutic. Connect Tissue Res 55: 13-19.
- Hotary KB, Allen ED, Brooks PC, Datta NS, Long MW, et al. (2003) Membrane type I matrix metalloproteinase usurps tumour growth control imposed by the three dimensional extracellular matrix. Cell 114: 33-45.
- Zaman MH, Trapani LM, Sieminski AL, Mackellar D, Gong H, et al. (2006) Migration of tumour cells in 3D matrices is governed by matrix stiffness along with cell-matrix adhesion and proteolysis. Proceedings of the National Academy of Sciences of the United States of America 103: 10889-10894.
- Dunn G, Heath J (1976) A new hypothesis of contact guidance in tissue cells. Exp Cell Res 101: 1-14.
- Manwaring ME, Walsh JF, Tresco PA (2004) Contact guidance induced organization of extracellular matrix. Biomaterials 25: 3631-3638.
- Testa JE (1992) Loss of the metastatic phenotype by a human epidermoid carcinoma cell line, HEp-3, is accompanied by increased expression of tissue inhibitor of metalloproteinase 2. Cancer research 52: 5597-603.
- Hogea CS, Murray BT, Sethian JA (2006) Simulating complex tumour dynamics from avascular to vascular growth using a general level-set method. Journal of mathematical biology 53: 86-134.
- Aldredge RC (2010) Semi-Lagrangian advection-propagation (SLAP) scheme for three-dimensional interface tracking. Journal of Computational Physics 229: 4683-4702.
- Leonard BP (1988) Universal limiter for transient interpolation modeling of the advective transport equations: the ULTIMATE conservative difference scheme. NASA Technical Memorandum 100916.
- Leonard BP, MacVean MK, Lock AP (1993) Positivity-preserving numerical schemes for multidimensional advection, NASA Technical Memorandum 106055.
- Zhilin L, Kazufumi I (2006) The Immersed Interface Method--Numerical Solutions of PDEs Involving Interfaces and Irregular Domains. SIAM Frontiers in Applied mathematics.
- Peng DP, Merriman B, Osher S, Zhao HK, Kang MJ (1999) A PDE-based fast local level set method. J Comput Phys 155: 410-438.
- Sethian JA, Smereka P (2003) Level set methods for fluid interfaces. Annu Rev Fluid Mech 35: 341-372.
- Macklin P, Lowengrub J (2005) Evolving interfaces via gradients of geometry-dependent interior Poisson problems: application to tumour growth, J Comput Phys 203: 191-220.
- Macklin P, Lowengrub J (2006) An improved geometry-aware curvature discretization for level set methods: Application to tumour growth, J Comput Phys 215: 392-401.
- Nargis NN, Aldredge RC (2016) The Influence of haptotaxis on avascular tumour growth and morphology. Computational and Mathematical Methods in Medicine submitted 1-218.
- Bearer EL, Lowengrub JS, Frieboes HB, Chuang YL, Jin F, et al. (2009) Multiparameter computational modeling of tumour invasion, Cancer research 69: 4493-4501.

Select your language of interest to view the total content in your interested language

- Anaerobic Biodegradation
- Aptamers
- Array biosensor
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Bio-electrochemistry
- Bioactuators
- Bioassay
- Biochips
- Biodegradable Balloons
- Biodegradable Confetti
- Biodegradable Diapers
- Biodegradable Plastics
- Biodegradable Sunscreen
- Biodegradation
- Bioelectronics
- Bioengineering Application
- Biomaterial Science
- Biomedical Equipment
- Biomedical Instrumentation
- Biomedical Services
- Biomedical science
- Biometrics
- Biomimetics
- Bioprocessing
- Bioreactors
- Bioremediation Bacteria
- Bioremediation Oil Spills
- Bioremediation Plants
- Bioremediation Products
- Biosensor Devices
- Biosensor applications
- Biosensor elements
- Biosensor packaging and assembly
- Biosensors
- Biosensors in drug delivery
- Biosensorâs clinical validation
- Biostatistics methods
- Biotechnology Engineering
- Cell Apoptosis
- Chemical sensor
- Clinical Trail
- Cross-Covariance and Cross-Correlation
- DNA-Based Biosensors
- Ex Situ Bioremediation
- Fermentation Science
- Gene Expressions
- Genetic Biochip
- Genetic Linkage
- Genetics
- Heavy Metal Bioremediation
- Hypothesis Testing
- Immuno sensors
- In Situ Bioremediation
- Integrated nanoscale devices
- Internal Medicine
- Label-Free Biosensor Cell Assays
- Large-scale Survey Data
- Matrix
- Medication
- Microarray Studies
- Microassay
- Microfluidic Biochips
- Microfluidics biosensors
- Microlithography
- Molecular and Medicine science
- Molecular recognition biomolecules
- Multivariate-Normal Model
- Mycoremediation
- Nanochips
- Nanorods
- Nanosensor
- Neural sensor
- Non Biodegradable
- Non rigid Image Registration
- Organic Electronics
- Photonic sensing
- Phytoremediation
- Protein Microarrays
- Regenerative Medication
- Regressions
- Rehabilitation Technology
- Robust Method
- Sewage Water Treatment
- Soft biometrics
- Soil Bioremediation
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Tissue Engineering
- Tissue Engineering Development
- Types of Upwelling
- Waste Degredation
- Wearable Biosensors
- Xenobiotics

- 5th Annual Congress and Medicare Expo on
**Trauma & Critical Care**March 05-06, 2019 Amsterdam, Netherlands

March 05-06, 2019 Amsterdam, Netherlands - 2nd Annual Summit on
**Cell Therapy**and**Regenerative Medicine**November 9-10, 2018 Atlanta, Georgia, USA

November 9-10, 2018 Atlanta, USA

- Total views:
**2058** - [From(publication date):

December-2016 - Jan 19, 2019] - Breakdown by view type
- HTML page views :
**1960** - PDF downloads :
**98**

Peer Reviewed Journals

International Conferences 2019-20