CDKN2A/p16 is inactivated in most melanoma cell lines

The CDKN2A gene maps to chromosome 9p21-22 and is responsible for melanoma susceptibility in some families. Its product, p16, binds specifically to CDK4 and CDK6 in vitro and in vivo, inhibiting their kinase activity. CDKN2A is homozygously deleted or mutated in a large proportion of tumor cell lines and some primary tumors, including melanomas. The aim of this study was to investigate the involvement of CDKN2A and elucidate the mechanisms of p16 inactivation in a panel of 60 cell lines derived from sporadic melanomas. Twenty-six (43%) of the melanoma lines were homozygously deleted for CDKN2A, and an additional 15 (25%) lines carried missense, nonsense, or frameshift mutations. All but one of the latter group were shown by microsatellite analysis to be hemizygous for the region of 9p surrounding CDKN2A. p16 was detected by Western blotting in only five of the cell lines carrying mutations. Immunoprecipitation of p16 in these lines, followed by Western blotting to detect the coprecipitation of CDK4 and CDK6, revealed that p16 was functionally compromised in all cell lines but the one that carried a heterozygous CDKN2A mutation. In the remaining 19 lines that carried wild-type CDKN2A alleles, Western blot analysis and immunoprecipitation indicated that 11 cell lines expressed a wild-type protein. Northern blotting was performed on the remaining eight cell lines and revealed that one cell line carried an aberrantly sized RNA transcript, and two other cell lines failed to express RNA. The promoter was found to be methylated in five cell lines that expressed CDKN2A transcript but not p16. Presumably, the message seen by Northern blotting in these cell lines is the result of cross-hybridization of the total cDNA probe with the exon 1beta transcript. Microsatellite analysis revealed that the majority of these cell lines were hemi/homozygous for the region surrounding CDKN2A, indicating that the wild-type allele had been lost. In the 11 cell lines that expressed functional p16, microsatellite analysis revealed loss of heterozygosity at the markers immediately surrounding CDKN2A in five cases, and the previously characterized R24C mutation of CDK4 was identified in one of the remaining 6 lines. These data indicate that 55 of 60 (92%) melanoma cell lines demonstrated some aberration of CDKN2A or CDK4, thus suggesting that this pathway is a primary genetic target in melanoma development.

CDKN2A is not the principal target of deletions on the short arm of chromosome 9 in neuroendocrine (Merkel cell) carcinoma of the skin

The majority of small-cell lung cancers (SCLCs) express p16 but not pRb. Given our previous study showing loss of pRb in Merkel cell carcinoma (MCC)/neuroendocrine carcinoma of the skin and the clinicopathological similarities between SCLC and MCC, we wished to determine if this was also the case in MCC. Twenty-nine MCC specimens from 23 patients were examined for deletions at 10 loci on 9p and 1 on 9q. No loss of heterozygosity (LOH) was seen in 9 patients including 2 for which tumour and cell line DNAs were examined. Four patients had LOH for all informative loci on 9p. Ten tumours showed more limited regions of loss on 9p, and from these 2 common regions of deletion were determined. Half of all informative cases had LOH at D9S168, the most telomeric marker examined, and 3 specimens showed loss of only D9S168. A second region (IFNA-D9S126) showed LOH in 10 (44%) cases, and case MCC26 showed LOH for only D9S126, implicating genes centromeric of the CDKN2A locus. No mutations in the coding regions of p16 were seen in 7 cell lines tested, and reactivity to anti-p16 antibody was seen in all 11 tumour specimens examined and in 6 of 7 cell lines from 6 patients. Furthermore, all cell lines examined reacted with anti-p14(ARF) antibody. These results suggest that neither transcript of the CDKN2A locus is the target of deletions on 9p in MCC and imply the existence of tumour-suppressor genes mapping both centromeric and telomeric of this locus.

Mutations in exon 3 of the beta-catenin gene are rare in melanoma cell lines

Mutations in exon 3 of the CTNNB1 gene encoding beta-catenin have been reported in colorectal cancer cell lines and tumours. Although one study reported mutations or deletions affecting beta-catenin in 20% of melanoma cell lines, subsequent reports detected a much lower frequency of aberrations in uncultured melanomas. To determine whether this difference in mutation frequency reflected an in vitro culturing artefact, exon 3 of CTNNB1 was screened in a panel of 62 melanoma cell lines. In addition, reverse transcription-polymerase chain reaction (RT-PCR) was performed to detect intragenic deletions affecting exon 3. One out of 62 (1.6%) cell lines was found to carry a mutation, indicating that aberration of the Wnt-1/wingless pathway through activation of beta-catenin is a rare event, even in melanoma cell lines.

