A BAYESIAN APPROACH TO DOSE-RESPONSE ASSESSMENT AND DRUG-DRUG INTERACTION ANALYSIS: APPLICATION TO IN VITRO STUDIES

The considerable search for synergistic agents in cancer research is motivated by the therapeutic benefits achieved by combining anti-cancer agents. Synergistic agents make it possible to reduce dosage while maintaining or enhancing a desired effect. Other favorable outcomes of synergistic agents include reduction in toxicity and minimizing or delaying drug resistance. Dose-response assessment and drug-drug interaction analysis play an important part in the drug discovery process, however analysis are often poorly done. This dissertation is an effort to notably improve dose-response assessment and drug-drug interaction analysis. The most commonly used method in published analysis is the Median-Effect Principle/Combination Index method (Chou and Talalay, 1984). The Median-Effect Principle/Combination Index method leads to inefficiency by ignoring important sources of variation inherent in dose-response data and discarding data points that do not fit the Median-Effect Principle. Previous work has shown that the conventional method yields a high rate of false positives (Boik, Boik, Newman, 2008; Hennessey, Rosner, Bast, Chen, 2010) and, in some cases, low power to detect synergy. There is a great need for improving the current methodology. We developed a Bayesian framework for dose-response modeling and drug-drug interaction analysis. First, we developed a hierarchical meta-regression dose-response model that accounts for various sources of variation and uncertainty and allows one to incorporate knowledge from prior studies into the current analysis, thus offering a more efficient and reliable inference. Second, in the case that parametric dose-response models do not fit the data, we developed a practical and flexible nonparametric regression method for meta-analysis of independently repeated dose-response experiments. Third, and lastly, we developed a method, based on Loewe additivity that allows one to quantitatively assess interaction between two agents combined at a fixed dose ratio. The proposed method makes a comprehensive and honest account of uncertainty within drug interaction assessment. Extensive simulation studies show that the novel methodology improves the screening process of effective/synergistic agents and reduces the incidence of type I error. We consider an ovarian cancer cell line study that investigates the combined effect of DNA methylation inhibitors and histone deacetylation inhibitors in human ovarian cancer cell lines. The hypothesis is that the combination of DNA methylation inhibitors and histone deacetylation inhibitors will enhance antiproliferative activity in human ovarian cancer cell lines compared to treatment with each inhibitor alone. By applying the proposed Bayesian methodology, in vitro synergy was declared for DNA methylation inhibitor, 5-AZA-2'-deoxycytidine combined with one histone deacetylation inhibitor, suberoylanilide hydroxamic acid or trichostatin A in the cell lines HEY and SKOV3. This suggests potential new epigenetic therapies in cell growth inhibition of ovarian cancer cells.

Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

THE USE OF MULTIVARIATE ANALYSIS OF COVARIANCE IN EVALUATING THE COMPARABILITY OF LABORATORY DETERMINATIONS IN COOPERATIVE STUDIES

The role of clinical chemistry has traditionally been to evaluate acutely ill or hospitalized patients. Traditional statistical methods have serious drawbacks in that they use univariate techniques. To demonstrate alternative methodology, a multivariate analysis of covariance model was developed and applied to the data from the Cooperative Study of Sickle Cell Disease.^ The purpose of developing the model for the laboratory data from the CSSCD was to evaluate the comparability of the results from the different clinics. Several variables were incorporated into the model in order to control for possible differences among the clinics that might confound any real laboratory differences.^ Differences for LDH, alkaline phosphatase and SGOT were identified which will necessitate adjustments by clinic whenever these data are used. In addition, aberrant clinic values for LDH, creatinine and BUN were also identified.^ The use of any statistical technique including multivariate analysis without thoughtful consideration may lead to spurious conclusions that may not be corrected for some time, if ever. However, the advantages of multivariate analysis far outweigh its potential problems. If its use increases as it should, the applicability to the analysis of laboratory data in prospective patient monitoring, quality control programs, and interpretation of data from cooperative studies could well have a major impact on the health and well being of a large number of individuals. ^

Analysis of interval -censored events in hypertension studies: Evaluation of treatment of coronary heart disease

Many statistical studies feature data with both exact-time and interval-censored events. While a number of methods currently exist to handle interval-censored events and multivariate exact-time events separately, few techniques exist to deal with their combination. This thesis develops a theoretical framework for analyzing a multivariate endpoint comprised of a single interval-censored event plus an arbitrary number of exact-time events. The approach fuses the exact-time events, modeled using the marginal method of Wei, Lin, and Weissfeld, with a piecewise-exponential interval-censored component. The resulting model incorporates more of the information in the data and also removes some of the biases associated with the exclusion of interval-censored events. A simulation study demonstrates that our approach produces reliable estimates for the model parameters and their variance-covariance matrix. As a real-world data example, we apply this technique to the Systolic Hypertension in the Elderly Program (SHEP) clinical trial, which features three correlated events: clinical non-fatal myocardial infarction, fatal myocardial infarction (two exact-time events), and silent myocardial infarction (one interval-censored event). ^

