Response of Tumor Spheroids to Radiation: Modeling and Parameter Estimation

We propose a spatially distributed continuous model for the spheroid response to radiation, in which the oxygen distribution is represented by means of a diffusion-consumption equation and the radiosensitivity parameters depend on the oxygen concentration. The induction of lethally damaged cells by a pulse of radiation, their death, and the degradation of dead cells are included. The compartments of lethally damaged cells and of dead cells are subdivided into different subcompartments to simulate the delays that occur in cell death and cell degradation, with a gain in model flexibility. It is shown that, for a single irradiation and under the hypothesis of a sufficiently small spheroid radius, the model can be reformulated as a linear stationary ordinary differential equation system. For this system, the parameter identifiability has been investigated, showing that the set of unknown parameters can be univocally identified by exploiting the response of the model to at least two different radiation doses. Experimental data from spheroids originated from different cell lines are used to identify the unknown parameters and to test the predictive capability of the model with satisfactory results.     

Influence of Polyols on the Stability and Kinetic Parameters of Invertase from Candida utilis: Correlation with the Conformational Stability and Activity

Invertase (β-d-fructofuranoside fructohydrolase-E.C. 3.2.1.26) is a sucrose hydrolyzing enzyme found in microbial, plant and animal sources. Invertase from Candida utilis is a dimeric glycoprotein composed of two identical monomer subunits with an apparent molecular mass of 150 kDa. We investigated the mechanism of stabilization of invertase with polyols (glycerol, xylitol, and sorbitol). Activity, thermodynamic and kinetic measurements of invertase were performed as a function of polyol concentration and showed that polyols act as very effective stabilizing agents. The result from the solvent-invertase interaction shows preferential exclusion of the polyols from the protein domain leading to preferential hydration of protein. Apparent thermal denaturation temperature of the protein (T _m ) rose from 75 °C to a maximum of 85 °C in 30% glycerol. The stabilization has been attributed to the preferential hydration of the enzyme.     

Modeling the Dynamics of Viral Evolution Considering Competition Within Individual Hosts and at Population Level: The Effects of Treatment

We consider two viral strains competing against each other within individual hosts (at cellular level) and at population level (for infecting hosts) by studying two cases. In the first case, the strains do not mutate into each other. In this case, we found that each individual in the population can be infected by only one strain and that co-existence in the population is possible only when the strain that has the greater basic intracellular reproduction number, R _0c , has the smaller population number R _0p . Treatment against the one strain shifts the population equilibrium toward the other strain in a complicated way (see Appendix B). In the second case, we assume that the strain that has the greater intracellular number R _0c can mutate into the other strain. In this case, individual hosts can be simultaneously infected by both strains (co-existence within the host). Treatment shifts the prevalence of the two strains within the hosts, depending on the mortality induced by the treatment, which is, in turn, dependent upon the doses given to each individual. The relative proportions of the strains at the population level, under treatment, depend both on the relative proportions within the hosts (which is determined by the dosage of treatment) and on the number of individuals treated per unit time, that is, the rate of treatment. Implications for cases of real diseases are briefly discussed.     

An Estimation of the Manganese Requirement for Broilers from 1 to 21 Days of Age

A dose–response experiment was conducted to find the sensitive and consistent biomarker for the estimation of dietary manganese (Mn) requirement and establish the optimal Mn level for broilers fed a practical corn–soybean meal diet from 1 to 21 days of age post-hatching. A total of 480 1-day-old Arbor Acres male chicks were randomly allotted to one of eight treatments with five replicates of 12 birds each and fed diets supplemented with 0, 20, 40, 60, 80, 100, 120, or 140 mg Mn/kg from reagent grade Mn sulfate. Tissue Mn concentrations, manganese-containing superoxide dismutase (MnSOD) activity, and MnSOD mRNA concentration within heart tissue were analyzed at 7, 14, and 21 days of age. Tissue Mn concentrations and heart MnSOD activity showed significant quadratic responses, and heart MnSOD mRNA concentration showed an asymptotic response to dietary supplemental Mn level, respectively. The estimate of dietary Mn for chicks from 1 to 21 days of age was 122–128 for heart Mn concentration, 141–159 for pancreas Mn concentration, 127–138 for liver Mn concentration, and 135–156 mg/kg for heart MnSOD activity, respectively. Heart MnSOD mRNA concentration was a consistent index for the estimation of the Mn requirement of broilers. Based on this index, the estimate of dietary Mn requirement for broilers from 1 to 21 days of age post-hatching was about 130 mg/kg, which was a little more than two times of the current NRC (1994) requirement.     

