Monte Carlo simulations of the HP model (the "Ising model" of protein folding)

Using Wang-Landau sampling with suitable Monte Carlo trial moves (pull moves and bond-rebridging moves combined) we have determined the density of states and thermodynamic properties for a short sequence of the HP protein model. For free chains these proteins are known to first undergo a collapse "transition" to a globule state followed by a second "transition" into a native state. When placed in the proximity of an attractive surface, there is a competition between surface adsorption and folding that leads to an intriguing sequence of "transitions". These transitions depend upon the relative interaction strengths and are largely inaccessible to "standard" Monte Carlo methods.Comment: 6 pages, 6 figures. Article in press. To be published in Computer Physics Communications (2011

Simple Models of the Protein Folding Problem

The protein folding problem has attracted an increasing attention from physicists. The problem has a flavor of statistical mechanics, but possesses the most common feature of most biological problems -- the profound effects of evolution. I will give an introduction to the problem, and then focus on some recent work concerning the so-called ``designability principle''. The designability of a structure is measured by the number of sequences that have that structure as their unique ground state. Structures differ drastically in terms of their designability; highly designable structures emerge with a number of associated sequences much larger than the average. These highly designable structures 1) possess ``proteinlike'' secondary structures and motifs, 2) are thermodynamically more stable, and 3) fold faster than other structures. These results suggest that protein structures are selected in nature because they are readily designed and stable against mutations, and that such selection simultaneously leads to thermodynamic stability and foldability. According to this picture, a key to the protein folding problem is to understand the emergence and the properties of the highly designable structures.Comment: 21 pages, 14 figures. Invited talk at Dynamics Days Asian Pacific, Hong Kong, July 13-16, 1999. To appear in Physica

Glassy transition in a disordered model for the RNA secondary structure

We numerically study a disordered model for the RNA secondary structure and we find that it undergoes a phase transition, with a breaking of the replica symmetry in the low temperature region (like in spin glasses). Our results are based on the exact evaluation of the partition function.Comment: 4 pages, 3 figure

Steric constraints in model proteins

A simple lattice model for proteins that allows for distinct sizes of the amino acids is presented. The model is found to lead to a significant number of conformations that are the unique ground state of one or more sequences or encodable. Furthermore, several of the encodable structures are highly designable and are the non-degenerate ground state of several sequences. Even though the native state conformations are typically compact, not all compact conformations are encodable. The incorporation of the hydrophobic and polar nature of amino acids further enhances the attractive features of the model.Comment: RevTex, 5 pages, 3 postscript figure

Human lymphocytes mobilized with exercise have an anti-tumor transcriptomic profile and exert enhanced graft-versus-leukemia effects in xenogeneic mice

