Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Computational selection of inhibitors of Abeta aggregation and neuronal toxicity.

Alzheimer's disease (AD) is characterized by the cerebral accumulation of misfolded and aggregated amyloid-beta protein (Abeta). Disease symptoms can be alleviated, in vitro and in vivo, by 'beta-sheet breaker' pentapeptides that reduce plaque load. However the peptide nature of these compounds, made them biologically unstable and unable to penetrate membranes with high efficiency. The main goal of this study was to use computational methods to identify small molecule mimetics with better drug-like properties. For this purpose, the docked conformations of the active peptides were used to identify compounds with similar activities. A series of related beta-sheet breaker peptides were docked to solid state NMR structures of a fibrillar form of Abeta. The lowest energy conformations of the active peptides were used to design three dimensional (3D)-pharmacophores, suitable for screening the NCI database with Unity. Small molecular weight compounds with physicochemical features and a conformation similar to the active peptides were selected, ranked by docking and biochemical parameters. Of 16 diverse compounds selected for experimental screening, 2 prevented and reversed Abeta aggregation at 2-3microM concentration, as measured by Thioflavin T (ThT) fluorescence and ELISA assays. They also prevented the toxic effects of aggregated Abeta on neuroblastoma cells. Their low molecular weight and aqueous solubility makes them promising lead compounds for treating AD.

Expression and Regulation of Human Cytochrome P450 4F Isoforms in Tissue Samples and Under TNF-Alpha Challenges

CYP4F (Cytochrome P4504F) enzymes metabolize endogenous molecules including leukotrienes, prostaglandins and arachidonic acid. The involvement of these endogenous compounds in inflammation has led to the hypothesis that changes in the inflamed tissue environment may affect the expression of CYP4Fs during the pro-inflammatory state, which in turn may modulate inflammatory conditions during the anti-inflammatory state. We demonstrated that inflamed tissues have different levels of CYP4F isoform expression profiles in a number of human samples when compared to the average population. The CYP4F isoform expression levels change with the degree of inflammation present in tissue. Further investigation in cell culture studies revealed that inflammatory cytokines, in particular TNF-α, play a role in regulating the expression of the CYP4F family. One of the isoforms, CYP4F11, had different characteristics than that of the other five CYP4F family members. CYP4F11 metabolizes xenobiotics while the other isoforms metabolize endogenous compounds with higher affinity. CYP4F11 also was expressed at high quantities in the brain, and was up-regulated by TNF-α, while the other isoforms were not expressed at high quantities in the brain and were down-regulated by TNF-α. We identified the AP-1 protein of the JNK pathway as the signaling protein that causes significant increase in CYP4F11 expression. Since TNF-α stimulation causes a simultaneous activation of both JNK pathway and NF-κB signaling, we investigated further the role that NF-κB plays on expression of the CYP4F11 gene. We concluded that although there is a significant increase in CYP4F11 expression in the presence of TNF-α, the activation of NF-κB signaling inhibits CYP4F11 expression in a time dependent manner. The expression of CYP4F11 is only significantly increased after 24 hours of treatment with TNF-α; at shorter time points NF-κB signaling overpowers the JNK pathway activation. We believe that these findings may in the future lead to improved drug design for modulating inflammation.

THE HISTONE DEACETYLASE INHIBITOR, MS-275, SENSITIZES METASTATIC OSTEOSARCOMA TO FASL-INDUCED CELL DEATH: A ROLE FOR C-FLIP