Pressure and Temperature Dependence of Growth and Morphology of?Escherichia coli: Experiments and Stochastic Model

We have investigated the growth of Escherichia coli, a mesophilic bacterium, as a function of pressure (P) and temperature (T). Escherichia coli can grow and divide in a wide range of pressure (1–400 atm) and temperature (23–40°C). For T > 30°C, the doubling time of E. coli increases exponentially with pressure and exhibits a departure from exponential behavior at pressures between 250 and 400 atm for all the temperatures studied in our experiments. The sharp change in doubling time is followed by a sharp change in phenotypic transition of E. coli at high pressures where bacterial cells switch to an elongating cell type. We propose a model that this phenotypic change in bacteria at high pressures is an irreversible stochastic process, whereas the switching probability to elongating cell type increases with increasing pressure. The model fits well the experimental data. We discuss our experimental results in the light of structural and thus functional changes in proteins and membranes.     

Kinetics of Extracellular ATP from Goldfish Hepatocytes: A Lesson from Mathematical Modeling

In goldfish hepatocytes, hypotonic exposure leads to cell swelling, followed by a compensatory shrinkage termed RVD. It has been previously shown that ATP is accumulated in the extracellular medium of swollen cells in a non-linear fashion, and that extracellular ATP (ATPe) is an essential intermediate to trigger RVD. Thus, to understand how RVD proceeds in goldfish hepatocytes, we developed two mathematical models accounting for the experimental ATPe kinetics reported recently by Pafundo et al. in Am. J. Physiol. 294, R220–R233, 2008. Four different equations for ATPe fluxes were built to account for the release of ATP by lytic (J _L ) and nonlytic mechanisms (J _NL ), ATPe diffusion (J _D ), and ATPe consumption by ectonucleotidases (J _V ). Particular focus was given to J _NL , defined as the product of a time function (J _R ) and a positive feedback mechanism whereby ATPe amplifies J _NL . Several J _R functions (Constant, Step, Impulse, Gaussian, and Lognormal) were studied. Models were tested without (model 1) or with (model 2) diffusion of ATPe. Mathematical analysis allowed us to get a general expression for each of the models. Subsequently, by using model dependent fit (simulations) as well as model analysis at infinite time, we observed that:

Native-State Heterogeneity of β2-Microglobulin as Revealed by Kinetic Folding and Real-Time NMR Experiments

The kinetic folding of β₂-microglobulin from the acid-denatured state was investigated by interrupted-unfolding and interrupted-refolding experiments using stopped-flow double-jump techniques. In the interrupted unfolding, we first unfolded the protein by a pH jump from pH 7.5 to pH 2.0, and the kinetic refolding assay was carried out by the reverse pH jump by monitoring tryptophan fluorescence. Similarly, in the interrupted refolding, we first refolded the protein by a pH jump from pH 2.0 to pH 7.5 and used a guanidine hydrochloride (GdnHCl) concentration jump as well as the reverse pH jump as unfolding assays. Based on these experiments, the folding is represented by a parallel-pathway model, in which the molecule with the correct Pro32 cis isomer refolds rapidly with a rate constant of 5–6 s^? 1, while the molecule with the Pro32 trans isomer refolds more slowly (pH 7.5 and 25 °C). At the last step of folding, the native-like trans conformer produced on the latter pathway isomerizes very slowly (0.001–0.002 s^? 1) into the native cis conformer. In the GdnHCl-induced unfolding assays in the interrupted refolding, the native-like trans conformer unfolded remarkably faster than the native cis conformer, and the direct GdnHCl-induced unfolding was also biphasic, indicating that the native-like trans conformer is populated at a significant level under the native condition. The one-dimensional NMR and the real-time NMR experiments of refolding further indicated that the population of the trans conformer increases up to 7–9% under a more physiological condition (pH 7.5 and 37 °C).     

