Proteogenomic characterization of small cell lung cancer identifies biological insights and subtype-specific therapeutic strategies.

We performed comprehensive proteogenomic characterization of small cell lung cancer (SCLC) using paired tumors and adjacent lung tissues from 112 treatment-naive patients who underwent surgical resection. Integrated multi-omics analysis illustrated cancer biology downstream of genetic aberrations and highlighted oncogenic roles of FAT1 mutation, RB1 deletion, and chromosome 5q loss. Two prognostic biomarkers, HMGB3 and CASP10, were identified. Overexpression of HMGB3 promoted SCLC cell migration via transcriptional regulation of cell junction-related genes. Immune landscape characterization revealed an association between ZFHX3 mutation and high immune infiltration and underscored a potential immunosuppressive role of elevated DNA damage response activity via inhibition of the cGAS-STING pathway. Multi-omics clustering identified four subtypes with subtype-specific therapeutic vulnerabilities. Cell line and patient-derived xenograft-based drug tests validated the specific therapeutic responses predicted by multi-omics subtyping. This study provides a valuable resource as well as insights to better understand SCLC biology and improve clinical practice.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Proteogenomic characterization of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

Proteogenomics of Non-smoking Lung Cancer in East Asia Delineates Molecular Signatures of Pathogenesis and Progression.

Lung cancer in East Asia is characterized by a high percentage of never-smokers, early onset and predominant EGFR mutations. To illuminate the molecular phenotype of this demographically distinct disease, we performed a deep comprehensive proteogenomic study on a prospectively collected cohort in Taiwan, representing early stage, predominantly female, non-smoking lung adenocarcinoma. Integrated genomic, proteomic, and phosphoproteomic analysis delineated the demographically distinct molecular attributes and hallmarks of tumor progression. Mutational signature analysis revealed age- and gender-related mutagenesis mechanisms, characterized by high prevalence of APOBEC mutational signature in younger females and over-representation of environmental carcinogen-like mutational signatures in older females. A proteomics-informed classification distinguished the clinical characteristics of early stage patients with EGFR mutations. Furthermore, integrated protein network analysis revealed the cellular remodeling underpinning clinical trajectories and nominated candidate biomarkers for patient stratification and therapeutic intervention. This multi-omic molecular architecture may help develop strategies for management of early stage never-smoker lung adenocarcinoma.

Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy.

The integration of mass spectrometry-based proteomics with next-generation DNA and RNA sequencing profiles tumors more comprehensively. Here this "proteogenomics" approach was applied to 122 treatment-naive primary breast cancers accrued to preserve post-translational modifications, including protein phosphorylation and acetylation. Proteogenomics challenged standard breast cancer diagnoses, provided detailed analysis of the ERBB2 amplicon, defined tumor subsets that could benefit from immune checkpoint therapy, and allowed more accurate assessment of Rb status for prediction of CDK4/6 inhibitor responsiveness. Phosphoproteomics profiles uncovered novel associations between tumor suppressor loss and targetable kinases. Acetylproteome analysis highlighted acetylation on key nuclear proteins involved in the DNA damage response and revealed cross-talk between cytoplasmic and mitochondrial acetylation and metabolism. Our results underscore the potential of proteogenomics for clinical investigation of breast cancer through more accurate annotation of targetable pathways and biological features of this remarkably heterogeneous malignancy.

Integrated Proteogenomic Characterization of HBV-Related Hepatocellular Carcinoma.

We performed the first proteogenomic characterization of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) using paired tumor and adjacent liver tissues from 159 patients. Integrated proteogenomic analyses revealed consistency and discordance among multi-omics, activation status of key signaling pathways, and liver-specific metabolic reprogramming in HBV-related HCC. Proteomic profiling identified three subgroups associated with clinical and molecular attributes including patient survival, tumor thrombus, genetic profile, and the liver-specific proteome. These proteomic subgroups have distinct features in metabolic reprogramming, microenvironment dysregulation, cell proliferation, and potential therapeutics. Two prognostic biomarkers, PYCR2 and ADH1A, related to proteomic subgrouping and involved in HCC metabolic reprogramming, were identified. CTNNB1 and TP53 mutation-associated signaling and metabolic profiles were revealed, among which mutated CTNNB1-associated ALDOA phosphorylation was validated to promote glycolysis and cell proliferation. Our study provides a valuable resource that significantly expands the knowledge of HBV-related HCC and may eventually benefit clinical practice.