Direct intrathecal implantation of mesenchymal stromal cells leads to enhanced neuroprotection via an NFkappaB-mediated increase in interleukin-6 production.

Mesenchymal stromal cell (MSC) therapy has shown promise for the treatment of traumatic brain injury (TBI). Although the mechanism(s) by which MSCs offer protection is unclear, initial in vivo work has suggested that modulation of the locoregional inflammatory response could explain the observed benefit. We hypothesize that the direct implantation of MSCs into the injured brain activates resident neuronal stem cell (NSC) niches altering the intracerebral milieu. To test our hypothesis, we conducted initial in vivo studies, followed by a sequence of in vitro studies. In vivo: Sprague-Dawley rats received a controlled cortical impact (CCI) injury with implantation of 1 million MSCs 6 h after injury. Brain tissue supernatant was harvested for analysis of the proinflammatory cytokine profile. In vitro: NSCs were transfected with a firefly luciferase reporter for NFkappaB and placed in contact culture and transwell culture. Additionally, multiplex, quantitative PCR, caspase 3, and EDU assays were completed to evaluate NSC cytokine production, apoptosis, and proliferation, respectively. In vivo: Brain supernatant analysis showed an increase in the proinflammatory cytokines IL-1alpha, IL-1beta, and IL-6. In vitro: NSC NFkappaB activity increased only when in contact culture with MSCs. When in contact with MSCs, NSCs show an increase in IL-6 production as well as a decrease in apoptosis. Direct implantation of MSCs enhances neuroprotection via activation of resident NSC NFkappaB activity (independent of PI3 kinase/AKT pathway) leading to an increase in IL-6 production and decrease in apoptosis. In addition, the observed NFkappaB activity depends on direct cell contact.

The Role of IL-6 in Adenosine-Mediated Pulmonary Fibrosis

Adenosine is a purinergic signaling molecule that regulates various aspects of inflammation and has been implicated in the pathogenesis of chronic lung diseases. Previous studies have demonstrated that adenosine up-regulates IL-6 production through the engagement of the A2B adenosine receptor in various cell types, including alveolar macrophages. IL-6 is elevated in mouse models and humans with chronic lung disease, suggesting a potential role in disease progression. Furthermore, chronic elevation of adenosine in the lungs of adenosine deaminase deficient (Ada-/-) mice leads to the development of pulmonary inflammation, alveolar destruction, and fibrosis, in conjunction with IL-6 elevation. Thus, it was hypothesized that IL-6 contributes to pulmonary inflammation and fibrosis in this model. To test this hypothesis, Ada/IL-6 double knockout mice (Ada/IL-6-/-) were generated to assess the consequences of genetically removing IL-6 on adenosine-dependent pulmonary injury. Ada/IL-6-/- mice exhibited a significant reduction in inflammation, alveolar destruction, and pulmonary fibrosis. Next, Ada-/- mice were treated systematically with IL-6 neutralizing antibodies to test the efficacy of blocking IL-6 on chronic lung disease. These treatments were associated with decreased pulmonary inflammation, alveolar destruction, and fibrosis. To determine the role of IL-6 in a second model of pulmonary fibrosis, wild type mice and IL-6-/- mice were subjected to intraperitoneal injections of bleomycin twice a week for four weeks. Results demonstrated that IL-6-/- mice developed reduced pulmonary fibrosis. To examine a therapeutic approach in this model, wild type mice exposed to bleomycin were treated with IL-6 neutralizing antibodies. Similar results were observed as with Ada-/- mice, namely diminished pulmonary inflammation and fibrosis. In both models, elevations in IL-6 were associated with increased phosphorylated STAT-3 in the nuclei of numerous cell types in the airways, including type II alveolar epithelial cells (AEC). Genetic removal and neutralization of IL-6 in both models was associated with decreased STAT-3 activation in type II AEC. The mechanism of activation in these cells that lack the membrane bound IL-6Ra suggests IL-6 trans-signaling may play a role in regulating fibrosis. Characterization of this mechanism demonstrated that the soluble IL-6Ra (sIL-6Ra) is upregulated in both models during chronic conditions. In vitro studies in MLE-12 alveolar epithelial cells confirmed that IL-6, in combination with the sIL-6Ra, activates STAT-3 and TWIST in association with enhancement of epithelial-to-mesenchymal transition, which can contribute to fibrosis. Similarly, patients with idiopathic pulmonary fibrosis demonstrated a similar pattern of increased IL-6 expression, STAT-3 activation, and sIL-6Ra increases. These findings demonstrate that adenosine-dependent elevations in IL-6 contribute to the development and progression of pulmonary inflammation and fibrosis. The implications from these studies are that adenosine and/or IL-6 neutralizing agents represent novel therapeutic targets for the treatment of pulmonary disorders where fibrosis is a detrimental component.

