Preclinical modeling of multi-drug cancer therapies

Anticancer drugs typically are administered in the clinic in the form of mixtures, sometimes called combinations. Only in rare cases, however, are mixtures approved as drugs. Rather, research on mixtures tends to occur after single drugs have been approved. The goal of this research project was to develop modeling approaches that would encourage rational preclinical mixture design. To this end, a series of models were developed. First, several QSAR classification models were constructed to predict the cytotoxicity, oral clearance, and acute systemic toxicity of drugs. The QSAR models were applied to a set of over 115,000 natural compounds in order to identify promising ones for testing in mixtures. Second, an improved method was developed to assess synergistic, antagonistic, and additive effects between drugs in a mixture. This method, dubbed the MixLow method, is similar to the Median-Effect method, the de facto standard for assessing drug interactions. The primary difference between the two is that the MixLow method uses a nonlinear mixed-effects model to estimate parameters of concentration-effect curves, rather than an ordinary least squares procedure. Parameter estimators produced by the MixLow method were more precise than those produced by the Median-Effect Method, and coverage of Loewe index confidence intervals was superior. Third, a model was developed to predict drug interactions based on scores obtained from virtual docking experiments. This represents a novel approach for modeling drug mixtures and was more useful for the data modeled here than competing approaches. The model was applied to cytotoxicity data for 45 mixtures, each composed of up to 10 selected drugs. One drug, doxorubicin, was a standard chemotherapy agent and the others were well-known natural compounds including curcumin, EGCG, quercetin, and rhein. Predictions of synergism/antagonism were made for all possible fixed-ratio mixtures, cytotoxicities of the 10 best-scoring mixtures were tested, and drug interactions were assessed. Predicted and observed responses were highly correlated (r2 = 0.83). Results suggested that some mixtures allowed up to an 11-fold reduction of doxorubicin concentrations without sacrificing efficacy. Taken together, the models developed in this project present a general approach to rational design of mixtures during preclinical drug development. ^

Catalytic antibody response to the Staphylococcus aureus virulence factor Efb

Staphylococcus aureus is a globally prevalent pathogen that can cause a wide variety of acute and chronic diseases in both adults and children, in both immune susceptible populations and healthy individuals. Its ability to cause persistent infections has been linked to multiple immune evasion strategies, including Efb-mediated complement inhibition. As new multi-drug-resistant strains emerge, therapeutic alternatives to traditional antibiotics must be developed. These experiments assessed the ability of healthy patient immunoglobulin to cleave Efb and disable the complement-inhibitory properties of Efb in vitro. Levels of immunoglobulin-mediated Efb catalysis varied both between immunoglobulin isoform/isotype and between individuals. Serum IgG showed the strongest catalytic activity of the immunoglobulin isotypes tested. Additionally, IgG hydrolyzed the virulence factor in a way that enabled only minimal binding to the complement component C3b, effectively blocking Efb-mediated inhibition of complement lysis. Salivary IgA and serum IgM did not block Efb-mediated inhibition of complement. Catalytic IgG selectively cleaved Efb and showed no cleavage of a variety of other proteins tested. Catalytic activity of IgG was inhibited by serine protease inhibitors, but not by other protease inhibitors, suggesting a serine-protease mechanism of catalysis. It is proposed that varying concentrations and activity levels of catalytic IgG between healthy individuals and those with current or recurrent S. aureus infections in both adult and pediatric populations be studied in order to assess the potential effectiveness of passive immunization therapy with catalytic monoclonal IgG. ^

The application of latent variable models to the assessment of determinants of HIV risk behavior

Studies on the relationship between psychosocial determinants and HIV risk behaviors have produced little evidence to support hypotheses based on theoretical relationships. One limitation inherent in many articles in the literature is the method of measurement of the determinants and the analytic approach selected. ^ To reduce the misclassification associated with unit scaling of measures specific to internalized homonegativity, I evaluated the psychometric properties of the Reactions to Homosexuality scale in a confirmatory factor analytic framework. In addition, I assessed the measurement invariance of the scale across racial/ethnic classifications in a sample of men who have sex with men. The resulting measure contained eight items loading on three first-order factors. Invariance assessment identified metric and partial strong invariance between racial/ethnic groups in the sample. ^ Application of the updated measure to a structural model allowed for the exploration of direct and indirect effects of internalized homonegativity on unprotected anal intercourse. Pathways identified in the model show that drug and alcohol use at last sexual encounter, the number of sexual partners in the previous three months and sexual compulsivity all contribute directly to risk behavior. Internalized homonegativity reduced the likelihood of exposure to drugs, alcohol or higher numbers of partners. For men who developed compulsive sexual behavior as a coping strategy for internalized homonegativity, there was an increase in the prevalence odds of risk behavior. ^ In the final stage of the analysis, I conducted a latent profile analysis of the items in the updated Reactions to Homosexuality scale. This analysis identified five distinct profiles, which suggested that the construct was not homogeneous in samples of men who have sex with men. Lack of prior consideration of these distinct manifestations of internalized homonegativity may have contributed to the analytic difficulty in identifying a relationship between the trait and high-risk sexual practices. ^