Kinetic properties and substrate inhibition of α-galactosidase from Aspergillus niger

The recombinant AglB produced by Pichia pastoris exhibited substrate inhibition behavior for the hydrolysis of p-nitrophenyl α-galactoside, whereas it hydrolyzed the natural substrates, including galactomanno-oligosaccharides and raffinose family oligosaccharides, according to the Michaelian kinetics. These contrasting kinetic behaviors can be attributed to the difference in the dissociation constant of second substrate from the enzyme and/or to the ability of the leaving group of the substrates. The enzyme displays the grater k_cat/K_m values for hydrolysis of the branched α-galactoside in galactomanno-oligosaccharides than that of raffinose and stachyose. A sequence comparison suggested that AglB had a shallow active-site pocket, and it can allow to hydrolyze the branched α-galactosides, but not linear raffinose family oligosaccharides.     

Improved Horvitz–Thompson Estimation of Model Parameters from Two-phase Stratified Samples: Applications in Epidemiology

The case-cohort study involves two-phase samplings: simple random sampling from an infinite superpopulation at phase one and stratified random sampling from a finite cohort at phase two. Standard analyses of case-cohort data involve solution of inverse probability weighted (IPW) estimating equations, with weights determined by the known phase two sampling fractions. The variance of parameter estimates in (semi)parametric models, including the Cox model, is the sum of two terms: (i) the model-based variance of the usual estimates that would be calculated if full data were available for the entire cohort; and (ii) the design-based variance from IPW estimation of the unknown cohort total of the efficient influence function (IF) contributions. This second variance component may be reduced by adjusting the sampling weights, either by calibration to known cohort totals of auxiliary variables correlated with the IF contributions or by their estimation using these same auxiliary variables. Both adjustment methods are implemented in the R survey package. We derive the limit laws of coefficients estimated using adjusted weights. The asymptotic results suggest practical methods for construction of auxiliary variables that are evaluated by simulation of case-cohort samples from the National Wilms Tumor Study and by log-linear modeling of case-cohort data from the Atherosclerosis Risk in Communities Study. Although not semiparametric efficient, estimators based on adjusted weights may come close to achieving full efficiency within the class of augmented IPW estimators.     

Estimation of Interaction Between Human Keratinocytes and Xenogenic Collagen in vitro

This paper describes the interaction observed between human keratinocytes and xenogenic collagen in vitro modified by HCl. Human keratinocytes were cultivated for 3–10 days, on modified and control support. Their growth, morphology and interaction with support were analyzed. It was found that on both control and experimental (modified) collagen cells proliferated in a similar way. Within 3–10 days, the culture became multilayered and mature and differentiation of cells was visible. Using electron microscope elements of basal membrane interacting with support were seen. On modified support processes of cells penetrating the support are occasionally seen. By use of the immunofluorescent, cytochemical techniques was found the presence of: BP-180 (antigen), β₄ integrin, laminin 5 and collagen IV, VII, VIIc. On the modified support the above listed elements appeared between 3 and 7 days of culture, whereas on the control between 7th and 10th days. On 10th day of culture, the presence of elements of basal membranes became less evident. Results give some hope for using xenogenic, modified collagen as support of keratinocytes culture in process of human skin engineering.     

