Immuno-transcriptomic profiling of extracranial pediatric solid malignancies.

We perform an immunogenomics analysis utilizing whole-transcriptome sequencing of 657 pediatric extracranial solid cancer samples representing 14 diagnoses, and additionally utilize transcriptomes of 131 pediatric cancer cell lines and 147 normal tissue samples for comparison. We describe patterns of infiltrating immune cells, T cell receptor (TCR) clonal expansion, and translationally relevant immune checkpoints. We find that tumor-infiltrating lymphocytes and TCR counts vary widely across cancer types and within each diagnosis, and notably are significantly predictive of survival in osteosarcoma patients. We identify potential cancer-specific immunotherapeutic targets for adoptive cell therapies including cell-surface proteins, tumor germline antigens, and lineage-specific transcription factors. Using an orthogonal immunopeptidomics approach, we find several potential immunotherapeutic targets in osteosarcoma and Ewing sarcoma and validated PRAME as a bona fide multi-pediatric cancer target. Importantly, this work provides a critical framework for immune targeting of extracranial solid tumors using parallel immuno-transcriptomic and -peptidomic approaches.

The somatic autosomal mutation matrix in cancer genomes.

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

Can structural features of kinase receptors provide clues on selectivity and inhibition? A molecular modeling study.

Cancer is a complex disease resulting from the uncontrolled proliferation of cell signaling events. Protein kinases have been identified as central molecules that participate overwhelmingly in oncogenic events, thus becoming key targets for anticancer drugs. A majority of studies converged on the idea that ligand-binding pockets of kinases retain clues to the inhibiting abilities and cross-reacting tendencies of inhibitor drugs. Even though these ideas are critical for drug discovery, validating them using experiments is not only difficult, but also in some cases infeasible. To overcome these limitations and to test these ideas at the molecular level, we present here the results of receptor-focused in-silico docking of nine marketed drugs to 19 different wild-type and mutated kinases chosen from a wide range of families. This investigation highlights the need for using relevant models to explain the correct inhibition trends and the results are used to make predictions that might be able to influence future experiments. Our simulation studies are able to correctly predict the primary targets for each drug studied in majority of cases and our results agree with the existing findings. Our study shows that the conformations a given receptor acquires during kinase activation, and their micro-environment, defines the ligand partners. Type II drugs display high compatibility and selectivity for DFG-out kinase conformations. On the other hand Type I drugs are less selective and show binding preferences for both the open and closed forms of selected kinases. Using this receptor-focused approach, it is possible to capture the observed fold change in binding affinities between the wild-type and disease-centric mutations in ABL kinase for Imatinib and the second-generation ABL drugs. The effects of mutation are also investigated for two other systems, EGFR and B-Raf. Finally, by including pathway information in the design it is possible to model kinase inhibitors with potentially fewer side-effects.

Guanine holes are prominent targets for mutation in cancer and inherited disease.

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

The role of methylation in the intrinsic dynamics of B- and Z-DNA.

Methylation of cytosine at the 5-carbon position (5 mC) is observed in both prokaryotes and eukaryotes. In humans, DNA methylation at CpG sites plays an important role in gene regulation and has been implicated in development, gene silencing, and cancer. In addition, the CpG dinucleotide is a known hot spot for pathologic mutations genome-wide. CpG tracts may adopt left-handed Z-DNA conformations, which have also been implicated in gene regulation and genomic instability. Methylation facilitates this B-Z transition but the underlying mechanism remains unclear. Herein, four structural models of the dinucleotide d(GC)(5) repeat sequence in B-, methylated B-, Z-, and methylated Z-DNA forms were constructed and an aggregate 100 nanoseconds of molecular dynamics simulations in explicit solvent under physiological conditions was performed for each model. Both unmethylated and methylated B-DNA were found to be more flexible than Z-DNA. However, methylation significantly destabilized the BII, relative to the BI, state through the Gp5mC steps. In addition, methylation decreased the free energy difference between B- and Z-DNA. Comparisons of α/γ backbone torsional angles showed that torsional states changed marginally upon methylation for B-DNA, and Z-DNA. Methylation-induced conformational changes and lower energy differences may contribute to the transition to Z-DNA by methylated, over unmethylated, B-DNA and may be a contributing factor to biological function.

Non-B DNA-forming sequences and WRN deficiency independently increase the frequency of base substitution in human cells.

