Statistical analysis and optimal classification of blood cell populations using Gaussian distributions

This dissertation develops a new figure of merit to measure the similarity (or dissimilarity) of Gaussian distributions through a novel concept that relates the Fisher distance to the percentage of data overlap. The derivations are expanded to provide a generalized mathematical platform for determining an optimal separating boundary of Gaussian distributions in multiple dimensions. Real-world data used for implementation and in carrying out feasibility studies were provided by Beckman-Coulter. It is noted that although the data used is flow cytometric in nature, the mathematics are general in their derivation to include other types of data as long as their statistical behavior approximate Gaussian distributions. ^ Because this new figure of merit is heavily based on the statistical nature of the data, a new filtering technique is introduced to accommodate for the accumulation process involved with histogram data. When data is accumulated into a frequency histogram, the data is inherently smoothed in a linear fashion, since an averaging effect is taking place as the histogram is generated. This new filtering scheme addresses data that is accumulated in the uneven resolution of the channels of the frequency histogram. ^ The qualitative interpretation of flow cytometric data is currently a time consuming and imprecise method for evaluating histogram data. This method offers a broader spectrum of capabilities in the analysis of histograms, since the figure of merit derived in this dissertation integrates within its mathematics both a measure of similarity and the percentage of overlap between the distributions under analysis. ^

Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

Chemical dictionary, American and British usage. Containing the words generally used in chemistry, and many of the terms used in the related sciences of physics, astrophysics, mineralogy, pharmacy, agriculture, biology, medicine, engineering, etc.; based on recent chemical literature.

Re-issue, with changes and additions of the 3d ed. pub. in 1944.

Diffusion tensor imaging of rat spinal cord in vivo

Magnetic resonance imaging, with its exquisite soft tissue contrast, is an ideal modality for investigating spinal cord pathology. While conventional MRI techniques are very sensitive for spinal cord pathology, their specificity is somewhat limited. Diffusion MRI is an advanced technique which is a very sensitive and specific indicator of the integrity of white matter tracts. Diffusion imaging has been shown to detect early ischemic changes in white matter, while conventional imaging demonstrates no change. By acquiring the complete apparent diffusion tensor (ADT), tissue diffusion properties can be expressed in terms of quantitative and rotationally invariant parameters. ^ Systematic study of SCI in vivo requires controlled animal models such as the popular rat model. To date, studies of spinal cord using ADT imaging have been performed exclusively in fixed, excised spinal cords, introducing inevitable artifacts and losing the benefits of MRI's noninvasive nature. In vivo imaging reflects the actual in vivo tissue properties, and allows each animal to be imaged at multiple time points, greatly reducing the number of animals required to achieve statistical significance. Because the spinal cord is very small, the available signal-to-noise ratio (SNR) is very low. Prior spin-echo based ADT studies of rat spinal cord have relied on high magnetic field strengths and long imaging times—on the order of 10 hours—for adequate SNR. Such long imaging times are incompatible with in vivo imaging, and are not relevant for imaging the early phases following SCI. Echo planar imaging (EPI) is one of the fastest imaging methods, and is popular for diffusion imaging. However, EPI further lowers the image SNR, and is very sensitive to small imperfections in the magnetic field, such as those introduced by the bony spine. Additionally, The small field-of-view (FOV) needed for spinal cord imaging requires large imaging gradients which generate EPI artifacts. The addition of diffusion gradients introduces yet further artifacts. ^ This work develops a method for rapid EPI-based in vivo diffusion imaging of rat spinal cord. The method involves improving the SNR using an implantable coil; reducing magnetic field inhomogeneities by means of an autoshim, and correcting EPI artifacts by post-processing. New EPI artifacts due to diffusion gradients described, and post-processing correction techniques are developed. ^ These techniques were used to obtain rotationally invariant diffusion parameters from 9 animals in vivo, and were validated using the gold-standard, but slow, spinecho based diffusion sequence. These are the first reported measurements of the ADT in spinal cord in vivo . ^ Many of the techniques described are equally applicable toward imaging of human spinal cord. We anticipate that these techniques will aid in evaluating and optimizing potential therapies, and will lead to improved patient care. ^

4HPR-X radiation mechanisms of interactions in non-small cell lung cancer (NSCLC) cells in vitro

