Modeling AIDS vaccines: the cellular level

This paper discuses current strategies for the development of AIDS vaccines wich allow immunzation to disturb the natural course of HIV at different detailed stages of its life cycle. Mathematical models describing the main biological phenomena (i.e. virus and vaccine induced T4 cell growth; virus and vaccine induced activation latently infected T4 cells; incremental changes immune response as infection progress; antibody dependent enhancement and neutralization of infection) and allowing for different vaccination strategies serve as a backgroud for computer simulations. The mathematical models reproduce updated information on the behavior of immune cells, antibody concentrations and free viruses. The results point to some controversial outcomes of an AIDS vaccine such as an early increase in virus concentration among vaccinated when compared to nonvaccinated individuals.

Neuropeptide Y at the cellular level – Studies on PreproNPY L7P Polymorphism and the Mitochondrial Form of NPY

Neuropeptide Y (NPY) is a widely expressed neurotransmitter in the central and peripheral nervous systems. Thymidine 1128 to cytocine substitution in the signal sequence of the preproNPY results in a single amino acid change where leucine is changed to proline. This L7P change leads to a conformational change of the signal sequence which can have an effect on the intracellular processing of NPY. The L7P polymorphism was originally associated with higher total and LDL cholesterol levels in obese subjects. It has also been associated with several other physiological and pathophysiological responses such as atherosclerosis and T2 diabetes. However, the changes on the cellular level due to the preproNPY signal sequence L7P polymorphism were not known. The aims of the current thesis were to study the effects of the [p.L7]+[p.L7] and the [p.L7]+[p.P7] genotypes in primary cultured and genotyped human umbilical vein endothelial cells (HUVEC), in neuroblastoma (SK-N-BE(2)) cells and in fibroblast (CHO-K1) cells. Also, the putative effects of the L7P polymorphism on proliferation, apoptosis and LDL and nitric oxide metabolism were investigated. In the course of the studies a fragment of NPY targeted to mitochondria was found. With the putative mitochondrial NPY fragment the aim was to study the translational preferences and the mobility of the protein. The intracellular distribution of NPY between the [p.L7]+[p.L7] and the [p.L7]+[p.P7] genotypes was found to be different. NPY immunoreactivity was prominent in the [p.L7]+[p.P7] cells while the proNPY immunoreactivity was prominent in the [p.L7]+[p.L7] genotype cells. In the proliferation experiments there was a difference in the [p.L7]+[p.L7] genotype cells between early and late passage (aged) cells; the proliferation was raised in the aged cells. NPY increased the growth of the cells with the [p.L7]+[p.P7] genotype. Apoptosis did not seem to differ between the genotypes, but in the aged cells with the [p.L7]+[p.L7] genotype, LDL uptake was found to be elevated. Furthermore, the genotype seemed to have a strong effect on the nitric oxide metabolism. The results indicated that the mobility of NPY protein inside the cells was increased within the P7 containing constructs. The existence of the mitochondria targeted NPY fragment was verified, and translational preferences were proved to be due to the origin of the cells. Cell of neuronal origin preferred the translation of mature NPY (NPY1-36) when compared to the non neuronal cells that translated both, NPY and the mitochondrial fragment of NPY. The mobility of the mitochondrial fragment was found to be minimal. The functionality of the mitochondrial NPY fragment remains to be investigated. L7P polymorphism in the preproNPY causes a series of intracellular changes. These changes may contribute to the state of cellular senescence, vascular tone and lead to endothelial dysfunction and even to increased susceptibility to diseases, like atherosclerosis and T2 diabetes.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

The photo-physiological response of a model cnidarian-dinoflagellate symbiosis to CO2-induced acidification at the cellular level

We measured the relationship between CO2-induced seawater acidification, photo-physiological performance and intracellular pH (pHi) in a model cnidarian-dinoflagellate symbiosis - the sea anemone Aiptasia sp. -under ambient (289.94 ± 12.54 µatm), intermediate (687.40 ± 25.10 µatm) and high (1459.92 ± 65.51 µatm) CO2 conditions. These treatments represented current CO2 levels, in addition to CO2 stabilisation scenarios IV and VI provided by the Intergovernmental Panel on Climate Change (IPCC). Anemones were exposed to each treatment for two months and sampled at regular intervals. At each time-point we measured a series of physiological responses: maximum dark-adapted fluorescent yield of PSII (Fv/Fm), gross photosynthetic rate, respiration rate, symbiont population density, and light-adapted pHi of both the dinoflagellate symbiont and isolated host anemone cell. We observed increases in all but one photo-physiological parameter (Pgross:R ratio). At the cellular level, increases in light-adapted symbiont pHi were observed under both intermediate and high CO2 treatments, relative to control conditions (pHi 7.35 and 7.46 versus pHi 7.25, respectively). The response of light-adapted host pHi was more complex, however, with no change observed under the intermediate CO2 treatment, but a 0.3 pH-unit increase under the high CO2 treatment (pHi 7.19 and 7.48, respectively). This difference is likely a result of a disproportionate increase in photosynthesis relative to respiration at the higher CO2 concentration. Our results suggest that, rather than causing cellular acidosis, the addition of CO2 will enhance photosynthetic performance, enabling both the symbiont and host cell to withstand predicted ocean acidification scenarios.

