Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease.

The genetics of complex disease produce alterations in the molecular interactions of cellular pathways whose collective effect may become clear through the organized structure of molecular networks. To characterize molecular systems associated with late-onset Alzheimer's disease (LOAD), we constructed gene-regulatory networks in 1,647 postmortem brain tissues from LOAD patients and nondemented subjects, and we demonstrate that LOAD reconfigures specific portions of the molecular interaction structure. Through an integrative network-based approach, we rank-ordered these network structures for relevance to LOAD pathology, highlighting an immune- and microglia-specific module that is dominated by genes involved in pathogen phagocytosis, contains TYROBP as a key regulator, and is upregulated in LOAD. Mouse microglia cells overexpressing intact or truncated TYROBP revealed expression changes that significantly overlapped the human brain TYROBP network. Thus the causal network structure is a useful predictor of response to gene perturbations and presents a framework to test models of disease mechanisms underlying LOAD.

A mechanistic biomarker investigation of fialuridine hepatotoxicity using the chimeric TK-NOG Hu-liver mouse model and in vitro micropatterned hepatocyte cocultures.

Fialuridine (FIAU) is a nucleoside-based drug that caused liver failure and deaths in a human clinical trial that were not predicted by nonclinical safety studies. A recent report concluded that a TK-NOG humanized liver (hu-liver) mouse model detected human-specific FIAU liver toxicity, and broader use of that model could improve drug safety testing. We further evaluated this model at similar dose levels to assess FIAU sensitivity and potential mechanistic biomarkers. Although we were unable to reproduce the marked acute liver toxicity with two separate studies (including one with a "sensitized" donor), we identified molecular biomarkers reflecting the early stages of FIAU mitochondrial toxicity, which were not seen with its stereoisomer (FIRU). Dose dependent FIAU-induced changes in hu-liver mice included more pronounced reductions in mitochondrial to nuclear DNA (mtDNA/nucDNA) ratios in human hepatocytes compared to mouse hepatocytes and kidneys of the same animals. FIAU treatment also triggered a p53 transcriptional response and opposing changes in transcripts of nuclear- and mitochondrial-encoded mitochondrial proteins. The time dependent accumulation of FIAU into mtDNA is consistent with the ≥9-week latency of liver toxicity observed for FIAU in the clinic. Similar changes were observed in an in vitro micro-patterned hepatocyte coculture system. In addition, FIAU-dependent mtDNA/nucDNA ratio and transcriptional alterations, especially reductions in mitochondrially encoded transcripts, were seen in livers of non-engrafted TK-NOG and CD-1 mice dosed for a shorter period. Conclusion: These mechanistic biomarker findings can be leveraged in an in vitro model and in a more routine preclinical model (CD-1 mice) to identify nucleosides with such a FIAU-like mitochondrial toxicity mechanistic liability potential. Further optimization of the TK-NOG hu-liver mouse model is necessary before broader adoption for drug safety testing.

A Collaborative Initiative to Establish Genomic Biomarkers for Assessing Tumorigenic Potential to Reduce Reliance on Conventional Rodent Carcinogenicity Studies.

There is growing recognition across broad sectors of the scientific community that use of genomic biomarkers has the potential to reduce the need for conventional rodent carcinogenicity studies of industrial chemicals, agrochemicals, and pharmaceuticals through a weight-of-evidence approach. These biomarkers fall into 2 major categories: (1) sets of gene transcripts that can identify distinct tumorigenic mechanisms of action; and (2) cancer driver gene mutations indicative of rapidly expanding growth-advantaged clonal cell populations. This call-to-action article describes a collaborative approach launched to develop and qualify biomarker gene expression panels that measure widely accepted molecular pathways linked to tumorigenesis and their activation levels to predict tumorigenic doses of chemicals from short-term exposures. Growing evidence suggests that application of such biomarker panels in short-term exposure rodent studies can identify both tumorigenic hazard and tumorigenic activation levels for chemical-induced carcinogenicity. In the future, this approach will be expanded to include methodologies examining mutations in key cancer driver gene mutation hotspots as biomarkers of both genotoxic and nongenotoxic chemical tumor risk. Analytical, technical, and biological validation studies of these complementary genomic tools are being undertaken by multisector and multidisciplinary collaborative teams within the Health and Environmental Sciences Institute. Success from these efforts will facilitate the transition from current heavy reliance on conventional 2-year rodent carcinogenicity studies to more rapid animal- and resource-sparing approaches for mechanism-based carcinogenicity evaluation supporting internal and regulatory decision-making.