Lung cancer is the leading cause of cancer death. However, poor survival using conventional therapies fuel the search for more rational interventions. The objective of this study was to design and implement a 4HPR-radiation interaction model in NSCLC, employing a traditional clinical modality (radiation), a relatively new, therapeutically unexplored agent (4HPR) and rationally combining them based on molecular mechanistic findings pertaining to their interactions. To test the hypothesis that 4HPR sensitizes cells to radiation-induced cell death via G2+M accumulation, we designed a working model consisting of H522 adenocarcinoma cells (p53, K-ras mutated) derived from an NSCLC patient; 4HPR at concentrations up to 10 μM; and X radiation up to 6 Gy generated by a patient-dedicated Phillips RT-250 X ray unit at 250 KV, 15 mA, 1.85 Gy/min. We found that 4HPR produced time- and dose-dependent morphological changes, growth inhibition, and DNA damage-inducing enhancement of reactive oxygen species. A transient G2+M accumulation of cells maximal at 24 h of continuous 4HPR exposure was used for irradiation time scheduling. Our data demonstrated enhanced cell death (both apoptotic and necrotic) in irradiated cells pre-treated with 4HPR versus those with either stressor alone. 4HPR's effect of increased NSCLC cells' radioresponse was confirmed by clonogenic assay. To explore these practical findings from a molecular mechanistic perspective, we further investigated and showed that levels of cyclin B1 and p34cdc2 kinase—both components of the mitosis promoting factor (MPF) regulating the G2/M transition—did not change following 4HPR treatment. Likewise, cdc25C phosphatase was not altered. However, enhanced p34cdc2 phosphorylation on its Thr14Tyr15 residues—indicative of its inactivation and increased expression of MPF negative regulators chk1 and wee1 kinases—were supportive of explaining 4HPR-treated cells' accumulation. Hence, p34cdc2 phosphorylation, chk1, and wee1 warrant further evaluation as potential molecular targets for 4HPR-X radiation combination. In summary, we (1) demonstrated that 4HPR not only induces cell death by itself, but also increases NSCLC cells' subsequent radioresponse, indicative of potential clinical applicability, and (2) for the first time, shed light on deciphering 4HPR-X radiation molecular mechanisms of interaction, including the finding of 4HPR's role as a p34cdc2 inactivator via Thr14Tyr15 phosphorylation. ^

Correction of aberrant FGFR1 RNA splicing through blocking intronic regulatory elements

Analysis of the human genome has revealed that more than 74% of human genes undergo alternative RNA splicing. Aberrations in alternative RNA splicing have been associated with several human disorders, including cancer. ^ We studied the aberrant expression of alternative RNA splicing isoforms of the Fibroblast Growth Factor Receptor 1 (FGFR1) gene in a human glioblastoma cancer model. Normal glial cells express the FGFR1α, which contains three extracellular domains. In tumors the most abundant isoform is the FGFR1β, which lacks the first extracellular domain due to the skipping of a single exon, termed alpha. The skipping of the α-exon is regulated by two intronic silencing sequences within the precursor mRNA. Since we observed no mutations on these elements in tumor cells, we hypothesized that the over-expression of regulatory proteins that recognize these sequences is responsible for the aberrant expression of splicing isoforms. Hence, we blocked the formation of protein complexes on the ISS using antisense RNA oligonucleotides in vitro. We also evaluated the impact of the ISS antisense oligonucleotides on the endogenous FGFR1 splicing, in a glioblastoma cell model. By targeting intronic regulatory elements we were able to increase the level of alpha exon inclusion up to 90% in glioblastoma cells. The effect was dose dependent, sequence specific and reproducible in glioblastoma and other cancer cells, which also exhibit an alpha exon skipping phenotype. Targeting FGFR1 endogenous ISS1 and ISS2 sequences did not have an additive or synergistic effect, which suggest a regulatory splicing mechanism that requires the interaction of complexes formed on these elements. An increase in the levels of the FGFR1α isoform resulted in a reduction in cell invasiveness. Also, a significant increase in the levels of caspase 3/7 activities, which is indicative of an elevation in apoptosis levels, suggests that expression of FGFR1β might be relevant for tumor survival. These studies demonstrate that it is possible to prevent aberrant expression of exon skipping events through the targeting of intronic regulatory elements, providing an important new therapeutic tool for the correction of human disease caused by alternative RNA splicing. ^

