In silico simulation of the effect of hypoxia on MCF-7 cell cycle kinetics under fractionated radiotherapy

The treatment outcome of a given fractionated radiotherapy scheme is affected by oxygen tension and cell cycle kinetics of the tumor population. Numerous experimental studies have supported the variability of radiosensitivity with cell cycle phase. Oxygen modulates the radiosensitivity through hypoxia-inducible factor (HIF) stabilization and oxygen fixation hypothesis (OFH) mechanism. In this study, an existing mathematical model describing cell cycle kinetics was modified to include the oxygen-dependent G1/S transition rate and radiation inactivation rate. The radiation inactivation rate used was derived from the linear-quadratic (LQ) model with dependence on oxygen enhancement ratio (OER), while the oxygen-dependent correction for the G1/S phase transition was obtained from numerically solving the ODE system of cyclin D-HIF dynamics at different oxygen tensions. The corresponding cell cycle phase fractions of aerated MCF-7 tumor population, and the resulting growth curve obtained from numerically solving the developed mathematical model were found to be comparable to experimental data. Two breast radiotherapy fractionation schemes were investigated using the mathematical model. Results show that hypoxia causes the tumor to be more predominated by the tumor subpopulation in the G1 phase and decrease the fractional contribution of the more radioresistant tumor cells in the S phase. However, the advantage provided by hypoxia in terms of cell cycle phase distribution is largely offset by the radioresistance developed through OFH. The delayed proliferation caused by severe hypoxia slightly improves the radiotherapy efficacy compared to that with mild hypoxia for a high overall treatment duration as demonstrated in the 40-Gy fractionation scheme.     

The deubiquitinating enzyme USP37 enhances CHK1 activity to promote the cellular response to replication stress

The deubiquitinating enzyme USP37 is known to contribute to timely onset of S phase and progression of mitosis. However, it is not clear if USP37 is required beyond S-phase entry despite expression and activity of USP37 peaking within S phase. We have utilized flow cytometry and microscopy to analyze populations of replicating cells labeled with thymidine analogs and monitored mitotic entry in synchronized cells to determine that USP37-depleted cells exhibited altered S-phase kinetics. Further analysis revealed that cells depleted of USP37 harbored increased levels of the replication stress and DNA damage markers γH2AX and 53BP1 in response to perturbed replication. Depletion of USP37 also reduced cellular proliferation and led to increased sensitivity to agents that induce replication stress. Underlying the increased sensitivity, we found that the checkpoint kinase 1 is destabilized in the absence of USP37, attenuating its function. We further demonstrated that USP37 deubiquitinates checkpoint kinase 1, promoting its stability. Together, our results establish that USP37 is required beyond S-phase entry to promote the efficiency and fidelity of replication. These data further define the role of USP37 in the regulation of cell proliferation and contribute to an evolving understanding of USP37 as a multifaceted regulator of genome stability.     

CpG island libraries from human Chromosomes 18 and 22: landmarks for novel genes

CpG islands are found at the 5′ end of approximately 60% of human genes and so are important genomic landmarks. They are concentrated in early-replicating, highly acetylated gene-rich regions. With respect to CpG island content, human Chrs 18 and 22 are very different from each other: Chr 18 appears to be CpG island poor, whereas Chr 22 appears to be CpG island rich. We have constructed and validated CpG island libraries from flow-sorted Chrs 18 and 22 and used these to estimate the difference in number of CpG islands found on these two chromosomes. These libraries contain normalized collections of sequences from the 5′ end of genes. Clones from the libraries were sequenced and compared with the sequence databases; one third matched ESTs, thus anchoring these ESTs at the 5′ end of their gene. However, it was striking that many clones either had no match or matched only existing CpG island clones. This suggests that a significant proportion of 5′ gene sequences are absent from databases, presumably either because they are difficult to clone or the gene is poorly expressed and/or has a restricted expression pattern. This point should be taken into consideration if the currently available libraries are those used for the elucidation of complete, as opposed to partial, gene sequences. The Chr 18 and 22 CpG island libraries are a sequence resource for the isolation of such 5′ gene sequences from specific human chromosomes. Received: 15 November 1999 / Accepted: 31 January 2000     

Genome-scale patterns in the loss of heterozygosity incidence in Saccharomyces cerevisiae

Former studies have established that loss of heterozygosity can be a key driver of sequence evolution in unicellular eukaryotes and tissues of metazoans. However, little is known about whether the distribution of loss of heterozygosity events is largely random or forms discernible patterns across genomes. To initiate our experiments, we introduced selectable markers to both arms of all chromosomes of the budding yeast. Subsequent extensive assays, repeated over several genetic backgrounds and environments, provided a wealth of information on the genetic and environmental determinants of loss of heterozygosity. Three findings stand out. First, the number of loss of heterozygosity events per unit time was more than 25 times higher for growing than starving cells. Second, loss of heterozygosity was most frequent when regions of homology around a recombination site were identical, about a half-% sequence divergence was sufficient to reduce its incidence. Finally, the density of loss of heterozygosity events was highly dependent on the genome’s physical architecture. It was several-fold higher on short chromosomal arms than on long ones. Comparably large differences were seen within a single arm where regions close to a centromere were visibly less affected than regions close, though usually not strictly adjacent, to a telomere. We suggest that the observed uneven distribution of loss of heterozygosity events could have been caused not only by an uneven density of initial DNA damages. Location-depended differences in the mode of DNA repair, or its effect on fitness, were likely to operate as well.     