Induction of Alzheimer's-like changes in brain of mice expressing mutant APP fed excess methionine.

Elevated plasma homocysteine, a risk factor for Alzheimer's disease, could result from increased production from methionine or by inefficient clearance by folate- and B-vitamin-dependent pathways. Understanding the relative contributions of these processes to pathogenesis is important for therapeutic strategies designed to lower homocysteine. To assess these alternatives, we elevated plasma homocysteine by feeding mutant amyloid precursor protein (APP)-expressing mice diets with either high methionine (HM) or deficient in B-vitamins and folate (B Def). Mutant APP mice fed HM demonstrated increased brain beta amyloid. Interestingly, this increase was not observed in mutant APP mice fed B Def diet, nor was it observed in C57Bl6 or YAC-APP mice fed HM. Furthermore, HM, but not B Def, produced a prolonged increase in brain homocysteine only in mutant APP mice but not wild-type mice. These changes were time-dependent over 10 weeks. Further, by 10 weeks HM increased brain cholesterol and phosphorylated tau in mutant APP mice. Transcriptional profiling experiments revealed robust differences in RNA expression between C57Bl6 and mutant APP mice. The HM diet in C57Bl6 mice transiently induced a transcriptional profile similar to mutant APP cortex, peaking at 2 weeks , following a time course comparable to brain homocysteine changes. Together, these data suggest a link between APP and methionine metabolism.

Understanding the effect of increased cell specific productivity on galactosylation of monoclonal antibodies produced using Chinese hamster ovary cells.

Achieving optimal productivity and desired product quality of the therapeutic monoclonal antibody (mAb) is one of the primary goals of process development. Across the various mAb programs at our company, we observed that increasing the specific productivity (qp) results in a decrease in the % galactosylation (%Gal) level on the protein. In order to gain further insight into this correlation, cells were cultured under different process conditions such as pH or media osmolality or in the presence of supplements such as sodium butyrate. A range of qp and N-glycan profiles were obtained with the greatest changes observed under high pH (lower qp, higher %Gal), higher osmolality (higher qp, lower %Gal) or sodium butyrate (moderately higher qp, moderately lower %Gal) conditions. Abundance of individual glycan species highlighted different bottlenecks in the N-glycosylation pathway depending on the treatment condition. Transcriptomics analysis was performed to identify changes in gene expression profiles that correlate with the inverse relationship between qp and %Gal. Results showed downregulation of Beta-1,4-galactosyltransferase 1 (B4GalT1), UDP-GlcNAc and Mn2+ transporter (slc35a3 and slc39a8 respectively) for the high osmolality conditions. Significant downregulation of slc39a8 (Mn2+ transporter) was observed for the sodium butyrate condition. No significant differences were observed for any of the genes in the N-glycosylation pathway under the high pH condition even though this condition showed highest %Gal. Together, data suggests that different treatments have distinct complex mechanisms by which the overall glycan levels of a mAb are influenced. Further studies based on these results will help build the knowledge necessary to design strategies to obtain the desired productivity and product quality of mAbs.

Quantitative Transcriptional Biomarkers of Xenobiotic Receptor Activation in Rat Liver for the Early Assessment of Drug Safety Liabilities.

