Estimation of the contribution of ionization and excitation to the lethal effect of ionizing radiation

Summary A simple theoretical model is proposed for estimating the differential contribution of ionization and excitation to the lethal effect of ionizing radiation. Numerical results were obtained on the basis of published experimental data on the ability of bacterial cellsEscherichia coli to undergo photoreactivation of radiation-induced damage. It was shown that inactivation by excitation may be highly significant for UV-hypersensitive cells capable of photoreactivation; inactivation by excitation increased with the energy of ionizing radiation and the volume of irradiated suspensions. The data are in qualitative agreement with the assumption of a possible contribution of the UV-component of erenkov radiation to the formation of excitations responsible for the lethal effect and the phenomenon of photoreactivation after ionizing radiation. Some predictions from the model are discussed.     

Estimation of the probability parameters for a model of radiation inactivation of cells from experimental survival curves

A simple method is proposed for estimation of a and b parameters of the probability model of radiation inactivation of cells with a reference to the experimental survival curves. The examples of such an estimation for bacteria, yeast, and mammalian cells are discussed.     

Cell inactivation by heavy charged particles

The inactivation of cells resulting in lethal or aberrant effects by charged particles is of growing interest. Charged particles at extremely high LET are capable of completely eliminating cell-type and cell-line differences in repair capacity. It is still not clear however whether the repair systems are inactivated, or merely that heavy-ion lesions are less repairable. Studies correlating the particle inactivation dose of radioresistant cells with intact DNA analyzed with pulse field gel electrophoresis and other techniques may be useful, but more experiments are also needed to assess the fidelity of repair. For particle irradiations between 40-100 keV/microns there is however evidence for particle-induced activation of specific genes in mammalian cells, and certain repair processes in bacteria. New data are available on the inactivation of developmental processes in several systems including seeds, and cells of the nematode C. elegans. Future experimental and theoretical modeling research emphasis should focus on exploring particle-induced inactivation of endpoints assessing functionality and not just lethality, and on analyzing molecular damage and genetic effects arising in damaged but non-inactivated survivors. The discrete nature of selective types of particle damage as a function of radiation quality indicates the value of accelerated ions as probes of normal and aberrant biological processes. Information obtained from molecular analyses of damage and repair must however be integrated into the context of cellular and tissue functions of the organism.     

DNA damage caused by ionizing radiation.

A survey is given of continuous-time Markov chain models for ionizing radiation damage to the genome of mammalian cells. In such models, immediate damage induced by the radiation is regarded as a batch-Poisson arrival process of DNA double-strand breaks (DSBs). Enzymatic modification of the immediate damage is modeled as a Markov process similar to those described by the master equation of stochastic chemical kinetics. An illustrative example is the restitution/complete-exchange model. The model postulates that, after being induced by radiation, DSBs subsequently either undergo enzymatically mediated restitution (repair) or participate pairwise in chromosome exchanges. Some of the exchanges make irremediable lesions such as dicentric chromosome aberrations. One may have rapid irradiation followed by enzymatic DSB processing or have prolonged irradiation with both DSB arrival and enzymatic DSB processing continuing throughout the irradiation period. Methods for analyzing the Markov chains include using an approximate model for expected values, the discrete-time Markov chain embedded at transitions, partial differential equations for generating functions, normal perturbation theory, singular perturbation theory with scaling, numerical computations, and certain matrix methods that combine Perron-Frobenius theory with variational estimates. Applications to experimental results on expected values, variances, and statistical distributions of DNA lesions are briefly outlined. Continuous-time Markov chains are the most systematic of those radiation damage models that treat DSB-DSB interactions within the cell nucleus as homogeneous (e.g., ignore diffusion limitations). They contain virtually all other relevant homogeneous models and semiempirical summaries as special cases, limiting cases, or approximations. However, the Markov models do not seem to be well suited for studying spatial dependence of DSB interactions, which is known to be important in some situations.     

Dose-response models for the radiation-induction of skin tumours in mice

Extensive data on radiation-induced skin tumours in mice were examined using 8 models, all based on the concept that incidences of radiation-induced tumours depend on a combination of two radiation effects: a tumour induction process and the loss of reproductive integrity by the potential tumour cells. Models with and without a threshold were used, in spite of theoretical objections to threshold models. No model fitted well both the epidermal and the dermal tumour data and models which proved to be statistically satisfactory for some of the data were rejected for biological reasons. It is concluded that, for skin tumours, dose-response curves depending on a combination of cancer induction and loss of cellular reproductive integrity are distorted by some special, relatively radio-resistant, factor which we have previously postulated as being involved in radiation skin carcinogenesis.     