The purpose of this study was to determine the effects of the histone deacetylase inhibitor, MS-275, on the Fas signaling pathway and susceptibility of osteosarcoma (OS) to Fas ligand (FasL)-induced cell death. OS metastasizes almost exclusively to the lungs. We have shown that Fas expression in OS cells is inversely correlated with their metastatic potential. Fas+ cells are rapidly eliminated when they enter the lungs via interaction with FasL, which is constitutively expressed in the lungs. Fas- OS cells escape this FasL-induced apoptosis and survive in the lung microenvironment. Moreover, upregulation of Fas in established OS lung metastases results in tumor regression. Therefore, agents that upregulate Fas expression or activate the Fas signaling pathway may have therapeutic potential. Treatment of Fas- metastatic OS cell lines with 2 μM MS-275 sensitized cells to FasL-induced cell death in vitro. We found that MS-275 did not alter the expression of Fas on the cell surface; rather it resulted in increased levels of Fas within the membrane lipid rafts, as demonstrated by an increase in Fas expression in detergent insoluble lipid raft fractions. We further demonstrated that following MS-275 treatment, Fas colocalized with GM1+ lipid rafts and that there was a decrease in c-FLIP (cellular FLICE-inhibitory protein) mRNA and protein. Downregulation of c-FLIP correlated with caspase activation and apoptosis induction. Transfection of cells with shRNA to c-FLIP also resulted in the localization of Fas to lipid rafts. These studies indicate that MS-275 sensitizes OS cells to FasL by upregulating the expression of Fas in membrane lipid rafts, which correlated with the downregulation of c-FLIP. Treatment of nu/nu-mice with established OS lung metastases with oral MS-275 resulted in increased apoptosis, a significant inhibition of c-FLIP expression in tumors and tumor regression. Histopathological examination of mice showed no significant organ toxicity. Overall, these results suggest that the mechanism by which MS-275 sensitizes OS cells and lung metastases to FasL-induced cell death may be by a reduction in the expression of c-FLIP.

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.

CHARACTERIZATION OF SELECTIVE INHIBITION OF ADENYLYL CYCLASE ACTIVITY BY SMALL MOLECULE INHIBITOR NKY80

The nine membrane-bound isoforms of adenylyl cyclase (AC), via synthesis of the signaling molecule cyclic AMP (cAMP), are involved in many isoform specific physiological functions. Decreasing AC5 activity has been shown to have potential therapeutic benefit, including reduced stress on the heart, pain relief, and attenuation of morphine dependence and withdrawal behaviors. However, AC structure is well conserved, and there are currently no isoform selective AC inhibitors in clinical use. P-site inhibitors inhibit AC directly at the catalytic site, but with an uncompetitive or noncompetitive mechanism. Due to this mechanism and nanomolar potency in cell-free systems, attempts at ligand-based drug design of novel AC inhibitors frequently use P-site inhibitors as a starting template. One small molecule inhibitor designed through this process, NKY80, is described as an AC5 selective inhibitor with low micromolar potency in vitro. P-site inhibitors reveal important ligand binding “pockets” in the AC catalytic site, but specific interactions that give NKY80 selectivity are unclear. Identifying and characterizing unique interactions between NKY80 and AC isoforms would significantly aid the development of isoform selective AC inhibitors. I hypothesized that NKY80’s selective inhibition is conferred by AC isoform specific interactions with the compound within the catalytic site. A structure-based virtual screen of the AC catalytic site was used to identify novel small molecule AC inhibitors. Identified novel inhibitors are isoform selective, supporting the catalytic site as a region capable of more potent isoform selective inhibition. Although NKY80 is touted commercially as an AC5 selective inhibitor, its characterization suggests strong inhibition of both AC5 and the closely related AC6. NKY80 was also virtually docked to AC to determine how NKY80 binds to the catalytic site. My results show a difference between NKY80 binding and the conformation of classic P-site inhibitors. The selectivity and notable differences in NKY80 binding to the AC catalytic site suggest a catalytic subregion more flexible in AC5 and AC6 that can be targeted by selective small molecule inhibitors.

Development of HIF-1α/HIF-1β heterodimerization inhibitors using a novel bioluminescence reporter assay system for in vitro high throughput screening and in vivo imaging