The robust transcriptional plasticity of liver mediated through xenobiotic receptors underlies its ability to respond rapidly and effectively to diverse chemical stressors. Thus, drug-induced gene expression changes in liver serve not only as biomarkers of liver injury, but also as mechanistic sentinels of adaptation in metabolism, detoxification, and tissue protection from chemicals. Modern RNA sequencing methods offer an unmatched opportunity to quantitatively monitor these processes in parallel and to contextualize the spectrum of dose-dependent stress, adaptation, protection, and injury responses induced in liver by drug treatments. Using this approach, we profiled the transcriptional changes in rat liver following daily oral administration of 120 different compounds, many of which are known to be associated with clinical risk for drug-induced liver injury by diverse mechanisms. Clustering, correlation, and linear modeling analyses were used to identify and optimize coexpressed gene signatures modulated by drug treatment. Here, we specifically focused on prioritizing 9 key signatures for their pragmatic utility for routine monitoring in initial rat tolerability studies just prior to entering drug development. These signatures are associated with 5 canonical xenobiotic nuclear receptors (AHR, CAR, PXR, PPARα, ER), 3 mediators of reactive metabolite-mediated stress responses (NRF2, NRF1, P53), and 1 liver response following activation of the innate immune response. Comparing paradigm chemical inducers of each receptor to the other compounds surveyed enabled us to identify sets of optimized gene expression panels and associated scoring algorithms proposed as quantitative mechanistic biomarkers with high sensitivity, specificity, and quantitative accuracy. These findings were further qualified using public datasets, Open TG-GATEs and DrugMatrix, and internal development compounds. With broader collaboration and additional qualification, the quantitative toxicogenomic framework described here could inform candidate selection prior to committing to drug development, as well as complement and provide a deeper understanding of the conventional toxicology study endpoints used later in drug development.

Development and Application of a Transcriptomic Signature of Bioactivation in an Advanced In Vitro Liver Model to Reduce Drug-induced Liver Injury Risk Early in the Pharmaceutical Pipeline.

Early risk assessment of drug-induced liver injury (DILI) potential for drug candidates remains a major challenge for pharmaceutical development. We have previously developed a set of rat liver transcriptional biomarkers in short-term toxicity studies to inform the potential of drug candidates to generate a high burden of chemically reactive metabolites that presents higher risk for human DILI. Here, we describe translation of those NRF1-/NRF2-mediated liver tissue biomarkers to an in vitro assay using an advanced micropatterned coculture system (HEPATOPAC) with primary hepatocytes from male Wistar Han rats. A 9-day, resource-sparing and higher throughput approach designed to identify new chemical entities with lower reactive metabolite-forming potential was qualified for internal decision making using 93 DILI-positive and -negative drugs. This assay provides 81% sensitivity and 90% specificity in detecting hepatotoxicants when a positive test outcome is defined as the bioactivation signature score of a test drug exceeding the threshold value at an in vitro test concentration that falls within 3-fold of the estimated maximum drug concentration at the human liver inlet following highest recommended clinical dose administrations. Using paired examples of compounds from distinct chemical series and close structural analogs, we demonstrate that this assay can differentiate drugs with lower DILI risk. The utility of this in vitro transcriptomic approach was also examined using human HEPATOPAC from a single donor, yielding 68% sensitivity and 86% specificity when the aforementioned criteria are applied to the same 93-drug test set. Routine use of the rat model has been adopted with deployment of the human model as warranted on a case-by-case basis. This in vitro transcriptomic signature-based strategy can be used early in drug discovery to derisk DILI potential from chemically reactive metabolites by guiding structure-activity relationship hypotheses and candidate selection.

Author Correction: The Transcriptomic Signature Of Disease Development And Progression Of Nonalcoholic Fatty Liver Disease.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Application of a Rat Liver Drug Bioactivation Transcriptional Response Assay Early in Drug Development That Informs Chemically Reactive Metabolite Formation and Potential for Drug-induced Liver Injury.

Drug-induced liver injury is a major reason for drug candidate attrition from development, denied commercialization, market withdrawal, and restricted prescribing of pharmaceuticals. The metabolic bioactivation of drugs to chemically reactive metabolites (CRMs) contribute to liver-associated adverse drug reactions in humans that often goes undetected in conventional animal toxicology studies. A challenge for pharmaceutical drug discovery has been reliably selecting drug candidates with a low liability of forming CRM and reduced drug-induced liver injury potential, at projected therapeutic doses, without falsely restricting the development of safe drugs. We have developed an in vivo rat liver transcriptional signature biomarker reflecting the cellular response to drug bioactivation. Measurement of transcriptional activation of integrated nuclear factor erythroid 2-related factor 2 (NRF2)/Kelch-like ECH-associated protein 1 (KEAP1) electrophilic stress, and nuclear factor erythroid 2-related factor 1 (NRF1) proteasomal endoplasmic reticulum (ER) stress responses, is described for discerning estimated clinical doses of drugs with potential for bioactivation-mediated hepatotoxicity. The approach was established using well benchmarked CRM forming test agents from our company. This was subsequently tested using curated lists of commercial drugs and internal compounds, anchored in the clinical experience with human hepatotoxicity, while agnostic to mechanism. Based on results with 116 compounds in short-term rat studies, with consideration of the maximum recommended daily clinical dose, this CRM mechanism-based approach yielded 32% sensitivity and 92% specificity for discriminating safe from hepatotoxic drugs. The approach adds new information for guiding early candidate selection and informs structure activity relationships (SAR) thus enabling lead optimization and mechanistic problem solving. Additional refinement of the model is ongoing. Case examples are provided describing the strengths and limitations of the approach.