The Cellular Level of PR500, a Protein Complex Related to the 19S Regulatory Particle of the Proteasome, Is Regulated in Response to Stresses in Plants

In Arabidopsis seedlings and cauliflower florets, Rpn6 (a proteasome non-ATPase regulatory subunit) was found in two distinct protein complexes of ∼800 and 500 kDa, respectively. The large complex likely represents the proteasome 19S regulator particle (RP) because it displays the expected subunit composition and all characteristics. The small complex, designated PR500, shares at least three subunits with the “lid” subcomplex of 19S RP and is loosely associated with an hsp70 protein. In Arabidopsis COP9 signalosome mutants, PR500 was specifically absent or reduced to an extent that correlates with the severity of the mutations. Furthermore, PR500 was also diminished in response to potential protein-misfolding stresses caused by the heat shock and canavanine treatment. Immunofluorescence studies suggest that PR500 has a distinct localization pattern and is enriched in specific nuclear foci. We propose that PR500 may be evolved in higher plants to cope with the frequently encountered environmental stresses.

Low-Level Laser Intensity Application in Masseter Muscle for Treatment Purposes

Objective: This study evaluated with histochemical analysis how the number of laser applications can affect the masseter muscle. Background: In dentistry today, the laser is used in patients with temporomandibular disorders (TMDs), mainly for radiating pain in the masticatory muscles, whose origins may be associated with malocclusion, although the laser effects are not well understood on the cellular level. Materials and Methods: Thirty mice (HRS/J lineage) were randomly distributed into groups according to the number of laser applications (three, six, and 10). For each group of laser applications (experimental, n = 5), it was considered the control group (n = 5), which was not irradiated. All animals inhaled halothane (2-bromo-2-chloro-1, 1, 1-trifluoroethane, minimum 99%, Sigma Aldrich, India) before each laser irradiation performed on the left masseter muscle region, on alternate days with 20 J/cm(2), 40mW, for 20 sec. The muscle samples were collected for histochemical analysis with succinate dehydrogenase (SDH) enzyme 72 h after the last application. Results: (a) A decrease in area of light fibers type (35.91% +/- 6.9%; 32.08% +/- 6.3%, and 27.88% +/- 6.3%), according to the increase of laser applications (p < 0.05); (b) significant increase (p < 0.05) in the area of intermediate fibers, with an increase of laser application (11.08% +/- 3.9%; 16.52% +/- 5.7%, and 15.96% +/- 3.9%), although the increase with 10 applications was small; (c) area increase of dark fibers in the group with three laser applications (0.16% +/- 0.3%) (p < 0.05), and in groups with six and 10 laser applications, respectively (9.68% +/- 6.0% and 9.60% +/- 4.0%). Conclusions: The SDH enzyme activity revealed that the number of laser applications increases the metabolic pattern of the muscle fibers. A minimal difference in metabolic activity between six and 10 applications of a laser suggests that further analyses should be done to confirm that six applications are enough to produce the same clinical effects, thereby contributing data to professionals from different fields in regard to the cost-benefit ratio of this therapy.

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.

Nano-electro-mechanical in dynamical cellular processes

Report for the scientific sojourn carried out at the Department of Structure and Constituents of Matter during 2007.The main focus of the work was on phenomena related to nano-electromechanical processes that take place on a cellular level. Additionally, it has also been performed independent work related to charge and energy transfer in bio molecules, energy transfer in coupled spin systems as well as electrodynamics of nonlinear metamaterials.

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.

System-level characterization of Th2 cell development and immune cell responses to ZnO and TiO2 nanoparticles