Tumor growth often outpaces its vascularization, leading to development of a hypoxic tumor microenvironment. In response, an intracellular hypoxia survival pathway is initiated by heterodimerization of hypoxia-inducible factor (HIF)-1α and HIF-1β, which subsequently upregulates the expression of several hypoxia-inducible genes, promotes cell survival and stimulates angiogenesis in the oxygen-deprived environment. Hypoxic tumor regions are often associated with resistance to various classes of radio- or chemotherapeutic agents. Therefore, development of HIF-1α/β heterodimerization inhibitors may provide a novel approach to anti-cancer therapy. To this end, a novel approach for imaging HIF-1α/β heterodimerization in vitro and in vivo was developed in this study. Using this screening platform, we identified a promising lead candidate and further chemically derivatized the lead candidate to assess the structure-activity relationship (SAR). The most effective first generation drug inhibitors were selected and their pharmacodynamics and anti-tumor efficacy in vivo were verified by bioluminescence imaging (BLI) of HIF-1α/β heterodimerization in the xenograft tumor model. Furthermore, the first generation drug inhibitors, M-TMCP and D-TMCP, demonstrated efficacy as monotherapies, resulting in tumor growth inhibition via disruption of HIF-1 signaling-mediated tumor stromal neoangiogenesis.

Design, Synthesis and Development of Transporter Targeting Agents for Image-guided Therapy and Drug Delivery

The purpose of this study was to design, synthesize and develop novel transporter targeting agents for image-guided therapy and drug delivery. Two novel agents, N4-guanine (N4amG) and glycopeptide (GP) were synthesized for tumor cell proliferation assessment and cancer theranostic platform, respectively. N4amG and GP were synthesized and radiolabeled with 99mTc and 68Ga. The chemical and radiochemical purities as well as radiochemical stabilities of radiolabeled N4amG and GP were tested. In vitro stability assessment showed both 99mTc-N4amG and 99mTc-GP were stable up to 6 hours, whereas 68Ga-GP was stable up to 2 hours. Cell culture studies confirmed radiolabeled N4amG and GP could penetrate the cell membrane through nucleoside transporters and amino acid transporters, respectively. Up to 40% of intracellular 99mTc-N4amG and 99mTc-GP was found within cell nucleus following 2 hours of incubation. Flow cytometry analysis revealed 99mTc-N4amG was a cell cycle S phase-specific agent. There was a significant difference of the uptake of 99mTc-GP between pre- and post- paclitaxel-treated cells, which suggests that 99mTc-GP may be useful in chemotherapy treatment monitoring. Moreover, radiolabeled N4amG and GP were tested in vivo using tumor-bearing animal models. 99mTc-N4amG showed an increase in tumor-to-muscle count density ratios up to 5 at 4 hour imaging. Both 99mTc-labeled agents showed decreased tumor uptake after paclitaxel treatment. Immunohistochemistry analysis demonstrated that the uptake of 99mTc-N4amG was correlated with Ki-67 expression. Both 99mTc-N4amG and 99mTc-GP could differentiate between tumor and inflammation in animal studies. Furthermore, 68Ga-GP was compared to 18F-FDG in rabbit PET imaging studies. 68Ga-GP had lower tumor standardized uptake values (SUV), but similar uptake dynamics, and different biodistribution compared with 18F-FDG. Finally, to demonstrate that GP can be a potential drug carrier for cancer theranostics, several drugs, including doxorubicin, were selected to be conjugated to GP. Imaging studies demonstrated that tumor uptake of GP-drug conjugates was increased as a function of time. GP-doxorubicin (GP-DOX) showed a slow-release pattern in in vitro cytotoxicity assay and exhibited anti-cancer efficacy with reduced toxicity in in vivo tumor growth delay study. In conclusion, both N4amG and GP are transporter-based targeting agents. Radiolabeled N4amG can be used for tumor cell proliferation assessment. GP is a potential agent for image-guided therapy and drug delivery.

Design, synthesis and cellular and molecular pharmacology of 2$\sp\prime$-bromo-4$\sp\prime$-epidaunorubicin (WP401): A novel non-cross-resistant anthracycline antibiotic with the intrinsic ability to circumvent P -glycoprotein-mediated drug resistance