Reconstruction and stability of secondary structure elements in the context of protein structure prediction.

Efficient and accurate reconstruction of secondary structure elements in the context of protein structure prediction is the major focus of this work. We present a novel approach capable of reconstructing alpha-helices and beta-sheets in atomic detail. The method is based on Metropolis Monte Carlo simulations in a force field of empirical potentials that are designed to stabilize secondary structure elements in room-temperature simulations. Particular attention is paid to lateral side-chain interactions in beta-sheets and between the turns of alpha-helices, as well as backbone hydrogen bonding. The force constants are optimized using contrastive divergence, a novel machine learning technique, from a data set of known structures. Using this approach, we demonstrate the applicability of the framework to the problem of reconstructing the overall protein fold for a number of commonly studied small proteins, based on only predicted secondary structure and contact map. For protein G and chymotrypsin inhibitor 2, we are able to reconstruct the secondary structure elements in atomic detail and the overall protein folds with a root mean-square deviation of <10 A. For cold-shock protein and the SH3 domain, we accurately reproduce the secondary structure elements and the topology of the 5-stranded beta-sheets, but not the barrel structure. The importance of high-quality secondary structure and contact map prediction is discussed.

AhR Activation in Pharmaceutical Development: Applying Liver Gene Expression Biomarker Thresholds to Identify Doses Associated With Tumorigenic Risks in Rats.

Aryl hydrocarbon receptor (AhR) activation is associated with carcinogenicity of non-genotoxic AhR-activating carcinogens such as 2,3,7,8-tetrachlorodibenzodioxin (TCDD), and is often observed with drug candidate molecules in development and raises safety concerns. As downstream effectors of AhR signaling, the expression and activity of Cyp1a1 and Cyp1a2 genes are commonly monitored as evidence of AhR activation to inform carcinogenic risk of compounds in question. However, many marketed drugs and phytochemicals are reported to induce these Cyps modestly and are not associated with dioxin-like toxicity or carcinogenicity. We hypothesized that a threshold of AhR activation needs to be surpassed in a sustained manner in order for the dioxin-like toxicity to manifest, and a simple liver gene expression signature based on Cyp1a1 and Cyp1a2 from a short-term rat study could be used to assess AhR activation strength and differentiate tumorigenic dose levels from non-tumorigenic ones. To test this hypothesis, short-term studies were conducted in Wistar Han rats with 2 AhR-activating carcinogens (TCDD and PCB126) at minimally carcinogenic and noncarcinogenic dose levels, and 3 AhR-activating noncarcinogens (omeprazole, mexiletine, and canagliflozin) at the top doses used in their reported 2-year rat carcinogenicity studies. A threshold of AhR activation was identified in rat liver that separated a meaningful "tumorigenic-strength AhR signal" from a statistically significant AhR activation signal that was not associated with dioxin-like carcinogenicity. These studies also confirmed the importance of the sustainability of AhR activation for carcinogenic potential. A sustained activation of AhR above the threshold could thus be used in early pharmaceutical development to identify dose levels of drug candidates expected to exhibit dioxin-like carcinogenic potential.

Comparative transcriptomics of choroid plexus in Alzheimer's disease, frontotemporal dementia and Huntington's disease: implications for CSF homeostasis.

