Maintenance DFMO Increases Survival in High Risk Neuroblastoma

High risk neuroblastoma (HRNB) accounts for 15% of all pediatric cancer deaths. Despite aggressive therapy approximately half of patients will relapse, typically with only transient responses to second-line therapy. This study evaluated the ornithine decarboxylase inhibitor difluoromethylornithine (DFMO) as maintenance therapy to prevent relapse following completion of standard therapy (Stratum 1) or after salvage therapy for relapsed/refractory disease (Stratum 2). This Phase II single agent, single arm multicenter study enrolled from June 2012 to February 2016. Subjects received 2 years of oral DFMO (750 ± 250 mg/m² twice daily). Event free survival (EFS) and overall survival (OS) were determined on an intention-to-treat (ITT) basis. 101 subjects enrolled on Stratum 1 and 100 were eligible for ITT analysis; two-year EFS was 84% (±4%) and OS 97% (±2%). 39 subjects enrolled on Stratum 2, with a two-year EFS of 54% (±8%) and OS 84% (±6%). DFMO was well tolerated. The median survival time is not yet defined for either stratum. DFMO maintenance therapy for HRNB in remission is safe and associated with high EFS and OS. Targeting ODC represents a novel therapeutic mechanism that may provide a new strategy for preventing relapse in children with HRNB.

Abstract 3191: The MDM2 inhibitor CGM097 synergizes with the BET inhibitor OTX015 to induce cell death in neuroblastoma cells

Background: Neuroblastoma (NB) is the most common extracranial solid tumor in infants and children, with amplification of the oncogene MYCN being a hallmark of high-risk disease and poor prognosis. Most NB tumors initially respond to chemotherapy; however, over half of these cases relapse, resulting in chemoresistant disease. In this stage of NB there is evidence of p53 inactivation, which suggests a role for p53 in chemo-sensitivity. Although highly mutated in other cancer types, p53 mutations in NB are rare. Previous data have shown that amplification or deregulation of MDM2, a negative regulator of p53, may be the more common mechanism of p53 suppression and chemoresistance in NB. MYCN and MDM2 have been shown to interact and contribute to NB growth and disease progression. In vitro treatment of NB cells with MDM2 inhibitors has shown promise in increasing the expression of p53, leading to a decrease in proliferation, and increasing apoptosis. BET (Bromodomain and Extra-Terminal domain) inhibitors have also been shown to be effective in treating NB cells in vitro, decreasing MYCN expression, and resulting in increased apoptosis and differentiation. We hypothesize that the combination treatment of a MDM2 inhibitor (CGM097) with a BET inhibitor (OTX015) will result in greater p53 activation and a subsequent synergistic increase in NB cell death.

Methods: Single-agent IC50 values for both CGM097 and OTX015 were obtained for the established NB cell lines SMS-KCNR, SH-SY5Y, BE(2)-C, and CHLA-90 using Calcein AM and Cell Titer-Glo cell viability assays after 48 hours of treatment. Drug combination studies were performed to evaluate synergism between CGM097 and OTX015 in the treatment of NB cells. IncuCyte ZOOM live cell imaging was used for kinetic monitoring of apoptosis with single-agent and combination drug treatment in NB cells. Western blot analysis was used to determine protein levels of apoptosis (Caspase-3) and MYCN/p53 pathway targets (MYCN, MDM2 and p53).

Results: IC50 values for the four established NB cell lines ranged from 0.384 µM-13.5 µM and 7.25 µM-42.79 µM for CGM097 and OTX015, respectively. The combination treatment exhibited synergistic effects in wild-type p53 NB cells, with supporting data obtained from IncuCyte imaging tracking real-time cell death. IncuCyte showed induction of caspase-mediated apoptosis confirmed by Western blot analysis of Cleaved Caspase 3. Pathway analysis shows the importance of wild-type p53 for treatment effect.

Conclusion: This study indicates that the combination of CGM097 and OTX015 synergistically decreases viability in NB cells with wild-type p53 expression. Further in vitro and in vivo work will be necessary to elucidate the mechanisms behind this interaction and in vivo efficacy.

Citation Format: Tyler P. Maser, Joseph W. Zagorski, Austin J. Goodyke, Elizabeth A. VanSickle, Jeffrey P. Bond, Giselle L. Saulnier Sholler. The MDM2 inhibitor CGM097 synergizes with the BET inhibitor OTX015 to induce cell death in neuroblastoma cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3191.

Abstract 68: Neuroblastoma drug response profiles are associated with gene expression profiles

Background: High risk neuroblastoma patients have a poor response to initial chemotherapy, with 20% unable to achieve a complete response. Genomic studies have identified 6 targeted agents whose pathways have been identified in some high risk neuroblastoma tumors. The low frequency of established driver mutations, along with the variety of molecular mechanisms for cancer pathway dysregulation, suggest development of gene expression-based predictive biomarkers in addition to sequence-based biomarkers. Our platform for developing gene expression-based biomarkers includes primary neuroblastoma cell lines, high throughput quantification of cell survival, exome sequencing, and genome-wide expression profiling.

Methods: Cell survival was quantified for each of thirty primary neuroblastoma cell lines in the presence of each of six drugs (bortezomib, crizotinib, dasatinib, lapatinib, sorafenib, and vorinostat) using high throughput cell survival measurements. Titration curves were based on treatment with vehicle as well as sixteen concentrations of each drug (replicated within and between plates). Summary descriptions of the response of each cell line to each drug were obtained using the four parameter logistic model as well as the area under the curve. Exome sequences were obtained using Illumina HiSeq technology. Genome-wide expression profiles were obtained using Affymetrix GeneChip technology. Inference was based on parametric univariate and non-parametric multivariate tests.

Results: For three of the drugs (bortezomib, dasatinib, and vorinostat) the variation in response between cell lines was statistically significant (p < 0.05; F ≥3) and assay results were reproducible in the context of the differences between cell lines (for the other three drugs the estimated standard deviation of log10 EC50 ≤ 0.2). Genome-wide expression profiles were associated with variation in drug response between cell lines, for example, high PDGFRB expression indicates sensitivity to dasatinib (p < 0.001).

Conclusion: High throughput cell survival measurements provided for identification of predictive biomarkers of neuroblastoma drug sensitivity. These biomarkers are currently being studied in the PEDS-PLAN (Pediatric Precision Laboratory Advanced Neuroblastoma Therapy) clinical trial.

Citation Format: Jeffrey P. Bond, Elizabeth VanSickle, Thomas S. Dexhemier, Richard R. Neubig, Giselle L. Sholler. Neuroblastoma drug response profiles are associated with gene expression profiles. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 68.

RAD51 variant proteins from human lung and kidney tumors exhibit DNA strand exchange defects

In human cells, error-free repair of DNA double-strand breaks requires the DNA pairing and strand exchange activities of RAD51 recombinase. Activation of RAD51 recombination activities requires the assembly of RAD51 presynaptic filaments on the single-stranded DNA that forms at resected DSB ends. Mutations in proteins that control presynaptic filament assembly, such as BRCA2, and in RAD51 itself, are associated with human breast cancer. Here we describe the properties of two mutations in RAD51 protein that derive from human lung and kidney tumors, respectively. Sequence variants Q268P and Q272L both map to the DNA binding loop 2 (L2) region of RAD51, a motif that is involved in DNA binding and in the allosteric activation of ATP hydrolysis and DNA strand exchange activities. Both mutations alter the thermal stability, DNA binding, and ATPase properties of RAD51, however both variants retain intrinsic DNA strand exchange activity towards oligonucleotide substrates under optimized conditions. In contrast, both Q268P and Q272L variants exhibit drastically reduced DNA strand exchange activity in reaction mixtures containing long homologous ssDNA and dsDNA substrates and human RPA protein. Mixtures of wild-type and variant proteins also exhibit reduced DNA strand exchange activity, suggesting that heterozygous mutations could negatively affect DNA recombination and repair processes in vivo. Together, the findings of this study suggest that hypomorphic missense mutations in RAD51 protein could be drivers of genomic instability in cancer cells, and thereby contribute to the etiology of metastatic disease.

Targeting ornithine decarboxylase reverses the LIN28/Let-7 axis and inhibits glycolytic metabolism in neuroblastoma

LIN28 has emerged as an oncogenic driver in a number of cancers, including neuroblastoma (NB). Overexpression of LIN28 correlates with poor outcome in NB, therefore drugs that impact the LIN28/Let-7 pathway could be beneficial in treating NB patients. The LIN28/Let-7 pathway affects many cellular processes including the regulation of cancer stem cells and glycolytic metabolism. Polyamines, regulated by ornithine decarboxylase (ODC) modulate eIF-5A which is a direct regulator of the LIN28/Let-7 axis. We propose that therapy inhibiting ODC will restore balance to the LIN28/Let-7 axis, suppress glycolytic metabolism, and decrease MYCN protein expression in NB. Difluoromethylornithine (DFMO) is an inhibitor of ODC in clinical trials for children with NB. In vitro experiments using NB cell lines, BE(2)-C, SMS-KCNR, and CHLA90 show that DFMO treatment reduced LIN28B and MYCN protein levels and increased Let-7 miRNA and decreased neurosphere formation. Glycolytic metabolic activity decreased with DFMO treatment in vivo. Additionally, sensitivity to DFMO treatment correlated with LIN28B overexpression (BE(2)-C>SMS-KCNR>CHLA90). This is the first study to demonstrate that DFMO treatment restores balance to the LIN28/Let-7 axis and inhibits glycolytic metabolism and neurosphere formation in NB and that PET scans may be a meaningful imaging tool to evaluate the therapeutic effects of DFMO treatment.