PTEN inactivation is rare in melanoma tumours but occurs frequently in melanoma cell lines

Deletions detected in cytogenetic and loss of heterozygosity (LOH) studies indicate that at least one tumour suppressor gene maps to the long arm of chromosome 10. Previous deletion mapping studies have observed LOH on 10q in about 30% of melanomas analysed. The PTEN gene, mapping to chromosome band 10q23.3, encodes a protein with both lipid and protein phosphatase activity. Somatic mutations and deletions in have been detected in a variety of cell lines and tumours, including melanoma samples. We performed mutation analyses and extensive allelic loss studies to investigate the role this gene plays in melanoma pathogenesis. We found that a total of 34 out of 57 (60%) melanoma cell lines carried hemizygous deletions of chromosome 10q encompassing the PTEN locus. A further three cell lines carried smaller deletions excluding PTEN. Inactivation of both PTEN alleles by exon-specific homozygous deletion or mutation was observed in 13 out of 57 (23%) melanoma cell lines. The mutation spectrum observed does not indicate an important role for ultraviolet radiation in the genesis of these mutations, and evidence from three cell lines supports the acquisition of PTEN aberrations in culture. Ten out of 49 (20%) matched melanoma tumour/normal samples harboured hemizygous deletions of either the whole chromosome or most of the long arm. Mutations within were detected in only one of the 10 tumours demonstrating LOH at 10q23 that were analysed. These results suggest that PTEN inactivation may be important for the propagation of melanoma cells in culture, and that another chromosome 10 tumour suppressor gene may be important for melanoma pathogenesis.

p53 prevents progression of nevi to melanoma predominantly through cell cycle regulation

p53 is the central member of a critical tumor suppressor pathway in virtually all tumor types, where it is silenced mainly by missense mutations. In melanoma, p53 predominantly remains wild type, thus its role has been neglected. To study the effect of p53 on melanocyte function and melanomagenesis, we crossed the 'high-p53'Mdm4+/- mouse to the well-established TP-ras0/+ murine melanoma progression model. After treatment with the carcinogen dimethylbenzanthracene (DMBA), TP-ras0/+ mice on the Mdm4+/- background developed fewer tumors with a delay in the age of onset of melanomas compared to TP-ras0/+ mice. Furthermore, we observed a dramatic decrease in tumor growth, lack of metastasis with increased survival of TP-ras0/+: Mdm4+/- mice. Thus, p53 effectively prevented the conversion of small benign tumors to malignant and metastatic melanoma. p53 activation in cultured primary melanocyte and melanoma cell lines using Nutlin-3, a specific Mdm2 antagonist, supported these findings. Moreover, global gene expression and network analysis of Nutlin-3-treated primary human melanocytes indicated that cell cycle regulation through the p21WAF1/CIP1 signaling network may be the key anti-melanomagenic activity of p53.

Stochastic modelling of T cell homeostasis for two competing clonotypes via the master equation

Stochastic models for competing clonotypes of T cells by multivariate, continuous-time, discrete state, Markov processes have been proposed in the literature by Stirk, Molina-París and van den Berg (2008). A stochastic modelling framework is important because of rare events associated with small populations of some critical cell types. Usually, computational methods for these problems employ a trajectory-based approach, based on Monte Carlo simulation. This is partly because the complementary, probability density function (PDF) approaches can be expensive but here we describe some efficient PDF approaches by directly solving the governing equations, known as the Master Equation. These computations are made very efficient through an approximation of the state space by the Finite State Projection and through the use of Krylov subspace methods when evolving the matrix exponential. These computational methods allow us to explore the evolution of the PDFs associated with these stochastic models, and bimodal distributions arise in some parameter regimes. Time-dependent propensities naturally arise in immunological processes due to, for example, age-dependent effects. Incorporating time-dependent propensities into the framework of the Master Equation significantly complicates the corresponding computational methods but here we describe an efficient approach via Magnus formulas. Although this contribution focuses on the example of competing clonotypes, the general principles are relevant to multivariate Markov processes and provide fundamental techniques for computational immunology.

Stochastic simulation for spatial modelling of dynamic process in a living cell

One of the fundamental motivations underlying computational cell biology is to gain insight into the complicated dynamical processes taking place, for example, on the plasma membrane or in the cytosol of a cell. These processes are often so complicated that purely temporal mathematical models cannot adequately capture the complex chemical kinetics and transport processes of, for example, proteins or vesicles. On the other hand, spatial models such as Monte Carlo approaches can have very large computational overheads. This chapter gives an overview of the state of the art in the development of stochastic simulation techniques for the spatial modelling of dynamic processes in a living cell.

