2Department of Haematology, Western General Hospital, Edinburgh, EH4 2XU, United Kingdom
3Department of Chemistry and Biochemistry, Medical University of Sofia, 1431 Sofia, Bulgaria
Abstract
The clinical course of B-CLL is heterogeneous. This heterogeneity leads to a clinical dilemma: can we identify those patients who will benefit from early treatment and predict the survival? In recent years, mathematical modelling has contributed significantly in understanding the complexity of diseases. In order to build a mathematical model for determining prognosis of B-CLL one has to identify, characterise and quantify key molecules involved in the disease. Here we discuss the need and role of mathematical modelling in predicting B-CLL disease pathogenesis and suggest a new systems biology approach for a personalised therapy of B-CLL patients.
Citation: Tummala H, Goltsov A, Khalil HS, Sproul A, Scott F, Mitev V et al. Advocating the need of a systems biology approach for personalised prognosis and treatment of B-CLL patients. Biodiscovery 2012; 6: 4; DOI: 10.7750/BioDiscovery.2012.6.4
Copyright: © 2012 Tummala et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, provided the original authors and source are credited.
Received: July 14, 2012; Accepted: December 23, 2012; Available online /Published: December 30, 2012
Keywords: B-CLL, cancer, systems biology
*Corresponding Author: Nikolai Zhelev: e-mail: N.Zhelev@abertay.ac.uk
Conflict of Interests: No potential conflict of interest was disclosed and there is not any competing financial interest in relation to the work described.
1. Introduction
B-cell chronic lymphocytic leukemia (B-CLL), is the most common leukaemia in the western world. It is not curable [1] and the ability of radiation and chemotherapeutic agents to induce clinically significant regression of the disease in patients with CLL has not improved for the last decades [2]. The development of treatment approaches for B-CLL has lagged in comparison to other hematologic malignancies for various reasons. B-CLL disease seems especially prone to familial occurrence being nearly 3 times higher than that expected for the general population [3]. Historical reports exists for identical twins [4], mother and son [5], grandfather, son, and grandson [6] affected due to familial aggregation of chronic lymphocytic leukaemia pathogenesis associated with immune defects. In 1980 Conley et al [7] reported an increasing frequency of B-CLL disease associated with autoimmune diseases, concluding that genetic factors in these families are the prime cause of deregulation in the immune system. Many reports from the Swedish family cancer database [8] evidently reported that there is significantly increased familial risk involved in B-CLL disease pathogenesis. CLL risks were similar in parents, siblings, and offspring cases, and in male and female relatives. They were not affected by the case’s age at diagnosis.
Although prognostic markers have been identified, there is as yet no proven indications for initiating treatment in patients with asymptomatic disease.. The lack of accuracy in predicting disease progression and survival on an individual basis has been a long asked research question in B-CLL. It’s a known fact that in B-CLL disease pathogenesis to date, no single biomarker has been useful in predicting the prognosis in CLL patients. This review focuses on our current understanding of clinical course and molecular mechanisms in B-CLL disease, and suggests the application of systems biology as an approach in personalized treatment for B-CLL patients.
2. Clinical paradigm
The clinical course of B-CLL disease is highly variable, with life expectancies ranging from months to decades. There is no apparent survival advantage for early treatment intervention, and some patients may never require treatment. At present the information available is only used to counsel patients and to help inform clinicians regarding the frequency of monitoring. The currently available clinical staging systems (e.g. Rai or Binets test) for B-CLL are simple and inexpensive but lack accuracy to predict disease progression and survival on an individual basis. Prediction was not significant when analyzed using life table methods performed by Goldin et al. [8] and it is thereby concluded that the familial component of CLL is shared with other lympho proliferative malignancies, suggesting common genetic pathways. As infection is a major cause of morbidity and mortality in B-CLL, previous studies on infection predictors in CLL were based on patients receiving chemotherapy, or on those diagnosed with infections during the ‘watch and wait’ period and during therapy [9-12]. Treatment strategies in preventing recurrent infection procedures in these patients, however, have not yet been devised [13].
B-CLL disease is characterized by an accumulation of mature, non-proliferating B lymphocytes in the blood, spleen, lymph nodes and bone marrow. The accumulation of B cells is the result of a proliferative defect, failing B cells to undergo apoptosis and leading to large numbers of cells being blocked in the G0/G1 phase of the cell cycle [14]. Not surprisingly, many studies have shown that B-CLL cells are resistant to various drug induced apoptosis [15]. The reasons for this are unclear, but may include a combination of various factors like reduced expression of several apoptosis inducing proteins that include receptors, kinases and enzymes.
The diagnosis of B-CLL is carried out using a number of classic and novel sensitive techniques, which allow diagnosing and differentiating the disease from other chronic lympho proliferative disorders [16-18]. Morphology and immuno phenotype are the initial diagnostic investigations, alongside molecular genetics and/or histology to exclude other B cell disorders. To date immuno phenotyping is the only method that can determine or confirm the diagnosis, as B-CLL lymphocytes carry a distinct immuno phenotypic signature (more than 5 x 109 B cells/litre). A scoring system compounding the results with a set of markers allows a refined confirmation of the diagnosis as described by Matutes et al. [18], but this scoring system cannot be applied globally to patient management, as the B-CLL disease heterogeneity exists irrespective of mutated and unmutated pattern of the IgVH region.
B-CLL disease has been presented as two distinct types. The presence of somatic mutations in the immunoglobulin (IgH) heavy chain gene define a group of patients exhibiting stable or slowly progressive disease that requires late or no treatment. By contrast, the absence of mutations in the IgH genes of B-CLL cells define a group, which exhibit a progressive clinical course requiring early treatment [19]. Other immunological markers such as surface antigens (for example, CD38 and ZAP-70) have an important prognostic impact but still lack accurate prediction in patient survival rates and overall treatment response rates. Moreover, optimal drug administration schedules are required in order to improve clinical response in CLL patients.
3. Cell cycle, proliferation and apoptosis
B-CLL lymphocytes lack the fundamental functions of regulating the cell cycle and apoptosis. That is why the cellular factors that govern the entry of resting B-cells into the cell cycle and promote their progression through the different phases of cell cycle, are considered to be optimal targets for drug therapy. Normally the matured naive lymphocytes circulating in blood remain quiescent in the G0 phase and reside in lymphoid tissues such as the spleen until they encounter an appropriate mitogenic stimulus. The progression of the cell cycle in naive B cells depends on the PI3K activation of AKT activation, and inducing NF-kB activation and cyclin D2 activity by deregulating FOXO transcription rate [20-21]. The rate at which B lymphocytes enter the cell cycle is tightly regulated by a very complex process involving multiple components of the cell cycle machinery, apoptosis, and the antibody-generating apparatus and by downstream effectors of Bcl-2 pathway (e.g. caspases).
The proliferation rate and resistance to drug-induced apoptosis in B cells are recognized as important factors in the outcome of treatment in CLL disease. Mutated and unmutated CLL patients clearly differ in terms of prognosis. The genomic aberrations occurring in B-CLL patients have their individual effects due to distinct expression profiles, as reported by Klein et al [22]. This is termed as gene dosage effect, which is correlative with the variable mutational status of IgVh [23]. CLL cells are dependent on the continuous expression of intrinsically short lived anti-apoptotic proteins for survival. This biological context has been a promising strategy to induce apoptosis by the inhibitors of transcription and translation. High levels of Mcl-1 and Bcl-2 mRNA [24] and protein [25] have been found in B-CLL. Mcl-1 is associated with the failure in response of B-CLL patients to initial therapy, due to its high expression [26, 27] whereas low expression of Mcl-1 in CLL cells is related to prolonged patient survival [24]. The presence of Bcl-2 (B-cell B-CLL/lymphoma 2), an anti-apoptotic protein that is highly expressed in B-CLL, and which cannot be neutralised with the levels of pro-apoptotic proteins, has been a characteristic feature for B-CLL cells when compared to normal B cells [28]. B-cell antigen receptor (BCR) is the key signalling molecule of B-lymphocytes, triggering pathways involved in B-cell proliferation, survival, differentiation, and apoptosis. In B-CLL, the status of the BCR is of important prognostic, biological, and therapeutic value. This is reflected in the great variety of emerging novel agents targeting the BCR signalling cascade.
To date large components of BCR signalling pathways have been elucidated in determining the prognosis of B-CLL, but failed to predict accurately. Pharmacological inhibition of BCR signalling in CLL can be achieved by targeting components of the BCR signalling cascade. Several groups have reported an increase in MAPK (ERK2) phosphorylation and NFAT transcriptional activity and a corresponding lack of AKT phosphorylation due to constitutive activation of distinct BCR signalling pathways in a subset of particular CLL cases [29]. The aberrant signalling mechanism in BCR signalling cascade includes the activation of PKC and PI3K pathways, which regulate B cells survival. The signal transduction inhibitors in these pathways have resulted promising results in vitro, but could not be used as effective single agents for B-CLL therapy. Existence of molecular heterogeneity in CLL that is dependent on the expression of genes defining BCR signalling has been pivotal, but still lacks in attaining accurate prediction of patient survival in clinical course.
Increased aggregation of B lymphocytes due to a lack of apoptotic signalling has been a well known cause in B- CLL disease pathogenesis, suggesting that B-CLL cells lack the fine tuning of a balance between pro-apoptotic and anti apoptotic factors in the microenvironment. The anti-apoptotic protein Bcl-2 is highly expressed in B-CLL patients and has been well characterized, but paradoxically, high levels of pro-apoptotic molecules have also been observed [30, 31]. Deregulation of the cell cycle regulators (CCND3, CCND1, and MYC) and apoptosis regulators (NPAT, CUL5 and PPP2R1B), and a role for these genes, has been implicated in the pathogenesis of B-CLL disease [32]. Moreover there is a rapidly growing number of the modulators of BCR components and the cell death machinery. These can be targeted to disrupt the turgid balance between pro apoptotic and anti apoptotic factors, and restore healthy balance, and are considered to be of optimal therapeutic value in the treatment of CLL [33, 34].
4. Microenvironment: Cross talk and interplay
Microenvironment regulating signals are crucial for the process of signaling in leukemic cells that induce proliferation and lead to the survival and accumulation of leukemic cells within lymphoid organs. The in vivo accumulation of leukaemic lymphocytes is potentiated by interactions of CLL cells with other cells such as mesenchymal stromal cells (MSCs), nurse-like cells (NLCs), T cells and other soluble factors that include Interleukins (IL-4) and cytokines such as CCL22 and CCL17 [35]. CLL cells grow in the microenvironment of a solid tissue where they initiate as monoclonal disease of autoreactive B cells, due to receptor-ligand interactions resembling an autoantigen and neoantigen. A process called auto antigen drive has been shown to be essential in selecting susceptible B cells out of the normal clone and driving them into the disease state [36].
Common genetic aberrations that occur in CLL disease and different stimuli originating from the microenvironment, cooperate in the selection and expansion of the malignant clone. CLL cells relentlessly accumulate in vivo but rapidly undergo spontaneous apoptosis in vitro. This implies that their apoptosis resistance, rather than being an intrinsic feature of leukemic stage, depends on external signals for survival since the maturation stages of B cells are highly dependent on microenvironment signals [35].
Exemplified crosstalk between CLL cells and accessory cells occurs in the marrow and/or lymphoid tissue microenvironments. Contact between CLL cells and NLCs or MSCs is established and maintained by chemokine receptors and adhesion molecules. NLCs express the chemokines CXCL12 and CXCL13, whereas MSCs predominantly express CXCL12. NLCs and MSCs attract CLL cells via the G protein coupled chemokine receptors CXCR4 and CXCR5, which are expressed at high levels on CLL cells [35]. CD38 expression allows CLL cells to interact with CD31, the ligand for CD38, expressed by stromal and NLCs. Ligation of CD38 activates ZAP-70 and downstream survival pathways. Self and/or environmental antigens (Ag) are considered a key factor in the stimulation and expansion of the CLL clone. BCR stimulation and co-culture with NLCs also induce CLL cells to secrete high levels potent T-cell attracting chemokines such as CCL3 and CCL4.
Through this mechanism, CLL cells can actively recruit T cells for cognate T-cell interactions with CLL cells. CD40L T-cells are preferentially found in CLL proliferation centers and can interact with CLL cells via CD40. Cytokines secreted by T-cells or CLL cells such as TNF-α or IL-4 are considered important regulators of CLL cell survival. Moreover, downregulation of proximal and intermediate T-cell receptor signalling cascades and globally reduced cytokine secretion, contributes to significant immunodeficiency in non-leukemic indolent B-cell lymphomas [37]. Collectively, this crosstalk between CLL cells and accessory cells results in activation of survival and drug resistance pathways [36]. Focusing on the therapeutic aspect of regulating the microenvironment, this includes disruption of targets in the complex interplay between CXCR’s and adhesion molecules, which is the precise definition of systemic T-cell dysfunction. This might eventually provide a rationale for appropriately tailored molecular therapies targeting not only B-CLL cells but also the microenvironment, to ensure improved treatment.
5. Micro RNAs
Micro RNA’s (miRNAs/miR’s) are a recently discovered class of molecules that regulate gene expression at the posttranscriptional level. The interactions between miRNAs, target genes, and pathways in CLL are clearly complex, as are the links between genotype and phenotype miRNA machinery, which include miRNAs involved in many cellular processes such as proliferation and apoptosis. Some microRNAs are referred to as oncomiR’s and exhibit differential expression level in cancer disease by acting as either oncogenes or tumour suppressor genes [38, 39]. Recent in vivo studies demonstrated that tumours became addicted to oncomiR’s. Their inactivation resulted in complete tumour regression, demonstrating the importance of miRNAs in therapeutic pharmacological intervention (for example miR-21). A role for miR-21 in cancer has been implicated strongly by the fact that it is overexpressed in haematological malignancies, especially B-CLL [38]. Several studies have proposed the use of microRNA-based classifications as signature profiles for B-CLL disease prediction and survival [41]. miR’s studies are likely to have an increasing influence in the diagnosis, prognosis and treatment of human cancers, including CLL [42].
Cimmino et al. [43] proposed that miR-15a and miR-16-1 function as tumor suppressor genes by modulating Bcl-2. The same miR signature was also associated with the presence or absence of disease progression, since Bcl-2 mRNA stability is not affected by the overexpression of these miR’s but instead was regulated by Bcl-2 expression at a posttranscriptional level. miR-34a expression has been found to be very low in B-CLL patients. B-CLL cells with reduced miR-34a expression showed increased viability after DNA damage, independently of 17p status. Low expression of miR-34a in CLL is therefore associated with p53 inactivation as well as chemotherapy refractory disease, impaired DNA damage response, and apoptosis resistance irrespective of 17p deletion/p53 mutation [44]. This study identified the role of miR-34a in B-CLL pathogenesis where they speculated that it took over the control of apoptotic machinery in both p53 dependent and independent pathway by controlling the expression of Bcl-2 and p21 [45].
The demonstration of down- regulation of cyclin-E2, CDK6, E2F5, cyclin D1 (CCND1) and Bcl-2 protein levels after induction of miR-34a expression has been noted to be important, suggesting that miR-34a in CLL mediates its functions within and potentially outside the p53 pathway [46]. A recent study report by Asslaber et al. [47] showed that a low level of miR-34a expression profile predicted a shorter time for treatment of disease in B-CLL patients, alternatively an over expression of mIR34-a induces apoptosis through p53 attenuation by negative feedback loop regulatory mechanism. This study specifically indicates that the p53 tumour-suppressor pathways and associated miR-34a expression may be a promising avenue of research into the pathogenesis and prognostic assessment of CLL, and may prove more useful in aiding early-stage identification of patients with fludarabine refractory CLL. Moreover karyotype specific signature profiles have been found to be associated with different forms of B-CLL such as indolent and aggressive [41].
Analysis of miRNA expression values between the predefined subgroups showed that the miR-223, the miR-29b, the miR-29c, and the miR-181 family are down regulated in 17p-aggressive cases showing the highest level of Tcl-1 expression compared with 17p-indolent cases [48]. The findings in this study suggest that interactions between Tcl-1 and miRs play an important role and suggest them as new markers to assess the disease course in CLL patients. Although such findings are exciting, larger studies demonstrating reproducibility are still required. Targeted miRNA therapeutics remain at the early stages of development and are limited primarily to in vitro and murine models of disease. Given the fact that miRNA’s expression profiles may prove to be novel surrogate markers for common cytogenetic lesions, the prognostic indicators like karyotype specific signature in B-CLL disease is highly difficult to achieve through miRNA based therapy. To conclude, micro RNA studies indicate the importance of considering more than one parameter to assess the tumour burden in B-CLL, but could not advocate a replacement of standard clinical and molecular markers. The research findings suggest that the addition of the miR’s expression is critical for the evaluation and more accurate management of CLL patient care.
6. Drug targets and resistance
Adequate progress has been made in recent times for therapeutic intervention and management B-CLL disease. New insights into the molecular pathology of B-CLL have generated a plethora of biological markers that predict the prognosis and influence therapeutic decisions. These markers include historical Rai and Binet staging systems, loss of p53 and ATM functions, the unmutated IgVH gene, and high expression of ZAP-70 or CD38 [49]. These factors in B-CLL cells have been identified as unfavourable prognosis to standard chemotherapy or are indicative of aggressive disease. Several strategies that induce apoptosis and circumvent resistance by a combination of inhibitors of transcription and translation such as CDK inhibitors, together with other approaches that interfere with the function of anti-apoptotic proteins that initiate synergistic killing in B-CLL (flavopiridol, R-Roscovitine, Danciclib and SNS-032) have been employed [50-53].
A recent review by García-Escobar et al. [54] stated that preclinical evidence for the efficacy of signal transduction inhibitors such as alkylating agents [55] and purine analogues [56] in B-CLL is particularly encouraging, but the results of clinical trials on using these as single agents have not been satisfactory. Combinational therapy on the other hand has been proved efficacious and is currently being explored in pre-clinical studies and early-phase clinical trials. These have not been used routinely in the clinic because they exhibit differential responses due to poor diagnostic and prognostic marker relevance. To date drugs like fludarabine, bendamustine and the two monoclonal antibodies, alemtuzumab and rituximab have been approved by European and/or American regulatory agencies. Monoclonal antibody Alemtuzumab holds a promising combination of therapeutic modality with fludarabine as a first line treatment for refractory CLL patients. Recently haematopoietic stem cell transplantation such as low toxic non-myeloablative allogeneic transplant has opened new perspectives in the management of CLL patients, and appears to be a promising therapy for elderly patients.
Additional monoclonal antibodies targeting CD20, CD23, CD37, CD38 or CD40, as well as drugs designed to interfere with proteins regulating the cell cycle, apoptotic machinery or leukemic microenvironment (e.g., flavopiridol, oblimersen, ABT-263 or lenalidomide) are being investigated in clinical trials. To date, single-agent clinical trials have indicated that the major clinical outcome is the stabilization of disease states. Lyse [57] showed that pre-treatment of TRAIL resistant B-CLL cells with histone deacetylase inhibitors such as Oxam, rendered cells to be susceptible for subsequent killing by recombinant TRAIL protein. This therapy preferentially increased surface expression of TRAIL-R1 and TRAIL-R2, there by leading to caspase 8 and caspase 3 activation. Importantly, this is a combination therapy that enables preferentially sensitized transformed cells to apoptosis. Recently Niedermeier et al. [58] and Herman et al [59] demonstrated that the isoform-selective PI3K inhibitors induced B-CLL cell apoptosis in vitro in MSC co-cultures synergized with fludarabine induce B-CLL cell apoptosis. Thus, combination therapies using specific and recombinant inhibitors in signal transduction pathway might prove to be efficacious for the treatment of B-CLL.
The main focus of therapeutic strategies in B-CLL is still cytotoxic chemotherapeutically available alkylating agents or purine analogues that trigger DNA damage response via p53 leading to a prominent cell death. Most of B-CLL diseased patients carry defects in the p53 pathway and therefore it has been challenging to overcome the resistance through p53 independent cell death pathways [60, 61]. Drug- induced resistance to B-CLL treatment enhances by activating survival signaling pathways (for example NF-kB) through activation of PI3K-Akt and therefore counteracting the effect of the treatment. Microenvironment plays major roles in B-CLL tumour resistance to therapies. Most of the factors and pathways discussed above that are responsible for sustained activation of survival pathways in B-CLL are provided by contact with the microenvironment.
7. Systems biology approach
Although clinical stages such as Rai and Binet tests are the required basis for the prognosis in the B-CLL patients, many other biological markers mentioned in the above section have been offering important insights to prognostic information. Nonetheless these prognostic factors have not been fully validated or standardized in large clinical trials. Moreover, it is practically impossible to validate all possible biomarkers in clinical trials due to the fact that response upon drug dosage in B-CLL patients is varied, and therefore could not sustain in clinical evaluation. Existence of conflicting data in research studies has added more complexity and confusion in the evaluation of B-CLL disease pathogenesis. For example Klein et al [22] reported that a large number of apoptotic genes deregulated in particular with 17p deletion, which reflected the distinctly aggressive biology of a subgroup cohort study in B-CLL patients and contrasted with reports published by Stankovic et al [62] and Schaffner et al [63], showing an indistinguishable native gene expression pattern in terms of quantitative folds when comparing wild type and TP53-mutant in B-CLL cells acquired from patients.
Notably, Klein et al [22] used highly sensitive quantitative RТPCR approach and a greater number of case studies in evaluating the data, when compared to the Stankovic et al. [64] microarray study. Moreover, low Bcl-2 levels in 17pdel, the subgroup with the worst clinical course was reported by Klein et al [22] due to the known anti-apoptotic function of Bcl-2 and the association of high Bcl-2/BAX ratios with aggressive disease, which argues against a relevant role of Bcl-2 overexpression or elevated Bcl-2/BAX ratios in the pathogenesis of the B-CLL or in B-CLL prognosis. The Kim et al. [65] study supported the use of a whole genome sequencing approach for comprehensively decoding the B-CLL genome in order to better understand the genetic defects caused by insertion, deletion, base change, and restriction site polymorphism present in both coding and non-coding regions in hundreds of loci within the genomes, which indicate the wide presence of small genomic aberrations in B-CLL. B-CLL cases display recurrent genetic aberrations including trisomy 12 and monoallelic or biallelic deletion/inactivation of chromosomal regions 17p, 11q and 13q14 [66]. Presently these deletions and the mutational status of the IgVH gene provide prognostic information and may determine the type of therapy [67].
The above mentioned molecular signatures may have a profound impact on prognosis for B-CLL disease pathogenesis, but the lack of routine clinical methods for diagnosis is a major challenge which should be accounted for in future research. Taking into consideration our current understanding of B-CLL disease pathogenesis, it is now convincing that application of systems biology should be employed as an approach for predicting patient response and designing individualized therapies. An iterative interplay between biological experiments producing quantitative data followed by suitable modeling strategies should enable us to interpret the biology in a global system perspective, thereby increasing our current understanding of the signaling pathways involved in the pathogenesis of a complex disease such as B-CLL. This approach should also include the extracellular aspects of the disease mechanism such as microenvironment. Due to the fact of these differential study conclusions, it is now time to consider prognosis and treatment on an individual basis rather than cohort studies involving large numbers of patients. Moreover it is now evident that B-CLL is not a single target regulated disease, but rather has multiple targets with two distinct disease phenomenae (indolent and aggressive forms).
Building network models by applying systems biology is an approach that will enhance our understanding of key molecular players in the B-CLL disease pathogenesis. To generate a computational model for prognosis of B-CLL patients in silico, one has to overcome a number of challenges such as working with a large number of chemically diverse molecules that are often rapidly changing their concentrations in terms of their expression and regulation [68]. In order to overcome these challenges, various sensitive detection systems or gene reporter systems should be implemented to detect the optimal expression in vivo [69] before its execution in building a model. It is clear from the many experiments and studies carried out to date that there are no single genetic or molecular aberrations which alone lead to the B-CLL phenotype.
Disease progression can vary considerably from one patient to another suggesting that future therapies will need to be designed on an individual basis. To do this, it requires an in-depth analysis of two major groups of variables. Firstly by encompassing the molecular interactions between the various intra-cellular proteins which implement the cell’s regulatory pathways, in particular survival, anti-growth, proliferation and apoptosis (Figure 1). Secondly, as these interactions are highly complex at the intra-cellular level where the loss of homeostasis has to be analysed, a model reporting on the molecular activity, and the imbalance between B-cell proliferation and apoptosis should be addressed. Due to this complexity, some form of holistic approach to the problem is required and we believe that one way forward is to describe the dynamics of the various molecular and cellular interactions by using a mathematical model. In this way the mechanistic links, which underpin the B-CLL disease phenotype can be defined and, subject to calibration at the individual patient level, these models would enable a better understanding of disease progression, as well as providing a methodology for targeted drug interventions.
The key aspects in cellular signaling that involve the cell cycle, apoptosis and various intracellular pathways such as DNA damage repair and cytokine receptor activated network activating PI3K/AKT, are considered to be crucial components that are frequently deregulated in B-CLL disease. Each pathway is subject to an intensive study by the computational systems biology approach: modelling of the cell cycle [65-73], apoptosis [74-76], MAPK and PI3K/AKT pathways [77, 78]. As these pathways are involved in cell proliferation, survival and cell death, computational models of these pathways are considered promising tools in cancer research for prediction of cancer disease progression, development of biomarkers, and drug therapy efficacy [79-82]. Henceforth, understanding the complex regulation of B-CLL disease pathogenesis, requires a strategic knowledge on the models of the cell cycle, and apoptosis that report on intrinsic and extrinsic cell signaling network. Systems computational approach to the complex regulation in B-CLL pathogenesis is now at the initial stage of development. We discuss below recent results from computational modeling of these pathways and address some specific results relevant to B-CLL disease.
The mathematical approach to modelling the cell cycle was developed in a series of works by Novak and Tyson [70, 83, 84], which demonstrated that the cell cycle could be modelled as an irreversible sequence of changes between successive steady states, driven through changes in cyclin levels and loss of their stability at cell cycle phase transition. Further development of this approach was made in a series of models which successively described key features of cell cycle dynamics, and predicted main control mechanisms of the transitions in the cell cycle (see a review of cell cycle models in [85]. The general result of cell cycle modelling showed the extreme robustness of cell cycle function at a wide variation of model parameters such as protein-protein interaction constants, rate constants of protein phosphorylations, protein expression and degradation rate constants [70, 80, 81]. These findings are in accordance with experimental observation that the cell cycle is one of the highly controlled processes in the cell [84]. The decision making for cell division, its delay, arrest or apoptosis is tightly orchestrated by integrating both external cellular signals, such as growth factors and hormones, and internal signaling molecules. Cell cycle modelling demonstrated that despite extensive control mechanisms in regulating phase transitions, it can lose its control and become a fragile system when variations in specific parameters of the model that affect robustness are introduced. From a clinical point of view, this fragility in the system can be referred to as cell cycle malfunction, observed as uncontrolled cell proliferation, a common phenomenon underlying many cancers [82].
7.1 Identification of fragile nodes
In this context a new area in systems biology is being developed that aims at not only the description of the normal functioning of cellular control systems, but also the modelling of the mechanics of their malfunction in pathological conditions including cancer [74, 81, 86]. Novel approaches were designed to analyze the existing models of the cell cycle to study the abnormality in its control mechanisms leading to uncontrolled cell proliferation and survival, causing diseases including solid and hematological cancers [80, 82]. One of the challenges in systems biology is in determining the fragile, sensitive nodes of signaling pathways, including the cell cycle regulation, which cause the disruption of the balance between cell proliferation and death. These fragile points are suggested as control hubs within cellular signaling networks, which are extremely sensitive to mutations or the aberrant protein expression observed in cancers [82]. These fragile hubs, possessing increased sensitivity to internal perturbations, are assumed to also be sensitive to external perturbations such as drug action. Fragile points in signaling networks therefore represent promising targets for drug intervention [80-82].
One of the powerful methods for determining fragile nodes within cellular signaling networks is sensitivity analysis of the computational models [79, 82]. Sensitivity analysis of the model consists of the calculation of sensitivity coefficients: ‘S’ which represents the response of a network output, ‘ΔY’ to a variation of model parameters, ‘Δp’ (protein-protein interaction constants, rate constants of protein phosphorylation, protein expression, degradation rate constants, others): S=ΔY/Δp. Sensitivity analysis permits ranking sensitivity of different outputs of the model to variations of the model parameters and determining a network locus which is most sensitive to these perturbations. Points in the network that exhibit extreme sensitivity to small perturbations of the parameters are often referred to as “fragile” (the converse of “robust”) [74, 80].
The results of sensitivity analysis of the cell cycle models obtained by Nayak S et al. [82], screened for fragile mechanisms in the existing cell cycle models using sensitivity analysis and predicted that such mechanisms are implicated in various solid and hematological cancers including B-CLL. The sensitivity coefficients calculated in sensitivity analysis [82] of the G1/S model of Qu et al. [71] and the G2/M-DNA damage model of Aguda [72] are shown in figure 2. Sensitivity coefficients with high ranks showed that the following reactions can be classified as the most fragile steps in the G1/S checkpoint: 1) dephosphorylation of CDC25, 2) expression of cyclin E, 3) degradation of the cyclin E-CDK2 complex, and 4) level of the transcription factor E2F. Experimental and clinical data supported the computational identification of the fragile reactions in the cell cycle which represent promising targets for drug therapy. The first ranked reaction is the activation of CycE-Cdk2 by CDC25 phosphatases in the G1/S network. The significant role of CDC25 phosphatase was observed in many solid and haematological cancer progressions, and several CDC25 inhibitors showed promising results in cancer treatment [82].
The second ranked fragile mechanism in the G1/S network obtained in sensitivity analysis is cyclin E-CDK2 complex. Inhibition of the active cyclin E-CDK2 is being considered as a treatment strategy in different types of cancer including B-CLL. For example, a synthetic flavone, flavopiridol induces cell cycle arrest by inhibiting multiple CDKs, and leads to p53-independent apoptosis in CLL cells [87]. The first ranked fragile points in the G2/M-DNA damage network is the activation of preMPF (cyclin B–CDK1 complex) catalyzed by CDC25 [82]. This finding correlates with the CLL treatment by bryostatin-1, a protein kinase C (PKC) inhibitor and antagonist of cyclin B-CDK1 complex [88]. The other types of fragile mechanisms obtained in sensitivity analysis [82] are related to protein synthesis and degradation (see sensitivity coefficients for concentrations of E2F, synthesis of CycE, CDC25, and other proteins in figure 2). Targeting these mechanisms now becomes an active area of therapeutic strategy by developing the manipulation methods of translational machinery and programmed proteolysis [82].
7.2 Signalling simulation models
The signaling pathway relevant to B-CLL disease is PI3K/AKT pathway, which plays a key role in cell proliferation, differentiation and survival [58, 59]. Constitutive activation of this pathway by chemokine and cytokine receptors in CLL cells results from interaction of leukaemia cells with their environment, mainly due to chemokines secreted by stromal cells, which activate CXCR4 receptors and downstream pro-survival PI3K/AKT signaling [89]. Computational modeling of the PI3K/AKT pathway was a subject of many theoretical studies, which focused on different aspects of the complexity in this signaling pathway [77, 78, 90-92]. Development of large-scale models of the PI3K/AKT network, reported the response of the signaling network upon different cytokine stimuli [91] and elucidated the control mechanisms in AKT signaling cross talk with other pathways [92]. The application of the computational modeling of PI3K/AKT signaling in cancer research [93] is aimed at the description of activation signaling in different cancer cell lines [78], analysis of the roles of cancer-driven mutations in this pathway [78,94], modeling of the response of the pro-apoptotic signal to anticancer drugs [78, 92], elucidation of mechanisms of drug resistance [95], identification of drug targets, [78, 96], optimization of anticancer therapy, and also the development of combination therapy [78, 97].
Sensitivity analysis of these models being a powerful method for analysing PI3K/AKT signaling in cancer was carried out in many works to identify fragile points in this pathway [78, 91, 96, 97], to determine sensitivity of the pathway to different mutations involved [97], and to find promising drug targets [78, 96, 97]. Note the different methods of sensitivity analysis used in different works [78, 91, 96, 97]. Below we discuss the results of sensitivity analysis of PI3K/AKT pathway model carried in [96, 97], which can be used to analyse drug targets within this network.
The results of sensitivity analysis of the PI3K/AKT pathway model [97] are given in figure 3. The comparative analysis of sensitivity coefficient ranking for different network modules revealed that one of the most sensitive loci within this network is the PI3K/PTEN cycle which forms a regulatory hub within cellular signaling. It includes PI3K and phosphatase PTEN, which jointly control the pool of the second lipid messenger, phosphatidylinositol-3,4,5-trisphosphate (PIP3). This conclusion is in agreement with the results of an integrative genomic and proteomic analysis of solid [98] and haematological cancers [99] which revealed that this regulation hub is vulnerable to mutations in cancer. PI3K is considered an attractive treatment for targeting several types of cancers including B-CLL [89]. The recent development of specific PI3K inhibitors targeting different PI3K catalytic subunits enable the inhibition of the p110δ isoform of PI3K, which is reported to be specific to hematopoietic cells [58, 89]. Clinical testing of CAL-101, a selective inhibitor of the PI3K p110δ isoform showed a suppression effect on stromal cell-derived survival of CLL cells [89, 100].
7.3 Apoptosis simulation models
The high sensitivity of PI3K/AKT signalling to AKT expression level revealed in sensitivity analysis (see figure 3B) showed its significant impact to anti-apoptotic signal. Activated AKT inhibits apoptosis by phosphorylation of Bad protein, procaspase-9, and involves in NF-kB activation, so AKT inhibition is an attractive target for CLL drug therapy. Preclinical trials showed that AKT inhibitor, A-443654 induces apoptosis in B-CLL cells [101, 102]. The cell cycle and PI3K/AKT pathways discussed above are directly connected to the apoptotic signalling pathway which plays a significant role in the pathogenesis in B-CLL, and is considered an attractive target in the development of novel strategies in CLL therapy [89].
A large number of theoretical studies are devoted to computational modelling of both the intrinsic and extrinsic apoptotic pathways [74, 75, 103-110]. The theoretical modeling of intrinsic apoptosis pathways are focused in part on apoptosis induced by DNA damage signals, and includes the description of p53-Mdm2 pathway with ATM/CHK2, ATR/CHK1 and p21 components, which regulate DNA repair, cell cycle arrest, and apoptosis (see figure 1) [80, 104, 106, 108]. Sensitivity analysis of p53-Mdm2 pathway models showed that the dynamics of the DNA damage response are sensitive in part to the level of p53 which dynamically changes in different cells [106, 108]. The fragility of the DNA damage pathway leading to promoted cell survival and genome instability, can play a particularly significant role in pathogenesis of B-CLL [89]. TP53 mutations are considered to be biomarkers of the overall survival, treatment-free survival, and relapse in CLL patients who showed the high prognostic value independent from del[17p]. 4.5% of untreated patients with CLL have TP53 mutations in the absence of del[17p]. Mutations in TP53 are detected in about 10%-15% of patients with CLL at diagnosis, and 30 % on combination with del[17p] [83]. p53 deregulation alone or concomitant with del[17p] is reported to be a significant factor in the effective treatment option for these patients [89]. Additionally, the ATM components of the DNA damage pathway was observed to be deregulated in B-CLL cells and might underlie refractory disease [63, 89].
The extrinsic apoptosis is triggered by binding the tumour necrosis factor receptors (TNF-Rs) with members of TNF family ligands (TNF-α, FasL, and TRAIL) [74]. Computational modeling of the extrinsic apoptotic pathway revealed key regulators of TNF-mediated apoptosis. Particularly important role plays the positive feedback in the bistability dynamics of apoptotic signals (i.e. all-or-none behaviour) and irreversibility in caspase activation, as well as multi-factorial controlling type I versus II apoptosis pathways [74, 75, 105, 106, 109-112]. Sensitivity analysis of TNF-mediated apoptosis models aimed at identifying control points in apoptotic signalling [75, 103, 107, 109-111, 113]. The results of sensitivity analysis of the computational models of apoptosis induced by TRAIL [114] and FasL [75] are given in figure 4A and B, respectively.
The high level of sensitivity of the model output (caspase-3 activation) showed the following sensitive points of this pathway: 1. kinetic parameters of TRAIL ligand-TRAIL receptor interaction and concentrations of membrane receptors; 2. activation of pro-apoptotic protein Bak/Bax negatively regulated by Bcl-2 protein; and 3. Smac/XIAP inhibition of caspase-3 and caspase-9 resulting in positive feedback in apoptosis. The sensitivity of caspase-3 activation to Bcl-2 concentration was obtained at upregulation of Bcl-2 as opposed to insensitivity to Bcl-2 at downregulation of Bcl-2 when Bcl-2 does not inhibit the type II apoptotic pathway significantly (see the representation of the sensitivity coefficients at the decrease and increase of protein expression levels in figure 4B) [75]. This finding is in agreement with clinical data on the significant role of Bcl-2 in B-CLL [89]. The B-CLL cells, acquire survival benefit against apoptotic signals, and resistance to therapy by upregulation of Bcl-2 [87]. Overexpression of Bcl-2 is considered as one of the biomarkers of CLL and a promising therapeutic target in CLL therapy as alone and in combination. Inhibition of Bcl-2 protein family activates programmed cell death through initiation of the mitochondrial pathway of apoptosis. Now several therapeutic agents targeting Bcl-2 family of proteins are in clinical development. Oblimersen (Bcl-2 antisense oligonucleotide) is at phase III trial and small-molecular mimetics of the BH3 peptide domain of Bcl-2 (obatoclax, navitoclax) are being clinically evaluated in patients with repulsed/refractory CLL [89, 115]. The different sensitivity of mitochondrial apoptotic pathway towards down-regulation and up-regulation of Bcl-2 revealed in sensitivity analysis showed that Bcl-2 inhibitors can possess selectivity for CLL cells versus healthy cells with a normal expression of Bcl-2 [75].
High sensitivity of caspase-3 activation with the affinity of TRAIL ligand to TRAIL receptor and their concentrations obtained in sensitivity analysis (figure 4), correlates with a significant role of the death receptor in initiation of apoptosis in B-CLL cells [101]. It was observed that CLL cells manifest resistance to TRAIL-induced apoptosis due to a low surface expression of death receptors, or over-expression of inhibitory proteins such as FLIP, IAP2, and XIAP [101, 116]. Sensitising B-CLL cells to TRAIL-induced apoptosis is possible by up-regulating surface expression of death receptors or down-regulating the inhibitory proteins [101]. Targeting TRAIL receptors by anti-TRAIL-R1 antibody, mapatumumab also showed in vitro efficacy in different haematological malignancies [101]. Apart from pro-apoptotic signaling, the TNF-family receptors activate pro-survival signaling cascades, mainly through the NF-kB pathway [89, 117].
The decision making in cell survival/apoptosis is determined by the balance between these competing signals. Two ligands that belong to the TNF superfamily, BAFF (B-cell–activating factor), APRIL (a proliferation-inducing ligand), and their receptors, were reported to be key regulators of B-CLL cell survival [89]. Receptors of BAFF and APRIL are expressed on B-CLL cells, and when activated can promote B-CLL cell survival. The release of BAFF and APRIL through an autocrine mechanism is reported to cause the constitutive activation of both the canonical and noncanonical NF-κB pathways [117, 118]. Computational modelling of NF-κB pathway was carried out to elucidate specific mechanisms governing NF-κB-induced gene expression, as in part-mechanisms of oscillatory dynamics in this pathway [104, 113, 119-121]. Sensitivity analysis of NF-κB pathway models revealed that formation and function of the inhibitor IκB kinase complex (IKK), significantly impact on the NF-κB signalling: [104, 119, 120]. The inhibitors of IKKb that inhibit canonical NF-κB pathway are expected to be attractive targets for B-CLL therapy [117]. In in vitro testing of the action of IKKb inhibitor, UTC on B-CLL cells showed suppression of the protective effects of BAFF on B-CLL cell survival [117]. The other therapeutic approach to inactivate pro-survival NF-κB pathway in B-CLL cells is developing, and in part designing, a fusion protein targeting both BAFF and APRIL ligands [89, 122]. In addition, novel compounds with pro-apoptotic activity should also be considered as potential for potential treatment of B-CLL patients [123, 124].
8. Concluding remarks and perspectives
Building regulatory network model is based on establishing the dynamics of interaction networks showing regulatory pathways involved in gene regulation during the cell cycle and apoptosis. These network models only form a very small part in terms of calibrating of the models, this can be done using patient samples to produce data in terms of the quantitative measurement of key regulatory proteins and cytokine profiling. Once calibrated on an individual patient basis, the models would then enable monitoring of that patient to determine significant changes in the molecular phenotype and, in terms of therapy, permit intervention therapies to be designed using a broad range of bio-therapeutic drugs. A multi-target approach would be possible by identifying key regulators that over come drug resistance trigger apoptosis. Moreover by determining sensitivity nodes or fragile points to induce apoptosis in B-CLL cells, can be useful as a robust drug sensitivity reporter model. These models can result in achieving maximum benefit in treatment regimes for B-CLL patients. Due to the inherent complexity however, particularly in regard to the number of regulatory feedback loops, we believe that a successful approach can only be achieved using mathematical methods.
Modelling the cell cycle has been shown to be extremely robust in terms of the parameter variations. Experimental and clinical data support the hypothesis that computationally identified fragile interactions in the cell cycle represent promising targets for drug therapy [89]. Modulation of translational efficiency and the manipulation of programmed proteolysis are the key components of fragile mechanisms across all the models, and therefore these areas should be considered in model network development for prediction scores in B-CLL disease in order to reach an appropriate therapeutic decision in patient treatment. Once the model satisfies the parameters such as sensitivity, robustness and fragile point determination, this can be applied for prediction of selective drug targets and treatment outcome. It is known that B-CLL is a slowly progressing disease. This allows drug selectivity of single or combination of therapeutic agents to be predicted using a model and then drug combinations to be tested on B-CLL ex vivo culture obtained form patient. As well as being useful for developing new therapies, the described models might be useful in establishing the causes of B-CLL. For example it may be that an antigen, either foreign or perceived foreign is involved, and the models could be used to further investigate molecular mechanism of the disease. It is also possible that the B-CLL phenotype covers a number of disease sub-types and hence the modelling approach would enable some of these sub-types to be described, thus leading to improved diagnosis and personalised therapies.
Acknowledgments
The authors would like to acknowledge financial support from the Northwood Charitable Trust and the Scottish Informatics and Computer Science Alliance (SICSA).
References
- Rozman C, Montserrat E. Chronic lymphocytic leukemia, New Engl J Med 1995; 333: 1052–1057.
REFERENCE LINK - Bergsagel DE. The chronic leukemias: a review of disease manifestations and the aims of therapy. Can Med Assoc J 1967; 96: 1615-1620.
- Gunz FW, Fitzgerald PH, Adams A. An abnormal chromosome in chroninc lymphocytic leukaemia. Brit Med J 1962; 2: 1097-1099.
REFERENCE LINK - Dameshek W, Savitz HA, Arbor B. Chronic lymphocytic leukemia in twin brothers aged fifty-six. JAMA 1929; 92: 1348-1349.
REFERENCE LINK - Branda RF, Ackerman SK, Handwerger BS, Howe RB, Douglas SD. Lymphocyte studies in familial chronic lymphatic leukemia. Am J Med 1978; 64: 508-514.
REFERENCE LINK - Furbetta D, Solinas P. Hereditary chronic lymphatic leukemia. Proc 2nd Int Cong Hum Genet. 1963; 2:1078-1079.
- Conley CL, Misiti J, Laster AJ. Genetic factors predisposing to chronic lymphocytic leukemia and to autoimmune disease. Medicine 1980; 59:323-334.
REFERENCE LINK - Goldin LR, Pfeiffer RM, Li X, Hemminki K. Familial risk of lymphoproliferative tumors in families of patients with chronic lymphocytic leukemia: results from the Swedish Family-Cancer Database. Blood 2004; 104: 1850-1854.
REFERENCE LINK - Anaissie EJ, Kontoyiannis DP, O’Brien S, Kantarjian H, Robertson L, Lerner S et al. Infections in patients with chronic lymphocytic leukemia treated with fludarabine. Ann Intern Med 1998; 129: 559-566.
- Hensel M, Kornacker M, Yammeni S, Egerer G, Ho AD. Disease activity and pretreatment, rather than hypogammaglobuli- naemia, are major risk factors for infectious complications in patients with chronic lymphocytic leukaemia. Brit J Haematol 2003; 122: 600-606.
REFERENCE LINK - Perkins JG, Flynn JM, Howard RS, Byrd JC. Frequency and type of serious infections in fludarabine-refractory B-cell chronic lymphocytic leukemia and small lymphocytic lymphoma: implications for clinical trials in this patient population. Cancer 2002; 94: 2033-2039.
REFERENCE LINK - Francis S, Karanth M, Pratt G, Starczynski J, Hooper L, Fegan C et al. The effect of immunoglobulin VH gene mutation status and other prognostic factors on the incidence of major infections in patientswith chronic lymphocytic leukemia. Cancer 2006; 107: 1023–1033.
REFERENCE LINK - Rossi D, De Paoli L, Rossi FM, Cerri M, Deambrogi C, Rasi S et al. Early stage chronic lymphocytic leukaemia carrying unmutated IGHV genes is at risk of recurrent infections during watch and wait. Br J Haematol 2008; 141: 734-736.
REFERENCE LINK - Hamblin TJ, Oscier DG. Chronic lymphocytic leukaemia: the nature of the leukaemic cell. Blood Rev 1997; 11: 119-128.
REFERENCE LINK - MacFarlane M, Harper N, Snowden TR, Dyer MJ, Barnett GA, Pringle JH et al. Mechanisms of resistance to TRAIL-induced apoptosis in primary B cell chronic lymphocytic leukaemia. Oncogene 2002; 21: 6809-6818.
REFERENCE LINK - Cheson BD, Bennett JM, Grever M, Kay N, Keating MJ, O’Brien S et al. National Cancer Institute-sponsored Working Group guidelines for chronic lymphocytic leukemia: revised guidelines for diagnosis and treatment. Blood 1996; 87(12): 4990-4997.
- Todorov IT, Philipova RN, Zhelev NZ, Hadjiolov AA. Monoclonal antibody to a nucleolar antigen of human B-lymphoblastoid cells. Cell Biol Int Rep 1987; 11(3): 181-187.
REFERENCE LINK - Matutes E, Attygalle A, Wotherspoon A, Catovsky D. Diagnostic issues in chronic lymphocytic leukaemia (CLL). Best Pract Res Clin Haematol 2010; 23(1): 3-20.
REFERENCE LINK - Staudt LM. Molecular diagnosis of the hematologic cancers. New Eng J Med 2003; 348: 1777-1785.
REFERENCE LINK - Grumont RJ, Strasser A, Gerondakis S. B cell growth is controlled by phosphatidylinosotol 3-kinase dependent induction of Rel/NF-kappaB regulated c-myc transcription. Mol Cell 2002; 10(6): 1283-1294.
REFERENCE LINK - Yusuf I, Zhu X, Kharas MG, Chen J, Fruman DA. Optimal B-cell proliferation requires phosphoinositide 3-kinase-dependent inactivation of FOXO transcription factors. Blood 2004; 104(3): 784-787.
REFERENCE LINK - Klein U, Dalla-Favera R. New insights into the phenotype and cell derivation of B cell chronic lymphocytic leukemia. Curr Top Microbiol Immunol 2005; 294 31-49.
REFERENCE LINK - Haslinger C, Schweifer N, Stilgenbauer S, Döhner H, Lichter P, Kraut N et al. Microarray gene expression profiling of B-cell chronic lymphocytic leukemia subgroups defined by genomic aberrations and VH mutation status. J Clin Oncol 2004; 22: 3937-3949.
- Gottardi D, Alfarano A, De Leo AM, Stacchini A, Aragno M, Rigo A. et al. In leukaemic CD5 þ B cells the expression of BCL-2 gene family is shifted toward protection from apoptosis. Br J Haematol 1996; 94: 612-618.
REFERENCE LINK - Kitada S, Andersen J, Akar S, Zapata JM, Takayama S, Krajewski S et al. Expression of apoptosis-regulating proteins in chronic lymphocytic leukemia correlations with in vitro and in vivo chemoresponses. Blood 1998; 91: 3379-3389.
- Pepper C, Lin TT, Pratt G, Hewamana S, Brennan P, Hiller L et al. Mcl-1 expression has in vitro and in vivo significance in chronic lymphocytic leukemia and is associated with other poor prognostic markers. Blood 2008; 112: 3807-3817.
REFERENCE LINK - Veronese L, Tournilhac O, Verrelle P, Davi F, Dighiero G, Chautard E et al. Low MCL-1 mRNA expression correlates with prolonged survival in B-cell chronic lymphocytic leukemia. Leukemia 2008; 22:1291-1293.
REFERENCE LINK - Chen R, Plunkett W. Strategy to induce apoptosis and circumvent resistance in chronic lymphocytic leukaemia. Best Pract Res Clin Haematol 2010; 23: 155-166.
REFERENCE LINK - Muzio M, Apollonio B, Scielzo C, Frenquelli M, Vandoni I, Boussiotis V et al. Constitutive activation of distinct BCR-signaling pathways in a subset of CLL patients: a molecular signature of energy. Blood 2008; 112:188-195.
REFERENCE LINK - Moore VDG, Brown JR, Certo M, Love TM, Novina CD, Letai A. Chronic lymphocytic leukemia requires BCL2 to sequester prodeath BIM, explaining sensitivity to BCL2 antagonist ABT-737. J Clin Invest 2007; 117: 112-121.
REFERENCE LINK - Del GM, Letai A. Rational design of therapeutics targeting the BCL-2 family: are some cancer cells primed for death but waiting for a final push? Adv Exp Med Biol 2008; 615: 159-175.
REFERENCE LINK - Kalla C, Scheuermann MO, Kube I, Schlotter M, Mertens D, Döhner H et al. Analysis of 11q22-q23 deletion target genes in B-cell chronic lymphocytic leukaemia: evidence for a pathogenic role of NPAT, CUL5, and PPP2R1B. Eur J Cancer 2007; 43: 1328-1335.
- Tumilasci VF, Oliere S, Nguyen TL, Shamy A, Bell J, Hiscott J. Targeting the apoptotic pathway with BCL-2 inhibitors sensitizes primary chronic lymphocytic leukemia cells to vesicular stomatitis virus induced oncolysis. J Virol 2008; 82(17): 8487-8499.
REFERENCE LINK - Zhong F, Harr MW, Bultynck G, Monaco G, Parys JB, De Smedt H et al. Induction of Ca(2)+-driven apoptosis in chronic lymphocytic leukemia cells by peptide-mediated disruption of Bcl-2-IP3 receptor interaction. Blood 2011; 117: 2924-2934.
REFERENCE LINK - Burger JA, Ghia P, Rosenwald A, Caligaris-Cappio F. The microenvironment in mature B-cell malignancies: a target for new treatment strategies. Blood 2009; 114: 3367-3375.
REFERENCE LINK - Ghia P, Chiorazzi N, Stamatopoulos K. Microenvironmental influences in chronic lymphocytic leukaemia: the role of antigen stimulation. J Intern Med 2008; 264: 549-562.
REFERENCE LINK - Christopoulos P, Pfeifer D, Bartholome K, Follo M, Timmer J, Fisch P et al. Definition and characterization of the systemic T-cell dysregulation in untreated indolent B-cell lymphoma and very early CLL. Blood 2011; 117: 3836-3846.
REFERENCE LINK - Esquela-Kerscher A, Slack FJ. Oncomirs – microRNAs with a role in cancer. Nat Rev Cancer 2006; 6: 259-269.
REFERENCE LINK - Slack FJ, Weidhaas JB. MicroRNAs as a potential magic bullet in cancer. Future Oncol 2006; 2: 73-82.
REFERENCE LINK - Fulci V, Chiaretti S, Goldoni M, Azzalin G, Carucci N, Tavolaro S, et al. Quantitative technologies establish a novel microRNA profile of chronic lymphocytic leukemia. Blood 2007; 109: 4944-4951.
REFERENCE LINK - Visone R, Croce CM. MiRNAs and cancer. Am J Pathol 2009; 174: 1131-1138.
REFERENCE LINK - Ward BP, Tsongalis GJ, Kaur P. MicroRNAs in chronic lymphocytic leukemia. Exp Mol Pathol 2011; 90: 173-178.
REFERENCE LINK - Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M et al. miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc Natl Acad Sci USA 2005; 102: 13944-13949.
REFERENCE LINK - Zenz T, Mohr J, Eldering E, Kater AP, Bühler A, Kienle D et al. miR-34a as part of the resistance network in chronic lymphocytic leukemia. Blood 2009; 113: 3801-3808.
- Bommer GT, Gerin I, Feng Y, Kaczorowski AJ, Kuick R, Love RE et al. p53-mediated activation of miRNA34 candidate tumor suppressor genes. Curr Biol 2007; 17: 1298-1307.
REFERENCE LINK - Sun F, Fu H, Liu Q, Tie Y, Zhu J, Xing R et al. Downregulation of CCND1 and CDK6 by miR-34a induces cell cycle arrest. FEBS Lett 2008; 582:1564-1568.
REFERENCE LINK - Asslaber D, Pinon JD, Seyfried I, Desch P, Stocher M, Tinhofer I et al. MicroRNA-34aexpression correlates with MDM2 SNP309 polymorphism and treatment-free survival in chronic lymphocytic leukemia. Blood 2010; 115 4191-4197.
REFERENCE LINK - Pekarsky Y, Santanam U, Cimmino A, Palamarchuk A, Efanov A, Maximov V et al. Tcl1 expression in chronic lymphocytic leukemia is regulated by miR-29 and miR-181. Cancer Res 2006; 66 11590- 11593.
REFERENCE LINK - Ibrahim S, Keating M, Do KA, O’Brien S, Huh YO, Jilani I et al. CD38 expression as an important prognostic factor in B-cell chronic lymphocytic leukemia. Blood 2001; 98: 181–186.
REFERENCE LINK - Węsierska-Gądek J, Krame MP, The impact of CDK inhibition in human malignancies associated with pronounced defects in apoptosis: advantages of multi-targeting small molecules. Future Med Chem 2012; 4(4): 395-424.
REFERENCE LINK - McClue SJ, Blake D, Clarke R, Cowan A, Cummings L, Fischer PM et al. In vitro and in vivo antitumor properties of the cyclin dependent kinase inhibitor CYC202 (R-roscovitine). Int J Cancer 2002; 102(5): 463-468.
REFERENCE LINK - Whittaker S, Walton M, Kelland L, Garrett M, Zhelev N, Workman P. RB phosphorylation as a pharmacodynamic marker of roscovitine (CYC202) activity in vitro and in vivo. P Am Assoc Canc Res 2001; 42: 926.
- McClue S, Fischer PM, Blake D, Clarke R, Duff S, Krauss E et al. Studies on the mechanism of action of CYC202 (R-roscovitine). P Am Assoc Canc Res 2002; 43: 666.
- García-Escobar I, Sepúlveda J, Castellano D, Cortés-Funes H. Therapeutic management of chronic lymphocytic leukaemia State of the art and future perspectives. Crit Rev Oncol Hematol 2011; 80: 100-113.
REFERENCE LINK - Robak T. Recent progress in the management of chronic lymphocytic leukemia, Cancer Treat Rev 2007; 33: 710–728.
REFERENCE LINK - Ricci F, Tedeschi A, Morra E, Montillo M. Fludarabine in the treatment of chronic lymphocytic leukemia: a review. Ther Clin Risk Manag 2009; 5: 187-207.
- Norian LA, Kucaba TA, Earel JK, Knutson T, Vanoosten RL, Griffith TS. Synergistic Induction of Apoptosis in Primary B-CLL Cells after Treatment withRecombinant Tumor Necrosis Factor-Related Apoptosis-Inducing Ligand and Histone Deacetylase Inhibitors. J Oncol 2009; 2009: 408038.
- Niedermeier M, Hennessy BT, Knight ZA, Henneberg M, Hu J, Kurtova AV et al. Isoform- selective phosphoinositide3′-kinase inhibitors inhibit CXCR4 signaling and overcome stromal cell-mediated drug resistance in chronic lymphocytic leukemia: a novel therapeutic approach. Blood 2009; 113: 5549-5557.
REFERENCE LINK - Herman SE, Gordon AL, Wagner AJ, Heerema NA, Zhao W, Flynn JM et al. The phosphatidylinositol 3-kinase-delta inhibitor CAL-101 demonstrates promising preclinical activity in chronic lymphocytic leukemia by antagonizing intrinsic and extrinsic cellular survival signals. Blood 2010; 116: 2078-2088.
REFERENCE LINK - Pleyer L, Egle A, Hartmann TN, Greil R. Molecular and cellular mechanisms of CLL: novel therapeutic approaches. Nat Rev Clin Oncol 2009; 6: 405-418.
REFERENCE LINK - Hartmann TN, Pleyer L, Desch P, Egle A, Greil R. Novel therapeutics approaches to chronic lymphocytic leukemia based on recent biological insights. Discov Med 2009; 8: 157-164.
- Stankovic T, Weber P, Stewart G, Bedenham T, Murray J, Byrd PJ. et al. Inactivation of ataxia telangiectasia mutated gene in B-cell chronic lymphocytic leukaemia. Lancet 1999; 353: 26–29.
REFERENCE LINK - Schaffner C, Idler I, Stilgenbauer S, Dohner H, Lichter P. Mantle cell lymphoma is characterized by inactivation of the ATM gene. Proc Natl Acad Sci USA 2000; 97: 2773–2778.
REFERENCE LINK - Stankovic T, Hubank M, Cronin D, Stewart GS, Fletcher D, Bignell CR et al. Microarray analysis reveals that TP53- and ATM- mutant B-CLLs share a defect in activating pro- apoptotic responses after DNA damage but are distinguished by major differences in activating prosurvival responses. Blood 2004; 103: 291-300.
REFERENCE LINK - Kim YC, Jung YC, Chen J, Alhasan AH, Kaewsaard P, Zhang Y. Evidences showing wide presence of small genomic aberrations in chronic lymphocytic leukemia. BMC Res Notes 2010; 3: 341.
REFERENCE LINK - Faratian D, Clyde RG, Crawford JW, Harrison DJ. Systems pathology – taking molecular pathology into a new dimension. Nat Rev Clin Oncol 2009; 6: 455-464.
REFERENCE LINK - Dohner H, Stilgenbauer S, Benner A, Leupolt E, Kröber A, Bullinger L et al. Genomic aberrations and survival in chronic lymphocytic leukemia. New Eng J Med 2000; 343: 1910-1916.
- Schrattenholz A, Groebe K, Soskic V. Systems biology approaches and tools for analysis of interactomes and multi-target drugs. Methods Mol Biol 2010; 662: 29-58.
REFERENCE LINK - Zhang GJ, Safran M, Wei W, Sorensen E, Lassota P, Zhelev N et al. Bioluminescent imaging of Cdk2 inhibition in vivo. Nat Med 2004; 10: 643-648.
REFERENCE LINK - Novak M, Tyson JJ. A model for restriction point control of the mammalian cell cycle. J Theor Biol 2004; 230: 563–579.
REFERENCE LINK - Qu Z, Weiss JN, MacLellan WR. Regulation of the mammalian cell cycle: a model of the G1-to-S transition. Am J Physiol-Cell Physiol 2003; 284: C349-C364
- Aguda BD. A qualitative analysis of the kinetics of the G2 DNA damage checkpoint system. Proc Natl Acad Sci USA 1999; 96: 11352–11357.
REFERENCE LINK - Clyde RG, Tummala H, Khalil HS, Goszcz K, Lucka I, Tupone MG et al. A novel quantitative systems biology approach to cancer research and treatment. Curr Opin Biotech 2011; 22(S1): S58.
REFERENCE LINK - Spencer SL, Sorger PK. Measuring and modeling apoptosis in single cells. Cell 2011; 144: 926-939.
REFERENCE LINK - Hua F, Cornejo MG, Cardone MH, Stokes CL, Lauffenburger DA. Effects of Bcl-2 Levels on Fas Signaling-Induced Caspase-3 Activation: Molecular Genetic Tests of Computational Model Predictions. J Immunol 2005; 175: 985-995.
- Shoemaker JE, Doyle III FJ. Identifying Fragilities in Biochemical Networks: Robust Performance Analysis of Fas Signaling-Induced Apoptosis. Biophys J 2008; 95: 2610–2623.
REFERENCE LINK - Birtwistle MR, Hatakeyama M, Yumoto N, Ogunnaike BA, Hoek JB, Kholodenko BN. Ligand-dependent responses of the ErbB signaling network: experimental and modeling analyses. Mol Syst Biol 2007; 3: 144.
REFERENCE LINK - Schoeberl B, Pace EA, Fitzgerald JB, Harms BD, Xu L, Nie L et al. Therapeutically Targeting ErbB3: A Key Node in Ligand-Induced Activation of the ErbB Receptor–PI3K Axis. Sci Signal 2009; 2(77):) ra31.
- Schoeberl B, Faber AC, Li D, Liang MC, Crosby K, M. Onsum et al. An ErbB3 antibody, MM-121, is active in cancers with ligand-dependent activation. Cancer Res 2010; 70: 2485-2494.
REFERENCE LINK - Kitano H. Cancer as a robust system: implications for anticancer therapy. Nat Rev Cancer 2004; 4: 227-235.
REFERENCE LINK - Kitano H. A robustness based approach to systems-oriented drug design. Nat Rev Drug Discovery 2007; 6: 202-210.
REFERENCE LINK - Nayak S, Salim S, Luan D, Zai M, Varner JD. A test of highly optimized tolerance reveals fragile cell-cycle mechanisms are molecular targets in clinical cancer trials. PLoS One 2008; 3: e2016.
REFERENCE LINK - Novak B, Tyson J. Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos. J Cell Sci 1993; 106: 1153-1168.
- Tyson J, Chen K, Novak B. Network dynamics and cell physiology. Nat Rev Mol Cell Biol 2001; 2: 908-916.
REFERENCE LINK - Clyde RG, Bown JL, Hupp TR, Zhelev N, Crawford JW. The role of modelling in identifying drug targets for diseases of the cell cycle. J R Soc Interface 2006; 3: 617-627.
REFERENCE LINK - Shiraishi T, Matsuyama S, Kitano H. Large-scale analysis of network bistability for human cancers. PLoS Comput Biol 2010; 6: e1000851.
REFERENCE LINK - Badoux XC, Keating MJ, Wierda WG. What is the best frontline therapy for patients with CLL and 17p deletion? Curr Hematol Malig Rep 2011; 6: 36-46.
REFERENCE LINK - Varterasian ML, Mohammad RM, Eilender DS, Hulburd K, Rodriguez DH, P.A. Pemberton PA et al. Phase I study of bryostatin 1 in patients with relapsed non-Hodgkin’s lymphoma and chronic lymphocytic leukemia. J Clin Oncol 1998; 16: 56–62.
- Wierda WG, Chiorazzi N, Dearden C, Brown JR, Montserrat E, Shpall E et al. Chronic lymphocytic leukemia: new concepts for future therapy. Clin Lymphoma Myeloma Leuk 2010; 10:369-378.
REFERENCE LINK - Kholodenko BN, Demin OV, Moehren G, Hoek JB. Quantification of short term signaling by the epidermal growth factor receptor. J Biol Chem 1999; 274: 30169-30181.
REFERENCE LINK - Chen WW, Schoeberl B, Jasper PJ, Niepel M, Nielsen UB, Lauffenburger DA et al. Input–output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic data. Mol Syst Biol 2009; 5 239.
REFERENCE LINK - Wang CC, Cirit M, Haugh JM. PI3K-dependent cross-talk interactions converge with Ras as quantifiable inputs integrated by Erk. Mol Syst Biol 2009; 5 246.
REFERENCE LINK - Kreeger PK, Lauffenburger DA. Cancer systems biology: a network modeling perspective. Carcinogenesis 2010; 31: 2-8.
REFERENCE LINK - Orton RJ, Adriaens ME, Gormand A, Sturm OE, Kolch W, Gilbert DR, Computational modelling of cancerous mutations in the EGFR/ERK signalling pathway. BMC Syst Biol 2009; 3 100.
REFERENCE LINK - Faratian D, Goltsov A, Lebedeva G, Moodie S, Mullen P, Kay C et al. Systems biology reveals new strategies for personalising cancer medicine and confirms PTENs role in resistance to trastuzumab. Cancer Res 2009; 69(16): 6713-6720.
REFERENCE LINK - Lebedeva G, Sorokin A, Faratian D, Mullen P, Goltsov A, Langdon SP et al. Model-based global sensitivity analysis as applied to identification of anti-cancer drug targets and biomarkers of drug resistance in the ErbB2/3 network. Eur J Pharm Sci 2011; 46(4): 244-258.
REFERENCE LINK - Goltsov A, Faratian D, Langdon SP, Harrison DJ, Bown J. Features of the reversible sensitivity-resistance transition in ERK/PI3K/PTEN/AKT signalling network at HER2 inhibition. Cell Signal 2012; 24: 493-504.
REFERENCE LINK - Yuan TL, Cantley LC. PI3K pathway alterations in cancer: variations on a theme, Oncogene 2008; 27: 5497-5510.
REFERENCE LINK - Martelli AM, Evangelisti C, Chiarini F, Grimaldi C, Cappellini A, Ognibene A et al. The emerging role of the phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin signaling network in normal myelopoiesis and leukemogenesis. Biochim Biophys Acta 2010; 1803: 991-1002.
- Niedermeier M, Hennessy BT, Knight ZA, Henneberg M, Hu J, Kurtova AV et al. Burger, Isoform-selective hosphoinositide 3′-kinase inhibitors inhibit CXCR4 signaling and overcome stromal cellmediated drug resistance in chronic lymphocytic leukemia: a novel therapeutic approach. Blood 2009; 113: 5549-5557.
REFERENCE LINK - Masood A, Sher T, Paulus A, Miller KC, Chitta S, Chanan-Khan A. Targeted treatment for chronic lymphocytic leukemia. Onco Targets Ther 2011; 4: 169-183.
- De Frias M, Iglesias-Serret D, Cosialls AM, Coll-Mulet L, Santidrián AF, González-Gironès DM et al. Akt inhibitors induce apoptosis in chronic lymphocytic leukemia cells. Haematologica 2009; 94: 1698-1707.
- Bentele M, Lavrik I, Ulrich M, Stösser S, Heermann DW, Kalthoff H, Mathematical modeling reveals threshold mechanism in CD95-induced apoptosis. J Cell Biol 2004; 166: 839-851.
REFERENCE LINK - Toettcher JE, Loewer A, Ostheimer GJ, Yaffe MB, Tidor B, Lahav G. Distinct mechanisms act in concert to mediate cell cycle arrest. Proc Natl Acad Sci USA 2009; 106: 785-790.
REFERENCE LINK - Aldridge BB, Gaudet S, Lauffenburger DA, Sorger PK. Lyapunov exponents and phase diagrams reveal multi-factorial control over TRAIL-induced apoptosis, Mol Syst Biol 2011; 7: 553.
REFERENCE LINK - Komorowski M, Costa MJ, Rand DA, Stumpf MP. Sensitivity, robustness, and identifiability in stochastic chemical kinetics models. Proc Natl Acad Sci USA 2011; 108: 8645-8650.
REFERENCE LINK - Koh G, Lee DY. Mathematical modeling and sensitivity analysis of the integrated TNFα-mediated apoptotic pathway for identifying key regulators. Comput Biol Med 2011; 41: 512-528.
REFERENCE LINK - Sun T, Yang W, Liu J, Shen P. Modeling the basal dynamics of p53 system. PLoS One 2011; 6: e27882.
REFERENCE LINK - Harrington HA, Ho KL, Ghosh S, Tung KC. Construction and analysis of a modular model of caspase activation in apoptosis. Theor Biol Med Model 2008; 5: 26.
REFERENCE LINK - F Hua , S. Hautaniemi, R. Yokoo, D.A. Lauffenburger, Integrated mechanistic and data-driven modelling for multivariate analysis of signalling pathways. J R Soc Interface 2006; 3: 515-526.
REFERENCE LINK - Perumal TM, Gunawan R. Understanding dynamics using sensitivity analysis: caveat and solution. BMC Syst Biol 2011; 5 41.
REFERENCE LINK - Legewie S, Bluthgen N, Herzel H. Mathematical modeling identifies inhibitors of apoptosis as mediators of positive feedback and bistability. PLoS Comput Biol 2006; 2: e120.
REFERENCE LINK - Cho KH, Shin SY, Kolch W, Wolkenhauer O. Experimental design in systems biology, based on parameter sensitivity analysis using Maonte Carlo method: a case study for the TNF-mediated NF-B signal transduction pathway. Simulation 2003; 79: 726.
REFERENCE LINK - Zhang T, Wu M, Chen Q, Sun Z. Investigation into the regulation mechanisms of TRAIL apoptosis pathway by mathematical modelling. Acta Biochim. Biophys Sin (Shanghai) 2010; 42: 98-108.
REFERENCE LINK - Roberts AW, Seymour JF, Brown JR, Wierda WG, Kipps TJ, Khaw SL. Substantial susceptibility of chronic lymphocytic leukemia to BCL2 inhibition: results of a phase I study of navitoclax in patients with relapsed or refractory disease. J Clin Oncol 2012; 30: 488-496.
REFERENCE LINK - MacFarlane M, Harper N, Snowden RT, Dyer MJ, Barnett GA, Pringle JH et al. Mechanisms of resistance to TRAIL-induced apoptosis in primary B cell chronic lymphocytic leukaemia. Oncogene 2002; 21: 6809-6818.
REFERENCE LINK - Endo T, Nishio M, Enzler T, Cottam HB, Fukuda T, James DF et al. BAFF and APRIL support chronic lymphocytic leukemia B-cell survival through activation of the canonical NF-kB pathway. Blood 2007; 109: 703-710.
REFERENCE LINK - M. Nishikori, Classical and alternative NF-kB activation pathways and their roles in lymphoid malignancies J. Clin. Exp. Hematopathol. 45 (2005) 15-24.
REFERENCE LINK - Ihekwaba AE, Broomhead DS, Grimley RL, Benson N, Kell DB. Sensitivity analysis of parameters controlling oscillatory signalling in the NF-kappaB pathway: the roles of IKK and IkappaBalpha. Syst Biol 2001; 1(1): 93-103.
REFERENCE LINK - Nikolov S, Vera J, Rath O, Kolch W, Wolkenhauer O. Role of inhibitory proteins as modulators of oscillations in NFkB signalling. IET Syst Biol 2009; 3: 59-76.
REFERENCE LINK - Tay S, Hughey JJ, Lee TK, Lipniacki T, Quake SR, Covert MW. Single-cell NF-kB dynamics reveal digital activation and analogue information processing, Nature 2010; 466: 267-271.
REFERENCE LINK - Kofler DM, Gawlik BB, Elter T, Gianella-Borradori A, Wendtner CM, Hallek M. Phase 1b trial of atacicept a recombinant protein binding BLyS and APRIL, in patients with chronic lymphocytic leukemia. Leukemia 2012; 26(4) 841-844.
REFERENCE LINK - Neychev VK, Nikolova E, Zhelev N, Mitev VI. Saponins from Tribulus terrestris L are less toxic for normal human fibroblasts than for many cancer lines: influence on apoptosis and proliferation. Exp Biol Med (Maywood) 2007; 232(1): 126-133.
- Sarek J, Klinot J, Dzubák P, Klinotová E, Nosková V, Krecek V, New lupane derived compounds with pro-apoptotic activity in cancer cells: synthesis and structure-activity relationships. J Med Chem 2003; 46(25): 5402-5415.
REFERENCE LINK