BACKGROUND: In Alzheimer's disease, there are striking changes in CSF composition that relate to altered choroid plexus (CP) function. Studying CP tissue gene expression at the blood-cerebrospinal fluid barrier could provide further insight into the epithelial and stromal responses to neurodegenerative disease states. METHODS: Transcriptome-wide Affymetrix microarrays were used to determine disease-related changes in gene expression in human CP. RNA from post-mortem samples of the entire lateral ventricular choroid plexus was extracted from 6 healthy controls (Ctrl), 7 patients with advanced (Braak and Braak stage III-VI) Alzheimer's disease (AD), 4 with frontotemporal dementia (FTD) and 3 with Huntington's disease (HuD). Statistics and agglomerative clustering were accomplished with MathWorks, MatLab; and gene set annotations by comparing input sets to GeneGo ( http://www.genego.com ) and Ingenuity ( http://www.ingenuity.com ) pathway sets. Bonferroni-corrected hypergeometric p-values of < 0.1 were considered a significant overlap between sets. RESULTS: Pronounced differences in gene expression occurred in CP of advanced AD patients vs. Ctrls. Metabolic and immune-related pathways including acute phase response, cytokine, cell adhesion, interferons, and JAK-STAT as well as mTOR were significantly enriched among the genes upregulated. Methionine degradation, claudin-5 and protein translation genes were downregulated. Many gene expression changes in AD patients were observed in FTD and HuD (e.g., claudin-5, tight junction downregulation), but there were significant differences between the disease groups. In AD and HuD (but not FTD), several neuroimmune-modulating interferons were significantly enriched (e.g., in AD: IFI-TM1, IFN-AR1, IFN-AR2, and IFN-GR2). AD-associated expression changes, but not those in HuD and FTD, were enriched for upregulation of VEGF signaling and immune response proteins, e.g., interleukins. HuD and FTD patients distinctively displayed upregulated cadherin-mediated adhesion. CONCLUSIONS: Our transcript data for human CP tissue provides genomic and mechanistic insight for differential expression in AD vs. FTD vs. HuD for stromal as well as epithelial components. These choroidal transcriptome characterizations elucidate immune activation, tissue functional resiliency, and CSF metabolic homeostasis. The BCSFB undergoes harmful, but also important functional and adaptive changes in neurodegenerative diseases; accordingly, the enriched JAK-STAT and mTOR pathways, respectively, likely help the CP in adaptive transcription and epithelial repair and/or replacement when harmed by neurodegeneration pathophysiology. We anticipate that these precise CP translational data will facilitate pharmacologic/transgenic therapies to alleviate dementia.

Learning about protein hydrogen bonding by minimizing contrastive divergence.

Defining the strength and geometry of hydrogen bonds in protein structures has been a challenging task since early days of structural biology. In this article, we apply a novel statistical machine learning technique, known as contrastive divergence, to efficiently estimate both the hydrogen bond strength and the geometric characteristics of strong interpeptide backbone hydrogen bonds, from a dataset of structures representing a variety of different protein folds. Despite the simplifying assumptions of the interatomic energy terms used, we determine the strength of these hydrogen bonds to be between 1.1 and 1.5 kcal/mol, in good agreement with earlier experimental estimates. The geometry of these strong backbone hydrogen bonds features an almost linear arrangement of all four atoms involved in hydrogen bond formation. We estimate that about a quarter of all hydrogen bond donors and acceptors participate in these strong interpeptide hydrogen bonds.

Transcriptional Profiling of Cholinergic Neurons From Basal Forebrain Identifies Changes in Expression of Genes Between Sleep and Wake.

Study objective: To assess differences in gene expression in cholinergic basal forebrain cells between sleeping and sleep-deprived mice sacrificed at the same time of day. Methods: Tg(ChAT-eGFP)86Gsat mice expressing enhanced green fluorescent protein (eGFP) under control of the choline acetyltransferase (Chat) promoter were utilized to guide laser capture of cholinergic cells in basal forebrain. Messenger RNA expression levels in these cells were profiled using microarrays. Gene expression in eGFP(+) neurons was compared (1) to that in eGFP(-) neurons and to adjacent white matter, (2) between 7:00 am (lights on) and 7:00 pm (lights off), (3) between sleep-deprived and sleeping animals at 0, 3, 6, and 9 hours from lights on. Results: There was a marked enrichment of ChAT and other markers of cholinergic neurons in eGFP(+) cells. Comparison of gene expression in these eGFP(+) neurons between 7:00 am and 7:00 pm revealed expected differences in the expression of clock genes (Arntl2, Per1, Per2, Dbp, Nr1d1) as well as mGluR3. Comparison of expression between spontaneous sleep and sleep-deprived groups sacrificed at the same time of day revealed a number of transcripts (n = 55) that had higher expression in sleep deprivation compared to sleep. Genes upregulated in sleep deprivation predominantly were from the protein folding pathway (25 transcripts, including chaperones). Among 42 transcripts upregulated in sleep was the cold-inducible RNA-binding protein. Conclusions: Cholinergic cell signatures were characterized. Whether the identified genes are changing as a consequence of differences in behavioral state or as part of the molecular regulatory mechanism remains to be determined.