Cell cycle modelling

Models able to describe the events of cellular growth and division and the dynamics of cell populations are useful for the understanding of functional control mechanisms and for the theoretical support for automated analysis of flow cytometric data and of cell volume distributions. This paper reports on models that we have developed with this aim for different kinds of cells. The models are composed by two subsystems: one describes the growth dynamics of RNA and protein, and the second accounts for DNA replication and cell division, and describe in a rather unitary frame the cell cycle of eukaryotic cells, like mammalian cells and yeast, and of prokaryotic cells. The model is also used to study the effects of various sources of variability on the statistical properties of cell populations, and we find that in microbial cells the main source of variability appears to be an inaccuracy of the molecular mechanism that monitors cell size. In normal mammalian cells another source of variability, that depends upon the interaction with growth factors which give competence, is apparent. An extended version of the model, which comprises also this additional variability, is presented and used to describe the properties of mammalian cell growth.     

Radiation mapping as one of main methods of creating high resolution maps of human and animal genomes

Radiation hybrid mapping (RH mapping) is considered as one of the main methods of constructing physical maps of mammalian genomes. In introduction, theoretical prerequisites of developing of the RH mapping and statistical methods of data analysis are discussed. Comparative characteristics of universal commercial panels of the radiation hybrid somatic cells (RH panels) are shown. In experimental part of the work, RH mapping is used to localise nucleotide sequences adjacent to NotI sites of human chromosome 3 with the aim to integrate contig map of NotI clones to comprehensive maps of human genome. Five nucleotide sequences adjacent to the sites of integration of papilloma virus in human genome and expressed in the cells of cervical cancer were localised. It was demonstrated that the region 13q14.3-q21.1 was enriched with nucleotide sequences involved in the processes of oncogenesis. RH mapping can be considered as one of the most perspective applications of the modern radiation biology in the field of molecular genetics, that is, in constructing physical maps of mammalian genomes with high resolution level.     

A systemic model of the reproductive death of mammalian cells. The effect of radiation dose rate and fractionation