Respiratory rhythm generation: A modeling study of the respiratory central pattern generator and the phrenic motor neuron

The respiratory central pattern generator is a collection of medullary neurons that generates the rhythm of respiration. The respiratory central pattern generator feeds phrenic motor neurons, which, in turn, drive the main muscle of respiration, the diaphragm. The purpose of this thesis is to understand the neural control of respiration through mathematical models of the respiratory central pattern generator and phrenic motor neurons. ^ We first designed and validated a Hodgkin-Huxley type model that mimics the behavior of phrenic motor neurons under a wide range of electrical and pharmacological perturbations. This model was constrained physiological data from the literature. Next, we designed and validated a model of the respiratory central pattern generator by connecting four Hodgkin-Huxley type models of medullary respiratory neurons in a mutually inhibitory network. This network was in turn driven by a simple model of an endogenously bursting neuron, which acted as the pacemaker for the respiratory central pattern generator. Finally, the respiratory central pattern generator and phrenic motor neuron models were connected and their interactions studied. ^ Our study of the models has provided a number of insights into the behavior of the respiratory central pattern generator and phrenic motor neurons. These include the suggestion of a role for the T-type and N-type calcium channels during single spikes and repetitive firing in phrenic motor neurons, as well as a better understanding of network properties underlying respiratory rhythm generation. We also utilized an existing model of lung mechanics to study the interactions between the respiratory central pattern generator and ventilation. ^

p21 activated kinase 5 links multiple GTPase signaling pathways at mitochondria

The p21-activated kinase 5 (PAK5) is a serine/threonine protein kinase associated with the group 2 subfamily of PAKs. Although our understanding about PAK5 is very limited, it is receiving increasing interest due to its tissue specific expression pattern and important signaling properties. PAK5 is highly expressed in brain. Its overexpression induces neurite outgrowth in neuroblastoma cells and promotes survival in fibroblasts. ^ The serine/threonine protein kinase Raf-1 is an essential mediator of Ras-dependent signaling that controls the ERK/MAPK pathway. In contrast to PAK5, Raf-1 has been the subject of intensive investigation. However due to the complexity of its activation mechanism, the biological inputs controlling Raf-1 activation are not fully understood. ^ PAKs 1-3 are the known kinases responsible for phosphorylation of Raf-1 on serine 338, which is a crucial phosphorylation site for Raf-1 activation. However, dominant negative versions of these kinases do not block EGF-induced Raf-1 activation, indicating that other kinases may regulate the phosphorylation of Raf-1 on serine 338. ^ This thesis work was initiated to test whether the group 2 PAKs 4, 5 and 6 are responsible for EGF-induced Raf-1 activation. We found that PAK5, and to a lesser extent PAK4, can activate Raf-1 in cells. Our studies thereafter focused on PAK5. With the progress of our study we found that PAK5 does not significantly stimulate serine 338 phosphorylation of Triton X-100 soluble Raf-1. PAK5, however, constitutively and specifically associates with Raf-1 and targets it to a Triton X-100 insoluble, mitochondrial compartment, where PAK5 phosphorylates serine 338 of Raf-1. We further demonstrated that endogenous PAK5 and Raf-1 colocalize in Hela cells at the mitochondrial outer membrane. In addition, we found that the mitochondria-targeting of PAK5 is determined by its C-terminal kinase domain plus the upstream proximal region, and facilitated by the N-terminal p21 binding domain. We also demonstrated that Rho GTPases Cdc42 and RhoD associate with and regulate the subcellular localization of PAK5. Taken together, this work suggests that the mitochondria-targeting of PAK5 may link Ras and Rho GTPase-mediated signaling pathways, and sheds light on aspects of PAK5 signaling that may be important for regulating neuronal homeostasis. ^

The role of AKT-mediated phosphorylation in suppression of MAD1 and the development of new approach in protein complex detection