STUDIES INTO THE MOLECULAR MECHANISMS REGULATING MUSCULAR ATROPHY RESULTING FROM HINDLIMB IMMOBILIZATION IN THE RAT (PROTEIN, SYNTHESIS, RNA)

The purpose of the work performed in this dissertation was to examine some of the possible regulatory mechanisms involved in the initiation of muscular atrophy during periods of decreased muscle utilization resulting from hindlimb immobilization in the rat. A 37% decrease in the rate of total muscle protein synthesis which has been observed to occur in the first 6 h of immobilization contributes significantly to the observed loss of protein during immobilization.^ The rates of cytochrome c and actin synthesis were determined in adult rat red vastus lateralis and gastrocnemius muscles, respectively, by the constant infusion and incorporation of ('3)H-tyrosine into protein. The fractional synthesis rates of both actin and cytochrome c were significantly decreased (P < 0.05) in the 6th h of hindlimb immobilization.^ RHA was extracted from adult rat gastrocnemius muscle by modification of the phenol: chloroform: SDS extraction procedures commonly used for preparation of RNA for hybridization analysis from other mammalian tissues. RNA content of rat gastrocnemius muscle, as determined by this method of extraction and its subsequent quantification by UV absorbance and orcinol assay, was significantly greater than the RNA content previously determined for adult rat gastrocnemius by other commonly employed methods.^ RNA extracted by this method from gastrocnemius muscles of control and 6h immobilized rats was subjected to "dot blot" hybridization to ('32)P-labelled probe from plasmid p749, containing a cDNA sequence complementary to (alpha)-actin mRNA and from rat skeletal muscle. (alpha)-Actin specific mRNA content as estimated by this procedure is not significantly decreased in rat gastrocnemius following 6h or hindlimb immobilization. However, (alpha)-actin specific mRNA content is significantly decreased (P < 0.05) in adult rat gastrocnemius (alpha)-actin specific mRNA is not decreased in adult rat gastrocnemius muscle following 6h of immobilization, a time when actin synthesis is significantly decreased, it is concluded that a change in (alpha)-actin specific mRNA content is not the initiating event responsible for the early decrease in actin synthesis observed in the 6th h of immobilization. ^

Rules to limit snacking in front of the TV and childhood BMI status: Findings from the "Fun Families Study"

Childhood obesity is increasing at epidemic rates, and thus there is a need to target appropriate childhood behaviors that contribute to obesity. Many factors contribute to childhood weight status. The aim of this study was to look at relationships between parental rules to limit snacking while watching television and childhood weight status. The study looked at the presence of the behavior of snacking while watching television yesterday, congruence between child- and parent-reported perception of the presence of rules to limit snacking while watching television, and parent-reported frequency of children following rules to limit snacking while watching television. The outcomes were examined in a multi-ethnic population of children ages 6 to 9 years in Southeast Texas.^ This study was a cross-sectional secondary data analysis of the pilot program, Fun Families. This study examined baseline data from 202 parent-child dyads, which included both the control ( N= 101) and intervention groups (N= 101). Data were gathered using validated questions that were administered to 6-9 year old children and their primary caregiver (referred to as parent in the rest of the discussion) in Southeast Texas, between 2006 and 2008. The main study outcome was childhood weight status based on CDC BMI-for-age categories. The independent variables are (1) the presence of parental rules to limit snacking while watching television, (2) the congruence between child and parent about the presence of rules to limit snacking while watching television, and (3) the parent-reported frequency of the child following the rules to limit snacking while watching television. Chi-Square analyses were used to determine if weight status was different for (1) children who reported rules to limit snacking yesterday, (2) children who reported snacking, (3) children whose parents reported rules were present, and (4) those who had rule congruence with the parents not. Chi-Square analyses also examined if there was a difference in the presence of snacking behavior for children who reported rules, for children whose parents reported rules, and for those children who had congruence about rules. Linear regressions were used to determine if any of the studied variables predicted increased weight status or reported snacking while watching television yesterday.^ This study found that child-reported snacking yesterday was significantly different for children who reported rules (4.12, p= 0.04). Child-reported rules was significantly associated with (p= -0.14, α= 0.04) and predicted child-reported snacking yesterday (R 2 0.021, p= 0.04, t= -2.04, 95% CI -0.31, -0.01). There was statistical significant incongruence between child and parent perception about the presence of rules to limit snacking yesterday (15.06, p= 0.00). For this population, parent education level was significantly associated with child-reported rules (r= -0.16, p= 0.02), child-reported snacking yesterday (r= -0.15, p= 0.04), and parent-reported frequency of child following rules to limit snacking (r= 0.29, p= -0.01). Parent-reported speaking another language besides English at home was significantly associated with parent-reported rules (r= 0.17, p= 0.02).^ Although the studied variables did not show any significant associations or predictors for childhood weight status, the significant discord between parent and child perception about the presence of rules provides valuable information to future interventions that aim to reduce childhood weight status. Including the creation and enforcement of parental rules in interventions to reduce childhood weight status will be beneficial for future studies.^