The Transcriptomic Signature Of Disease Development And Progression Of Nonalcoholic Fatty Liver Disease.

A longitudinal molecular model of the development and progression of nonalcoholic fatty liver disease (NAFLD) over time is lacking. We have recently validated a high fat/sugar water-induced animal (an isogenic strain of C57BL/6 J:129S1/SvImJ mice) model of NAFLD that closely mimics most aspects of human disease. The hepatic transcriptome of such mice with fatty liver (8 weeks), steatohepatitis with early fibrosis (16-24 weeks) and advanced fibrosis (52 weeks) after initiation of the diet was evaluated and compared to mice on chow diet. Fatty liver development was associated with transcriptional activation of lipogenesis, FXR-RXR, PPAR-α mediated lipid oxidation and oxidative stress pathways. With progression to steatohepatitis, metabolic pathway activation persisted with additional activation of IL-1/inhibition of RXR, granulocyte diapedesis/adhesion, Fc macrophage activation, prothrombin activation and hepatic stellate cell activation. Progression to advanced fibrosis was associated with dampening of metabolic, oxidative stress and cell stress related pathway activation but with further Fc macrophage activation, cell death and turnover and activation of cancer-related networks. The molecular progression of NAFLD involves a metabolic perturbation which triggers subsequent cell stress and inflammation driving cell death and turnover. Over time, inflammation and fibrogenic pathways become dominant while in advanced disease an inflammatory-oncogenic profile dominates.

Bayesian segmental models with multiple sequence alignment profiles for protein secondary structure and contact map prediction.

In this paper, we develop a segmental semi-Markov model (SSMM) for protein secondary structure prediction which incorporates multiple sequence alignment profiles with the purpose of improving the predictive performance. The segmental model is a generalization of the hidden Markov model where a hidden state generates segments of various length and secondary structure type. A novel parameterized model is proposed for the likelihood function that explicitly represents multiple sequence alignment profiles to capture the segmental conformation. Numerical results on benchmark data sets show that incorporating the profiles results in substantial improvements and the generalization performance is promising. By incorporating the information from long range interactions in beta-sheets, this model is also capable of carrying out inference on contact maps. This is an important advantage of probabilistic generative models over the traditional discriminative approach to protein secondary structure prediction. The Web server of our algorithm and supplementary materials are available at http://public.kgi.edu/-wild/bsm.html.

Discovering Biology in Periodic Data through Phase Set Enrichment Analysis (PSEA).

Several tools use prior biological knowledge to interpret gene expression data. However, existing enrichment tools assume that variables are monotonic and incorrectly measure the distance between periodic phases. As a result, these tools are poorly suited for the analysis of the cell cycle, circadian clock, or other periodic systems. Here, we develop Phase Set Enrichment Analysis (PSEA) to incorporate prior knowledge into the analysis of periodic data. PSEA identifies biologically related gene sets showing temporally coordinated expression. Using synthetic gene sets of various sizes generated from von Mises (circular normal) distributions, we benchmarked PSEA alongside existing methods. PSEA offered enhanced sensitivity over a broad range of von Mises distributions and gene set sizes. Importantly, and unlike existing tools, the sensitivity of PSEA is independent of the mean expression phase of the set. We applied PSEA to 4 published datasets. Application of PSEA to the mouse circadian atlas revealed that several pathways, including those regulating immune and cell-cycle function, demonstrate temporal orchestration across multiple tissues. We then applied PSEA to the phase shifts following a restricted feeding paradigm. We found that this perturbation disrupts intraorgan metabolic synchrony in the liver, altering the timing between anabolic and catabolic pathways. Reanalysis of expression data using custom gene sets derived from recent ChIP-seq results revealed circadian transcriptional targets bound exclusively by CLOCK, independently of BMAL1, differ from other exclusive circadian output genes and have well-synchronized phases. Finally, we used PSEA to compare 2 cell-cycle datasets. PSEA increased the apparent biological overlap while also revealing evidence of cell-cycle dysregulation in these cancer cells. To encourage its use by the community, we have implemented PSEA as a Java application. In sum, PSEA offers a powerful new tool to investigate large-scale, periodic data for biological insight.