We designed and synthesized a novel daunorubicin (DNR) analogue that effectively circumvents P-glycoprotein (P-gp)-mediated drug resistance. The fully protected carbohydrate intermediate 1,2-dibromoacosamine was prepared from acosamine and effectively coupled to daunomycinone in high yield. Deprotection under alkaline conditions yielded 2$\sp\prime$-bromo-4$\sp\prime$-epidaunorubicin (WP401). The in vitro cytotoxicity and cellular and molecular pharmacology of WP401 were compared with those of DNR in a panel of wild-type cell lines (KB-3-1, P388S, and HL60S) and their multidrug-resistant (MDR) counterparts (KB-V1, P388/DOX, and HL60/DOX). Fluorescent spectrophotometry, flow cytometry, and confocal laser scanning microscopy were used to measure intracellular accumulation, retention, and subcellular distribution of these agents. All MDR cell lines exhibited reduced DNR uptake that was restored, upon incubation with either verapamil (VER) or cyclosporin A (CSA), to the level found in sensitive cell lines. In contrast, the uptake of WP401 was essentially the same in the absence or presence of VER or CSA in all tested cell lines. The in vitro cytotoxicity of WP401 was similar to that of DNR in the sensitive cell lines but significantly higher in resistant cell lines (resistance index (RI) of 2-6 for WP401 vs 75-85 for DNR). To ascertain whether drug-mediated cytotoxicity and retention were accompanied by DNA strand breaks, DNA single- and double-strand breaks were assessed by alkaline elution. High levels of such breaks were obtained using 0.1-2 $\mu$g/mL of WP401 in both sensitive and resistant cells. In contrast, DNR caused strand breaks only in sensitive cells and not much in resistant cells. We also compared drug-induced DNA fragmentation similar to that induced by DNR. However, in P-gp-positive cells, WP401 induced 2- to 5-fold more DNA fragmentation than DNR. This increased DNA strand breakage by WP401 was correlated with its increased uptake and cytotoxicity in these cell lines. Overall these results indicate that WP401 is more cytotoxic than DNR in MDR cells and that this phenomenon might be related to the reduced basicity of the amino group and increased lipophilicity of WP401. ^

An Innovative Phase I Trial Design Allowing for the Identification of Multiple Potential Maximum Tolerated Doses with Combination Therapy of Targeted Agents

Treatment for cancer often involves combination therapies used both in medical practice and clinical trials. Korn and Simon listed three reasons for the utility of combinations: 1) biochemical synergism, 2) differential susceptibility of tumor cells to different agents, and 3) higher achievable dose intensity by exploiting non-overlapping toxicities to the host. Even if the toxicity profile of each agent of a given combination is known, the toxicity profile of the agents used in combination must be established. Thus, caution is required when designing and evaluating trials with combination therapies. Traditional clinical design is based on the consideration of a single drug. However, a trial of drugs in combination requires a dose-selection procedure that is vastly different than that needed for a single-drug trial. When two drugs are combined in a phase I trial, an important trial objective is to determine the maximum tolerated dose (MTD). The MTD is defined as the dose level below the dose at which two of six patients experience drug-related dose-limiting toxicity (DLT). In phase I trials that combine two agents, more than one MTD generally exists, although all are rarely determined. For example, there may be an MTD that includes high doses of drug A with lower doses of drug B, another one for high doses of drug B with lower doses of drug A, and yet another for intermediate doses of both drugs administered together. With classic phase I trial designs, only one MTD is identified. Our new trial design allows identification of more than one MTD efficiently, within the context of a single protocol. The two drugs combined in our phase I trial are temsirolimus and bevacizumab. Bevacizumab is a monoclonal antibody targeting the vascular endothelial growth factor (VEGF) pathway which is fundamental for tumor growth and metastasis. One mechanism of tumor resistance to antiangiogenic therapy is upregulation of hypoxia inducible factor 1α (HIF-1α) which mediates responses to hypoxic conditions. Temsirolimus has resulted in reduced levels of HIF-1α making this an ideal combination therapy. Dr. Donald Berry developed a trial design schema for evaluating low, intermediate and high dose levels of two drugs given in combination as illustrated in a recently published paper in Biometrics entitled “A Parallel Phase I/II Clinical Trial Design for Combination Therapies.” His trial design utilized cytotoxic chemotherapy. We adapted this design schema by incorporating greater numbers of dose levels for each drug. Additional dose levels are being examined because it has been the experience of phase I trials that targeted agents, when given in combination, are often effective at dosing levels lower than the FDA-approved dose of said drugs. A total of thirteen dose levels including representative high, intermediate and low dose levels of temsirolimus with representative high, intermediate, and low dose levels of bevacizumab will be evaluated. We hypothesize that our new trial design will facilitate identification of more than one MTD, if they exist, efficiently and within the context of a single protocol. Doses gleaned from this approach could potentially allow for a more personalized approach in dose selection from among the MTDs obtained that can be based upon a patient’s specific co-morbid conditions or anticipated toxicities.