Although alternative DNA secondary structures (non-B DNA) can induce genomic rearrangements, their associated mutational spectra remain largely unknown. The helicase activity of WRN, which is absent in the human progeroid Werner syndrome, is thought to counteract this genomic instability. We determined non-B DNA-induced mutation frequencies and spectra in human U2OS osteosarcoma cells and assessed the role of WRN in isogenic knockdown (WRN-KD) cells using a supF gene mutation reporter system flanked by triplex- or Z-DNA-forming sequences. Although both non-B DNA and WRN-KD served to increase the mutation frequency, the increase afforded by WRN-KD was independent of DNA structure despite the fact that purified WRN helicase was found to resolve these structures in vitro. In U2OS cells, ∼70% of mutations comprised single-base substitutions, mostly at G·C base-pairs, with the remaining ∼30% being microdeletions. The number of mutations at G·C base-pairs in the context of NGNN/NNCN sequences correlated well with predicted free energies of base stacking and ionization potentials, suggesting a possible origin via oxidation reactions involving electron loss and subsequent electron transfer (hole migration) between neighboring bases. A set of ∼40,000 somatic mutations at G·C base pairs identified in a lung cancer genome exhibited similar correlations, implying that hole migration may also be involved. We conclude that alternative DNA conformations, WRN deficiency and lung tumorigenesis may all serve to increase the mutation rate by promoting, through diverse pathways, oxidation reactions that perturb the electron orbitals of neighboring bases. It follows that such "hole migration" is likely to play a much more widespread role in mutagenesis than previously anticipated.

Abundance and length of simple repeats in vertebrate genomes are determined by their structural properties.

Microsatellites are abundant in vertebrate genomes, but their sequence representation and length distributions vary greatly within each family of repeats (e.g., tetranucleotides). Biophysical studies of 82 synthetic single-stranded oligonucleotides comprising all tetra- and trinucleotide repeats revealed an inverse correlation between the stability of folded-back hairpin and quadruplex structures and the sequence representation for repeats > or =30 bp in length in nine vertebrate genomes. Alternatively, the predicted energies of base-stacking interactions correlated directly with the longest length distributions in vertebrate genomes. Genome-wide analyses indicated that unstable sequences, such as CAG:CTG and CCG:CGG, were over-represented in coding regions and that micro/minisatellites were recruited in genes involved in transcription and signaling pathways, particularly in the nervous system. Microsatellite instability (MSI) is a hallmark of cancer, and length polymorphism within genes can confer susceptibility to inherited disease. Sequences that manifest the highest MSI values also displayed the strongest base-stacking interactions; analyses of 62 tri- and tetranucleotide repeat-containing genes associated with human genetic disease revealed enrichments similar to those noted for micro/minisatellite-containing genes. We conclude that DNA structure and base-stacking determined the number and length distributions of microsatellite repeats in vertebrate genomes over evolutionary time and that micro/minisatellites have been recruited to participate in both gene and protein function.

Proteomic analysis identifies oxidative stress induction by adaphostin.

PURPOSE: Activities distinct from inhibition of Bcr/abl have led to adaphostin (NSC 680410) being described as "a drug in search of a mechanism." In this study, proteomic analysis of adaphostin-treated myeloid leukemia cell lines was used to further elucidate a mechanism of action. EXPERIMENTAL DESIGN: HL60 and K562 cells treated with adaphostin for 6, 12, or 24 h were analyzed using two-dimensional PAGE. Differentially expressed spots were excised, digested with trypsin, and analyzed by liquid chromatography-tandem mass spectrometry. The contribution of the redox-active hydroquinone group in adaphostin was also examined by carrying out proteomic analysis of HL60 cells treated with a simple hydroquinone (1,4-dihydroxybenzene) or H(2)O(2). RESULTS: Analysis of adaphostin-treated cells identified 49 differentially expressed proteins, the majority being implicated in the response to oxidative stress (e.g., CALM, ERP29, GSTP1, PDIA1) or induction of apoptosis (e.g., LAMA, FLNA, TPR, GDIS). Interestingly, modulation of these proteins was almost fully prevented by inclusion of an antioxidant, N-acetylcysteine. Validation of the proteomic data confirmed GSTP1 as an adaphostin resistance gene. Subsequent analysis of HL60 cells treated with 1,4-dihydroxybenzene or H(2)O(2) showed similar increases in intracellular peroxides and an almost identical proteomic profiles to that of adaphostin treatment. Western blotting of a panel of cell lines identified Cu/Zn superoxide dismutase (SOD) as correlating with adaphostin resistance. The role of SOD as a second adaphostin resistance gene was confirmed by demonstrating that inhibition of SOD using diethyldithiocarbamate increased adaphostin sensitivity, whereas transfection of SOD I attenuated toxicity. Importantly, treatment with 1,4-dihydroxybenzene or H(2)O(2) replicated adaphostin-induced Bcr/abl polypeptide degradation, suggesting that kinase inhibition is a ROS-dependent phenomenon. CONCLUSION: Adaphostin should be classified as a redox-active-substituted dihydroquinone.

