Mechanistic insights into the pathogenesis of microtubule-targeting agent-induced peripheral neuropathy from pharmacogenetic and functional studies.

Chemotherapy-induced peripheral neuropathy (CIPN) is a common dose-limiting toxicity that affects 30%-40% of patients undergoing cancer treatment. Although multiple mechanisms of chemotherapy-induced neurotoxicity have been described in preclinical models, these have not been translated into widely effective strategies for the prevention or treatment of CIPN. Predictive biomarkers to inform therapeutic approaches are also lacking. Recent studies have examined genetic risk factors associated with CIPN susceptibility. This review provides an overview of the clinical and pathologic features of CIPN and summarizes efforts to identify target pathways through genetic and functional studies. Structurally and mechanistically diverse chemotherapeutics are associated with CIPN; however, the current review is focused on microtubule-targeting agents since these are the focus of most pharmacogenetic association and functional studies of CIPN. Genome-wide pharmacogenetic association studies are useful tools to identify not only causative genes and genetic variants but also genetic networks implicated in drug response or toxicity and have been increasingly applied to investigations of CIPN. Induced pluripotent stem cell-derived models of human sensory neurons are especially useful to understand the mechanistic significance of genomic findings. Combined genetic and functional genomic efforts to understand CIPN hold great promise for developing therapeutic approaches for its prevention and treatment.

Human Induced Pluripotent Stem Cell Derived Sensory Neurons are Sensitive to the Neurotoxic Effects of Paclitaxel.

Chemotherapy-induced peripheral neuropathy (CIPN) is a dose-limiting adverse event associated with treatment with paclitaxel and other chemotherapeutic agents. The prevention and treatment of CIPN are limited by a lack of understanding of the molecular mechanisms underlying this toxicity. In the current study, a human induced pluripotent stem cell-derived sensory neuron (iPSC-SN) model was developed for the study of chemotherapy-induced neurotoxicity. The iPSC-SNs express proteins characteristic of nociceptor, mechanoreceptor, and proprioceptor sensory neurons and show Ca2+ influx in response to capsaicin, α,ß-meATP, and glutamate. The iPSC-SNs are relatively resistant to the cytotoxic effects of paclitaxel, with half-maximal inhibitory concentration (IC50 ) values of 38.1 µM (95% confidence interval (CI) 22.9-70.9 µM) for 48-hour exposure and 9.3 µM (95% CI 5.7-16.5 µM) for 72-hour treatment. Paclitaxel causes dose-dependent and time-dependent changes in neurite network complexity detected by ßIII-tubulin staining and high content imaging. The IC50 for paclitaxel reduction of neurite area was 1.4 µM (95% CI 0.3-16.9 µM) for 48-hour exposure and 0.6 µM (95% CI 0.09-9.9 µM) for 72-hour exposure. Decreased mitochondrial membrane potential, slower movement of mitochondria down the neurites, and changes in glutamate-induced neuronal excitability were also observed with paclitaxel exposure. The iPSC-SNs were also sensitive to docetaxel, vincristine, and bortezomib. Collectively, these data support the use of iPSC-SNs for detailed mechanistic investigations of genes and pathways implicated in chemotherapy-induced neurotoxicity and the identification of novel therapeutic approaches for its prevention and treatment.

Genomewide Meta-Analysis Validates a Role for S1PR1 in Microtubule Targeting Agent-Induced Sensory Peripheral Neuropathy.

Microtubule targeting agents (MTAs) are anticancer therapies commonly prescribed for breast cancer and other solid tumors. Sensory peripheral neuropathy (PN) is the major dose-limiting toxicity for MTAs and can limit clinical efficacy. The current pharmacogenomic study aimed to identify genetic variations that explain patient susceptibility and drive mechanisms underlying development of MTA-induced PN. A meta-analysis of genomewide association studies (GWAS) from two clinical cohorts treated with MTAs (Cancer and Leukemia Group B (CALGB) 40502 and CALGB 40101) was conducted using a Cox regression model with cumulative dose to first instance of grade 2 or higher PN. Summary statistics from a GWAS of European subjects (n = 469) in CALGB 40502 that estimated cause-specific risk of PN were meta-analyzed with those from a previously published GWAS of European ancestry (n = 855) from CALGB 40101 that estimated the risk of PN. Novel single nucleotide polymorphisms in an enhancer region downstream of sphingosine-1-phosphate receptor 1 (S1PR1 encoding S1PR1 ; e.g., rs74497159, ßCALGB 40101 per allele log hazard ratio (95% confidence interval (CI)) = 0.591 (0.254-0.928), ßCALGB 40502 per allele log hazard ratio (95% CI) = 0.693 (0.334-1.053); PMETA  = 3.62 × 10-7 ) were the most highly ranked associations based on P values with risk of developing grade 2 and higher PN. In silico functional analysis identified multiple regulatory elements and potential enhancer activity for S1PR1 within this genomic region. Inhibition of S1PR1 function in induced pluripotent stem cell-derived human sensory neurons shows partial protection against paclitaxel-induced neurite damage. These pharmacogenetic findings further support ongoing clinical evaluations to target S1PR1 as a therapeutic strategy for prevention and/or treatment of MTA-induced neuropathy.