BackgroundEvery bout of exercise mobilizes and redistributes large numbers of effector lymphocytes with a cytotoxic and tissue migration phenotype. The frequent redistribution of these cells is purported to increase immune surveillance and play a mechanistic role in reducing cancer risk and slowing tumor progression in physically active cancer survivors. Our aim was to provide the first detailed single cell transcriptomic analysis of exercise-mobilized lymphocytes and test their effectiveness as a donor lymphocyte infusion (DLI) in xenogeneic mice engrafted with human leukemia.MethodsPeripheral blood mononuclear cells (PBMCs) were collected from healthy volunteers at rest and at the end of an acute bout of cycling exercise. Flow cytometry and single-cell RNA sequencing was performed to identify phenotypic and transcriptomic differences between resting and exercise-mobilized cells using a targeted gene expression panel curated for human immunology. PBMCs were injected into the tail vein of xenogeneic NSG-IL-15 mice and subsequently challenged with a luciferase tagged chronic myelogenous leukemia cell line (K562). Tumor growth (bioluminescence) and xenogeneic graft-versus-host disease (GvHD) were monitored bi-weekly for 40-days.ResultsExercise preferentially mobilized NK-cell, CD8+ T-cell and monocyte subtypes with a differentiated and effector phenotype, without significantly mobilizing CD4+ regulatory T-cells. Mobilized effector lymphocytes, particularly effector-memory CD8+ T-cells and NK-cells, displayed differentially expressed genes and enriched gene sets associated with anti-tumor activity, including cytotoxicity, migration/chemotaxis, antigen binding, cytokine responsiveness and alloreactivity (e.g. graft-versus-host/leukemia). Mice receiving exercise-mobilized PBMCs had lower tumor burden and higher overall survival (4.14E+08 photons/s and 47%, respectively) at day 40 compared to mice receiving resting PBMCs (12.1E+08 photons/s and 22%, respectively) from the same donors (p<0.05). Human immune cell engraftment was similar for resting and exercise-mobilized DLI. However, when compared to non-tumor bearing mice, K562 increased the expansion of NK-cell and CD3+/CD4-/CD8- T-cells in mice receiving exercise-mobilized but not resting lymphocytes, 1-2 weeks after DLI. No differences in GvHD or GvHD-free survival was observed between groups either with or without K562 challenge.ConclusionExercise in humans mobilizes effector lymphocytes with an anti-tumor transcriptomic profile and their use as DLI extends survival and enhances the graft-versus-leukemia (GvL) effect without exacerbating GvHD in human leukemia bearing xenogeneic mice. Exercise may serve as an effective and economical adjuvant to increase the GvL effects of allogeneic cell therapies without intensifying GvHD

Rare Tumor-Normal Matched Whole Exome Sequencing Identifies Novel Genomic Pathogenic Germline and Somatic Aberrations

Whole exome sequencing (WES) of matched tumor-normal pairs in rare tumors has the potential to identify genome-wide mutations and copy number alterations (CNAs). We evaluated 27 rare cancer patients with tumor-normal matching by WES and tumor-only next generation sequencing (NGS) as a comparator. Our goal was to: 1) identify known and novel variants and CNAs in rare cancers with comparison to common cancers; 2) examine differences between germline and somatic variants and how that functionally impacts rare tumors; 3) detect and characterize alleles in biologically relevant genes-pathways that may be of clinical importance but not represented in classical cancer genes. We identified 3343 germline single nucleotide variants (SNVs) and small indel variants-1670 in oncogenes and 1673 in tumor suppressor genes-generating an average of 124 germline variants/case. The number of somatic SNVs and small indels detected in all cases was 523:306 in oncogenes and 217 in tumor suppressor genes. Of the germline variants, six were identified to be pathogenic or likely pathogenic. In the 27 analyzed rare cancer cases, CNAs are variable depending on tumor type, germline pathogenic variants are more common. Cell fate pathway mutations (e.g., Hippo, Notch, Wnt) dominate pathogenesis and double hit (mutation + CNV) represent similar to 18% cases.Open access journalThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Leukocyte expression profiles reveal gene sets with prognostic value for seizure-free outcome following stereotactic laser amygdalohippocampotomy

Among patients with intractable epilepsy, the most commonly performed surgical procedure is craniotomy for amygdalohippocampectomy (AH). Stereotactic laser amygdalohippocampotomy (SLAH) has also been recently employed as a minimally invasive treatment for intractable temporal lobe epilepsy (TLE). Among patients treated with AH and SLAH approximately 65% and 54% of patients become seizure-free, respectively. Therefore, selection criteria for surgical candidates with improved prognostic value for post-operative seizure-free outcome are greatly needed. In this study, we perform RNA sequencing (RNA-Seq) on whole blood leukocyte samples taken from 16 patients with intractable TLE prior to SLAH to test the hypothesis that pre-operative leukocyte RNA expression profiles are prognostic for post-operative seizure outcome. Multidimensional scaling analysis of the RNA expression data indicated separate clustering of patients with seizure free (SF) and non-seizurefree (NSF) outcomes. Differential expression (DE) analysis performed on SF versus NSF groups revealed 24 significantly differentially expressed genes (>= 2.0-fold change, p-value < 0.05, FDR < 0.05). Network and pathway analyses identified differential activation of pathways involved in lipid metabolism, morphology of oligodendrocytes, inflammatory response, and development of astrocytes. These results suggest that pre-operative leukocyte expression profiles have prognostic value for seizure outcome following SLAH.NIH from Temple University [RO1MH065151, HCS 04-42]; Department of Surgery 2016 Faculty Seed Grant, University of Arizona College of Medicine; 2016 University of Arizona Health Sciences Clinical Research Pilot Program Award (UAHS-CRPA)Open access journal.This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Predictive network analysis identifies JMJD6 and other potential key drivers in Alzheimer’s disease