Characterization of Quinone Outside Inhibitor Fungicide Resistance in Cercospora sojina and Development of Diagnostic Tools for its Identification

Frogeye leaf spot of soybean, caused by the fungus Cercospora sojina, reduces soybean yields in most of the top-producing countries around the world. Control strategies for frogeye leaf spot can rely heavily on quinone outside inhibitor (QoI) fungicides. In 2010, QoI fungicide-resistant C. sojina isolates were identified in Tennessee for the first time. As the target of QoI fungicides, the cytochrome b gene present in fungal mitochondria has played a key role in the development of resistance to this fungicide class. The cytochrome b genes from three QoI-sensitive and three QoI-resistant C. sojina isolates were cloned and sequenced. The complete coding sequence of the cytochrome b gene was identified and found to encode 396 amino acids. The QoI-resistant C. sojina isolates contained the G143A mutation in the cytochrome b gene, a guanidine to cytosine transversion at the second position in codon 143 that causes an amino acid substitution of alanine for glycine. C. sojina-specific polymerase chain reaction primer sets and TaqMan probes were developed to efficiently discriminate QoI-resistant and -sensitive isolates. The molecular basis of QoI fungicide resistance in field isolates of C. sojina was identified as the G143A mutation, and specific molecular approaches were developed to discriminate and to track QoI-resistant and -sensitive isolates of C. sojina.

Feasibility of implementing molecular-guided therapy for the treatment of patients with relapsed or refractory neuroblastoma

The primary objective of the study was to evaluate the feasibility and safety of a process which would utilize genome-wide expression data from tumor biopsies to support individualized treatment decisions. Current treatment options for recurrent neuroblastoma are limited and ineffective, with a survival rate of <10%. Molecular profiling may provide data which will enable the practitioner to select the most appropriate therapeutic option for individual patients, thus improving outcomes. Sixteen patients with neuroblastoma were enrolled of which fourteen were eligible for this study. Feasibility was defined as completion of tumor biopsy, pathological evaluation, RNA quality control, gene expression profiling, bioinformatics analysis, generation of a drug prediction report, molecular tumor board yielding a treatment plan, independent medical monitor review, and treatment initiation within a 21 day period. All eligible biopsies passed histopathology and RNA quality control. Expression profiling by microarray and RNA sequencing were mutually validated. The average time from biopsy to report generation was 5.9 days and from biopsy to initiation of treatment was 12.4 days. No serious adverse events were observed and all adverse events were expected. Clinical benefit was seen in 64% of patients as stabilization of disease for at least one cycle of therapy or partial response. The overall response rate was 7% and the progression free survival was 59 days. This study demonstrates the feasibility and safety of performing real-time genomic profiling to guide treatment decision making for pediatric neuroblastoma patients.

Histamine H(3) receptor integrates peripheral inflammatory signals in the neurogenic control of immune responses and autoimmune disease susceptibility

Histamine H₃ receptor (Hrh3/H₃R) is primarily expressed by neurons in the central nervous system (CNS) where it functions as a presynaptic inhibitory autoreceptor and heteroreceptor. Previously, we identified an H₃R-mediated central component in susceptibility to experimental allergic encephalomyelitis (EAE), the principal autoimmune model of multiple sclerosis (MS), related to neurogenic control of blood brain barrier permeability and peripheral T cell effector responses. Furthermore, we identified Hrh3 as a positional candidate for the EAE susceptibility locus Eae8. Here, we characterize Hrh3 polymorphisms between EAE-susceptible and resistant SJL and B10.S mice, respectively, and show that Hrh3 isoform expression in the CNS is differentially regulated by acute peripheral inflammatory stimuli in an allele-specific fashion. Next, we show that Hrh3 is not expressed in any subpopulations of the immune compartment, and that secondary lymphoid tissue is anatomically poised to be regulated by central H₃R signaling. Accordingly, using transcriptome analysis, we show that, inflammatory stimuli elicit unique transcriptional profiles in the lymph nodes of H₃RKO mice compared to WT mice, which is indicative of negative regulation of peripheral immune responses by central H₃R signaling. These results further support a functional link between the neurogenic control of T cell responses and susceptibility to CNS autoimmune disease coincident with acute and/or chronic peripheral inflammation. Pharmacological targeting of H₃R may therefore be useful in preventing the development and formation of new lesions in MS, thereby limiting disease progression.

Milk yield responses to changes in milking frequency during early lactation are associated with coordinated and persistent changes in mammary gene expression

Background

The lactating mammary gland responds to changes in milking frequency by modulating milk production. This response is locally regulated and, in dairy cows, the udder is particularly sensitive during early lactation. Relative to cows milked twice-daily throughout lactation, those milked four-times-daily for just the first 3 weeks of lactation produce more milk throughout that lactation. We hypothesized that the milk yield response would be associated with increased mammary cell turnover and changes in gene expression during frequent milking and persisting thereafter. Cows were assigned to unilateral frequent milking (UFM; left udder halves milked twice-daily; right udder halves milked four-times daily) on days 1 to 21 of lactation, followed by twice-daily milking for the remainder of lactation. Relative to udder halves milked twice-daily, those milked four-times produced more milk during UFM; the difference in milk yield declined acutely upon cessation of UFM after day 21, but remained significantly elevated thereafter. We obtained mammary biopsies from both udder halves on days 21, 23, and 40 of lactation.

Results

Mammary cell proliferation and apoptosis were not affected by milking frequency. We identified 75 genes that were differentially expressed between paired udder halves on day 21 but exhibited a reversal of differential expression on day 23. Among those genes, we identified four clusters characterized by similar temporal patterns of differential expression. Two clusters (11 genes) were positively correlated with changes in milk yield and were differentially expressed on day 21 of lactation only, indicating involvement in the initial milk yield response. Two other clusters (64 genes) were negatively correlated with changes in milk yield. Twenty-nine of the 75 genes were also differentially expressed on day 40 of lactation.

Conclusions

Changes in milking frequency during early lactation did not alter mammary cell population dynamics, but were associated with coordinated changes in mammary expression of at least 75 genes. Twenty-nine of those genes were differentially expressed 19 days after cessation of treatment, implicating them in the persistent milk yield response. We conclude that we have identified a novel transcriptional signature that may mediate the adaptive response to changes in milking frequency.

PCI-24781 (abexinostat), a novel histone deacetylase inhibitor, induces reactive oxygen species-dependent apoptosis and is synergistic with bortezomib in neuroblastoma

In this study, we investigated the cytotoxic effects of a broad-spectrum histone deacetylase (HDAC) inhibitor, PCI-24781, alone and in combination with the proteasome inhibitor bortezomib in neuroblastoma cell lines. The combination was shown to induce synergistic cytotoxity involving the formation of reactive oxygen species. The cleavage of caspase-3 and PARP, as determined by western blotting, indicated that cell death was primarily due to apoptosis. Xenograft mouse models indicated increased survival among animals treated with this combination. The Notch signaling pathway and MYCN gene expression were quantified by reverse transcription-polymerase chain reaction (PCR) in cells treated with PCI-24781 and bortezomib, alone and in combination. Notch pathway expression increased in response to an HDAC inhibitor. NFKB1 and MYCN were both significantly down regulated. Our results suggest that PCI-24781 and bortezomib are synergistic in neuroblastoma cell lines and may be a new therapeutic strategy for this disease.

Identification of Orch3, a locus controlling dominant resistance to autoimmune orchitis, as kinesin family member 1C

Experimental autoimmune orchitis (EAO), the principal model of non-infectious testicular inflammatory disease, can be induced in susceptible mouse strains by immunization with autologous testicular homogenate and appropriate adjuvants. As previously established, the genome of DBA/2J mice encodes genes that are capable of conferring dominant resistance to EAO, while the genome of BALB/cByJ mice does not and they are therefore susceptible to EAO. In a genome scan, we previously identified Orch3 as the major quantitative trait locus controlling dominant resistance to EAO and mapped it to chromosome 11. Here, by utilizing a forward genetic approach, we identified kinesin family member 1C (Kif1c) as a positional candidate for Orch3 and, using a transgenic approach, demonstrated that Kif1c is Orch3. Mechanistically, we showed that the resistant Kif1c^D2 allele leads to a reduced antigen-specific T cell proliferative response as a consequence of decreased MHC class II expression by antigen presenting cells, and that the L⁵⁷⁸→P⁵⁷⁸ and S¹⁰²⁷→P¹⁰²⁷ polymorphisms distinguishing the BALB/cByJ and DBA/2J alleles, respectively, can play a role in transcriptional regulation. These findings may provide mechanistic insight into how polymorphism in other kinesins such as KIF21B and KIF5A influence susceptibility and resistance to human autoimmune diseases.