Heterosexual Transmission of Subtype C HIV-1 Selects Consensus-Like Variants without Increased Replicative Capacity or Interferon-α Resistance

Heterosexual transmission of HIV-1 is characterized by a genetic bottleneck that selects a single viral variant, the transmitted/founder (TF), during most transmission events. To assess viral characteristics influencing HIV-1 transmission, we sequenced 167 near full-length viral genomes and generated 40 infectious molecular clones (IMC) including TF variants and multiple non-transmitted (NT) HIV-1 subtype C variants from six linked heterosexual transmission pairs near the time of transmission. Consensus-like genomes sensitive to donor antibodies were selected for during transmission in these six transmission pairs. However, TF variants did not demonstrate increased viral fitness in terms of particle infectivity or viral replicative capacity in activated peripheral blood mononuclear cells (PBMC) and monocyte-derived dendritic cells (MDDC). In addition, resistance of the TF variant to the antiviral effects of interferon-α (IFN-α) was not significantly different from that of non-transmitted variants from the same transmission pair. Thus neither in vitro viral replicative capacity nor IFN-α resistance discriminated the transmission potential of viruses in the quasispecies of these chronically infected individuals. However, our findings support the hypothesis that within-host evolution of HIV-1 in response to adaptive immune responses reduces viral transmission potential.     

197 Modeling the maturation pathway of the HIV-1 5’-UTR RNA dimer

HIV-1 genomic RNA is packaged as a dimer into the virions. The initial metastable RNA dimer is believed to be formed by virtue of “kissing interactions” between two copies of the palindromic apical loops of stem-loop SL1 of the 5’-untranslated region (5’-UTR) of the genomic RNA. Viral nucleocapsid protein NCp7 promotes maturation of the RNA dimer into more stable form, which involves extended or linear form of SL1 dimer (reviewed in Paillart et al., 2004; Moore & Hu, 2009; Lu et al., 2011). In vitro experiments have shown that this conversion occurs at stoichiometric amounts of NCp7 without breaking interactions between the two copies of the SL1 apical loops (Mujeeb et al., 2007). We have proposed a hypothetical pathway and calculated models of the intermediate structures for the SL1 stem-loop dimer maturation that does not require simultaneous dissociation of all base pairs in SL1 stems; this pathway involves formation of an RNA analog of the Holliday junction intermediate between the two stems of the SL1 dimer and a following branch migration towards the palindromic duplex (Ulyanov et al., 2011). Here, we extend these models to the dimer of the 1–344 fragment of HIV-1 RNA, which includes all of the 5’-UTR and the gag start AUG codon region, and show that the branch-migration mechanism of the dimer maturation is also feasible for the full 5’-UTR RNA. All RNA models have been calculated with the miniCarlo program (Zhurkin et al., 1991).     

Serological identification and molecular characterization of B (A) 02 subtype in patients and blood donors from Eastern China

This study was designed to identify the rare?type?ABO?blood?groups, B(A) 02, from Eastern China. Three samples with discordant serological results during routine blood type identification and four samples from one sample’ family were selected. All of them were detected by serological method. The exon 6 and 7 of the ABO genes were amplified by PCR and sequenced. They were typed as AsubB by serology and as BO by genotype. In AsubB samples, nt 700C>G mutation was detected in B gene, which was previously defined as B(A)02 alleles. In these seven samples, six showed B(A)02/O01 and one showed B(A)02/O02.B(A)02 allele was found to be more common in this study than B(A)04 which is considered to be more frequent than B(A)02. The careful identification of rare blood types is important for the safety of clinical blood transfusion.     

Growth kinetics of 1–2 mm and 3–4 mm colonies of Nostoc sphaeroides (Cyanophyta) in outdoor culture