Despite decades of genetic studies on late-onset Alzheimer’s disease, the underlying molecular mechanisms remain unclear. To better comprehend its complex etiology, we use an integrative approach to build robust predictive (causal) network models using two large human multi-omics datasets. We delineate bulk-tissue gene expression into single cell-type gene expression and integrate clinical and pathologic traits, single nucleotide variation, and deconvoluted gene expression for the construction of cell type-specific predictive network models. Here, we focus on neuron-specific network models and prioritize 19 predicted key drivers modulating Alzheimer’s pathology, which we then validate by knockdown in human induced pluripotent stem cell-derived neurons. We find that neuronal knockdown of 10 of the 19 targets significantly modulates levels of amyloid-beta and/or phosphorylated tau peptides, most notably JMJD6. We also confirm our network structure by RNA sequencing in the neurons following knockdown of each of the 10 targets, which additionally predicts that they are upstream regulators of REST and VGF. Our work thus identifies robust neuronal key drivers of the Alzheimer’s-associated network state which may represent therapeutic targets with relevance to both amyloid and tau pathology in Alzheimer’s disease

Rare variants of small effect size in neuronal excitability genes influence clinical outcome in Japanese cases of <i>SCN1A</i> truncation-positive Dravet syndrome

<div><p>Dravet syndrome (DS) is a rare, devastating form of childhood epilepsy that is often associated with mutations in the voltage-gated sodium channel gene, <i>SCN1A</i>. There is considerable variability in expressivity within families, as well as among individuals carrying the same primary mutation, suggesting that clinical outcome is modulated by variants at other genes. To identify modifier gene variants that contribute to clinical outcome, we sequenced the exomes of 22 individuals at both ends of a phenotype distribution (i.e., mild and severe cognitive condition). We controlled for variation associated with different mutation types by limiting inclusion to individuals with a <i>de novo</i> truncation mutation resulting in <i>SCN1A</i> haploinsufficiency. We performed tests aimed at identifying 1) single common variants that are enriched in either phenotypic group, 2) sets of common or rare variants aggregated in and around genes associated with clinical outcome, and 3) rare variants in 237 candidate genes associated with neuronal excitability. While our power to identify enrichment of a common variant in either phenotypic group is limited as a result of the rarity of mild phenotypes in individuals with <i>SCN1A</i> truncation variants, our top candidates did not map to functional regions of genes, or in genes that are known to be associated with neurological pathways. In contrast, we found a statistically-significant excess of rare variants predicted to be damaging and of small effect size in genes associated with neuronal excitability in severely affected individuals. A <i>KCNQ2</i> variant previously associated with benign neonatal seizures is present in 3 of 12 individuals in the severe category. To compare our results with the healthy population, we performed a similar analysis on whole exome sequencing data from 70 Japanese individuals in the 1000 genomes project. Interestingly, the frequency of rare damaging variants in the same set of neuronal excitability genes in healthy individuals is nearly as high as in severely affected individuals. Rather than a single common gene/variant modifying clinical outcome in <i>SCN1A</i>-related epilepsies, our results point to the cumulative effect of rare variants with little to no measurable phenotypic effect (i.e., typical genetic background) unless present in combination with a disease-causing truncation mutation in <i>SCN1A</i>.</p></div