Liver glycogen phosphorylase is upregulated in glioblastoma and provides a metabolic vulnerability to high dose radiation

Channelling of glucose via glycogen, known as the glycogen shunt, may play an important role in the metabolism of brain tumours, especially in hypoxic conditions. We aimed to dissect the role of glycogen degradation in glioblastoma (GBM) response to ionising radiation (IR). Knockdown of the glycogen phosphorylase liver isoform (PYGL), but not the brain isoform (PYGB), decreased clonogenic growth and survival of GBM cell lines and sensitised them to IR doses of 10–12 Gy. Two to five days after IR exposure of PYGL knockdown GBM cells, mitotic catastrophy and a giant multinucleated cell morphology with senescence-like phenotype developed. The basal levels of the lysosomal enzyme alpha-acid glucosidase (GAA), essential for autolysosomal glycogen degradation, and the lipidated forms of gamma-aminobutyric acid receptor-associated protein-like (GABARAPL1 and GABARAPL2) increased in shPYGL U87MG cells, suggesting a compensatory mechanism of glycogen degradation. In response to IR, dysregulation of autophagy was shown by accumulation of the p62 and the lipidated form of GABARAPL1 and GABARAPL2 in shPYGL U87MG cells. IR increased the mitochondrial mass and the colocalisation of mitochondria with lysosomes in shPYGL cells, thereby indicating reduced mitophagy. These changes coincided with increased phosphorylation of AMP-activated protein kinase and acetyl-CoA carboxylase 2, slower ATP generation in response to glucose loading and progressive loss of oxidative phosphorylation. The resulting metabolic deficiencies affected the availability of ATP required for mitosis, resulting in the mitotic catastrophy observed in shPYGL cells following IR. PYGL mRNA and protein levels were higher in human GBM than in normal human brain tissues and high PYGL mRNA expression in GBM correlated with poor patient survival. In conclusion, we show a major new role for glycogen metabolism in GBM cancer. Inhibition of glycogen degradation sensitises GBM cells to high-dose IR indicating that PYGL is a potential novel target for the treatment of GBMs.Subject terms: Cancer metabolism, CNS cancer     

The bacterial yjdF riboswitch regulates translation through its tRNA-like fold

Unlike most riboswitches, which have one cognate effector, the bacterial yjdF riboswitch binds to diverse azaaromatic compounds, only a subset of which cause it to activate translation. We examined the yjdF aptamer domain by small-angle X-ray scattering and found that in the presence of activating ligands, the RNA adopts an overall shape similar to that of tRNA. Sequence analyses suggested that the yjdF aptamer is a homolog of tRNA^Lys, and that two of the conserved loops of the riboswitch are equivalent to the D-loop and T-loop of tRNA, associating to form an elbow-like tertiary interaction. Chemical probing indicated that this association is promoted by activating ligands such as chelerythrine and harmine. In its native mRNA context, activator ligands stabilize the tRNA-like fold of the yjdF aptamer, outcompeting the attenuated state in which its T-loop base pairs to the Shine–Dalgarno element of the mRNA. Moreover, we demonstrate that the liganded aptamer itself activates translation, as authentic tRNAs, when grafted into mRNA, can potently activate translation. Taken together, our data demonstrate the ability of tRNA to function as a small-molecule responsive cis regulatory element.     

Characterizing the Propagation of Uterine Electrophysiological Signals Recorded with a Multi-Sensor Abdominal Array in Term Pregnancies

The objective of this study was to quantify the number of segments that have contractile activity and determine the propagation speed from uterine electrophysiological signals recorded over the abdomen. The uterine magnetomyographic (MMG) signals were recorded with a 151 channel SARA (SQUID Array for Reproductive Assessment) system from 36 pregnant women between 37 and 40 weeks of gestational age. The MMG signals were scored and segments were classified based on presence of uterine contractile burst activity. The sensor space was then split into four quadrants and in each quadrant signal strength at each sample was calculated using center-of-gravity (COG). To this end, the cross-correlation analysis of the COG was performed to calculate the delay between pairwise combinations of quadrants. The relationship in propagation across the quadrants was quantified and propagation speeds were calculated from the delays. MMG recordings were successfully processed from 25 subjects and the average values of propagation speeds ranged from 1.3–9.5 cm/s, which was within the physiological range. The propagation was observed between both vertical and horizontal quadrants confirming multidirectional propagation. After the multiple pairwise test (99% CI), significant differences in speeds can be observed between certain vertical or horizontal combinations and the crossed pair combinations. The number of segments containing contractile activity in any given quadrant pair with a detectable delay was significantly higher in the lower abdominal pairwise combination as compared to all others. The quadrant-based approach using MMG signals provided us with high spatial-temporal information of the uterine contractile activity and will help us in the future to optimize abdominal electromyographic (EMG) recordings that are practical in a clinical setting.