S100A4 and its role in metastasis – simulations of knockout and amplification of epithelial growth factor receptor and matrix metalloproteinases

Interplay of mutations, alternate mechanisms, and treatment breaks in leukaemia: Understanding and implications studied with stochastic models

Bcr-Abl1 kinase domain mutations are the most prevalent cause of treatment resistance in chronic myeloid leukaemia (CML). Alternate resistance pathways nevertheless exist, and cell line experiments show certain patterns in the gain, and loss, of some of these alternate adaptations. These adaptations have clinical consequences when the tumour develops mechanisms that are beneficial to its growth under treatment, but slow down its growth when not treated. The results of temporarily halting treatment in CML have not been widely discussed in the clinic and there is no robust theoretical model that could suggest when such a pause in therapy can be tolerated. We constructed a dynamic model of how mechanisms such as Bcr-Abl1 overexpression and drug transporter upregulation evolve to produce resistance in cell lines, and investigate its behaviour subject to different treatment schedules, in particular when the treatment is paused ('drug holiday'). Our study results suggest that the presence of additional resistance mechanisms creates an environment which favours mutations that are either preexisting or occur late during treatment. Importantly, the results suggest the existence of tumour drug addiction, where cancer cells become dependent on the drug for (optimal) survival, which could be exploited through a treatment holiday. All simulation code is available at https://github.com/Sandalmoth/dual-adaptation.

The molecular mechanisms behind activation of FLT3 in acute myeloid leukemia and resistance to therapy by selective inhibitors

Acute myeloid leukemia is an aggressive cancer, which, in spite of increasingly better understanding of its genetic background remains difficult to treat. Mutations in the FLT3 gene are observed in ≈30% of the patients. Most of these mutations are internal tandem duplications (ITDs) of a sequence within the protein coding region, an activation mechanism that is almost non-existent with other genes and cancers. As patients each carry their own unique set of mutations, it is challenging to understand how ITDs activate the protein, and ascertain the risk for each individual patient. Available treatment options are limited due to development of drug resistance. Here, recent studies are reviewed that help to better understand the molecular mechanism behind activation of the FLT3 protein due to mutations. It is argued that difference in mutation sequences and especially location might be coupled to prognosis. When it comes to FLT3 inhibitors, key differences between them can be attributed to the mode of inhibition (type-1 and type-2 inhibitors), effective inhibitory coefficient in the blood plasma and off-target binding. Accounting for the position and length of insertions may in the future be used to predict prognosis and rationalise treatment. Development of new inhibitors must take into account the potential for resistance mutations. Inhibitors aimed at multiple specific targets are currently being developed. These, and as well as combination therapies will hopefully lead to longer periods during which targeted FLT3 therapy will remain effective.

Dexamethasone Induces Changes in Osteogenic Differentiation of Human Mesenchymal Stromal Cells via SOX9 and PPARG, but Not RUNX2

Despite the huge body of research on osteogenic differentiation and bone tissue engineering, the translation potential of in vitro results still does not match the effort employed. One reason might be that the protocols used for in vitro research have inherent pitfalls. The synthetic glucocorticoid dexamethasone is commonly used in protocols for trilineage differentiation of human bone marrow mesenchymal stromal cells (hBMSCs). However, in the case of osteogenic commitment, dexamethasone has the main pitfall of inhibiting terminal osteoblast differentiation, and its pro-adipogenic effect is well known. In this work, we aimed to clarify the role of dexamethasone in the osteogenesis of hBMSCs, with a particular focus on off-target differentiation. The results showed that dexamethasone does induce osteogenic differentiation by inhibiting SOX9 expression, but not directly through RUNX2 upregulation as it is commonly thought. Rather, PPARG is concomitantly and strongly upregulated, leading to the formation of adipocyte-like cells within osteogenic cultures. Limiting the exposure to dexamethasone to the first week of differentiation did not affect the mineralization potential. Gene expression levels of RUNX2, SOX9, and PPARG were simulated using approximate Bayesian computation based on a simplified theoretical model, which was able to reproduce the observed experimental trends but with a different range of responses, indicating that other factors should be integrated to fully understand how dexamethasone influences cell fate. In summary, this work provides evidence that current in vitro differentiation protocols based on dexamethasone do not represent a good model, and further research is warranted in this field.

Combating drug resistance in acute myeloid leukaemia by drug rotations: the effects of quizartinib and pexidartinib

BACKGROUND