Although the etiology of autoimmunity is not well known, a variety of studies have demonstrated that genetic predisposition is a major contributor to disease susceptibility and resistance. The major histocompatibility complex (MHC) is the primary genetic determinant of autoimmune disease susceptibility with multiple additional interacting loci required. However, the identification and characterization of non–MHC genes has been problematic, since most autoimmune diseases are polygenic with the individual genes exhibiting only partial or minimal penetrance. We previously identified Orch3 (mouse chromosome 11) as the most important immune-suppressive locus controlling dominant resistance to autoimmune orchitis, the principal animal model of non-infectious testicular inflammatory/autoimmune disease. Here, using congenic mapping, we identified kinesin family member 1C (Kif1c) as a positional candidate for Orch3. Furthermore, over-expression of the Kif1c resistant allele in susceptible mice rendered animals autoimmune orchitis resistant, demonstrating that Kif1c is Orch3. We propose that Kif1c plays an immunoregulatory role by controlling the levels of MHC class II in antigen presenting cells and consequently impacting autoreactive orchitogenic T cell responses. These finding are particularly relevant since polymorphism in other kinesins such as KIF21B and KIF5A have been associated with susceptibility and resistance to human autoimmune disease.

Classifying variants of CDKN2A using computational and laboratory studies

Variants in the CDKN2A tumor suppressor are associated with Familial Melanoma (FM), although for many variants the linkage is weak. The effects of missense variants on protein function and pathogenicity are often unclear. Multiple methods (e.g., laboratory, computational, epidemiological) have been developed to analyze whether a missense variant is pathogenic or not. It is not yet clear how to integrate these data types into a strategy for variant classification. We studied 51 CDKN2A missense variants using a cell cycle arrest assay. There was a continuum of results ranging from full wild-type effect through partial activity to complete loss of arrest. A reproducible decrease of 30% of cell cycle arrest activity correlated with FM association. We analyzed missense CDKN2A germline variants using a Bayesian method to combine multiple data types and derive a probability of pathogenicity. When equal to or more than two data types could be evaluated with this method, 22 of 25 FM-associated variants and 8 of 15 variants of uncertain significance were classified as likely pathogenic with >95% probability. The other 10 variants were classified as uncertain (probability 5-95%). For most variants, there were insufficient data to draw a conclusion. The Bayesian model appears to be a sound method of classifying missense variants in cancer susceptibility genes.

Mammary gland morphological and gene expression changes underlying pregnancy protection of breast cancer tumorigenesis

A full-term pregnancy early in life reduces lifetime risk of developing breast cancer, and the effect can be mimicked in rodents by full-term pregnancy or short-term treatment with exogenous estrogen and progesterone. To gain insight into the protective mechanism, 15 3-mo-old postpubertal virgin Lewis rats were randomly assigned to three groups: control (C), pregnancy (P), or hormone (H). The P group animals underwent a full-term pregnancy, and H group animals were implanted subcutaneously with silastic capsules filled with ethynyl estradiol and megesterol acetate for 21 days. C and P animals were implanted with sham capsules. On day 21 capsules were removed, which was followed by a 49-day involution period, euthanasia, and mammary tissue collection. Global gene expression was measured using Rat Genome 230.2 Arrays. Histological analysis revealed that P and H treatments induced sustained morphological changes in the mammary gland with significantly increased percentages of mammary parenchyma and stromal tissues and higher ratio of stroma to parenchyma. Transcriptome analysis showed that P and H treatments induced sustained global changes in gene expression in the mammary gland. Analysis of commonly up- and downregulated genes in P and H relative to C treatment showed increased expression of three matrix metallopeptidases (Mmp3, 8, and 12), more differentiated mammary phenotype, enhanced innate and adaptive immunity, and reduced cell proliferation and angiogenic signatures. The sustained morphological and global gene expression changes in mammary tissue after pregnancy and hormone treatment may function together to provide the protective effect against breast cancer.

Acute milk yield response to frequent milking during early lactation is mediated by genes transiently regulated by milk removal