Frozen tissue coring and layered histological analysis improves cell type-specific proteogenomic characterization of pancreatic adenocarcinoma.

BACKGROUND: Omics characterization of pancreatic adenocarcinoma tissue is complicated by the highly heterogeneous and mixed populations of cells. We evaluate the feasibility and potential benefit of using a coring method to enrich specific regions from bulk tissue and then perform proteogenomic analyses. METHODS: We used the Biopsy Trifecta Extraction (BioTExt) technique to isolate cores of epithelial-enriched and stroma-enriched tissue from pancreatic tumor and adjacent tissue blocks. Histology was assessed at multiple depths throughout each core. DNA sequencing, RNA sequencing, and proteomics were performed on the cored and bulk tissue samples. Supervised and unsupervised analyses were performed based on integrated molecular and histology data. RESULTS: Tissue cores had mixed cell composition at varying depths throughout. Average cell type percentages assessed by histology throughout the core were better associated with KRAS variant allele frequencies than standard histology assessment of the cut surface. Clustering based on serial histology data separated the cores into three groups with enrichment of neoplastic epithelium, stroma, and acinar cells, respectively. Using this classification, tumor overexpressed proteins identified in bulk tissue analysis were assigned into epithelial- or stroma-specific categories, which revealed novel epithelial-specific tumor overexpressed proteins. CONCLUSIONS: Our study demonstrates the feasibility of multi-omics data generation from tissue cores, the necessity of interval H&E stains in serial histology sections, and the utility of coring to improve analysis over bulk tissue data.

Identification of a specific exporter that enables high production of aconitic acid in Aspergillus pseudoterreus.

Aconitic acid is an unsaturated tricarboxylic acid that is attractive for its potential use in manufacturing biodegradable and biocompatible polymers, plasticizers, and surfactants. Previously Aspergillus pseudoterreus was engineered as a platform to produce aconitic acid by deleting the cadA (cis-aconitic acid decarboxylase) gene in the itaconic acid biosynthetic pathway. In this study, the aconitic acid transporter gene (aexA) was identified using comparative global discovery proteomics analysis between the wild-type and cadA deletion strains. The protein AexA belongs to the Major Facilitator Superfamily (MFS). Deletion of aexA almost abolished aconitic acid secretion, while its overexpression led to a significant increase in aconitic acid production. Transportation of aconitic acid across the plasma membrane is a key limiting step in its production. In vitro, proteoliposome transport assay further validated AexA's function and substrate specificity. This research provides new approaches to efficiently pinpoint and characterize exporters of fungal organic acids and accelerate metabolic engineering to improve secretion capability and lower the cost of bioproduction.

Exercise therapy improves pain and mouth opening in temporomandibular disorders: A systematic review with meta-analysis.

OBJECTIVE: To analyse the effectiveness of exercise therapy in improving pain and active or passive maximum mouth opening in patients with temporomandibular disorders. DATA SOURCES: PubMed Medline, Web of Science, Scopus, CINAHL Complete and Physiotherapy Evidence Database, until April 2022, in accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. REVIEW METHODS: We included randomized controlled trials evaluating the effect of exercise therapy on pain and on active and passive maximum mouth opening in patients with temporomandibular disorders. Effect size was calculated using Cohen's standardized mean difference (SMD) and their 95% confidence interval (95% CI) in a random-effects model. RESULTS: A total of 16 studies with 812 participants were included. Exercise therapy is effective in reducing pain (SMD: -0.58; 95% CI: -1.01 to -0.12) and increasing the pain pressure threshold (SMD: 0.45; 95% CI: 0.14-0.76), active and passive maximum mouth opening (SMD: 0.43; 95% CI: 0.14-0.71 and SMD: 0.4; 95% CI: 0.06-0.75, respectively). Subgroup analyses showed more effect of exercise therapy more splints versus splints on pain (SMD: -0.5; 95% CI: -0.73 to -0.26), active and passive maximum mouth opening (SMD: 1.14; 95% CI: 0.22-2.07 and SMD: 0.56; 95% CI: 0.06-1.06, respectively). On pain pressure threshold, exercise therapy was better than physiotherapy approach (manual therapy and electrotherapy) (SMD: 0.48; 95% CI: 0.09-0.87). CONCLUSIONS: Therapeutic exercise is an effective therapy to reduce pain and increase pain pressure threshold and active and passive maximum mouth opening in patients with temporomandibular disorders.

Neoplastic cell enrichment of tumor tissues using coring and laser microdissection for proteomic and genomic analyses of pancreatic ductal adenocarcinoma.