MAX dimerization protein 1 (MAD1) is a basic-helix-loop-helix transcription factors that recruits transcription repressor such as HDAC to suppress target genes transcription. It antagonizes to MYC because the promoter binding sites for MYC are usually also serve as the binding sites for MAD1 so they compete for it. However, the mechanism of the switch between MYC and MAD1 in turning on and off of genes' transcription is obscure. In this study, we demonstrated that AKT-mediated MAD1 phosphorylation inhibits MAD1 transcription repression function. The association between MAD1 and its target genes' promoter is reduced after been phosphorylated by AKT; therefore, consequently, allows MYC to occupy the binding site and activates transcription. Mutation of such phosphorylation site abrogates the inhibition from AKT. In addition, functional assays demonstrated that AKT suppressed MAD1-mediated transcription repression of its target genes hTERT and ODC. Cell cycle and cell growth were also been released from inhibition by MAD1 in the presents of AKT. Taken together, our study suggests that MAD1 is a novel substrate of AKT and AKT-mediated MAD1 phosphorylation inhibits MAD1function; therefore, activates MAD1 target genes expression. ^ Furthermore, analysis of protein-protein interaction is indispensable for current molecular biology research, but multiplex protein dynamics in cells is too complicated to be analyzed by using existing biochemical methods. To overcome the disadvantage, we have developed a single molecule level detection system with nanofluidic chip. Single molecule was analyzed based on their fluorescent profile and their profiles were plotted into 2 dimensional time co-incident photon burst diagram (2DTP). From this 2DTP, protein complexes were characterized. These results demonstrate that the nanochannel protein detection system is a promising tool for future molecular biology. ^

Immunologic characteristics of immunogenic variants generated by {\it in vitro\/} treatment of a methylcholanthrene-induced murine fibrosarcoma MCA-F with 1 methyl-3-nitro-1-nitrosoguanidine, 5 -aza-2$\sp\prime$-deoxycytidine, or ultraviolet radiation

This study addresses the questions of whether the frequency of generation and in vivo cross-reactivity of highly immunogenic tumor clones induced in a single parental murine fibrosarcoma cell line MCA-F is more closely related to the agent used to induce the Imm$\sp{+}$ clone or whether these characteristics are independent of the agents used. These questions were addressed by treating the parental tumor cell line MCA-F with UV-B radiation (UV-B), 1-methyl-3-nitro-1-nitrosoguanidine (MNNG), or 5-aza-2$\sp\prime$-deoxycytidine (5-azaCdR). The frequency of Imm$\sp{+}$ variant generation was similarly high for the three different agents, suggesting that the frequency of Imm$\sp{+}$ generation was related more closely to the cell line than to the inducing agent used. Cross-reactivity was tested with two Imm$\sp{+}$ clones from each treatment group in a modified immunoprotection assay that selectively engendered antivariant, but not antiparental immunity. Under these conditions each clone, except one, immunized against itself. The MNNG-induced clones engendered stronger antivariant immunity but a weaker variant cross-reactive immunity could also be detected.^ This study also characterized the lymphocyte populations responsible for antivariant and antiparental immunity in vivo. Using the local adoptive transfer assay (LATA) and antibody plus complement depletion of T-cell subsets, we showed that immunity induced by the Imm$\sp{+}$ variants against the parent MCA-F was transferred by the Thy1.2$\sp{+}$, L3T4a$\sp{+}$, Lyt2.1$\sp{-}$ (CD4$\sp{+}$) population, without an apparent contribution by Thy1.2$\sp{+}$, L3T4a$\sp{-}$, Lyt2.1$\sp{+}$ (CD8$\sp{+}$) cells. A role for Lyt2.1$\sp{+}$T lymphocytes in antivariant, but not antiparent immunity was supported by the results of LATA and CTL assays. Immunization with low numbers of viable Imm$\sp{+}$ cells, or with high numbers of non viable Imm$\sp{+}$ cells engendered only antivariant immunity without parental cross-protection. The associative recognition of parental antigens and variant neoantigens resulting in strong antiparent immunity was investigated using somatic cells hybrids of Imm$\sp{+}$ variants of MCA-F and an antigenically distinct tumor MCA-D. An unexpected result of these latter experiments was the expression of a unique tumor-specific antigen by the hybrid cells. These studies demonstrate that the parental tumor-specific antigen and the variant neoantigen must be coexpressed on the cell surface to engender parental cross-protective immunity. (Abstract shortened with permission of author.) ^