Efficient estimation of grouped survival models.

BACKGROUND: Time- and dose-to-event phenotypes used in basic science and translational studies are commonly measured imprecisely or incompletely due to limitations of the experimental design or data collection schema. For example, drug-induced toxicities are not reported by the actual time or dose triggering the event, but rather are inferred from the cycle or dose to which the event is attributed. This exemplifies a prevalent type of imprecise measurement called grouped failure time, where times or doses are restricted to discrete increments. Failure to appropriately account for the grouped nature of the data, when present, may lead to biased analyses. RESULTS: We present groupedSurv, an R package which implements a statistically rigorous and computationally efficient approach for conducting genome-wide analyses based on grouped failure time phenotypes. Our approach accommodates adjustments for baseline covariates, and analysis at the variant or gene level. We illustrate the statistical properties of the approach and computational performance of the package by simulation. We present the results of a reanalysis of a published genome-wide study to identify common germline variants associated with the risk of taxane-induced peripheral neuropathy in breast cancer patients. CONCLUSIONS: groupedSurv enables fast and rigorous genome-wide analysis on the basis of grouped failure time phenotypes at the variant, gene or pathway level. The package is freely available under a public license through the Comprehensive R Archive Network.

A Pharmacogenetic Prediction Model of Progression-Free Survival in Breast Cancer using Genome-Wide Genotyping Data from CALGB 40502 (Alliance).

Genome-wide genotyping data are increasingly available for pharmacogenetic association studies, but application of these data for development of prediction models is limited. Prediction methods, such as elastic net regularization, have recently been applied to genetic studies but only limitedly to pharmacogenetic outcomes. An elastic net was applied to a pharmacogenetic study of progression-free survival (PFS) of 468 patients with advanced breast cancer in a clinical trial of paclitaxel, nab-paclitaxel, and ixabepilone. A final model included 13 single nucleotide polymorphisms (SNPs) in addition to clinical covariates (prior taxane status, hormone receptor status, disease-free interval, and presence of visceral metastases) with an area under the curve (AUC) integrated over time of 0.81, an increase compared to an AUC of 0.64 for a model with clinical covariates alone. This model may be of value in predicting PFS with microtubule targeting agents and may inform reverse translational studies to understand differential response to these drugs.

Multiple structural states exist throughout the helical nucleation sequence of the intrinsically disordered protein stathmin, as reported by electron paramagnetic resonance spectroscopy.

The intrinsically disordered protein (IDP) stathmin plays an important regulatory role in cytoskeletal maintenance through its helical binding to tubulin and microtubules. However, it lacks a stable fold in the absence of its binding partner. Although stathmin has been a focus of research over the past two decades, the solution-phase conformational dynamics of this IDP are poorly understood. It has been reported that stathmin is purely monomeric in solution and that it bears a short helical region of persistent foldedness, which may act to nucleate helical folding in the C-terminal direction. Here we report a comprehensive study of the structural equilibria local to this region in stathmin that contradicts these two claims. Using the technique of electron paramagnetic resonance (EPR) spectroscopy on spin-labeled stathmin mutants in the solution-phase and when immobilized on Sepharose solid support, we show that all sites in the helical nucleation region of stathmin exhibit multiple spectral components that correspond to dynamic states of differing mobilities and stabilities. Importantly, a state with relatively low mobility dominates each spectrum with an average population greater than 50%, which we suggest corresponds to an oligomerized state of the protein. This is in contrast to a less populated, more mobile state, which likely represents a helically folded monomeric state of stathmin, and a highly mobile state, which we propose is the random coil conformer of the protein. Our interpretation of the EPR data is confirmed by further characterization of the protein using the techniques of native and SDS PAGE, gel filtration chromatography, and multiangle and dynamic light scattering, all of which show the presence of oligomeric stathmin in solution. Collectively, these data suggest that stathmin exists in a diverse equilibrium of states throughout the purported helical nucleation region and that this IDP exhibits a propensity toward oligomerization.