On the basis of the previously developed systemic model a study was made of the effect of dose rate on the survival of mammalian cells, RBE of small doses of heavy ions, and fractionation of radiation. There was a good agreement between theoretical and experimental results. The calculations showed that D10 (10% survival dose) is a function of dose rate P even for such ions as helium and boron which, however, exhibited an insignificant dependence of D10 on P (within the range from (10(-1) to 1 cGy/min). The influence of repair and the rate of cell division on RBE of radiation was determined.     

Mutation induction by inhaled radon progeny modeled at the tissue level

The observable responses of living systems to ionizing radiation depend on the level of biological organization studied. Understanding the relationships between the responses characteristic of the different levels of organization is of crucial importance. The main objective of the present study is to investigate how some cellular effects of radiation manifest at the tissue level by modeling mutation induction due to chronic exposure to inhaled radon progeny. For this purpose, a mathematical model of the bronchial epithelium was elaborated to quantify cell nucleus hits and cell doses. Mutagenesis was modeled considering endogenous as well as radiation-induced DNA damages and cell cycle shortening due to cell inactivation. The model parameters describing the cellular effects of radiation are obtained from experimental data. Cell nucleus hits, cell doses, and mutation induction were computed for the activity hot spots of the large bronchi at different exposures. Results demonstrate that the mutagenic effect of densely ionizing radiation is dominated by cell cycle shortening due to cell inactivation and not by DNA damages. This suggests that radiation burdens of non-progenitor cells play a significant role in mutagenesis in case of protracted exposures to densely ionizing radiation. Mutation rate as a function of dose rate exhibits a convex shape below a threshold. This threshold indicates the exhaustion of the tissue regeneration capacity of local progenitor cells. It is suggested that progenitor cell hyperplasia occurs beyond the threshold dose rate, giving a possible explanation of the inverse dose-rate effect observed in the epidemiology of lung cancer among uranium miners.     

Modeling temporal behavior of postnatal cat retinal ganglion cells

During development, mammalian retinal ganglion cells (RGCs) go through marked ontogenetic changes with respect to their excitable membrane properties. Voltage-clamp studies conducted in our laboratory have shown that the amplitude, voltage-dependence and kinetics of activation and inactivation (where present) of Na(+), K(+) and Ca(2+) conductances all exhibit developmental changes during a time when the firing patterns of mammalian ganglion cells shift from being transient to being predominantly sustained in nature. In order to better understand the contribution of each conductance to the generation of spikes and spiking patterns, we have developed a model based on our experimental data. For simplicity, we have initially used experimental data obtained from postnatal ganglion cells. At this age the ontogenetic changes observed in the characteristics of the various ionic currents are complete. Utilizing the methods adopted by Hodgkin and Huxley for the giant squid axon, we have determined rate equations for the activation and inactivation properties of the I(A), I(K dr), I(Na), I(Ca L), I(Ca N), and I(leak) currents in postnatal cat RGCs. Combining these with a simplified model of the calcium-activated potassium current (I(KCa)), we have solved and analysed the resulting differential equations. While spikes and spiking patterns resembling experimental data could be obtained from a model in which [Ca(2+)i] was averaged across the whole cell, more accurate simulations were obtained when the diffusion of intracellular Ca(2+) was modeled spatially. The resulting spatial calcium gradients were more effective in gating I(KCa), and our simulations more accurately matched the recorded amplitude and shape of individual spikes as well as the frequency of maintained discharges observed in mammalian postnatal RGCs.     

The role of the repair of double-stranded DNA breaks in the radioresistance of yeast cells

The role of DNA double-strand break (DSB) repair in radioresistance of Saccharomyces cerevisiae G1 cells is discussed. The contribution of rapid and slow DNA DSB repair to radioresistance of diploid yeast has been estimated. The contribution of the DNA DSB repair involving no homologous chromosome interaction is shown to be insignificant in comparison with the recombinational repair. The rapid DNA DSB repair efficiency calculation method based on the proposed yeast radiation inactivation model is given. The calculations are in a satisfactory agreement with the experimental data. Possible mechanisms of radiation induction of lethal sectoring in yeast are discussed. This phenomenon is supposed to be due to the DNA DSB processing during vegetative division of irradiated cells. A general scheme of radiation inactivation of yeast cells is proposed.     

Mathematical Modeling of the Phoenix Rising Pathway

Apoptosis is a tightly controlled process in mammalian cells. It is important for embryogenesis, tissue homoeostasis, and cancer treatment. Apoptosis not only induces cell death, but also leads to the release of signals that promote rapid proliferation of surrounding cells through the Phoenix Rising (PR) pathway. To quantitatively understand the kinetics of interactions of different molecules in this pathway, we developed a mathematical model to simulate the effects of various changes in the PR pathway on the secretion of prostaglandin E2 (PGE2), a key factor for promoting cell proliferation. These changes include activation of caspase 3 (C3), caspase 7 (C7), and nuclear factor κB (NFκB). In addition, we simulated the effects of cyclooxygenase-2 (COX2) inhibition and C3 knockout on the level of secreted PGE2. The model predictions on PGE2 in MEF and 4T1 cells at 48 hours after 10-Gray radiation were quantitatively consistent with the experimental data in the literature. Compared to C7, the model predicted that C3 activation was more critical for PGE2 production. The model also predicted that PGE2 production could be significantly reduced when COX2 expression was blocked via either NFκB inactivation or treatment of cells with exogenous COX2 inhibitors, which led to a decrease in the rate of conversion from arachidonic acid to prostaglandin H2 in the PR pathway. In conclusion, the mathematical model developed in this study yielded new insights into the process of tissue regrowth stimulated by signals from apoptotic cells. In future studies, the model can be used for experimental data analysis and assisting development of novel strategies/drugs for improving cancer treatment or normal tissue regeneration.     

Experimental and numerical study of heterogeneous pressure-temperature-induced lethal and sublethal injury of Lactococcus lactis in a medium scale high-pressure autoclave

The present contribution is dedicated to experimental and theoretical assessment of microbiological process heterogeneities of the high-pressure (HP) inactivation of Lactococcus lactis ssp. cremoris MG 1363. The inactivation kinetics are determined in dependence of pressure, process time, temperature and absence or presence of co-solutes in the buffer system namely 4 M sodium chloride and 1.5 M sucrose. The kinetic analysis is carried out in a 0.1-L autoclave in order to minimise thermal and convective effects. Upon these data, a deterministic inactivation model is formulated with the logistic equation. Its independent variables represent the counts of viable cells (viable but injured) and of the stress-resistant cells (viable and not injured). This model is then coupled to a thermo-fluiddynamical simulation method, high-pressure computer fluid dynamics technique (HP-CFD), which yields spatiotemporal temperature and flow fields occurring during the HP application inside any considered autoclave. Besides the thermo-fluiddynamic quantities, the coupled model predicts also the spatiotemporal distribution of both viable (VC) and stress-resistant cell counts (SRC). In order to assess the process non-uniformity of the microbial inactivation in a 3.3-L autoclave experimentally, microbial samples are placed at two distinct locations and are exposed to various process conditions. It can be shown with both, experimental and theoretical models that thermal heterogeneities induce process non-uniformities of more than one decimal power in the counts of the viable cells at the end of the treatment.     

Computational Modeling of Orientation Tuning Dynamics in Monkey Primary Visual Cortex

In the primate visual pathway, orientation tuning of neurons is first observed in the primary visual cortex. The LGN cells that comprise the thalamic input to V1 are not orientation tuned, but some V1 neurons are quite selective. Two main classes of theoretical models have been offered to explain orientation selectivity: feedforward models, in which inputs from spatially aligned LGN cells are summed together by one cortical neuron; and feedback models, in which an initial weak orientation bias due to convergent LGN input is sharpened and amplified by intracortical feedback. Recent data on the dynamics of orientation tuning, obtained by a cross-correlation technique, may help to distinguish between these classes of models. To test this possibility, we simulated the measurement of orientation tuning dynamics on various receptive field models, including a simple Hubel-Wiesel type feedforward model: a linear spatiotemporal filter followed by an integrate-and-fire spike generator. The computational study reveals that simple feedforward models may account for some aspects of the experimental data but fail to explain many salient features of orientation tuning dynamics in V1 cells. A simple feedback model of interacting cells is also considered. This model is successful in explaining the appearance of Mexican-hat orientation profiles, but other features of the data continue to be unexplained.     

The dose-response relationships for myeloid leukemia and malignant lymphoma in BC3F1 mice

The present analysis of data on the induction of lymphoma and myeloid leukemia in BC3F1 mice has indicated some new and interesting aspects regarding the shapes of the dose-effect curves. The incidence data can be interpreted by radiobiological models of the induction process coupled with cell inactivation. In particular, for malignant lymphoma the dose-response curve after X rays can be described assuming a quadratic model corrected for cell inactivation, while the incidence data after fission neutrons are best fitted by a linear model which also allows for cell inactivation. Myeloid leukemia has also been induced in BC3F1 mice. The bell-shaped dose-response curves observed after irradiation with either X rays or neutrons are explained by assuming simultaneous initial transforming events and cell inactivation with the data for cell inactivation at higher doses being in agreement with data reported for other strains of mice. A value for relative biological effectiveness of 4 is obtained at the lowest neutron dose used. The value of the inactivation parameters can be compared with those of the cell inactivation probability per unit dose for the bone marrow hematopoietic stem cells, which are believed to be the target cells for these tumors.     

Na channel kinetics: developing models from non-stationary current fluctuations by analytic methods

In a previous study, we analyzed Na current fluctuations in myelinated axons from Xenopus laevis under voltage clamp conditions. The statistical properties were analyzed in terms of covariance functions for consecutive time intervals of varying duration during the pulse step. The underlying channel kinetics was analyzed by performing stochastic simulations of published Na channel models and calculating corresponding covariance functions. None of the models explained the fluctuation results. We therefore developed a novel minimal Na channel model that satisfactorily described the results. In the present paper, we extend the analysis and specify the possible models explaining the experimental data by using analytical methods. We derive general relations between the experimental data, including the covariance functions, and the rate constants of specific one-open-state models. A general feature of these models is that they comprise an inactivation step from the first closed state and a relatively low backward rate from the open state. This is in accordance the minimal model inferred from numerical stochastic calculations in the previous study.     

Deriving the probabilities of water loss and ammonia loss for amino acids from tandem mass spectra

In protein identification through tandem mass spectrometry, it is critical to accurately predict the theoretical spectrum for a peptide sequence. The widely used prediction models, such as SEQUEST and MASCOT, ignore the intensity of the ions with important neutral losses, including water loss and ammonia loss. However, ignoring these neutral losses results in a significant deviation between the predicted theoretical spectrum and its experimental counterpart. Here, based on the "one peak, multiple explanations" observation, we proposed an expectation-maximization (EM) method to automatically learn the probabilities of water loss and ammonia loss for each amino acid. Then we employed these probabilities to design an improved statistical model for theoretical spectrum prediction. We implemented these methods and tested them on practical data. On a training set containing 1803 spectra, the experimental results show a good agreement with some known knowledge about neutral losses, such as the tendency of water loss from Asp, Glu, Ser, and Thr. Furthermore, on a testing set containing 941 spectra, the improved similarity between the experimental and predicted spectra demonstrates that this method can generate more reasonable predictions relative to the model that ignores neutral losses. As an application of the derived probabilities, we implemented a database searching method adopting the improved theoretical spectrum model with neutral loss ions estimated. Experimental results on Keller's data set demonstrate that this method can identify peptides more accurately than SEQUEST. In another application to validate SEQUEST's results, the reported peptide-spectrum pairs are reranked with respect to the similarity between experimental and predicted spectra. Experimental results on both LTQ and QSTAR data sets suggest that this reranking strategy can effectively distinguish the false negative predictions reported by SEQUEST.     

Stem cells and embryo-derived cell lines: tools for study of gene expression

Embryonal carcinoma (EC) cells obtained either from teratocarcinomas or directly from in vitro cultures of mouse embryos (EK) can be used as models for the early stages of normal mammalian development. A few known examples of experimental designs with such cells are reviewed: aggregates with normal embryos, promotion of parthenogenetic development by injection of EK cells into blastocysts, EK cells homozygous for a lethal gene, timing of expression in differentiating EC cells of a tissue-specific gene product, and X chromosome inactivation.     

Effects of radiation quality on interactions between oxidative stress, protein and DNA damage in Deinococcus radiodurans

Ionizing radiation damages DNA and also induces oxidative stress, which can affect the function of proteins involved in DNA repair, thereby causing repair of DNA damage to become less efficient. We previously developed a mathematical model of this potentially synergistic relationship and applied it to γ-ray exposure data on the radiation-resistant prokaryote Deinococcus radiodurans. Here, we investigate the effects of radiation quality on these processes by applying the model to data on exposures of D. radiodurans to heavy ions with linear energy transfer (LET) of 18.5–11,300 keV/μm. The model adequately describes these data using three parameters combinations: radiogenic DNA damage induction, repair protein inactivation and cellular repair capacity. Although statistical uncertainties around best-fit parameter estimates are substantial, the behaviors of model parameters are consistent with current knowledge of LET effects: inactivation cross-sections for both DNA and proteins increase with increasing LET; DNA damage yield per unit of radiation dose also increases with LET; protein damage per unit dose tends to decrease with LET; DNA and especially protein damage yields are reduced when cells are irradiated in the dry state. These results suggest that synergism between oxidative stress and DNA damage may play an important role not only during γ-ray exposure, but during high-LET radiation exposure as well.     

Inactivation and mutation of cultured mammalian cells by aluminium characteristic ultrasoft X-rays. II. Dose-responses of Chinese hamster and human diploid cells to aluminium X-rays and radiations of different LET.

The induction of inactivation and mutation to thioguanine-resistance of two types of cultured mammalian cells, V79 Chinese hamster and HF19 human diploid, was studied after irradiation with aluminium K characteristic ultrasoft X-rays, helium ion track intersections of different LET, 42 MeV d-Be neutrons, and hard X- or gamma-rays. The form of the dose-response curves was different for the two cell-types, and there was an overall difference in radiosensitivity, the human cells being the more sensitive to all radiations. However, for both inactivation and mutation-induction, the relative responses of both cell-types to these radiations was similar. Aluminium X-rays were considerably more effective than hard X- or gamma-rays and were at least as effective as helium ions of 20-28 keV micron-1, although aluminium X-rays produce tracks of very limited range (less than about 0.07 micron). Single track effects by aluminium X-rays cannot, therefore, extend beyond about 0.07 micron, and the subcellular sites involved in inactivation and mutation cannot be greater than this dimension or else the effectiveness of aluminium X-rays would be similar to that of low-LET radiations. This observation is in contradiction to models of radiation action which require relatively large sensitive sites; for example the 'theory of dual radiation action' requires a site diameter of about 0.4 micron to explain the shape of the dose-response curves for V79 hamster cells.