Transcriptional repression of serum amyloid A1 by AP -2 contributes to its liver -specific expression

Studies on the transcriptional regulation of serum amyloid A1 (SAA1) gene, a liver specific acute-phase gene, identified a regulatory element in its promoter that functioned to repress (SAA1) gene transcription in nonliver cells. This silencer element interacts with a nuclear protein that is detectable in HeLa cells, fibroblasts and placental tissues but not in liver or liver-derived cells. As the expression pattern of this repressor is consistent with its potential regulatory role in repressing SAA1 expression, and that many other liver gene promoters also contain this repressor binding site, we sought to investigate whether this repressor may have a broader functional role in repressing liver genes. ^ We have utilized protein purification, cell culture, transient and stable gene transfection, and molecular biology approaches to identify this protein and investigate its possible function in the regulation of (SAA1) and other liver genes. Analyses of amino acid sequence of the purified nuclear protein, and western blot and gel shift studies identified the repressor as transcription factor AP-2 or AP-2-like protein. Using transient transfection of DNA into cultured cells, we demonstrate that AP-2 can indeed function as a repressor to inhibit transcription of SAA1 gene promoter. This conclusion is supported by the following experimental results: (1) overexpression of AP-2 in hepatoma cells inhibits conditioned medium (CM)-induced expression of SAA1 promoter; (2) binding of AP-2 to the SAA1 promoter is required for AP-2 repression function; (3) one mechanism by which AP-2 inhibits SAA1 may be by antagonizing the activation function of the strong transactivator NFκB; (4) mutation of AP-2 binding sites results in derepression of SAM promoter in HeLa cells; and (5) inhibition of endogenous AP-2 activity by a dominant-negative mutant abolishes AP-2's inhibitory effect on SAM promoter in HeLa cells. In addition to the SAM promoter, AP-2 also can bind to the promoter regions of six other liver genes tested, suggesting that it may have a broad functional role in restricting the expression of many liver genes in nonliver cells. Consistent with this notion, ectopic expression of AP-2 also represses CM-mediated activation of human third component of complement 3 promoter. Finally, in AP-2-expressing stable hepatoma cell lines, AP-2 inhibits not only the expression of endogenous SAA, but also the expression of several other endogenous liver genes including albumin, α-fetoprotein. ^ Our findings that AP-2 has the ability to repress the expression of liver genes in nonliver cells opens a new avenue of investigation of negative regulation of gene transcription, and should improve our understanding of tissue-specific expression of liver genes. In summary, our data provide evidence suggesting a novel role of AP-2 as a repressor, inhibiting the expression of liver genes in nonliver cells. Thus, the tissue-specific expression of AP-2 may constitute an important mechanism contributing to the liver-specific expression of liver genes. ^

Cancer therapy combining the modalities of hyperthermia and chemotherapy: In vitro cellular response after rapid heat accumulation in the cancer cell

Hyperthermia is usually used at a sub-lethal level in cancer treatment to potentiate the effects of chemotherapy. The purpose of this study is to investigate the role of heating rate in achieving synergistic cell killing by chemotherapy and hyperthermia. For this purpose, in vitro cell culture experiments with a uterine cancer cell line (MES-SA) and its multidrug resistant (MDR) variant MES-SA/Dx5 were conducted. The cytotoxicity, mode of cell death, induction of thermal tolerance and P-gp mediated MDR following the two different modes of heating were studied. Doxorubicin (DOX) was used as the chemotherapy drug. Indocyanine green (ICG), which absorbs near infrared light at 808nm (ideal for tissue penetration), was chosen for achieving rapid rate hyperthermia. A slow rate hyperthermia was provided by a cell culture incubator. The results show that the potentiating effect of hyperthermia to chemotherapy can be maximized by increasing the rate of heating as evident by the results from the cytotoxicity assay. When delivered at the same thermal dose, a rapid increase in temperature from 37°C to 43°C caused more cell membrane damage than gradually heating the cells from 37°C to 43°C and thus allowed for more intracellular accumulation of the chemotherapeutic agents. Different modes of cell death are observed by the two hyperthermia delivery methods. The rapid rate laser-ICG hyperthermia @ 43°C caused cell necrosis whereas the slow rate incubator hyperthermia @ 43°C induced very mild apoptosis. At 43°C a positive correlation between thermal tolerance and the length of hyperthermia exposure is identified. This study shows that by increasing the rate of heating, less thermal dose is needed in order to overcome P-gp mediated MDR.