Interaction of noncompetitive inhibitors with an immobilized alpha3beta4 nicotinic acetylcholine receptor investigated by affinity chromatography, quantitative-structure activity relationship analysis, and molecular docking.

A large number of drug substances act as noncompetitive inhibitors (NCIs) of the nicotinic acetylcholine receptor (nAChR) by blocking the ion flux through the channel. An affinity chromatography technique has been developed for investigating the interactions between NCIs and the alpha3beta4 subtype of neuronal nAChR. The data obtained from the chromatographic study were used to construct QSAR models of the NCI-nAChR binding with both electronic and steric parameters observed as important descriptors. A molecular model of the transmembrane domain of the alpha3beta4 subtype of nAChR was constructed and used to simulate the docking of a series of NCIs. A key aspect of the model was the discovery of the cleft produced by the incorporation of the bulky phenylalanine moiety into the nonpolar section of the lumen by the beta4 subunit. Quantitatively, the results of docking simulations modeled the experimental affinity data better than QSAR results. The computational approach, combined with the modeling of NCI-nAChR interaction by affinity chromatography, can be used to predict possible toxicities and adverse interactions.

Theoretical studies on the hydrolysis of mono-phosphate and tri-phosphate in gas phase and aqueous solution.

Phosphate hydrolysis by GTPases plays an important role as a molecular switch in signal transduction and as an initiator of many other biological processes. Despite the centrality of this ubiquitous reaction, the mechanism is still poorly understood. As a first step to understand the mechanisms of this process, the nonenzymatic hydrolysis of mono-phosphate and tri-phosphate esters were systematically studied in gas phase and aqueous solution using hybrid density functional methods. The dielectric effect of the environment on the energetics of these processes was also explored. Theoretical results show that for mono-phosphate ester, the dissociative pathway is much more favorable than the associative pathway. However, the reaction barriers for the dissociative and associative pathways of tri-phosphate hydrolysis are very close in aqueous solution, though the dissociative pathway is more favorable in the gas phase. High dielectric solvents, such as water, significantly lower the activation barrier of the associative pathway due to the greater solvation energy of the associative transition states than that of the reactant complex. By contrast, the barrier of the dissociative pathway, with respect to the gas phase, is less sensitive to the surrounding dielectric. In the associative hydrolysis pathway of the tri-phosphate ester, negative charge is transferred from the gamma-phosphate to beta-phosphate through the bridging ester oxygen and results in Pgamma-O bond dissociation. No analogous charge transfer was observed in the dissociative pathway, where Pgamma-O bond dissociation resulted from proton transfer from the gamma-phosphate to the bridge oxygen. Finally, the active participation of local water molecules can significantly lower the activation energy of the dissociative pathway for both mono-phosphate and tri-phosphate.

Human topoisomerase I inhibition: docking camptothecin and derivatives into a structure-based active site model.

Human topoisomerase I (top1) is an important target for anti-cancer drugs, which include camptothecin (CPT) and its derivatives. To elucidate top1 inhibition in vitro, we made a series of duplex DNA substrates containing a deoxyadenosine stereospecifically modified by a covalent adduct of benzo[a]pyrene (BaP) diol epoxide [Pommier, Y., et al. (2000) Proc. Natl. Acad. Sci. U.S.A. 97, 10739-10744]. The known orientation of the hydrocarbon adduct in the DNA duplex relative to the top1 cleavage site, in combination with a top1/DNA crystal structure [Redinbo, M. R., et al. (1998) Science 279, 1504-1513], was used to construct a structure-based model to explain the in vitro top1 inhibition results obtained with adducted DNA duplexes. Here we experimentally determined that the lactone form of CPT was stabilized by an irreversible top1/DNA covalent complex. We removed the BaP moiety from the DNA in the published model, and docked the lactone forms of CPT and derivatives into the top1/DNA active site cavity. The docked ligands were minimized, and interaction energy scores between the ligands and the top1/DNA complex were determined. CPT docks perpendicular to the DNA backbone, projects outward from the major groove, and makes a network of potential H-bonds with the active site DNA and top1 residues, including Arg364, Lys532, and Asn722. The results are consistent with the known structure-activity relationships of CPT and derivatives. In addition, the model proposed a novel top1/N352A "resistance" mutation for 10-OH derivatives of CPT. The in vitro biochemical characterization of the top1/N352A mutant supported the model.