Milking dairy cows four times daily (4×) instead of twice daily (2×) during early lactation stimulates an increase in milk yield that partly persists through late lactation; however, the mechanisms behind this response are unknown. We hypothesized that the acute mammary response to regular milkings would be transient and would involve different genes from those that may be specifically regulated in response to 4×. Nine multiparous cows were assigned at parturition to unilateral frequent milking (UFM; 2× of the left udder half, 4× of the right udder half). Mammary biopsies were obtained from both rear quarters at 5 days in milk (DIM), immediately after 4× glands had been milked (experiment 1, n = 4 cows), or 2.5 h after both udder halves had last been milked (experiment 2, n = 5 cows). Affymetrix GeneChip Bovine Genome Arrays were used to measure gene expression. We found 855 genes were differentially expressed in mammary tissue between 2× vs. 4× glands of cows in experiment 1 (false discovery rate ≤ 0.05), whereas none were differentially expressed in experiment 2 using the same criterion. We conclude that there is an acute transcriptional response to milk removal, but 4× milking did not elicit differential expression of unique genes. Therefore, there does not appear to be a sustained transcriptional response to 4× milking on day 5 of lactation. Using a differential expression plot of data from both experiments, as well as qRT-PCR, we identified at least two genes (chitinase 3-like-1 and low-density lipoprotein-related protein-2 that may be responsive to both milk removal and to 4× milking. Therefore, the milk yield response to 4× milking may be mediated by genes that are acutely regulated by removal of milk from the mammary gland.

Using shifts in amino acid frequency and substitution rate to identify latent structural characters in base-excision repair enzymes

Protein evolution includes the birth and death of structural motifs. For example, a zinc finger or a salt bridge may be present in some, but not all, members of a protein family. We propose that such transitions are manifest in sequence phylogenies as concerted shifts in substitution rates of amino acids that are neighbors in a representative structure. First, we identified rate shifts in a quartet from the Fpg/Nei family of base excision repair enzymes using a method developed by Xun Gu and coworkers. We found the shifts to be spatially correlated, more precisely, associated with a flexible loop involved in bacterial Fpg substrate specificity. Consistent with our result, sequences and structures provide convincing evidence that this loop plays a very different role in other family members. Second, then, we developed a method for identifying latent protein structural characters (LSC) given a set of homologous sequences based on Gu's method and proximity in a high-resolution structure. Third, we identified LSC and assigned states of LSC to clades within the Fpg/Nei family of base excision repair enzymes. We describe seven LSC; an accompanying Proteopedia page (http://proteopedia.org/wiki/index.php/Fpg_Nei_Protein_Family) describes these in greater detail and facilitates 3D viewing. The LSC we found provided a surprisingly complete picture of the interaction of the protein with the DNA capturing familiar examples, such as a Zn finger, as well as more subtle interactions. Their preponderance is consistent with an important role as phylogenetic characters. Phylogenetic inference based on LSC provided convincing evidence of independent losses of Zn fingers. Structural motifs may serve as important phylogenetic characters and modeling transitions involving structural motifs may provide a much deeper understanding of protein evolution.

Abstract 5260: Bortezomib and HDAC inhibitor PCI-24781 combination therapy for neuroblastoma

Background: Neuroblastoma is an aggressive solid tumor with poor prognosis. Novel therapies are needed to help increase the survival and quality of life for affected children. It has been previously shown that bortezomib results in apoptosis in neuroblastoma both in vitro and in vivo. PCI-24781 is a novel pan-HDAC inhibitor which has been shown to be synergistic with bortezomib in lymphoma models. We hypothesize the combination of bortezomib with PCI-24781 would increase cytotoxicity to neuroblastoma.

Methods: Several established neuroblastoma cell lines as well as patient-derived primary neuroblastoma cultures were grown in 96-well plates and treated with bortezomib and PCI-24781, both separately and in combination for 48 hours. Cell viability was assessed by calcein AM assays and IC50's and drug synergy was determined at each cell line tested. mRNA was collected from drug treated cell lines and effects on cell signaling was evaluated using U133 mRNA expression arrays and Ingenuity analysis. Cell lysates prepared from drug treated cells were evaluated by western blot for caspase-3 and PARP cleavage. Nude mice were injected with 107 SMS-KCNR cells subcutaneously in the left flank. Mice were treated with daily doses of either 0.5 mg/kg bortezomib, 12.5 mg/kg PCI-24781, or the combination. Mice tumors were imaged using the IVIS lumina imager twice per week and caliper measurements were also obtained weekly until tumor max was reached.

Results: All neuroblastoma cell lines and patient cells tested showed sensitivity to bortezomib and PCI-24781 treatment in calcien AM cell viability assays with IC50's for bortezomib less than 50nM and IC50's for PCI-24781 less than 500nM. The combination of bortezomib and PCI-24781 was more cytotoxic than either drug alone. Cells treated with combined treatment showed an increase in caspase-3 and PARP activity confirming apoptosis. Expression analysis showed differential expression in the Notch and Wnt pathways between drug treated and vehicle treated cells. The xenograft models showed a significant decrease in tumor volume in mice treated with both bortezomib and PCI-24781 when compared to the single agent treatment groups or to the control group. Conclusions: Bortezomib and PCI-24781 inhibit neuroblastoma growth both in vitro and in vivo to a greater extent than either drug alone. The mechanism of action is currently being investigated further but certainly involves caspase-3 mediated apoptosis. This combination therapy appears to be effective and well tolerated in the mouse model and would be a novel therapy for neuroblastoma.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 5260.

Utilization of gene profiling and proteomics to determine mineral pathogenicity in a human mesothelial cell line (LP9/TERT-1)

Identifying and understanding the early molecular events that underscore mineral pathogenicity using in vitro screening tests is imperative, especially given the large number of synthetic and natural fibers and particles being introduced into the environment. The purpose of the work described here was to examine the ability of gene profiling (Affymetrix microarrays) to predict the pathogenicity of various materials in a human mesothelial cell line (LP9/TERT-1) exposed to equal surface area concentrations (15 x 10(6) or 75 x 10(6) microm(2)/cm(2)) of crocidolite asbestos, nonfibrous talc, fine titanium dioxide (TiO(2)), or glass beads for 8 or 24 h. Since crocidolite asbestos caused the greatest number of alterations in gene expression, multiplex analysis (Bio-Plex) of proteins released from LP9/TERT-1 cells exposed to crocidolite asbestos was also assessed to reveal if this approach might also be explored in future assays comparing various mineral types. To verify that LP9/TERT-1 cells were more sensitive than other cell types to asbestos, human ovarian epithelial cells (IOSE) were also utilized in microarray studies. Upon assessing changes in gene expression via microarrays, principal component analysis (PCA) of these data was used to identify patterns of differential gene expression. PCA of microarray data confirmed that LP9/TERT-1 cells were more responsive than IOSE cells to crocidolite asbestos or nonfibrous talc, and that crocidolite asbestos elicited greater responses in both cell types when compared to nonfibrous talc, TiO(2), or glass beads. Bio-Plex analysis demonstrated that asbestos caused an increase in interleukin-13 (IL-13), basic fibroblast growth factor (bFGF), granulocyte colony-stimulating factor (G-CSF), and vascular endothelial growth factor (VEGF). These responses were generally dose-dependent (bFGF and G-CSF only) and tumor necrosis factor (TNF)-alpha independent (except for G-CSF). Thus, microarray and Bio-Plex analyses are valuable in determining early molecular responses to fibers/particles and may directly contribute to understanding the etiology of diseases caused by them. The number and magnitude of changes in gene expression or "profiles" of secreted proteins may serve as valuable metrics for determining the potential pathogenicity of various mineral types. Hence, alterations in gene expression and cytokine/chemokine changes induced by crocidolite asbestos in LP9/TERT-1 cells may be indicative of its increased potential to cause mesothelioma in comparison to the other nonfibrous materials examined.

Alterations in gene expression in human mesothelial cells correlate with mineral pathogenicity

Human mesothelial cells (LP9/TERT-1) were exposed to low and high (15 and 75 microm(2)/cm(2) dish) equal surface area concentrations of crocidolite asbestos, nonfibrous talc, fine titanium dioxide (TiO2), or glass beads for 8 or 24 hours. RNA was then isolated for Affymetrix microarrays, GeneSifter analysis and QRT-PCR. Gene changes by asbestos were concentration- and time-dependent. At low nontoxic concentrations, asbestos caused significant changes in mRNA expression of 29 genes at 8 hours and of 205 genes at 24 hours, whereas changes in mRNA levels of 236 genes occurred in cells exposed to high concentrations of asbestos for 8 hours. Human primary pleural mesothelial cells also showed the same patterns of increased gene expression by asbestos. Nonfibrous talc at low concentrations in LP9/TERT-1 mesothelial cells caused increased expression of 1 gene Activating Transcription Factor 3 (ATF3) at 8 hours and no changes at 24 hours, whereas expression levels of 30 genes were elevated at 8 hours at high talc concentrations. Fine TiO2 or glass beads caused no changes in gene expression. In human ovarian epithelial (IOSE) cells, asbestos at high concentrations elevated expression of two genes (NR4A2, MIP2) at 8 hours and 16 genes at 24 hours that were distinct from those elevated in mesothelial cells. Since ATF3 was the most highly expressed gene by asbestos, its functional importance in cytokine production by LP9/TERT-1 cells was assessed using siRNA approaches. Results reveal that ATF3 modulates production of inflammatory cytokines (IL-1 beta, IL-13, G-CSF) and growth factors (VEGF and PDGF-BB) in human mesothelial cells.

Plant and fungal Fpg homologs are formamidopyrimidine DNA glycosylases but not 8-oxoguanine DNA glycosylases

Formamidopyrimidine DNA glycosylase (Fpg) and endonuclease VIII (Nei) share an overall common three-dimensional structure and primary amino acid sequence in conserved structural motifs but have different substrate specificities, with bacterial Fpg proteins recognizing formamidopyrimidines, 8-oxoguanine (8-oxoG) and its oxidation products guanidinohydantoin (Gh), and spiroiminodihydantoin (Sp) and bacterial Nei proteins recognizing primarily damaged pyrimidines. In addition to bacteria, Fpg has also been found in plants, while Nei is sparsely distributed among the prokaryotes and eukaryotes. Phylogenetic analysis of Fpg and Nei DNA glycosylases demonstrated, with 95% bootstrap support, a clade containing exclusively sequences from plants and fungi. Members of this clade exhibit sequence features closer to bacterial Fpg proteins than to any protein designated as Nei based on biochemical studies. The Candida albicans (Cal) Fpg DNA glycosylase and a previously studied Arabidopsis thaliana (Ath) Fpg DNA glycosylase were expressed, purified and characterized. In oligodeoxynucleotides, the preferred glycosylase substrates for both enzymes were Gh and Sp, the oxidation products of 8-oxoG, with the best substrate being a site of base loss. GC/MS analysis of bases released from gamma-irradiated DNA show FapyAde and FapyGua to be excellent substrates as well. Studies carried out with oligodeoxynucleotide substrates demonstrate that both enzymes discriminated against A opposite the base lesion, characteristic of Fpg glycosylases. Single turnover kinetics with oligodeoxynucleotides showed that the plant and fungal glycosylases were most active on Gh and Sp, less active on oxidized pyrimidines and exhibited very little or no activity on 8-oxoG. Surprisingly, the activity of AthFpg1 on an AP site opposite a G was extremely robust with a k(obs) of over 2500min(-1).

The induction of antibody production by IL-6 is indirectly mediated by IL-21 produced by CD4+ T cells

Interleukin (IL) 6 is a proinflammtory cytokine produced by antigen-presenting cells and nonhematopoietic cells in response to external stimuli. It was initially identified as a B cell growth factor and inducer of plasma cell differentiation in vitro and plays an important role in antibody production and class switching in vivo. However, it is not clear whether IL-6 directly affects B cells or acts through other mechanisms. We show that IL-6 is sufficient and necessary to induce IL-21 production by naive and memory CD4⁺ T cells upon T cell receptor stimulation. IL-21 production by CD4⁺ T cells is required for IL-6 to promote B cell antibody production in vitro. Moreover, administration of IL-6 with inactive influenza virus enhances virus-specific antibody production, and importantly, this effect is dependent on IL-21. Thus, IL-6 promotes antibody production by promoting the B cell helper capabilities of CD4⁺ T cells through increased IL-21 production. IL-6 could therefore be a potential coadjuvant to enhance humoral immunity.

Clostridium acetobutylicum 8-oxoguanine DNA glycosylase (Ogg) differs from eukaryotic Oggs with respect to opposite base discrimination

During repair of damaged DNA, the oxidized base 8-oxoguanine (8-oxoG) is removed by 8-oxoguanine-DNA glycosylase (Ogg) in eukaryotes and most archaea, whereas in most bacteria it is removed by formamidopyrimidine-DNA glycosylase (Fpg). We report the first characterization of a bacterial Ogg, Clostridium acetobutylicum Ogg (CacOgg). Like human OGG1 and Escherichia coli Fpg (EcoFpg), CacOgg excised 8-oxoguanine. However, unlike hOGG1 and EcoFpg, CacOgg showed little preference for the base opposite the damage during base excision and removed 8-oxoguanine from single-stranded DNA. Thus, our results showed unambiguous qualitative functional differences in vitro between CacOgg and both hOGG1 and EcoFpg. CacOgg differs in sequence from the eukaryotic enzymes at two sequence positions, M132 and F179, which align with amino acids (R154 and Y203) in human OGG1 (hOGG1) found to be involved in opposite base interaction. To address the sequence basis for functional differences with respect to opposite base interactions, we prepared three CacOgg variants, M132R, F179Y, and M132R/F179Y. All three variants showed a substantial increase in specificity for 8-oxoG.C relative to 8-oxoG.A. While we were unable to definitively associate these qualitative functional differences with differences in selective pressure between eukaryotes, Clostridia, and other bacteria, our results are consistent with the idea that evolution of Ogg function is based on kinetic control of repair.

Microarray analyses of pituitary adenylate cyclase activating polypeptide (PACAP)-regulated gene targets in sympathetic neurons

The high and preferential expression of the PAC(1)(short)HOP1 receptor in postganglionic sympathetic neurons facilitates microarray studies for mechanisms underlying PACAP-mediate neurotrophic signaling in a physiological context. Replicate primary sympathetic neuronal cultures were treated with 100 nM PACAP27 either acutely (9 h) or chronically (96 h) before RNA extraction and preparation for Affymetrix microarray analysis. Compared to untreated control cultures, acute PACAP treatment modulated significantly the expression of 147 transcripts of diverse functional groups, including peptides, growth factors/cytokines, transcriptional factors, receptors/signaling effectors and cell cycle regulators, that collectively appeared to facilitate neuronal plasticity, differentiation and/or regeneration processes. Some regulated transcripts, for example, were related to BDNF/TrkB, IL-6/Jak2/Socs2 and TGF/follistatin signaling; many transcripts affected bioactive peptide and polyamine biosynthesis. Although chronic PACAP treatments altered the expression of 109 sympathetic transcripts, only 43 transcripts were shared between the acute and chronic treatment data sets. The PACAP-mediated changes in transcript expression were corroborated independently by quantitative PCR measurement. The PACAP-regulated transcripts in sympathetic neurons did not bear strong resemblance to those in PACAP-treated pheochromocytoma cells. However, many PACAP-targeted sympathetic transcripts, especially those related to peptide plasticity and nerve regeneration processes, coincided significantly with genes altered after peripheral nerve injury. The ability for sympathetic PAC(1)(short)HOP1 receptors to engage multiple downstream signaling cascades appeared to be reflected in the number and diversity of genes targeted in a multifaceted strategy for comprehensive neurotrophic responses.

Estimation of DNA sequence context-dependent mutation rates using primate genomic sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.

Interpreting missense variants: comparing computational methods in human disease genes CDKN2A, MLH1, MSH2, MECP2, and tyrosinase (TYR)

The human genome contains frequent single-basepair variants that may or may not cause genetic disease. To characterize benign vs. pathogenic missense variants, numerous computational algorithms have been developed based on comparative sequence and/or protein structure analysis. We compared computational methods that use evolutionary conservation alone, amino acid (AA) change alone, and a combination of conservation and AA change in predicting the consequences of 254 missense variants in the CDKN2A (n = 92), MLH1 (n = 28), MSH2 (n = 14), MECP2 (n = 30), and tyrosinase (TYR) (n = 90) genes. Variants were validated as either neutral or deleterious by curated locus-specific mutation databases and published functional data. All methods that use evolutionary sequence analysis have comparable overall prediction accuracy (72.9-82.0%). Mutations at codons where the AA is absolutely conserved over a sufficient evolutionary distance (about one-third of variants) had a 91.6 to 96.8% likelihood of being deleterious. Three algorithms (SIFT, PolyPhen, and A-GVGD) that differentiate one variant from another at a given codon did not significantly improve predictive value over conservation score alone using the BLOSUM62 matrix. However, when all four methods were in agreement (62.7% of variants), predictive value improved to 88.1%. These results confirm a high predictive value for methods that use evolutionary sequence conservation, with or without considering protein structural change, to predict the clinical consequences of missense variants. The methods can be generalized across genes that cause different types of genetic disease. The results support the clinical use of computational methods as one tool to help interpret missense variants in genes associated with human genetic disease.

Methylation-controlled J protein promotes c-Jun degradation to prevent ABCB1 transporter expression

Methylation-controlled J protein (MCJ) is a newly identified member of the DnaJ family of cochaperones. Hypermethylation-mediated transcriptional silencing of the MCJ gene has been associated with increased chemotherapeutic resistance in ovarian cancer. However, the biology and function of MCJ remain unknown. Here we show that MCJ is a type II transmembrane cochaperone localized in the Golgi network and present only in vertebrates. MCJ is expressed in drug-sensitive breast cancer cells but not in multidrug-resistant cells. The inhibition of MCJ expression increases resistance to specific drugs by inducing expression of the ABCB1 drug transporter that prevents intracellular drug accumulation. The induction of ABCB1 gene expression is mediated by increased levels of c-Jun due to an impaired degradation of this transcription factor in the absence of MCJ. Thus, MCJ is required in these cells to prevent c-Jun-mediated expression of ABCB1 and maintain drug response.

Mechanism of presynaptic filament stabilization by the bacteriophage T4 UvsY recombination mediator protein

UvsY is the recombination mediator protein (RMP) of bacteriophage T4, which promotes homologous recombination by facilitating presynaptic filament assembly. The results of previous studies suggest that UvsY promotes the assembly of presynaptic filaments in part by stabilizing interactions between T4 UvsX recombinase and single-stranded DNA (ssDNA). To test this hypothesis, we studied the interactions of UvsX and UvsY with a fluorescein-derivatized oligonucleotide. This assay distinguishes between bipartite UvsX- or UvsY-ssDNA and tripartite UvsX-UvsY-ssDNA complex formation via differential fluorescence quenching effects. Salt stabilities of the three complexes were measured at equilibrium in the presence and absence of various nucleotide ligands of the UvsX protein and also under steady-state conditions for UvsX-catalyzed ssDNA-dependent ATP hydrolysis. The results demonstrate that UvsY globally stabilizes UvsX-ssDNA complexes, consistent with an increase in the apparent equilibrium binding affinity, K(ss)omega, of the UvsX-ssDNA interactions. The UvsY-mediated affinity increase is observed at equilibrium in the presence of ADP, ATPgammaS, or in the absence of the nucleotide and also at steady-state in the presence of ATP. Intriguingly, the stabilizing effects of UvsY and ATPgammaS on UvsX-ssDNA interactions are synergistic, indicating nonredundant mechanisms for UvsX-ssDNA complex stabilization by RMP versus nucleoside triphosphate effectors. Experiments with UvsY missense mutants defective in ssDNA binding demonstrate that UvsY-ssDNA interactions are of major importance in stabilizing UvsX-ssDNA complexes, whereas UvsY-UvsX protein-protein interactions provide residual stabilization energy. Together, the data is consistent with a mechanism in which UvsY stabilizes presynaptic filaments by organizing the ssDNA lattice into a structure that is favorable for UvsX-ssDNA interactions.

Engineering functional changes in Escherichia coli endonuclease III based on phylogenetic and structural analyses

Escherichia coli endonuclease III (EcoNth) plays an important cellular role by removing premutagenic pyrimidine damages produced by reactive oxygen species. EcoNth is a bifunctional enzyme that has DNA glycosylase and apurinic/apyrimidinic lyase activities. Using a phylogeny of natural sequences, we selected to study EcoNth serine 39, aspartate 44, and arginine 184, which are presumed to be in the vicinity of the damaged base in the glycosylase-substrate complex. These three amino acids are highly conserved among Nth orthologs, although not among homologous glycosylases, such as MutY, that have different base specificities and no lyase activity. To examine the role of these amino acids in catalysis, we constructed three mutants of EcoNth, in which Ser39 was replaced with leucine (S39L), Asp44 was replaced with valine (D44V), and Arg184 was replaced with alanine (R184A), which are the corresponding residues in EcoMutY. We showed that EcoNth S39L does not have significant glycosylase activity for oxidized pyrimidines, although it maintained AP lyase activity. In contrast, EcoNth D44V retained glycosylase activity against oxidized pyrimidines, but the apparent rate constant for the lyase activity of EcoNth D44V was significantly lower than that of EcoNth, indicating that Asp44 in EcoNth is required for beta-elimination. Finally, EcoNth R184A maintained lyase activity but exhibited glycosylase specificity different from that of EcoNth. The functional consequences of each of these three substitutions can be rationalized in the context of high resolution protein structures. Thus phylogeny-based scanning mutagenesis has allowed us to identify novel roles for amino acids in the substrate binding pocket of EcoNth in base recognition and/or catalysis.

The crystal structure of human endonuclease VIII-like 1 (NEIL1) reveals a zincless finger motif required for glycosylase activity

In prokaryotes, two DNA glycosylases recognize and excise oxidized pyrimidines: endonuclease III (Nth) and endonuclease VIII (Nei). The oxidized purine 8-oxoguanine, on the other hand, is recognized by Fpg (also known as MutM), a glycosylase that belongs to the same family as Nei. The recent availability of the human genome sequence allowed the identification of three human homologs of Escherichia coli Nei. We report here the crystal structure of a human Nei-like (NEIL) enzyme, NEIL1. The structure of NEIL1 exhibits the same overall fold as E. coli Nei, albeit with an unexpected twist. Sequence alignments had predicted that NEIL1 would lack a zinc finger, and it was therefore expected to use a different DNA-binding motif instead. Our structure revealed that, to the contrary, NEIL1 contains a structural motif composed of two antiparallel beta-strands that mimics the antiparallel beta-hairpin zinc finger found in other Fpg/Nei family members but lacks the loops that harbor the zinc-binding residues and, therefore, does not coordinate zinc. This "zincless finger" appears to be required for NEIL1 activity, because mutating a very highly conserved arginine within this motif greatly reduces the glycosylase activity of the enzyme.

The CDKN2A database: Integrating allelic variants with evolution, structure, function, and disease association

In this report, we introduce the CDKN2A Database, an online database of germline and somatic variants of the CDKN2A tumor suppressor gene recorded in human disease through the year 2002, annotated with evolutionary, structural, and functional information. The CDKN2A Database improves upon existing resources by: 1) including both somatic mutations and germline variants, thereby adding the perspective of somatic cell carcinogenesis to that of hereditary cancer predisposition; 2) including information that assists with the interpretation of allelic variants, such as other primary data (sequences, structures, alignments, functional measurements, and literature references) and annotations (extensive text, figures, and a tree-based phylogenetic classification); and 3) providing the information in a format that allows a user to either download the database or to easily manipulate it online. We describe the database structure, content, current uses, and potential implications (http://biodesktop.uvm.edu/perl/p16).

Response of Soybean Sudden Death Syndrome to Subsoil Tillage

The soilborne pathogen Fusarium solani f. sp. glycines causes sudden death syndrome (SDS) of soybean. Previous research indicated that soil compaction related directly to disease foliar symptoms. Therefore, we hypothesized that decreasing soil compaction would increase soil porosity and provide a more aerated root zone that would hinder root infection by the fungus and decrease SDS foliar symptom severity. Two experimental areas (110 by 120 m) were established to evaluate the relationship between soil variables and SDS. Across the experimental area, strips (9.14 m wide) were subsoiled perpendicular to soybean rows to a depth of 40 to 45 cm, which alternated with strips that were not tilled. In both 1999 and 2000, subsoiling dramatically reduced foliar symptoms of SDS. Compared with no-till plots, subsoiled plots had lower soil bulk density, greater soil porosity, and less soil moisture. In areas where SDS occurs and soil compaction exists, the use of subsoiling can be used to reduce severity of foliar symptoms of SDS.

Detailed computational study of p53 and p16: using evolutionary sequence analysis and disease-associated mutations to predict the functional consequences of allelic variants

Deciding whether a missense allelic variant affects protein function is important in many contexts. We previously demonstrated that a detailed analysis of p53 intragenic conservation correlates with somatic mutation hotspots. Here we refine these evolutionary studies and expand them to the p16/Ink4a gene. We calculated that in order for 'absolute conservation' of a codon across multiple species to achieve P<0.05, the evolutionary substitution database must contain at least 3(M) variants, where M equals the number of codons in the gene. Codons in p53 were divided into high (73% of codons), intermediate (29% of codons), and low (0 codons) likelihood of being mutation hotspots. From a database of 263 somatic missense p16 mutations, we identified only four codons that are mutational hotspots at P<0.05 (8 mutations). However, data on function, structure, and disease association support the conclusion that 11 other codons with > or =5 somatic mutations also likely indicate functionally critical residues, even though P0.05. We calculated p16 evolution using amino acid substitution matrices and nucleotide substitution distances. We looked for evolutionary parameters at each codon that would predict whether missense mutations were disease associated or disrupted function. The current p16 evolutionary substitution database is too small to determine whether observations of 'absolute conservation' are statistically significant. Increasing the number of sequences from three to seven significantly improved the predictive value of evolutionary computations. The sensitivity and specificity for conservation scores in predicting disease association of p16 codons is 70-80%. Despite the small p16 sequence database, our calculations of high conservation correctly predicted loss of cell cycle arrest function in 75% of tested codons, and low conservation correctly predicted wild-type function in 80-90% of codons. These data validate our hypothesis that detailed evolutionary analyses help predict the consequences of missense amino-acid variants.

A novel human DNA glycosylase that removes oxidative DNA damage and is homologous to Escherichia coli endonuclease VIII

Prokaryotes and lower eukaryotes possess redundant activities that remove the plethora of oxidative DNA base damages produced during normal oxidative metabolism and which have been associated with cancer and aging. Thus far, only one oxidized pyrimidine-specific DNA glycosylase has been identified in humans, hNthl. Here, we report the identification of three new putative human DNA glycosylases that are phylogenetically members of the Fpg/Nei family primarily found in the bacterial kingdom. We have characterized one of these, hNEI1, and show it to be functionally homologous to bacterial Nei, that is, its principal substrates are oxidized pyrimidines, it undergoes a lyase reaction by, beta,delta-elimination and traps a Schiff base with a substrate containing thymine glycol (Tg). Furthermore, inactivation of active site residues shown to be important in Escherichia coli Nei inactivate the human enzyme. The hNEI1 gene is located on the long arm of chromosome 15 that is frequently deleted in human cancers.

Expression and purification of the extracellular domain of human myelin protein zero

Myelin protein zero (P0), an adhesion protein of the immunoglobulin superfamily, is the major protein of peripheral nervous system myelin in higher vertebrates. Protein zero is required for the formation and maintenance of myelin structure in the internode, likely through homophilic interactions at both the extracellular and the intracellular domains. Mutations and deletions in the P0 gene correlate with hereditary peripheral neuropathies of varying severity. Comparisons between the human and rat isoforms, whose three-dimensional structure has been determined by X-ray crystallography, suggest that these disease-associated genetic alterations lead to structural changes in the protein that alter P0-P0 interactions and hence affect myelin functionality. Knowing the crystal structures of native and altered human P0 isoforms could help to elucidate the structural changes in myelin membrane packing that underlie the altered functionality. Alterations of P0 extracellular domain (P0-ED) are of additional interest as previous X-ray diffraction studies on myelin membrane packing suggest that P0-ED molecules can assume distinct adhesive arrangements. Here, we describe an improved method to express and purify human P0-ED (hP0-ED) suitable for crystallographic analysis. A fusion protein consisting of maltose binding protein fused to hP0-ED was secreted to the periplasm of Escherichia coli to allow an appropriate folding pathway. The fusion protein was extracted via osmotic shock and purified by affinity chromatography. Factor Xa was used to cleave the fusion protein, and a combination of affinity and ion-exchange chromatography was used to further purify hP0-ED. We document several significant improvements to previous protocols, including bacterial growth to approximately 15 OD using orbital shakers and the use of diafiltration, which result in yields of approximately 150 mg highly pure protein per liter of medium.

TP53 mutations in breast cancer associated with BRCA1 or BRCA2 germ-line mutations: distinctive spectrum and structural distribution

Several groups have studied the molecular pathology of inherited breast cancer. By combining several such studies, we show in this study that somatic TP53 abnormalities are more common in breast cancer associated with BRCA1 or BRCA2 germ-line mutations than in sporadic breast cancers (odds ratio, 2.8; P = 0.0003). Then, we compared the spectrum of TP53 mutations for breast cancers in the IARC TP53 mutation database with the 82 mutations reported in BRCA1/2-associated breast cancers. The spectrum differed significantly both in distribution (P < 1 x 10(-6)) and in base changes (P = 0.025). Mutations at A:T bp were more common in BRCA1/2-associated tumors and strand bias suggesting DNA repair abnormalities was found. Changes were common at TP53 codons that are not mutation hotspots. Structural modeling showed that most of these p53 non-hotspot amino acids characterized in breast tumors isolated from patients with deficient BRCA1/2 function are distributed in a region of the protein on the opposite side of the p53 DNA-binding surface. Our results suggest that BRCA1/2 mutations influence the type and distribution of TP53 mutations seen in breast cancer.

Distribution of plant-parasitic nematodes on sugarcane in louisiana and efficacy of nematicides

A survey conducted from May 1995 through August 1998 revealed diverse nematode communities in Louisiana sugarcane fields. High populations of Mesocriconema, Paratrichodorus, Pratylenchus, and Tylenchorhynchus were widespread in nine sugarcane production parishes. Comparisons of plant cane and ratoon sugarcane crops indicated that nematode community levels increase significantly in successive ratoon crops. Nematicide trials evaluated the efficacy of aldicarb, ethoprop, and phorate against indigenous nematode populations. Aldicarb consistently increased the number of millable stalks, cane tonnage, and yield of sucrose in soils with a high sand content. Yield increases were concomitant with reductions in the density of the nematode community shortly after planting and at harvest. In soils with a higher clay content, the chemicals were less effective in controlling nematode populations and, as a result, yield increases were minimal.

Proteobacterial histidine-biosynthetic pathways are paraphyletic

In Lactococcus lactis there is a protein, HisZ, in the histidine-biosynthetic operon that exhibits significant sequence identity with histidyl-tRNA synthetase (HisRS) but does not aminoacylate tRNA. HisRS homologs that, like HisZ, cannot aminoacylate tRNA are represented in a highly divergent set of bacteria (including an aquificale, cyanobacteria, firmicutes, and proteobacteria), yet are missing from other bacteria, including mycrobacteria and certain proteobacteria. Phylogenetic analysis of the HisRS and HisRS-like family suggests that the HisZ proteins form a monophyletic group that attaches outside the predominant bacterial HisRS clade. These observations are consistent with a model in which the absences of HisZ from bacteria are due to its loss during evolution. It has recently been shown that HisZ from L. lactis binds to the ATP-PRPP transferase (HisG) and that both HisZ and HisG are required for catalyzing the first reaction in histidine biosynthesis. Phylogenetic analysis of HisG sequences shows conclusively that proteobacterial HisG and histidinol dehydrogenase (HisD) sequences are paraphyletic and that the partition of the Proteobacteria associated with the presence/absence of HisZ corresponds to that based on HisG and HisD paraphyly. Our results suggest that horizontal gene transfer played an important role in the evolution of the regulation of histidine biosynthesis.

Defective nucleotide excision repair in xpc mutant mice and its association with cancer predisposition

Mice that are genetically engineered are becoming increasingly more powerful tools for understanding the molecular pathology of many human hereditary diseases, especially those that confer an increased predisposition to cancer. We have generated mouse strains defective in the Xpc gene, which is required for nucleotide excision repair (NER) of DNA. Homozygous mutant mice are highly prone to skin cancer following exposure to UVB radiation, and to liver and lung cancer following exposure to the chemical carcinogen acetylaminofluorene (AAF). Skin cancer predisposition is significantly augmented when mice are additionally defective in Trp53 (p53) gene function. We also present the results of studies with mice that are heterozygous mutant in the Apex (Hap1, Ref-1) gene required for base excision repair and with mice that are defective in the mismatch repair gene Msh2. Double and triple mutant mice mutated in multiple DNA repair genes have revealed several interesting overlapping roles of DNA repair pathways in the prevention of mutation and cancer.

Genotype-specific Trp53 mutational analysis in ultraviolet B radiation-induced skin cancers in Xpc and Xpc Trp53 mutant mice

We have examined the mutational spectrum in the Trp53 gene from UVB radiation-induced skin cancers in Trp53+/+ and Trp53+/- mutant mice of all three possible Xpc genotypes. Mutations were detected in exons 2-10 of the Trp53 coding region in approximately 90% of >80 different skin cancers examined. In contrast to Trp53+/+ mice in which most mutations in the Trp53 gene were located in exons 5-8, the majority of the mutations in Trp53+/- mice were at other exons. We observed a high predilection for C-->T transition mutations at a unique CpG site in codon 122 (exon 4) of the Trp53 gene in Xpc-/- Trp53+/- mice. This site is not part of a pyrimidine dinucleotide. Mutations at this codon, as well as in codons 124 and 210, were observed exclusively in Xpc-/- or Xpc+/- mice. Mutations at the corresponding codons (127 and 213) in the human p53 gene have been reported in skin tumors from human patients with xeroderma pigmentosum. Hence, mutations at codons 122 (125), 124 (127), and 210 (213) may constitute signatures for defective or deficient nucleotide excision repair in mice (humans). In Xpc-/- mice, the majority of mutations were located at C residues in CpG sites, in which the C is presumably methylated. A similar bias can be deduced from studies in human XP individuals.

Fabs specific for 8-oxoguanine: control of DNA binding

Free radicals produce a broad spectrum of DNA base modifications including 7,8-dihydro-8-oxoguanine (8-oxoG). Since free radicals have been implicated in many pathologies and in aging, 8-oxoG has become a benchmark for factors that influence free radical production. Fab g37 is a monoclonal antibody that was isolated by phage display in an effort to create a reagent for detecting 8-oxoG in DNA. Although this antibody exhibited a high degree of specificity for the 8-oxoG base, it did not appear to recognize 8-oxoG when present in DNA. Fab g37 was modified using HCDR1 and HCDR2 segment shuffling and light chain shuffling. Fab 166 and Fab 366 which bound to 8-oxoG in single-stranded DNA were isolated. Fab 166 binds more selectively to single-stranded oligonucleotides containing 8-oxoG versus control oligonucleotides than does Fab 366 which binds DNA with reduced dependency on 8-oxoG. Numerous other clones were also isolated and characterized that contained a spectrum of specificities for 8-oxoG and for DNA. Analysis of the primary sequences of these clones and comparison with their binding properties suggested the importance of different complementarity determining regions and residues in determining the observed binding phenotypes. Subsequent chain shuffling experiments demonstrated that mutation of SerH53 to ArgH53 in the Fab g37 heavy chain slightly decreased the Fab's affinity for 8-oxoG but significantly improved its binding to DNA in an 8-oxoG-dependent manner. The light chain shuffling experiments also demonstrated that numerous promiscuous light chains could enhance DNA binding when paired with either the Fab g37 or Fab 166 heavy chains; however, only the Fab 166 light chain did so in an additive manner when combined with the Fab 166 heavy chain that contains ArgH53. A three-point model for Fab 166 binding to oligonucleotides containing 8-oxoG is proposed. We describe a successful attempt to generate a desired antibody specificity, which was not present in the animal's original immune response.

Streptococcus parasanguis pepO encodes an endopeptidase with structure and activity similar to those of enzymes that modulate peptide receptor signaling in eukaryotic cells

Studies in our laboratory have identified two fimbria-associated adhesins, FimA and Fap1, of Streptococcus parasanguis FW213. In this study, we isolated and sequenced DNA fragments linked to fimA to determine if they contained additional factors associated with adherence, virulence, or survival in the host. An open reading frame just upstream and divergently transcribed from the fimA operon was identified and named pepO. Northern hybridization indicated that pepO is transcribed as a monocistronic message. pepO encodes a predicted 631-amino-acid protein with a molecular mass of approximately 70.6 kDa. PepO contains the essential motif HEXXH, typical of many zinc-dependent metalloproteases and metallopeptidases. PepO has significant sequence identity to mammalian metallopeptidases, including endothelin-converting enzyme, which converts a potent vasoconstrictor into its active form, and neutral endopeptidase (NEP), which is involved in terminating the activity of opioid peptides. The opioid peptide metenkephalin is a natural substrate of NEP. Cell extracts of FW213 cleaved metenkephalin at the same site as does NEP, while an extract from an insertionally inactivated pepO mutant did not. These results indicate that FW213 pepO encodes an enzyme with activity similar to that of known mammalian endopeptidases. Phylogenetic analysis of PepO and its homologues suggests lateral genetic exchange between bacteria and eukaryotes.

Evolutionary conservation and somatic mutation hotspot maps of p53: correlation with p53 protein structural and functional features

Missense mutations in p53 frequently occur at 'hotspot' amino acids which are highly conserved and represent regions of structural or functional importance. Using the p53 mutation database and the p53 DNA sequences for 11 species, we more precisely defined the relationships among conservation, mutation frequency and protein structure. We aligned the p53 sequences codon-by-codon and determined the degree of substitution among them. As a whole, p53 is evolving at an average rate for a mammalian protein-coding gene. As expected, the DNA binding domain is evolving more slowly than the carboxy and amino termini. A detailed map of evolutionary conservation shows that within the DNA binding domain there are repeating peaks and valleys of higher and lower evolutionary constraint. Mutation hotspots were identified by comparing the observed distribution of mutations to the pattern expected from a random multinomial distribution. Seventy-three hotspots were identified; these 19% of codons account for 88% of all reported p53 mutations. Both high evolutionary constraint and mutation hotspots are noted at amino acids close to the protein-DNA interface and at others more distant from DNA, often buried within the core of the folded protein but sometimes on its surface. The results indicate that targeting highly conserved regions for mutational and functional analysis may be efficient strategies for the study of cancer-related genes.

Enzymatic processing of uracil glycol, a major oxidative product of DNA cytosine

A major stable oxidation product of DNA cytosine is uracil glycol (Ug). Because of the potential of Ug to be a strong premutagenic lesion, it is important to assess whether it is a blocking lesion to DNA polymerase as is its structural counterpart, thymine glycol (Tg), and to evaluate its pairing properties. Here, a series of oligonucleotides containing Ug or Tg were prepared and used as templates for a model enzyme, Escherichia coli DNA polymerase I Klenow fragment (exo-). During translesion DNA synthesis, Ug was bypassed more efficiently than Tg in all sequence contexts examined. Furthermore, only dAMP was incorporated opposite template Ug and Tg and the kinetic parameters of incorporation showed that dAMP was inserted opposite Ug more efficiently than opposite Tg. Ug opposite G and A was also recognized and removed in vitro by the E. coli DNA repair glycosylases, endonuclease III (endo III), endonuclease VIII (endo VIII), and formamidopyrimidine DNA glycosylase. The steady state kinetic parameters indicated that Ug was a better substrate for endo III and formamidopyrimidine DNA glycosylase than Tg; for endonuclease VIII, however, Tg was a better substrate.

Uracil glycol deoxynucleoside triphosphate is a better substrate for DNA polymerase I Klenow fragment than thymine glycol deoxynucleoside triphosphate

A major stable oxidation product of DNA cytosine is 5,6-dihydroxy-5, 6-dihydrouracil (Ug). Ug can be formed directly in DNA or in the cellular nucleotide pools by deamination of the unstable primary product, cytosine glycol. Here, we synthesized dUgTP and showed that dUgTP was incorporated in place of dTTP and was a much better substrate for the model enzyme DNA polymerase I Klenow fragment lacking proofreading activity, Kf (exo-), than deoxythymidine glycol triphosphate (dTgTP). The relative efficiency for dUgTP insertion opposite A was 10 times higher than for dTgTP; however, the extension of a primer with 3' dUg was about 100 times more efficient than the extension of a primer with 3' dTg. At the insertion step, the differences in Vmax appeared to be responsible since the apparent Kms for dUgTP and dTgTP were about the same. In contrast, both the apparent Km and Vmax for elongation of dUg were markedly different from those of dTg. Molecular modeling was performed with both Tg and Ug and provides a rational structural explanation for these observations.

First Report of Erwinia carotovora subsp. carotovora Causing Bacterial Root Rot of Sweetpotato (Ipomoea batatas) in Louisiana

Bacterial root and stem rot of sweetpotato (Ipomoea batatas (L.) Lam.) was first fully characterized in the U.S. in 1977 (2). It was thought to be caused exclusively by Erwinia chrysanthemi. Although a previous report described sweetpotato as a host for E. carotovora subsp. carotovora, based on artificial inoculations, others have reported that neither E. carotovora subsp. carotovora nor E. carotovora subsp. atroseptica decay sweetpotato storage roots (1). In October 1995, storage roots of sweetpotato cv. Beauregard were received from St. Landry Parish, LA, that displayed typical bacterial root rot. Isolations from these roots yielded bacteria that showed a similarity of 0.945 to E. carotovora subsp. carotovora with the Biolog GN Bacterial Identification System (version 3.50). This isolate (Ecc-LH) also differed from isolates of E. chrysanthemi (Ech) from sweetpotato and other hosts in that it was insensitive to erythromycin, did not produce phosphatase or lecithinase, and did not produce gas from glucose. Ecc-LH differed from known strains of E. carotovora subsp. atroseptica in that it did not produce reducing substances from sucrose or acid from palatinose. When Beauregard storage roots were inoculated by inserting micropipette tips containing 50 μl of 1.0 × 10⁸ CFU/ml, both Ecc-LH and Ech-48 produced typical bacterial root rot symptoms. However, when they were compared by infectivity titrations at 28 to 32°C, Ecc-LH was less virulent than Ech-48. Ecc-LH had an ED₅₀ of approximately 1.0 × 10⁶ CFU/ml and did not cause appreciable disease below inoculum concentrations of 1.0 × 10⁵, whereas Ech-48 had an ED₅₀ of approximately 1.0 × 10⁸ and caused soft rot at the lowest concentration tested, 1.0 × 10³. Similar disease incidence was observed in infectivity titrations at 22 to 24°C, but Ech-48 caused less severe soft rot. E. carotovora subsp. carotovora was reisolated from inoculated storage roots and its identity was reconfirmed by Biolog. When terminal vine cuttings of Beauregard were dipped in 1.0 × 10⁸ CFU/ml and planted in a greenhouse, bacterial stem rot symptoms developed on plants inoculated with Ech-48 at about 4 weeks postinoculation, or when new growth began. However, no symptoms developed on plants inoculated with Ecc-LH. This is the first report of natural occurrence of E. carotovora subsp. carotovora causing bacterial root rot of sweetpotato in Louisiana. E. chrysanthemi remains the most important pathogen causing bacterial soft rot in sweetpotato since it is widely associated with sweetpotato, is more virulent on storage roots and also causes a stem rot. E. carotovora subsp. carotovora can cause root rot, but has been isolated in only one location to date, is less virulent on storage roots, and apparently does not cause stem rot on the predominant cultivar in U.S. sweetpotato production, Beauregard. References: (1) C. A. Clark and J. W. Moyer. 1988. Compendium of Sweet Potato Diseases. American Phytopathological Society, St. Paul, MN. (2) N. W. Schaad and D. Brenner. Phytopathology 67:302, 1977.

Thermodynamics of the interactions of lac repressor with variants of the symmetric lac operator: effects of converting a consensus site to a non-specific site

What are the thermodynamic consequences of the stepwise conversion of a highly specific (consensus) protein-DNA interface to one that is nonspecific? How do the magnitudes of key favorable contributions to complex stability (burial of hydrophobic surfaces and reduction of DNA phosphate charge density) change as the DNA sequence of the specific site is detuned? To address these questions we investigated the binding of lac repressor (LacI) to a series of 40 bp fragments carrying symmetric (consensus) and variant operator sequences over a range of temperatures and salt concentrations. Variant DNA sites contained symmetrical single and double base-pair substitutions at positions 4 and/or 5 [sequence: see text] in each 10 bp half site of the symmetric lac operator (Osym). Non-specific interactions were examined using a 40 bp non-operator DNA fragment. Disruption of the consensus interface by a single symmetrical substitution reduces the observed equilibrium association constant (K(obs)) for Osym by three to four orders of magnitude; double symmetrical substitutions approach the six orders in magnitude difference between specific and non-specific binding to a 40 bp fragment. At these adjacent positions in the consensus site, the free energy effects of multiple substitutions are non-additive: the first reduces /deltaG(obs)o/ by 3 to 5 kcal mol(-1), approximately halfway to the non-specific level, whereas the second is less deleterious, reducing /deltaG(obs)o/ by less than 3 kcal mol(-1). Variant-specific dependences of K(obs) on temperature and salt concentration characterize these LacI-operator interactions. In general, binding constants and standard free energies of binding both exhibit characteristic extrema near 290 K. As a consequence, both the enthalpic and entropic contributions to stability of Osym and variant complexes change from positive (i.e. entropy driven) at lower temperatures to negative (i.e. enthalpy driven) at higher temperatures, indicating that the heat capacity change upon binding, deltaC(obs)o, is large and negative. In general, /deltaC(obs)o/ decreases as the specificity and stability of the variant complex decreases. Stabilities of complexes of LacI with Osym and all variant operators are strongly [salt]-dependent. Binding constants for the variant complexes exhibit a power-dependence on [salt] that is larger in magnitude (i.e. more negative) than for Osym, but no obvious trend relates changes in contributions from the polyelectrolyte effect and the observed reductions in stability (delta deltaG(obs)o). These variant-specific thermodynamic signatures provide novel insights into the consequences of converting a consensus interface to a less specific one; such insights are not obtained from comparisons at the level of delta deltaG(obs)o. We propose that this variant-specific behavior arises from a strong effect of operator sequence on the extent of induced conformational changes in the protein (and possibly also in the DNA site) which accompany binding.

Spinal cord myelin is vulnerable to decompression

Spinal cord white matter is the major site of tissue damage resulting from decompression sickness (DCS or "the bends"). Damage is thought to result from bubble nucleation within the tissue. Why DCS occurs predominantly in the spinal cord and not in the brain is not known; neither is the exact pathological mechanism by which the spinal cord is damaged, nor how multiple sclerosis (MS)-like symptoms may ensue. To investigate the molecular basis of white matter damage, we subjected myelinated mouse tissues to varying durations of decompression, and then after recompression to one atmosphere, examined them for changes in myelin structure and composition. X-ray diffraction showed that the myelin period in spinal cord decreased by 4%, whereas those of optic and sciatic nerves were stable. The change in period was accompanied by a change in membrane bilayer profile--i.e., relative to control, the width of the bilayer decreased by approximately 6 A, whereas the interbilayer spaces each increased by approximately 3 A. The changes in electron density levels suggested a redistribution of matter from the interbilayer spaces into the lipid headgroup layers. By contrast with these structural changes, sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and high-performance thin layer chromatography (HPTLC) revealed no noticeable change in myelin composition--i.e., there was no release of myelin-specific proteins or lipids. Our findings indicate that spinal cord myelin has an inherent structural vulnerability that may facilitate the targeting of this tissue during pressure changes.

Concurrent factitious disorder and factitious disorder by proxy. Double jeopardy

Two forms of medical dissimulation-factitious disorder and factitious disorder by proxy-present enormous challenges to clinicians accustomed to receiving valid symptom reports from their patients. The consequences of such "disease forgery" are heightened when a patient simultaneously engages in both forms of deception. We discuss a 34-year-old nurse who simulated or induced a panoply of physical and psychological ailments in both herself and her daughter. The staff's insistence on access to outside information sources proved indispensible in establishing both diagnoses, facilitating ongoing treatment for the patient and ensuring appropriate protection of the child.