Acute myeloid leukaemia (AML) is an aggressive blood cancer. In approximately 30% of the cases, driver mutations in the FLT3 gene are identified. FLT3 inhibitors are used in treatment of such patients together with cytotoxic drugs or (in refractory AML) as single agents. Unfortunately, resistance to FLT3 inhibitors limits their efficacy. Resistance is often due to secondary mutations in the gene encoding the molecular target. The gatekeeper mutation F691L confers resistance to specific FLT3 inhibitors such as quizartinib, but pexidartinib is much less resistance to this mutation. Pexidartinib alone is however sensitive to many other resistance mutations. In chronic myeloid leukaemia (CML), it has been suggested that rotation between drugs with a different landscape of resistance mutations might postpone the emergence of resistance.

METHODS

We studied the effect of quizartinib and pexidartinib in AML cell lines that express FLT3 (MOLM-14 and MV4-11). Using a rotation protocol, we further examined whether the emergence of resistance could be postponed. Computational modelling was used to analyse the onset of resistance and suggest which mutations are most likely to occur in a quantitative fashion.

RESULTS

The cells were sensitive to both inhibitors but quickly developed resistance that could be inherited, suggesting a genetic origin. Rotation protocols were not useful to postpone the emergence of resistance, which implies that such protocols, or changing from pexidartinib to quizartinib (or vice-versa) should not be used in patients. The computational modelling led to similar conclusions and suggested that F691L is the most common mutation to occur with quizartinib, and also when both drugs are used in rotation.

CONCLUSIONS

AML patients are not likely to benefit from a quizartinib/pexidartinib rotation protocol. A combination of tyrosine kinase inhibitors (with different molecular targets) might be more useful in the future. Development of specific FLT3 inhibitors that are less sensitive to resistance mutations might also lead to a better outcome.

Inhibition of the invasion and metastasis of mammary carcinoma cells by NBD peptide targeting S100A4 via the suppression of the Sp1/MMP‑14 axis

A synthetic peptide that blocks the interaction between the metastasis‑enhancing calcium‑binding protein, S100A4, and its effector protein, methionine aminopeptidase 2 (MetAP2) (the NBD peptide), was previously demonstrated to inhibit the angiogenesis of endothelial cells, leading to the regression of human prostate cancer in a xenograft model. However, the effects of the NBD peptide on the malignant properties of cancer cells that express S100A4 remain to be elucidated. The present study demonstrates that the NBD peptide inhibits the invasiveness and metastasis of highly metastatic human mammary carcinoma cells. The introduction of the peptide into MDA‑MB‑231 variant cells resulted in the suppression of matrix degradation in a gelatin invadopodia assay and invasiveness in a Matrigel invasion assay. In line with these results, the peptide significantly downregulated the expression of matrix metalloproteinase (MMP)‑14 (MT1‑MMP). Mechanistic analysis of the downregulation of MMP‑14 revealed the suppression of the expression of the transcription factor, specificity protein 1 (Sp1), but not that of nuclear factor (NF)‑κB, early growth response 1 (EGR1) or ELK3, all of which were reported to be involved in transcriptional regulation of the MMP‑14 gene. At the same time, evidence suggested that the NBD peptide also suppressed Sp1 and MMP‑14 expression levels in MDA‑MB‑468 cells. Importantly, the intravenous administration of the NBD peptide encapsulated in liposomes inhibited pulmonary metastasis from mammary gland tumors in mice with xenograft tumors. These results indicate that the NBD peptide can suppress malignant tumor growth through the suppression of the Sp1/MMP‑14 axis. Taken together, these results reveal that the NBD peptide acts on not only endothelial cells, but also on tumor cells in an integrated manner, suggesting that the peptide may prove to be a promising cancer therapeutic peptide drug.

The effects of combination treatments on drug resistance in chronic myeloid leukaemia: an evaluation of the tyrosine kinase inhibitors axitinib and asciminib

BACKGROUND

Chronic myeloid leukaemia is in principle a treatable malignancy but drug resistance is lowering survival. Recent drug discoveries have opened up new options for drug combinations, which is a concept used in other areas for preventing drug resistance. Two of these are (I) Axitinib, which inhibits the T315I mutation of BCR-ABL1, a main source of drug resistance, and (II) Asciminib, which has been developed as an allosteric BCR-ABL1 inhibitor, targeting an entirely different binding site, and as such does not compete for binding with other drugs. These drugs offer new treatment options.