The critical role and regulation of transcription factor FoxM1 in gastric cancer development and progression

The mammalian Forkhead Box (Fox) transcription factor (FoxM1) is implicated in tumorgenesis. However, the role and regulation of FoxM1 in gastric cancer remain unknown.^ I examined FoxM1 expression in 86 cases of primary gastric cancer and 57 normal gastric tissue specimens. I found weak expression of FoxM1 protein in normal gastric mucosa, whereas I observed strong staining for FoxM1 in tumor-cell nuclei in various gastric tumors and lymph node metastases. The aberrant FoxM1 expression is associated with VEGF expression and increased angiogenesis in human gastric cancer. A Cox proportional hazards model revealed that FoxM1 expression was an independent prognostic factor in multivariate analysis. Furthermore, overexpression of FoxM1 by gene transfer significantly promoted the growth and metastasis of gastric cancer cells in orthotopic mouse models, whereas knockdown of FoxM1 expression by small interfering RNA did the opposite. Next, I observed that alteration of tumor growth and metastasis by elevated FoxM1 expression was directly correlated with alteration of VEGF expression and angiogenesis. In addition, promotion of gastric tumorigenesis by FoxM1 directly and significantly correlated with transactivation of vascular endothelial growth factor (VEGF) expression and elevation of angiogenesis. ^ To further investigate the underlying mechanisms that result in FoxM1 overexpression in gastric cancer, I investigated FoxM1 and Krüppel-like factor 4 (KLF4) expressions in primary gastric cancer and normal gastric tissue specimens. Concomitance of increased expression of FoxM1 protein and decreased expression of KLF4 protein was evident in human gastric cancer. Enforced KLF4 expression suppressed FoxM1 protein expression. Moreover, a region within the proximal FoxM1 promoter was identified to have KLF4-binding sites. Finally, I found an increased FoxM1 expression in gastric mucosa of villin-Cre -directed tissue specific Klf4-null mice.^ In summary, I offered both clinical and mechanistic evidence that dysregulated expression of FoxM1 play an important role in gastric cancer development and progression, while KLF4 mediates negative regulation of FoxM1 expression and its loss significantly contributes to FoxM1 dysregulation. ^

Ethnicity as a risk factor for diabetes incidence in a population of diagnosed hypertensives

Significant racial/ethnic differences exist in prevalence of hypertension (HTN) and non-insulin dependent diabetes mellitus (NIDDM). Hypertension is more common in diabetics than in non-diabetics, and an etiologic link between the two conditions has been proposed. Since there are few longitudinal studies of persons with both HTN and NIDDM, a retrospective cohort study was conducted to determine if ethnicity (Black, Hispanic (Mexican-American), and non-Hispanic White) was related to NIDDM incidence in a low-SES, multi-ethnic clinic population of diagnosed hypertensives. Two thousand nine hundred forty-one hypertensives free of NIDDM at baseline were followed for up to 10 years. Mean baseline age was 56 $\pm$ 12 years, M:F percent was 33:67, and Black:Hispanic:White percent was 63:17:20. There were 236 incident cases of NIDDM. In Cox proportional hazards analysis, the risk of developing NIDDM over 10 years was not related to ethnicity after controlling for significant covariates, including age, baseline blood glucose and body mass index (adjusted RR for Blacks compared to Whites =.82, 95 percent CI =.57-1.18; adjusted RR for Hispanics compared to Whites =.84, 95 percent CI =.51-1.38). This result contrasts with the increased risk of NIDDM among Blacks and Hispanics compared to Whites found in the general population. The study suggests that a diagnosis of hypertension equalizes the risk of developing NIDDM among the three ethnic groups. ^

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.