BACKGROUND: The identification of differentially expressed tumor-associated proteins and genomic alterations driving neoplasia is critical in the development of clinical assays to detect cancers and forms the foundation for understanding cancer biology. One of the challenges in the analysis of pancreatic ductal adenocarcinoma (PDAC) is the low neoplastic cellularity and heterogeneous composition of bulk tumors. To enrich neoplastic cells from bulk tumor tissue, coring, and laser microdissection (LMD) sampling techniques have been employed. In this study, we assessed the protein and KRAS mutation changes associated with samples obtained by these enrichment techniques and evaluated the fraction of neoplastic cells in PDAC for proteomic and genomic analyses. METHODS: Three fresh frozen PDAC tumors and their tumor-matched normal adjacent tissues (NATs) were obtained from three sampling techniques using bulk, coring, and LMD; and analyzed by TMT-based quantitative proteomics. The protein profiles and characterizations of differentially expressed proteins in three sampling groups were determined. These three PDACs and samples of five additional PDACs obtained by the same three sampling techniques were also subjected to genomic analysis to characterize KRAS mutations. RESULTS: The neoplastic cellularity of eight PDACs ranged from less than 10% to over 80% based on morphological review. Distinctive proteomic patterns and abundances of certain tumor-associated proteins were revealed when comparing the tumors and NATs by different sampling techniques. Coring and bulk tissues had comparable proteome profiles, while LMD samples had the most distinct proteome composition compared to bulk tissues. Further genomic analysis of bulk, cored, or LMD samples demonstrated that KRAS mutations were significantly enriched in LMD samples while coring was less effective in enriching for KRAS mutations when bulk tissues contained a relatively low neoplastic cellularity. CONCLUSIONS: In addition to bulk tissues, samples from LMD and coring techniques can be used for proteogenomic studies. The greatest enrichment of neoplastic cellularity is obtained with the LMD technique.

The association of salivary caffeine levels with serum concentrations in premature infants with apnea of prematurity.

The purpose of this paper is to verify whether the concentrations of caffeine in saliva are comparable to serum concentrations in preterm infants who are treated with caffeine for apnea of prematurity. This is a prospective observational study. Eligible participants were newborn infants < 37 weeks of gestational age treated with oral or intravenous caffeine for apnea of prematurity. Two paired samples of saliva and blood were collected per patient. Tube solid-phase microextraction coupled online to capillary liquid chromatography with diode array detection was used for analysis. A total of 47 infants with a median gestational age of 28 [26-30] weeks and a mean of 1.11 ± 0.4 kg of birth weight. Median postmenstrual age, when samples were collected, was 31 [29-33] weeks. Serum caffeine median levels of 19.30 µg/mL [1.9-53.90] and salivary caffeine median levels of 16.36 µg/mL [2.20-56.90] were obtained. There was a strong positive Pearson's correlation between the two variables r = 0.83 (p < 0.001). CONCLUSION: The measurement of salivary caffeine concentrations after intravenous or oral administration offers an alternative to serum caffeine monitoring in apnea of prematurity. Measurement of salivary concentration minimizes blood draws, improves blood conservation, and subsequently minimizes painful procedures in premature infants. WHAT IS KNOWN: • Salivary sampling may be useful when is applied to extremely low birth weight infant, in whom blood sampling must be severely restricted. WHAT IS NEW: • The measurement of caffeine salivary concentrations after intravenous or oral administration offers an alternative to serum caffeine monitoring in apnoea of prematurity. • Salivary sampling may be a valid non-invasive alternative that could be used to individualize and optimize caffeine dose.

Clinical Manifestations, Mutational Analysis, and Immunological Phenotype in Patients with RAG1/2 Mutations: First Cases Series from Mexico and Description of Two Novel Mutations.

Mutations in recombinase activating genes 1 and 2 (RAG1/2) result in human severe combined immunodeficiency (SCID). The products of these genes are essential for V(D)J rearrangement of the antigen receptors during lymphocyte development. Mutations resulting in null-recombination activity in RAG1 or RAG2 are associated with the most severe clinical and immunological phenotypes, whereas patients with hypomorphic mutations may develop leaky SCID, including Omenn syndrome (OS). A group of previously unrecognized clinical phenotypes associated with granulomata and/or autoimmunity have been described as a consequence of hypomorphic mutations. Here, we present six patients from unrelated families with missense variants in RAG1 or RAG2. Phenotypes observed in these patients ranged from OS to severe mycobacterial infections and granulomatous disease. Moreover, we report the first evidence of two variants that had not been associated with immunodeficiency. This study represents the first case series of RAG1- or RAG2-deficient patients from Mexico and Latin America.