Linear models for endocytic transformations from live cell imaging

Endocytosis is the process by which cells internalise molecules including nutrient proteins from the extracellular media. In one form, macropinocytosis, the membrane at the cell surface ruffles and folds over to give rise to an internalised vesicle. Negatively charged phospholipids within the membrane called phosphoinositides then undergo a series of transformations that are critical for the correct trafficking of the vesicle within the cell, and which are often pirated by pathogens such as Salmonella. Advanced fluorescent video microscopy imaging now allows the detailed observation and quantification of these events in live cells over time. Here we use these observations as a basis for building differential equation models of the transformations. An initial investigation of these interactions was modelled with reaction rates proportional to the sum of the concentrations of the individual constituents. A first order linear system for the concentrations results. The structure of the system enables analytical expressions to be obtained and the problem becomes one of determining the reaction rates which generate the observed data plots. We present results with reaction rates which capture the general behaviour of the reactions so that we now have a complete mathematical model of phosphoinositide transformations that fits the experimental observations. Some excellent fits are obtained with modulated exponential functions; however, these are not solutions of the linear system. The question arises as to how the model may be modified to obtain a system whose solution provides a more accurate fit.

Phenomenological modeling of cell-to-cell and beat-to-beat variability in isolated Guinea Pig ventricular myocytes

Experimental action potential (AP) recordings in isolated ventricular myoctes display significant temporal beat-to-beat variability in morphology and duration. Furthermore, significant cell-to-cell differences in AP also exist even for isolated cells originating from the same region of the same heart. However, current mathematical models of ventricular AP fail to replicate the temporal and cell-to-cell variability in AP observed experimentally. In this study, we propose a novel mathematical framework for the development of phenomenological AP models capable of capturing cell-to-cell and temporal variabilty in cardiac APs. A novel stochastic phenomenological model of the AP is developed, based on the deterministic Bueno-Orovio/Fentonmodel. Experimental recordings of AP are fit to the model to produce AP models of individual cells from the apex and the base of the guinea-pig ventricles. Our results show that the phenomenological model is able to capture the considerable differences in AP recorded from isolated cells originating from the location. We demonstrate the closeness of fit to the available experimental data which may be achieved using a phenomenological model, and also demonstrate the ability of the stochastic form of the model to capture the observed beat-to-beat variablity in action potential duration.

An analytical solution for diffusion and nonlinear uptake of oxygen in a spherical cell

We develop a new analytical solution for a reactive transport model that describes the steady-state distribution of oxygen subject to diffusive transport and nonlinear uptake in a sphere. This model was originally reported by Lin (Journal of Theoretical Biology, 1976 v60, pp449–457) to represent the distribution of oxygen inside a cell and has since been studied extensively by both the numerical analysis and formal analysis communities. Here we extend these previous studies by deriving an analytical solution to a generalized reaction-diffusion equation that encompasses Lin’s model as a particular case. We evaluate the solution for the parameter combinations presented by Lin and show that the new solutions are identical to a grid-independent numerical approximation.

Noble metal photocatalysts under visible light and UV light irradiation

One of the greatest challenges for the study of photocatalysts is to devise new catalysts that possess high activity under visible light illumination. This would allow the use of an abundant and green energy source, sunlight, to drive chemical reactions. Gold nanoparticles strongly absorb both visible light and UV light. It is therefore possible to drive chemical reactions utilising a significant fraction of full sunlight spectrum. Here we prepared gold nanoparticles supported on various oxide powders, and reported a new finding that gold nanoparticles on oxide supports exhibit significant activity for the oxidation of formaldehyde and methanol in the air at ambient temperature, when illuminated with visible light. We suggested that visible light can greatly enhance local electromagnetic fields and heat gold nanoparticles due to surface plasmon resonance effect which provides activation energy for the oxidation of organic molecules. Moreover, the nature of the oxide support has an important influence on the activity of the gold nanoparticles. The finding reveals the possibility to drive chemical reactions with sunlight on gold nanoparticles at ambient temperature, highlighting a new direction for research on visible light photocatalysts. Gold nanoparticles supported on oxides also exhibit significant dye oxidation activity under visible light irradiation in aqueous solution at ambient temperature. Turnover frequencies of the supported gold nanoparticles for the dye degradation are much higher than titania based photocatalysts under both visible and UV light. These gold photocatalysts can also catalyse phenol degradation as well as selective oxidation of benzyl alcohol under UV light. The reaction mechanism for these photocatalytic oxidations was studied. Gold nanoparticles exhibit photocatalytic activity due to visible light heating gold electrons in 6sp band, while the UV absorption results in electron holes in gold 5d band to oxidise organic molecules. Silver nanoparticles also exhibit considerable visible light and UV light absorption due to surface plasmon resonance effect and the interband transition of 4d electrons to the 5sp band, respectively. Therefore, silver nanoparticles are potentially photocatalysts that utilise the solar spectrum effectively. Here we reported that silver nanoparticles at room temperature can be used to drive chemical reactions when illuminated with light throughout the solar spectrum. The significant activities for dye degradation by silver nanoparticles on oxide supports are even better than those by semiconductor photocatalysts. Moreover, silver photocatalysts also can degrade phenol and drive the oxidation of benzyl alcohol to benzaldehyde under UV light. We suggested that surface plasmon resonance effect and interband transition of silver nanoparticles can activate organic molecule oxidations under light illumination.