METHODS

We measured the proliferation of KCL-22 cells exposed to imatinib-dasatinib, imatinib-asciminib and dasatinib-asciminib combinations and calculated combination index graphs for each case. Moreover, using the median-effect equation we calculated how much axitinib can reduce the growth advantage of T315I mutant clones in combination with available drugs. In addition, we calculated how much the total drug burden could be reduced by combinations using asciminib and other drugs, and evaluated which mutations such combinations might be sensitive to.

RESULTS

Asciminib had synergistic interactions with imatinib or dasatinib in KCL-22 cells at high degrees of inhibition. Interestingly, some antagonism between asciminib and the other drugs was present at lower degrees on inhibition. Simulations revealed that asciminib may allow for dose reductions, and its complementary resistance profile could reduce the risk of mutation based resistance. Axitinib, however, had only a minor effect on T315I growth advantage.

CONCLUSIONS

Given how asciminib combinations were synergistic in vitro, our modelling suggests that drug combinations involving asciminib should allow for lower total drug doses, and may result in a reduced spectrum of observed resistance mutations. On the other hand, a combination involving axitinib was not shown to be useful in countering drug resistance.

Reverse engineering directed gene regulatory networks from transcriptomics and proteomics data of biomining bacterial communities with approximate Bayesian computation and steady-state signalling simulations

BACKGROUND

Network inference is an important aim of systems biology. It enables the transformation of OMICs datasets into biological knowledge. It consists of reverse engineering gene regulatory networks from OMICs data, such as RNAseq or mass spectrometry-based proteomics data, through computational methods. This approach allows to identify signalling pathways involved in specific biological functions. The ability to infer causality in gene regulatory networks, in addition to correlation, is crucial for several modelling approaches and allows targeted control in biotechnology applications.

METHODS

We performed simulations according to the approximate Bayesian computation method, where the core model consisted of a steady-state simulation algorithm used to study gene regulatory networks in systems for which a limited level of details is available. The simulations outcome was compared to experimentally measured transcriptomics and proteomics data through approximate Bayesian computation.

RESULTS

The structure of small gene regulatory networks responsible for the regulation of biological functions involved in biomining were inferred from multi OMICs data of mixed bacterial cultures. Several causal inter- and intraspecies interactions were inferred between genes coding for proteins involved in the biomining process, such as heavy metal transport, DNA damage, replication and repair, and membrane biogenesis. The method also provided indications for the role of several uncharacterized proteins by the inferred connection in their network context.

CONCLUSIONS

The combination of fast algorithms with high-performance computing allowed the simulation of a multitude of gene regulatory networks and their comparison to experimentally measured OMICs data through approximate Bayesian computation, enabling the probabilistic inference of causality in gene regulatory networks of a multispecies bacterial system involved in biomining without need of single-cell or multiple perturbation experiments. This information can be used to influence biological functions and control specific processes in biotechnology applications.

Stochastic modelling of tyrosine kinase inhibitor rotation therapy in chronic myeloid leukaemia

Background

Resistance towards targeted cancer treatments caused by single nucleotide variations is a major issue in many malignancies. Currently, there are a number of available drugs for chronic myeloid leukaemia (CML), which are overcome by different sets of mutations. The main aim of this study was to explore if it can be possible to exploit this and create a treatment protocol that outperforms each drug on its own.

Methods

We present a computer program to test different treatment protocols against CML, based on available resistance mutation growth data. The evolution of a relatively stable pool of cancer stem cells is modelled as a stochastic process, with the growth of cells expressing a tumourigenic protein (here, Abl1) and any emerging mutants determined principally by the drugs used in the therapy.

Results

There can be some benefit to Bosutinib-Ponatinib rotation therapy even if the mutation status is unknown, whereas Imatinib-Nilotinib rotation is unlikely to improve the outcomes. Furthermore, an interplay between growth inhibition and selection effects generates a non-linear relationship between drug doses and the risk of developing resistance.

Conclusions

Drug rotation therapy might be able to delay the onset of resistance in CML patients without costly ongoing observation of mutation status. Moreover, the simulations give credence to the suggestion that lower drug concentrations may achieve better results following major molecular response in CML.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5690-5) contains supplementary material, which is available to authorized users.

New Insights into the Occurrence of Matrix Metalloproteases -2 and -9 in a Cohort of Breast Cancer Patients and Proteomic Correlations