Effects of Long-Term Odanacatib Treatment on Bone Gene Expression in Ovariectomized Adult Rhesus Monkeys: Differentiation From Alendronate.

Similar efficacy of the cathepsin K inhibitor odanacatib (ODN) and the bisphosphonate alendronate (ALN) in reducing bone turnover markers and increasing bone mineral density in spine and hip were previously demonstrated in ovariectomized (OVX)-monkeys treated for 20 months in prevention mode. Here, we profiled RNA from tibial metaphysis and diaphysis of the same study using Affymetrix microarrays, and selected 204 probe sets (p < 0.001, three-group ANOVA) that were differentially regulated by ODN or ALN versus vehicle. Both drugs produced strikingly different effects on known bone-related genes and pathways at the transcriptional level. Although ALN either reduced or had neutral effects on bone resorption-related genes, ODN significantly increased the expression of osteoclast genes (eg, APC5, TNFRSF11A, CTSK, ITGB3, and CALCR), consistent with previous findings on the effects of this agent in enhancing the number of nonresorbing osteoclasts. Conversely, ALN reduced the expression of known bone formation-related genes (eg, TGFBR1, SPP1, RUNX2, and PTH1R), whereas ODN either increased or had neutral effects on their expression. These differential effects of ODN versus ALN on bone resorption and formation were highly correlative to the changes in bone turnover markers, cathepsin K (Catk) target engagement marker serum C-terminal cross-linked telopeptide (1-CTP) and osteoclast marker tartrate resistant acid phosphatase isoform 5b (TRAP5b) in the same monkeys. Overall, the molecular profiling results are consistent with the known pharmacological actions of these agents on bone remodeling and clearly differentiate the molecular mechanisms of ODN from the bisphosphonates.

Network-driven plasma proteomics expose molecular changes in the Alzheimer's brain.

BACKGROUND: Biological pathways that significantly contribute to sporadic Alzheimer's disease are largely unknown and cannot be observed directly. Cognitive symptoms appear only decades after the molecular disease onset, further complicating analyses. As a consequence, molecular research is often restricted to late-stage post-mortem studies of brain tissue. However, the disease process is expected to trigger numerous cellular signaling pathways and modulate the local and systemic environment, and resulting changes in secreted signaling molecules carry information about otherwise inaccessible pathological processes. RESULTS: To access this information we probed relative levels of close to 600 secreted signaling proteins from patients' blood samples using antibody microarrays and mapped disease-specific molecular networks. Using these networks as seeds we then employed independent genome and transcriptome data sets to corroborate potential pathogenic pathways. CONCLUSIONS: We identified Growth-Differentiation Factor (GDF) signaling as a novel Alzheimer's disease-relevant pathway supported by in vivo and in vitro follow-up experiments, demonstrating the existence of a highly informative link between cellular pathology and changes in circulatory signaling proteins.

Exhaustive Metropolis Monte Carlo sampling and analysis of polyalanine conformations adopted under the influence of hydrogen bonds.

We propose a novel Metropolis Monte Carlo procedure for protein modeling and analyze the influence of hydrogen bonding on the distribution of polyalanine conformations. We use an atomistic model of the polyalanine chain with rigid and planar polypeptide bonds, and elastic alpha carbon valence geometry. We adopt a simplified energy function in which only hard-sphere repulsion and hydrogen bonding interactions between the atoms are considered. Our Metropolis Monte Carlo procedure utilizes local crankshaft moves and is combined with parallel tempering to exhaustively sample the conformations of 16-mer polyalanine. We confirm that Flory's isolated-pair hypothesis (the steric independence between the dihedral angles of individual amino acids) does not hold true in long polypeptide chains. In addition to 3(10)- and alpha-helices, we identify a kink stabilized by 2 hydrogen bonds with a shared acceptor as a common structural motif. Varying the strength of hydrogen bonds, we induce the helix-coil transition in the model polypeptide chain. We compare the propensities for various hydrogen bonding patterns and determine the degree of cooperativity of hydrogen bond formation in terms of the Hill coefficient. The observed helix-coil transition is also quantified according to Zimm-Bragg theory.