Performance control of distributed generation using digital estimation of signal parameters

The Queensland University of Technology (QUT) allows the presentation of a thesis for the Degree of Doctor of Philosophy in the format of published or submitted papers, where such papers have been published, accepted or submitted during the period of candidature. This thesis is composed of seven published/submitted papers, of which one has been published, three accepted for publication and the other three are under review. This project is financially supported by an Australian Research Council (ARC) Discovery Grant with the aim of proposing strategies for the performance control of Distributed Generation (DG) system with digital estimation of power system signal parameters. Distributed Generation (DG) has been recently introduced as a new concept for the generation of power and the enhancement of conventionally produced electricity. Global warming issue calls for renewable energy resources in electricity production. Distributed generation based on solar energy (photovoltaic and solar thermal), wind, biomass, mini-hydro along with use of fuel cell and micro turbine will gain substantial momentum in the near future. Technically, DG can be a viable solution for the issue of the integration of renewable or non-conventional energy resources. Basically, DG sources can be connected to local power system through power electronic devices, i.e. inverters or ac-ac converters. The interconnection of DG systems to power system as a compensator or a power source with high quality performance is the main aim of this study. Source and load unbalance, load non-linearity, interharmonic distortion, supply voltage distortion, distortion at the point of common coupling in weak source cases, source current power factor, and synchronism of generated currents or voltages are the issues of concern. The interconnection of DG sources shall be carried out by using power electronics switching devices that inject high frequency components rather than the desired current. Also, noise and harmonic distortions can impact the performance of the control strategies. To be able to mitigate the negative effect of high frequency and harmonic as well as noise distortion to achieve satisfactory performance of DG systems, new methods of signal parameter estimation have been proposed in this thesis. These methods are based on processing the digital samples of power system signals. Thus, proposing advanced techniques for the digital estimation of signal parameters and methods for the generation of DG reference currents using the estimates provided is the targeted scope of this thesis. An introduction to this research – including a description of the research problem, the literature review and an account of the research progress linking the research papers – is presented in Chapter 1. One of the main parameters of a power system signal is its frequency. Phasor Measurement (PM) technique is one of the renowned and advanced techniques used for the estimation of power system frequency. Chapter 2 focuses on an in-depth analysis conducted on the PM technique to reveal its strengths and drawbacks. The analysis will be followed by a new technique proposed to enhance the speed of the PM technique while the input signal is free of even-order harmonics. The other techniques proposed in this thesis as the novel ones will be compared with the PM technique comprehensively studied in Chapter 2. An algorithm based on the concept of Kalman filtering is proposed in Chapter 3. The algorithm is intended to estimate signal parameters like amplitude, frequency and phase angle in the online mode. The Kalman filter is modified to operate on the output signal of a Finite Impulse Response (FIR) filter designed by a plain summation. The frequency estimation unit is independent from the Kalman filter and uses the samples refined by the FIR filter. The frequency estimated is given to the Kalman filter to be used in building the transition matrices. The initial settings for the modified Kalman filter are obtained through a trial and error exercise. Another algorithm again based on the concept of Kalman filtering is proposed in Chapter 4 for the estimation of signal parameters. The Kalman filter is also modified to operate on the output signal of the same FIR filter explained above. Nevertheless, the frequency estimation unit, unlike the one proposed in Chapter 3, is not segregated and it interacts with the Kalman filter. The frequency estimated is given to the Kalman filter and other parameters such as the amplitudes and phase angles estimated by the Kalman filter is taken to the frequency estimation unit. Chapter 5 proposes another algorithm based on the concept of Kalman filtering. This time, the state parameters are obtained through matrix arrangements where the noise level is reduced on the sample vector. The purified state vector is used to obtain a new measurement vector for a basic Kalman filter applied. The Kalman filter used has similar structure to a basic Kalman filter except the initial settings are computed through an extensive math-work with regards to the matrix arrangement utilized. Chapter 6 proposes another algorithm based on the concept of Kalman filtering similar to that of Chapter 3. However, this time the initial settings required for the better performance of the modified Kalman filter are calculated instead of being guessed by trial and error exercises. The simulations results for the parameters of signal estimated are enhanced due to the correct settings applied. Moreover, an enhanced Least Error Square (LES) technique is proposed to take on the estimation when a critical transient is detected in the input signal. In fact, some large, sudden changes in the parameters of the signal at these critical transients are not very well tracked by Kalman filtering. However, the proposed LES technique is found to be much faster in tracking these changes. Therefore, an appropriate combination of the LES and modified Kalman filtering is proposed in Chapter 6. Also, this time the ability of the proposed algorithm is verified on the real data obtained from a prototype test object. Chapter 7 proposes the other algorithm based on the concept of Kalman filtering similar to those of Chapter 3 and 6. However, this time an optimal digital filter is designed instead of the simple summation FIR filter. New initial settings for the modified Kalman filter are calculated based on the coefficients of the digital filter applied. Also, the ability of the proposed algorithm is verified on the real data obtained from a prototype test object. Chapter 8 uses the estimation algorithm proposed in Chapter 7 for the interconnection scheme of a DG to power network. Robust estimates of the signal amplitudes and phase angles obtained by the estimation approach are used in the reference generation of the compensation scheme. Several simulation tests provided in this chapter show that the proposed scheme can very well handle the source and load unbalance, load non-linearity, interharmonic distortion, supply voltage distortion, and synchronism of generated currents or voltages. The purposed compensation scheme also prevents distortion in voltage at the point of common coupling in weak source cases, balances the source currents, and makes the supply side power factor a desired value.