Matrix metalloproteases (MMPs) are a family of well-known enzymes which operate prevalently in the extracellular domain, where they fulfil the function of remodeling the extracellular matrix (ECM). Within the 26 family members, encoded by 24 genes in humans, MMP-2 and MMP-9 have been regarded as primarily responsible for the basement membrane and peri-cellular ECM rearrangement. In cases of infiltrating carcinomas, which arise from the epithelial tissues of a gland or of an internal organ, a marked alteration of the expression and the activity levels of both MMPs is known to occur. The present investigation represents the continuation and upgrading of our previous studies, now focusing on the occurrence and intensity levels of MMP-2 and -9 and their proteomic correlations in a cohort of 80 breast cancer surgical tissues.

A computational study of hedgehog signalling involved in basal cell carcinoma reveals the potential and limitation of combination therapy

BACKGROUND

The smoothened (SMO) receptor is an essential component of the Sonic hedgehog (SHH) signalling, which is associated with the development of skin basal cell carcinoma (BCC). SMO inhibitors are indicated for BCC patients when surgical treatment or radiation therapy are not possible. Unfortunately, SMO inhibitors are not always well tolerated due to severe side effects, and their therapeutical success is limited by resistance mutations.

METHODS

We investigated how common are resistance-causing mutations in two genomic databases which are not linked to BCC or other cancers, namely 1000 Genomes and ExAC. To examine the potential for combination therapy or other treatments, we further performed knowledge-based simulations of SHH signalling, in the presence or absence of SMO and PI3K/Akt inhibitors.

RESULTS

The database analysis revealed that of 18 known mutations associated with Vismodegib-resistance, three were identified in the databases. Treatment of individuals carrying such mutations is thus liable to fail a priori. Analysis of the simulations suggested that a combined inhibition of SMO and the PI3K/Akt signalling pathway may provide an effective reduction in tumour proliferation. However, the inhibition dosage of SMO and PI3K/Akt depended on the activity of phosphodiesterases (PDEs). Under high PDEs activities, SMO became the most important control node of the network. By applying PDEs inhibition, the control potential of SMO decreased and PI3K appeared as a significant factor in controlling tumour proliferation.

CONCLUSIONS

Our systems biology approach employs knowledge-based computer simulations to help interpret the large amount of data available in public databases, and provides application-oriented solutions for improved cancer resistance treatments.

Computer simulations of the signalling network in FLT3 +-acute myeloid leukaemia - indications for an optimal dosage of inhibitors against FLT3 and CDK6

Background

Mutations in the FMS-like tyrosine kinase 3 (FLT3) are associated with uncontrolled cellular functions that contribute to the development of acute myeloid leukaemia (AML). We performed computer simulations of the FLT3-dependent signalling network in order to study the pathways that are involved in AML development and resistance to targeted therapies.

Results

Analysis of the simulations revealed the presence of alternative pathways through phosphoinositide 3 kinase (PI3K) and SH2-containing sequence proteins (SHC), that could overcome inhibition of FLT3. Inhibition of cyclin dependent kinase 6 (CDK6), a related molecular target, was also tested in the simulation but was not found to yield sufficient benefits alone.

Conclusions

The PI3K pathway provided a basis for resistance to treatments. Alternative signalling pathways could not, however, restore cancer growth signals (proliferation and loss of apoptosis) to the same levels as prior to treatment, which may explain why FLT3 resistance mutations are the most common resistance mechanism. Finally, sensitivity analysis suggested the existence of optimal doses of FLT3 and CDK6 inhibitors in terms of efficacy and toxicity.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2145-y) contains supplementary material, which is available to authorized users.

Increased S100A15 expression and decreased DNA methylation of its gene promoter are involved in high metastasis potential and poor outcome of lung adenocarcinoma

Purpose

This study aims to determine the functional role of S100A15 and its promoter DNA methylation patterns in lung cancer progression.

Experimental Design

We analyzed 178 formalin-fixed paraffin embedded specimens from lung cancer patients, including 24 early stage and 91 advanced stage adenocarcinoma. S100A15 protein expression was evaluated by immunohistochemistry stain, and its DNA methylation levels were measured by pyrosequencing.

Results