Extensive pollen-mediated gene flow across intensively managed landscapes in an insect-pollinated shrub native to semiarid habitats.

Our knowledge of the impact of landscape fragmentation on gene flow patterns is mainly drawn from tropical and temperate ecosystems, where landscape features, such as the distance of a tree to the forest edge, drive connectivity and mating patterns. Yet, the structure of arid and semiarid plant communities - with open canopies and a scattered distribution of trees - differs greatly from those that are well-characterized in the literature. As a result, we ignore whether the documented consequences of landscape fragmentation on plant mating and gene flow patterns also hold for native plant communities in arid and semiarid regions. We investigated the relative contribution of plant traits, pollinator activity, and individual neighbourhood in explaining variation in mating and gene flow patterns of an insect-pollinated semiarid arborescent shrub, Ziziphus lotus, at three sites embedded in highly altered agriculture landscapes. We used 14 SSRs, seed paternity analyses, and individual mixed effect mating models (MEMMi) to estimate the individual mating variables and the pollen dispersal kernel at each site. Individual spatial location, flower density, and floral visitation rate explained most of the variation of mating variables. Unexpectedly, individual correlated paternity was very low and shrubs surrounded by the most degraded matrix exhibited an increased fraction of pollen immigration and a high effective number of pollen donors per mother shrub. Overall, our results reveal that an active pollinator assemblage ensures highly efficient mating, and maintains pollen-mediated gene flow and notable connectivity levels, even in highly altered landscapes, potentially halting genetic isolation within and between distant sites.

EML4-ALK induces cellular senescence in mortal normal human cells and promotes anchorage-independent growth in hTERT-transduced normal human cells.

BACKGROUND: Chromosomal inversions involving anaplastic lymphoma kinase (ALK) and echinoderm microtubule associated protein like 4 (EML4) generate a fusion protein EML4-ALK in non-small cell lung cancer (NSCLC). The understanding of EML4-ALK function can be improved by a functional study using normal human cells. METHODS: Here we for the first time conduct such study to examine the effects of EML4-ALK on cell proliferation, cellular senescence, DNA damage, gene expression profiles and transformed phenotypes. RESULTS: The lentiviral expression of EML4-ALK in mortal, normal human fibroblasts caused, through its constitutive ALK kinase activity, an early induction of cellular senescence with accumulated DNA damage, upregulation of p16INK4A and p21WAF1, and senescence-associated ß-galactosidase (SA-ß-gal) activity. In contrast, when EML4-ALK was expressed in normal human fibroblasts transduced with telomerase reverse transcriptase (hTERT), which is activated in the vast majority of NSCLC, the cells showed accelerated proliferation and acquired anchorage-independent growth ability in soft-agar medium, without accumulated DNA damage, chromosome aberration, nor p53 mutation. EML4-ALK induced the phosphorylation of STAT3 in both mortal and hTERT-transduced cells, but RNA sequencing analysis suggested that the different signaling pathways contributed to the different phenotypic outcomes in these cells. While EML4-ALK also induced anchorage-independent growth in hTERT-immortalized human bronchial epithelial cells in vitro, the expression of EML4-ALK alone did not cause detectable in vivo tumorigenicity in immunodeficient mice. CONCLUSIONS: Our data indicate that the expression of hTERT is critical for EML4-ALK to manifest its in vitro transforming activity in human cells. This study provides the isogenic pairs of human cells with and without EML4-ALK expression.

Genetic, Immunological, and Clinical Features of the First Mexican Cohort of Patients with Chronic Granulomatous Disease.