Multifunctional nanoparticles in cancer: In vitro characterization, in vivo distribution, and cellular response after laser-NIR dye-induced heating

A novel biocompatible and biodegradable polymer, termed poly(Glycerol malate co-dodecanedioate) (PGMD), was prepared by thermal condensation method and used for fabrication of nanoparticles (NPs). PGMD NPs were prepared using the single oil emulsion technique and loaded with an imaging/hyperthermia agent (IR820) and a chemotherapeutic agent (doxorubicin, DOX). The size of the void PGMD NPs, IR820-PGMD NPs and DOX-IR820-PGMD NPs were approximately 90 nm, 110 nm, and 125 nm respectively. An acidic environment (pH=5.0) induced higher DOX and IR820 release compared to pH=7.4. DOX release was also enhanced by exposure to laser, which increased the temperature to 42°C. Cytotoxicity of DOX-IR820-PGMD NPs was comparable in MES-SA but was higher in Dx5 cells compared to free DOX plus IR820 (p<0.05). The combination of hyperthermia (HT) and chemotherapy improved cytotoxicity in both cell lines. We also explored the cellular response after rapid, short-term and low thermal dose (laser/Dye/NP) induced-heating, and compared it to slow, long-term and high thermal dose cell incubator heating by investigating the reactive oxygen species (ROS) level, hypoxia-inducible factor-1&agr; (HIF-1&agr;) and vascular endothelial growth factor (VEGF) expression. The cytotoxicity of IR820-PGMD NPs after laser/Dye/NP HT resulted in higher cancer cell killing compared to incubator HT. ROS level, HIF-1&agr; and VEGF expression were elevated under incubator HT, while maintained at the baseline level under the laser/Dye/NP HT. In vivo mouse studies showed that NP formulation significantly improved the plasma half-life of IR820 after tail vein injection. Significant lower IR820 content was observed in kidney in DOX-IR820-PGMD NP treatment as compared to free IR820 treatment in our biodistribution studies (p<0.05). In conclusion, both IR820-PGMD NPs and DOX-IR820-PGMD NPs were successfully developed and used for both imaging and therapeutic purposes. Rapid and short-term laser/Dye/NP HT, with a low thermal dose, did not up-regulate HIF-1&agr; and VEGF expression, whereas slow and long-term incubator HT, with a high thermal dose, can enhance expression of both HIF-1&agr; and VEGF.^

Pore size regulates cell and tissue interactions with PLGA-CaP scaffolds used for bone engineering

A common subject in bone tissue engineering is the need for porous scaffolds to support cell and tissue interactions aiming at repairing bone tissue. As poly(lactide-co-glycolide)calcium phosphate (PLGACaP) scaffolds can be manufactured with different pore sizes, the aim of this study was to evaluate the effect of pore diameter on osteoblastic cell responses and bone tissue formation. Scaffolds were prepared with 85% porosity, with pore diameters in the ranges 470590, 590850 and 8501200 mu m. Rat bone marrow stem cells differentiated into osteoblasts were cultured on the scaffolds for up to 10 days to evaluate cell growth, alkaline phosphatase (ALP) activity and the gene expression of the osteoblast markers RUNX2, OSX, COL, MSX2, ALP, OC and BSP by real-time PCR. Scaffolds were implanted in critical size rat calvarial defects for 2, 4, and 8 weeks for histomorphometric analysis. Cell growth and ALP activity were not affected by the pore size; however, there was an increase in the gene expression of osteoblastic markers with the increase in the pore sizes. At 2 weeks all scaffolds displayed a similar amount of bone and blood vessels formation. At 4 and 8 weeks much more bone formation and an increased number of blood vessels were observed in scaffolds with pores of 470590 mu m. These results show that PLGACaP is a promising biomaterial for bone engineering. However, ideally, combinations of larger (similar to 1000 mu m) and smaller (similar to 500 mu m) pores in a single scaffold would optimize cellular and tissue responses during bone healing. Copyright (C) 2011 John Wiley & Sons, Ltd.