Melanoma cell invasiveness is regulated by miR-211 suppression of the BRN2 transcription factor

To identify microRNAs potentially involved in melanomagenesis, we compared microRNA expression profiles between melanoma cell lines and cultured melanocytes. The most differentially expressed microRNA between the normal and tumor cell lines was miR-211. We focused on this pigment-cell-enriched miRNA as it is derived from the microphthalmia-associated transcription factor (MITF)-regulated gene, TRPM1 (melastatin). We find that miR-211 expression is greatly decreased in melanoma cells and melanoblasts compared to melanocytes. Bioinformatic analysis identified a large number of potential targets of miR-211, including POU3F2 (BRN2). Inhibition of miR-211 in normal melanocytes resulted in increased BRN2 protein, indicating that endogenous miR-211 represses BRN2 in differentiated cells. Over-expression of miR-211 in melanoma cell lines changed the invasive potential of the cells in vitro through directly targeting BRN2 translation. We propose a model for the apparent non-overlapping expression levels of BRN2 and MITF in melanoma, mediated by miR-211 expression.

Multifractal analysis of geomagnetic storm and solar flare indices and their class dependence

The multifractal properties of two indices of geomagnetic activity, D st (representative of low latitudes) and a p (representative of the global geomagnetic activity), with the solar X-ray brightness, X l , during the period from 1 March 1995 to 17 June 2003 are examined using multifractal detrended fluctuation analysis (MF-DFA). The h(q) curves of D st and a p in the MF-DFA are similar to each other, but they are different from that of X l , indicating that the scaling properties of X l are different from those of D st and a p . Hence, one should not predict the magnitude of magnetic storms directly from solar X-ray observations. However, a strong relationship exists between the classes of the solar X-ray irradiance (the classes being chosen to separate solar flares of class X-M, class C, and class B or less, including no flares) in hourly measurements and the geomagnetic disturbances (large to moderate, small, or quiet) seen in D st and a p during the active period. Each time series was converted into a symbolic sequence using three classes. The frequency, yielding the measure representations, of the substrings in the symbolic sequences then characterizes the pattern of space weather events. Using the MF-DFA method and traditional multifractal analysis, we calculate the h(q), D(q), and τ (q) curves of the measure representations. The τ (q) curves indicate that the measure representations of these three indices are multifractal. On the basis of this three-class clustering, we find that the h(q), D(q), and τ (q) curves of the measure representations of these three indices are similar to each other for positive values of q. Hence, a positive flare storm class dependence is reflected in the scaling exponents h(q) in the MF-DFA and the multifractal exponents D(q) and τ (q). This finding indicates that the use of the solar flare classes could improve the prediction of the D st classes.