PURPOSE: Chronic granulomatous disease (CGD) is a primary immunodeficiency characterized by an inability of phagocytes to produce reactive oxygen species, impairing their killing of various bacteria and fungi. We summarize here the 93 cases of CGD diagnosed in Mexico from 2011 to 2019. METHODS: Thirteen Mexican hospitals participated in this study. We describe the genetic, immunological, and clinical features of the 93 CGD patients from 78 unrelated kindreds. RESULTS: Eighty-two of the patients (88%) were male. All patients developed bacterial infections and 30% suffered from some kind of fungal infection. Fifty-four BCG-vaccinated patients (58%) presented infectious complications of BCG vaccine. Tuberculosis occurred in 29%. Granulomas were found in 56% of the patients. Autoimmune and inflammatory diseases were present in 15% of patients. A biological diagnosis of CGD was made in 89/93 patients, on the basis of NBT assay (n = 6), DHR (n = 27), and NBT plus DHR (n = 56). The deficiency was complete in all patients. The median age of biological diagnosis was 17 months (range, 0-186 months). A genetic diagnosis was made in 83/93 patients (when material was available), corresponding to CYBB (n = 64), NCF1 (n = 7), NCF2 (n = 7), and CYBA (n = 5) mutations. CONCLUSIONS: The clinical manifestations in these Mexican CGD patients were similar to those in patients elsewhere. This cohort is the largest in Latin America. Mycobacterial infections are an important cause of morbidity in Mexico, as in other countries in which tuberculosis is endemic and infants are vaccinated with BCG. X-linked CGD accounted for most of the cases in Mexico, as in other Latin American countries. However, a significant number of CYBA and NCF2 mutations were identified, expanding the spectrum of known causal mutations.

Abatacept in interstitial lung disease associated with rheumatoid arthritis: national multicenter study of 263 patients.

OBJECTIVE: To assess the efficacy of abatacept (ABA) in RA patients with interstitial lung disease (ILD) (RA-ILD). METHODS: This was an observational, multicentre study of RA-ILD patients treated with at least one dose of ABA. ILD was diagnosed by high-resolution CT (HRCT). We analysed the following variables at baseline (ABA initiation), 12 months and at the end of the follow-up: Modified Medical Research Council (MMRC) scale (1-point change), forced vital capacity (FVC) or diffusion lung capacity for carbon monoxide (DLCO) (improvement or worsening ≥10%), HRCT, DAS on 28 joints evaluated using the ESR (DAS28ESR) and CS-sparing effect. RESULTS: We studied 263 RA-ILD patients [150 women/113 men; mean (s.d.) age 64.6 (10) years]. At baseline, they had a median duration of ILD of 1 (interquartile range 0.25-3.44) years, moderate or severe degree of dyspnoea (MMRC grade 2, 3 or 4) (40.3%), FVC (% of the predicted) mean (s.d.) 85.9 (21.8)%, DLCO (% of the predicted) 65.7 (18.3) and DAS28ESR 4.5 (1.5). The ILD patterns were: usual interstitial pneumonia (UIP) (40.3%), non-specific interstitial pneumonia (NSIP) (31.9%) and others (27.8%). ABA was prescribed at standard dose, i.v. (25.5%) or s.c. (74.5%). After a median follow-up of 12 (6-36) months the following variables did not show worsening: dyspnoea (MMRC) (91.9%); FVC (87.7%); DLCO (90.6%); and chest HRCT (76.6%). A significant improvement of DAS28ESR from 4.5 (1.5) to 3.1 (1.3) at the end of follow-up (P < 0.001) and a CS-sparing effect from a median 7.5 (5-10) to 5 (2.5-7.5) mg/day at the end of follow-up (P < 0.001) was also observed. ABA was withdrawn in 62 (23.6%) patients due to adverse events (n = 30), articular inefficacy (n = 27), ILD worsening (n = 3) and other causes (n = 2). CONCLUSION: ABA may be an effective and safe treatment for patients with RA-ILD.

Topiramate pharmacokinetics in neonates undergoing therapeutic hypothermia and proposal of an optimised dosing schedule.

AIM: The adequate dosing of topiramate in neonates undergoing therapeutic hypothermia has not been established. The aim of this study was to design a dosing schedule capable of providing topiramate serum concentrations within the accepted therapeutic range. METHODS: Neonates (n = 52) with hypoxic ischaemic encephalopathy and subjected to therapeutic hypothermia were dosed with topiramate, 5 mg/kg on day one and 3 mg/kg on days two to five, to decrease seizure events. A total of 451 topiramate serum concentrations obtained in the patients were used to develop a population pharmacokinetic model using a non-linear mixed-effects modelling approach. RESULTS: A one-compartment model with first-order absorption and two different clearance terms, one for the cooling period and another for the post-warming period, were used to describe the concentration-time topiramate data. The probability of no-seizure events could not be related to topiramate concentrations, which was attributed to excessively low topiramate concentrations. A modified dosage schedule was designed with the aim of obtaining more than 90% of patients with topiramate concentrations within the therapeutic range after the first dose. CONCLUSION: The dosage schedule of topiramate in these patients should be modified with the aim of decreasing the frequency of seizure events.