Design, synthesis and biological evaluation of 5$\sp\prime$-mononucleotide prodrugs: Strategies to introduce nucleotides into cells

The neutral bis ((pivaloyloxy)methyl) (PIV$\sb2\rbrack$ derivatives of FdUMP, ddUMP, and AZTMP were synthesized as potential membrane-permeable prodrugs of FdUMP, ddUMP, and AZTMP. These compounds were designed to enter cells by passive diffusion and revert to the parent nucleotides after removal of the PIV groups by hydrolytic enzymes. These prodrugs were prepared by condensation of FUdR, ddU, and AZT with PIV$\sb2$ phosphate in the presence of triphenylphosphine and diethyl azodicarboxylate (the Mitsunobo reagent). PIV$\sb2$-FdUMP, PIV$\sb2$-ddUMP, and PIV$\sb2$-AZTMP were stable in the pH range 1.0-4.0 (t$\sb{1/2} = {>}$100 h). They were also fairly stable at pH 7.4 (t$\sb{1/2} = {>}$40 h). In 0.05 M NaOH solution, however, they were rapidly degraded (t$\sb{1/2} < 2$ min). In the presence hog liver carboxylate esterase, they were converted quantitatively to the corresponding phosphodiesters, PIV$\sb1$-FdUMP, PIV$\sb1$-ddUMP, and PIV$\sb1$-AZTMP; after 24 h incubation, only trace amounts of FdUMP, ddUMP, and AZTMP (1-5%) were observed indicating that the PIV$\sb1$ compounds were poor substrates for the enzyme. In human plasma, the PIV$\sb2$ compounds were rapidly degraded with half-lives of less than 5 min. The rate of degradation of the PIV$\sb2$ compounds in the presence of phosphodiesterase I was the same as that in buffer controls, indicating that they were not substrates for this enzyme. In the presence of phosphodiesterase I, PIV$\sb1$-FdUMP, PIV$\sb1$-ddUMP, and PIV$\sb1$-AZTMP were converted quantitatively to FdUMP, ddUMP, and AZTMP.^ PIV$\sb2$-ddUMP and PIV$\sb2$-AZTMP were effective at controlling HIV type 1 infection in MT-4 and CEM tk$\sp-$ cells in culture. Mechanistic studies demonstrated that PIV$\sb2$-ddUMP and PIV$\sb2$-AZTMP were taken up by the cells and converted to ddUTP and AZTTP, both potent inhibitors of HIV reverse transcriptase. However, a potential shortcoming of PIV$\sb2$-ddUMP and PIV$\sb2$-AZTMP as clinical therapeutic agents is that they are rapidly degraded (t$\sb{1/2}$ = approx. 4 minutes) in human plasma by carboxylate esterases. To circumvent this limitation, chemically-labile nucleotide prodrugs and liposome-encapsulated nucleotide prodrugs were investigated. In the former approach, the protective groups bis(N, N-(dimethyl)carbamoyloxymethyl) (DM$\sb2$) and bis (N-(piperidino)carbamoyloxymethyl) (DP$\sb2$) were used to synthesize DM$\sb2$-ddUMP and DP$\sb2$-ddUMP, respectively. In aqueous buffers (pH range 1.0-9.0) these compounds were degraded with half-lives of 3 to 4 h. They had similar half-lives in human plasma demonstrating that they were resistant to esterase-mediated cleavage. However, neither compound gave rise to significant concentrations of ddUMP in CEM or CEM tk$\sp-$ cells. In the liposome-encapsulated nucleotide prodrug approach, three different liposomal formulations of PIV$\sb2$-ddUMP (L-PIV$\sb2$-ddUMP) were investigated. The half-lifes of these L-PIV$\sb2$-ddUMP preparations in human plasma were 2 h compared with 4 min for the free drug. The preparations were more effective at controlling HIV-1 infection than free PIV$\sb2$-ddUMP in human T cells in culture. Collectively, these data indicate that PIV$\sb2$-FdUMP, PIV$\sb2$-ddUMP, and PIV$\sb2$-AZTMP are effective membrane-permeable prodrugs of FdUMP, ddUMP, and AZTMP. ^

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. ^

Club drug use and unprotected anal intercourse among HIV-seronegative men who have sex with men seeking HIV/STD testing at a local STD clinic in Houston

Objective. To explore (1) the association between "club drug" use and unprotected anal intercourse (UAI) and (2) the association between binge drug use and UAI among HIV seronegative men who have sex with men (MSM) seeking HIV/STD testing at a local clinic in Houston. ^ Study design. A sub-sample of 297 HIV seronegative MSM from a cross-sectional study of drug and sexual behavior in Houston was conducted in 2006. Patients who were seeking HIV/STD testing at a local MSM-identified STD clinic were recruited for an anonymous computer-assisted interview. Analysis of identified secondary data consisted of self-reported information about demographic characteristics, use of drugs, and sexual behaviors. ^ Results. With new and casual sex partners, there was a strong and statistically significant association between use of "club drugs" and UAI. No association between binge drug use and UAI was evident. Men aware of HIV seropositivity or unaware of the HIV serostatus of their primary partner were less likely to report UAI. ^ Conclusion. These data suggest that in the Houston area, HIV-negative MSM club drug users, particularly multiple drug users, are at higher risk of UAI than comparable MSMs who do not use club drugs. Episode-level data regarding binge use of these and other drugs, and UAI should be collected in future studies to explore their relationship. The 'new partner' category should be added to sex partner types to measure sex and drug use behaviors in future studies.^ Keywords. HIV-negative MSM; club drugs; unprotected anal intercourse; binge drug use. ^

Relevant drug to drug interactions between itraconazole and non-nucleoside reverse transcriptase inhibitors in HIV-infected patients on maintenance therapy for disseminated histoplasmosis

Objective. Itraconazole is recommended life-long for preventing relapse of disseminated histoplasmosis in HIV-infected patients. I sought to determine if serum itraconazole levels are affected by the type of Highly Active Anti-Retroviral Therapy (NNRTI or PI) being taken concomitantly to treat HIV. ^ Design. Retrospective cohort. ^ Methods. De-identified data were used from an IRB-approved parent study which identified patients on HAART and maintenance itraconazole for confirmed disseminated histoplasmosis between January 2003 and December 2006. Available itraconazole blood levels were abstracted as well as medications taken by each patient at the time of the blood tests. Mean itraconazole levels were compared using the student's t-test. ^ Results. 11 patients met study criteria. Patient characteristics were: median age 36, 91% men, 18% white, 18% black, 55% Hispanic and 9% Asians, median CD4 cell count 120 cells/mm3. 14 blood levels were available for analysis—8 on PI, 4 on NNRTI and 2 on both. 8/8 itraconazole levels obtained while taking concomitant PI were therapeutic (>0.4 μg/mL) in contrast to 0/4 obtained while taking NNRTI. Two patients switched from NNRTI to PI and reached therapeutic levels. Mean levels on NNRTI (0.05 μg/mL, s.d. 0.0) and on PI (2.45 μg/mL, s.d. 0.21) for these two patients were compared via a paired t-test (t = 16.00, d.f. = 1, P = 0.04). Remaining patient levels were compared using an unpaired t-test. Mean itraconazole on concomitant PI (n = 6) was 1.37 μg/mL (s.d. 0.74), while the mean on concomitant NNRTI was 0.05 μg/mL (s.d. 0.0), t = 2.39, d.f. = 6, P = 0.05. ^ Conclusions. Co-administration of NNRTI and itraconazole results in significant decreases in itraconazole blood levels, likely by inducing the CYP3A4 enzyme system. Itraconazole drug levels should be monitored in patients on concomitant NNRTI. PI-based HAART may be preferred over NNRTI-based HAART when using itraconazole to treat HIV-infected patients with disseminated histoplasmosis. ^