S100A15 nuclear staining was increased in lung adenocarcinoma patients with distant metastasis versus those without distant metastasis. There was reduced one/three-year overall survival in adenocarcinoma patients receiving first line target therapy and harboring high nuclear expressions of S100A15. Both DNA methylation levels over -423 and -248 CpG sites of the S100A15 gene promoter were decreased in adenocarcinoma patients with distant metastasis, and the former was associated with lower one-year overall survival. The highly invasive CL1-5 cell lines display decreased DNA methylation over −412/−248/−56 CpG sites of the S100A15 gene promoter and increased S100A15 gene/protein expressions as compared with the less invasive CL1-0 cell lines. Knockdown of S100A15 in CL1-5 cell line inhibited cell proliferation, migration, and invasion, while over-expression of S100A15 in CL1-0 cell line promoted cell proliferation, migration, and invasion. RNA sequencing analysis revealed potential biological effects of S100A15 over-expression and knock-down with CTNNB1, ZEB1, CDC42, HSP90AA1, BST2, and PCNA being the pivotal down-stream mediators.

Conclusions

Increased S100A15 expression and decreased DNA methylation of its gene promoter region were associated with high metastasis potential and poor outcome in lung adenocarcinoma, probably through triggering CTNNB1 -centered pathways.

Drug resistance in cancer: molecular evolution and compensatory proliferation

Targeted therapies have revolutionized cancer treatment. Unfortunately, their success is limited due to the development of drug resistance within the tumor, which is an evolutionary process. Understanding how drug resistance evolves is a prerequisite to a better success of targeted therapies. Resistance is usually explained as a response to evolutionary pressure imposed by treatment. Thus, evolutionary understanding can and should be used in the design and treatment of cancer. In this article, drug-resistance to targeted therapies is reviewed from an evolutionary standpoint. The concept of apoptosis-induced compensatory proliferation (AICP) is developed. It is shown that AICP helps to explain some of the phenomena that are observed experimentally in cancers. Finally, potential drug targets are suggested in light of AICP.

A DNA-directed covalent conjugation fluorescence probe for in vitro detection of functional matrix metalloproteinases

Matrix metalloproteinases (MMPs) have been considered to contribute to the progression of tumorigenesis and tumor invasion; MMP-9 in particular, has been regarded as a priority target in cancer treatment due to its massive up-regulation in malignant tissues and its ability to degrade type IV collagen. In this work, we employed a DNA-directed covalent conjugation method to design a fluorescence probe for in vitro detection of functional matrix metalloproteinases, by which a nitrilotriacetic acid (NTA)-modified DNA probe can combine with the Zn²⁺ in the active site of MMPs, and then a molecule beacon (MB) modified FITC and BHQ1 can open to bond with their complementary base, NTA-modified DNA. We can evaluate the amount of MMPs in the medium according to the fluorescence intensity. The detection procedure can be finished in 30 min with good selectivity, cheap reagents and easy preparation. All the results and the amount of secreted MMPs under three different cell culture conditions are in accordance with previous reports. Satisfactory results are obtained. Furthermore, owing to the importance of MMP-9, we designed an approach to achieve the desired selectivity and specificity of our work, using dual amplification for improving fluorescence intensity based on RCA to detect the amount of MMP-9.

Sensitivity Analysis of the NPM-ALK Signalling Network Reveals Important Pathways for Anaplastic Large Cell Lymphoma Combination Therapy

A large subset of anaplastic large cell lymphoma (ALCL) patients harbour a somatic aberration in which anaplastic lymphoma kinase (ALK) is fused to nucleophosmin (NPM) resulting in a constitutively active signalling fusion protein, NPM-ALK. We computationally simulated the signalling network which mediates pathological cell survival and proliferation through NPM-ALK to identify therapeutically targetable nodes through which it may be possible to regain control of the tumourigenic process. The simulations reveal the predominant role of the VAV1-CDC42 (cell division control protein 42) pathway in NPM-ALK-driven cellular proliferation and of the Ras / mitogen-activated ERK kinase (MEK) / extracellular signal-regulated kinase (ERK) cascade in controlling cell survival. Our results also highlight the importance of a group of interleukins together with the Janus kinase 3 (JAK3) / signal transducer and activator of transcription 3 (STAT3) signalling in the development of NPM-ALK derived ALCL. Depending on the activity of JAK3 and STAT3, the system may also be sensitive to activation of protein tyrosine phosphatase-1 (SHP1), which has an inhibitory effect on cell survival and proliferation. The identification of signalling pathways active in tumourigenic processes is of fundamental importance for effective therapies. The prediction of alternative pathways that circumvent classical therapeutic targets opens the way to preventive approaches for countering the emergence of cancer resistance.