Nostoc sphaeroides Kützing was cultivated in paddlewheel-driven raceway ponds and the growth kinetics of 1-2 mm and 3-4 mm colonies of N. sphaeroides was studied. The biomass productivities in 2.5 m(2) raceway ponds inoculated with 1-2 mm and 3-4 mm colonies were 5.2 and 0.25 g dry wt m(-2) d(-1), respectively. Furthermore, differently sized colonies showed different relative water content, total soluble carbohydrates, chlorophyll a content and density of filaments. This is the first report on mass culture of N. sphaeroides under outdoor conditions.     

Effect of Silt, Water and Periphyton Quality on Survival and Growth of the Mayfly Heptagenia Sulphurea

The herbivorous mayfly Heptagenia sulphurea is characteristic of rivers with stony bottoms. Records from the 20th century showed that this species had disappeared from the Common Meuse in the Netherlands, probably due to river regulation or changes in water quality. A field survey in 2003 showed that H. sulphurea was present in the Geul tributary, approximately 300 m upstream of its confluence with the Common Meuse. H. sulphurea has not recolonized the Common Meuse despite improvements in water quality over the last decades. Concentration of suspended sediments in the River Meuse, however, is still high, much higher than in the beginning of the 20th century. The presence of a silt layer may limit the return of H. sulphurea in this river by reducing availability and quality of its food. The prime objective of this study was to investigate the impact of silt on survival and growth of H. sulphurea in a laboratory experiment. Furthermore, the impact of water and periphyton quality from the Common Meuse on survival and growth of this mayfly was also investigated. Results showed that neither water quality nor cultured periphyton from the Common Meuse reduced survival and growth of H. sulphurea. The presence of a silt layer resulted in a significantly lower growth rate of the mayfly larvae. It is concluded that the silt layer reduces the accessibility of the food. Covering of food is possibly one of the main factors limiting the recolonization of H. sulphurea and probably other benthic grazers in the Common Meuse.     

Pre-steady-state Kinetic and Structural Analysis of Interaction of Methionine γ-Lyase from Citrobacter freundii with Inhibitors

Methionine γ-lyase (MGL) catalyzes the γ-elimination of l-methionine and its derivatives as well as the β-elimination of l-cysteine and its analogs. These reactions yield α-keto acids and thiols. The mechanism of chemical conversion of amino acids includes numerous reaction intermediates. The detailed analysis of MGL interaction with glycine, l-alanine, l-norvaline, and l-cycloserine was performed by pre-steady-state stopped-flow kinetics. The structure of side chains of the amino acids is important both for their binding with enzyme and for the stability of the external aldimine and ketimine intermediates. X-ray structure of the MGL·l-cycloserine complex has been solved at 1.6 Å resolution. The structure models the ketimine intermediate of physiological reaction. The results elucidate the mechanisms of the intermediate interconversion at the stages of external aldimine and ketimine formation.     

Modeling and Analysis of the Effect of Training on \dot{V}O_{2} Kinetics and Anaerobic Capacity

In this paper, we present an application of a number of tools and concepts for modeling and analyzing raw, unaveraged, and unedited breath-by-breath oxygen uptake data. A method for calculating anaerobic capacity is used together with a model, in the form of a set of coupled nonlinear ordinary differential equations to make predictions of the kinetics, the time to achieve a percentage of a certain constant oxygen demand, and the time limit to exhaustion at intensities other than those in which we have data. Speeded oxygen kinetics and increased time limit to exhaustion are also investigated using the eigenvalues of the fixed points of our model. We also use a way of analyzing the oxygen uptake kinetics using a plot of vs which allows one to observe both the fixed point solutions and also the presence of speeded oxygen kinetics following training. A method of plotting the eigenvalue versus oxygen demand is also used which allows one to observe where the maximum amplitude of the so-called slow component will be and also how training has changed the oxygen uptake kinetics by changing the strength of the attracting fixed point for a particular demand.     

Efficient Generation of CA3 Neurons from Human Pluripotent Stem Cells Enables Modeling of Hippocampal Connectivity In Vitro

Download video (9MB)Help with mp4 files