Asthma and allergy are common diseases and their prevalence is increasing. One of the hypotheses that explains this trend is exposure to inhalable chemicals such as traffi c-related air pollution. Epidemiological research supports this theory, as a correlation between environmental chemicals and allergic respiratory diseases has been found. In addition to ambient airborne particles, one may be exposed to engineered nanosized materials that are actively produced due to their favorable physico-chemical properties compared to their bulk size counterparts. On the cellular level, improper activity of T helper (Th) cells has been connected to allergic reactions. Th cells can differentiate into functionally different effector subsets, which are identifi ed according to their characteristic cytokine profi les resulting in specifi c ability to communicate with other cells. Th2 cells activate humoral immunity and stimulate eradication of extracellular pathogens. However, persistent predominance of Th2 cells is involved in a development of number of allergic diseases. The cytokine environment at the time of antigen recognition is the major factor determining the polarization of a naïve Th cell. Th2 cell differentiation is initiated by IL4, which signals via transcription factor STAT6. Although the importance of this pathway has been evaluated in the mouse studies, the signaling components involved have been largely unknown. The aim of this thesis was to identify molecules, which are under the control of IL4 and STAT6 in Th cells. This was done by using system-level analysis of STAT6 target genes at genome, mRNA and protein level resulting in identifi cation of various genes previously not connected to Th2 cell phenotype acquisition. In the study, STAT6-mediated primary and secondary target genes were dissection from each other and a detailed transcriptional kinetics of Th2 cell polarization of naïve human CD4+ T cells was collected. Integration of these data revealed the hierarchy of molecular events that mediates the differentiation towards Th2 cell phenotype. In addition, the results highlighted the importance of exploiting proteomics tools to complement the studies on STAT6 target genes identifi ed through transcriptional profi ling. In the last subproject, the effects of the exposure with ZnO and TiO2 nanoparticles was analyzed in Jurkat T cell line and in primary human monocyte-derived macrophages and dendritic cells to evaluate their toxicity and potential to cause infl ammation. Identifi cation of ZnO-derived gene expression showed that the same nanoparticles may elicit markedly distinctive responses in different cell types, thus underscoring the need for unbiased profi ling of target genes and pathways affected. The results gave additional proof that the cellular response to nanosized ZnO is due to leached Zn2+ ions. The approach used in ZnO and TiO2 nanoparticle study demonstrated the value of assessing nanoparticle responses through a toxicogenomics approach. The increased knowledge of Th2 cell signaling will hopefully reveal new therapeutic nodes and eventually improve our possibilities to prevent and tackle allergic infl ammatory diseases.

A 2-dimensional computational model to analyze the effects of cellular heterogeinity on cardiac pacemaking

The mechanical action of the heart is made possible in response to electrical events that involve the cardiac cells, a property that classifies the heart tissue between the excitable tissues. At the cellular level, the electrical event is the signal that triggers the mechanical contraction, inducing a transient increase in intracellular calcium which, in turn, carries the message of contraction to the contractile proteins of the cell. The primary goal of my project was to implement in CUDA (Compute Unified Device Architecture, an hardware architecture for parallel processing created by NVIDIA) a tissue model of the rabbit sinoatrial node to evaluate the heterogeneity of its structure and how that variability influences the behavior of the cells. In particular, each cell has an intrinsic discharge frequency, thus different from that of every other cell of the tissue and it is interesting to study the process of synchronization of the cells and look at the value of the last discharge frequency if they synchronized.

The role of the D2 dopamine receptor (D2R) in A2A adenosine receptor (A2AR)-mediated behavioral and cellular responses as revealed by A2A and D2 receptor knockout mice

The A2AR is largely coexpressed with D2Rs and enkephalin mRNA in the striatum where it modulates dopaminergic activity. Activation of the A2AR antagonizes D2R-mediated behavioral and neurochemical effects in the basal ganglia through a mechanism that may involve direct A2AR–D2R interaction. However, whether the D2R is required for the A2AR to exert its neural function is an open question. In this study, we examined the role of D2Rs in A2AR-induced behavioral and cellular responses, by using genetic knockout (KO) models (mice deficient in A2ARs or D2Rs or both). Behavioral analysis shows that the A2AR agonist 2–4-(2-carboxyethyl)phenethylamino-5′-N-ethylcarboxamidoadenosine reduced spontaneous as well as amphetamine-induced locomotion in both D2 KO and wild-type mice. Conversely, the nonselective adenosine antagonist caffeine and the A2AR antagonist 8-(3-chlorostyryl)caffeine produced motor stimulation in mice lacking the D2R, although the stimulation was significantly attentuated. At the cellular level, A2AR inactivation counteracted the increase in enkephalin expression in striatopallidal neurons caused by D2R deficiency. Consistent with the D2 KO phenotype, A2AR inactivation partially reversed both acute D2R antagonist (haloperidol)-induced catalepsy and chronic haloperidol-induced enkephalin mRNA expression. Together, these results demonstrate that A2ARs elicit behavioral and cellular responses despite either the genetic deficiency or pharmacological blockade of D2Rs. Thus, A2AR-mediated neural functions are partially independent of D2Rs. Moreover, endogenous adenosine acting at striatal A2ARs may be most accurately viewed as a facilitative modulator of striatal neuronal activity rather than simply as an inhibitory modulator of D2R neurotransmission.