RESUMO
Identifying host genetic factors modulating immune checkpoint inhibitor (ICI) efficacy has been experimentally challenging because of variations in both host and tumor genomes, differences in the microbiome, and patient life exposures. Utilizing the Collaborative Cross (CC) multi-parent mouse genetic resource population, we developed an approach that fixes the tumor genomic configuration while varying host genetics. With this approach, we discovered that response to anti-PD-1 (aPD1) immunotherapy was significantly heritable in four distinct murine tumor models (H2 between 0.18-0.40). For the MC38 colorectal carcinoma system (H2 = 0.40), we mapped four significant ICI response quantitative trait loci (QTL) localized to mouse chromosomes (mChr) 5, 9, 15 and 17, and identified significant epistatic interactions between specific QTL pairs. Differentially expressed genes within these QTL were highly enriched for immune genes and pathways mediating allograft rejection and graft vs host disease. Using a cross species analytical approach, we found a core network of 48 genes within the four QTLs that showed significant prognostic value for overall survival in aPD1 treated human cohorts that outperformed all other existing validated immunotherapy biomarkers, especially in human tumors of the previously defined immune subtype 4. Functional blockade of two top candidate immune targets within the 48 gene network, GM-CSF and high affinity IL-2/IL-15 signaling, completely abrogated the MC38 tumor transcriptional response to aPD1 therapy in vivo. Thus, we have established a powerful cross species in vivo platform capable of uncovering host genetic factors that establish the tumor immune microenvironment configuration propitious for ICI response.
RESUMO
Outcomes for adult patients with a high-grade glioma continue to be dismal and new treatment paradigms are urgently needed. To optimize the opportunity for discovery, we performed a phase 0/1 dose-escalation clinical trial that investigated tumor pharmacokinetics, pharmacodynamics, and single nucleus transcriptomics following combined ribociclib (CDK4/6 inhibitor) and everolimus (mTOR inhibitor) treatment in recurrent high-grade glioma. Patients with a recurrent high-grade glioma (n = 24) harboring 1) CDKN2A / B deletion or CDK4 / 6 amplification, 2) PTEN loss or PIK3CA mutations, and 3) wild-type retinoblastoma protein (Rb) were enrolled. Patients received neoadjuvant ribociclib and everolimus treatment and no dose-limiting toxicities were observed. The median unbound ribociclib concentrations in Gadolinium non-enhancing tumor regions were 170 nM (range, 65 - 1770 nM) and 634 nM (range, 68 - 2345 nM) in patients receiving 5 days treatment at the daily dose of 400 and 600 mg, respectively. Unbound everolimus concentrations were below the limit of detection (< 0.1 nM) in both enhancing and non-enhancing tumor regions at all dose levels. We identified a significant decrease in MIB1 positive cells suggesting ribociclib-associated cell cycle inhibition. Single nuclei RNAseq (snRNA) based comparisons of 17 IDH-wild-type on-trial recurrences to 31 IDH-wild-type standard of care treated recurrences data demonstrated a significantly lower fraction of cycling and neural progenitor-like (NPC-like) malignant cell populations. We validated the CDK4/6 inhibitor-directed malignant cell state shifts using three patient-derived cell lines. The presented clinical trial highlights the value of integrating pharmacokinetics, pharmacodynamics, and single nucleus transcriptomics to assess treatment effects in phase 0/1 surgical tissues, including malignant cell state shifts. ClinicalTrials.gov identifier: NCT03834740 .
RESUMO
For over a century, early researchers sought to study biological organisms in a laboratory setting, leading to the generation of both in vitro and in vivo model systems. Patient-derived models of cancer (PDMCs) have more recently come to the forefront of preclinical cancer models and are even finding their way into clinical practice as part of functional precision medicine programs. The PDMC Consortium, supported by the Division of Cancer Biology in the National Cancer Institute of the National Institutes of Health, seeks to understand the biological principles that govern the various PDMC behaviors, particularly in response to perturbagens, such as cancer therapeutics. Based on collective experience from the consortium groups, we provide insight regarding PDMCs established both in vitro and in vivo, with a focus on practical matters related to developing and maintaining key cancer models through a series of vignettes. Although every model has the potential to offer valuable insights, the choice of the right model should be guided by the research question. However, recognizing the inherent constraints in each model is crucial. Our objective here is to delineate the strengths and limitations of each model as established by individual vignettes. Further advances in PDMCs and the development of novel model systems will enable us to better understand human biology and improve the study of human pathology in the lab.
RESUMO
Endometriosis is characterized by the growth of endometrial-like tissue outside the uterus. It affects many women during their reproductive age, causing years of pelvic pain and potential infertility. Its pathophysiology remains largely unknown, which limits early diagnosis and treatment. We characterized peritoneal and ovarian lesions at single-cell transcriptome resolution and compared them to matched eutopic endometrium, unaffected endometrium and organoids derived from these tissues, generating data on over 122,000 cells across 14 individuals. We spatially localized many of the cell types using imaging mass cytometry. We identify a perivascular mural cell specific to the peritoneal lesions, with dual roles in angiogenesis promotion and immune cell trafficking. We define an immunotolerant peritoneal niche, fundamental differences in eutopic endometrium and between lesion microenvironments and an unreported progenitor-like epithelial cell subpopulation. Altogether, this study provides a holistic view of the endometriosis microenvironment that represents a comprehensive cell atlas of the disease in individuals undergoing hormonal treatment, providing essential information for future therapeutics and diagnostics.
Assuntos
Coristoma , Endometriose , Cistos Ovarianos , Neoplasias Ovarianas , Coristoma/complicações , Coristoma/genética , Coristoma/metabolismo , Endometriose/genética , Endometriose/metabolismo , Endométrio/metabolismo , Feminino , Humanos , Cistos Ovarianos/complicações , Cistos Ovarianos/metabolismo , Cistos Ovarianos/patologia , Neoplasias Ovarianas/patologia , Análise de Célula Única , Microambiente TumoralRESUMO
A major rate-limiting step in developing more effective immunotherapies for GBM is our inadequate understanding of the cellular complexity and the molecular heterogeneity of immune infiltrates in gliomas. Here, we report an integrated analysis of 201,986 human glioma, immune, and other stromal cells at the single cell level. In doing so, we discover extensive spatial and molecular heterogeneity in immune infiltrates. We identify molecular signatures for nine distinct myeloid cell subtypes, of which five are independent prognostic indicators of glioma patient survival. Furthermore, we identify S100A4 as a regulator of immune suppressive T and myeloid cells in GBM and demonstrate that deleting S100a4 in non-cancer cells is sufficient to reprogram the immune landscape and significantly improve survival. This study provides insights into spatial, molecular, and functional heterogeneity of glioma and glioma-associated immune cells and demonstrates the utility of this dataset for discovering therapeutic targets for this poorly immunogenic cancer.
Assuntos
Imunoterapia , Proteína A4 de Ligação a Cálcio da Família S100/isolamento & purificação , Análise de Célula Única/métodos , Animais , Neoplasias Encefálicas/imunologia , Feminino , Glioma/imunologia , Humanos , Masculino , Camundongos , Camundongos Endogâmicos C57BL , Células Mieloides , Prognóstico , Proteína A4 de Ligação a Cálcio da Família S100/genética , Microambiente Tumoral/imunologiaRESUMO
Insulin resistance is a pathological state often associated with obesity, representing a major risk factor for type 2 diabetes. Limited mechanism-based strategies exist to alleviate insulin resistance. Here, using single-cell transcriptomics, we identify a small, critically important, but previously unexamined cell population, p21Cip1 highly expressing (p21high) cells, which accumulate in adipose tissue with obesity. By leveraging a p21-Cre mouse model, we demonstrate that intermittent clearance of p21high cells can both prevent and alleviate insulin resistance in obese mice. Exclusive inactivation of the NF-κB pathway within p21high cells, without killing them, attenuates insulin resistance. Moreover, fat transplantation experiments establish that p21high cells within fat are sufficient to cause insulin resistance in vivo. Importantly, a senolytic cocktail, dasatinib plus quercetin, eliminates p21high cells in human fat ex vivo and mitigates insulin resistance following xenotransplantation into immuno-deficient mice. Our findings lay the foundation for pursuing the targeting of p21high cells as a new therapy to alleviate insulin resistance.
Assuntos
Diabetes Mellitus Tipo 2 , Resistência à Insulina , Tecido Adiposo/metabolismo , Animais , Senescência Celular/fisiologia , Diabetes Mellitus Tipo 2/metabolismo , Dieta Hiperlipídica , Camundongos , Camundongos Endogâmicos C57BL , Obesidade/metabolismoRESUMO
The lateral hypothalamic area (LHA) is a highly conserved brain region critical for maintaining physiological homeostasis and goal-directed behavior. LHA neurons that express melanin-concentrating hormone (MCH) are key regulators of arousal, energy balance, and motivated behavior. However, cellular and functional diversity among LHAMCH neurons is not well understood. Previous anatomic and molecular data suggest that LHAMCH neurons may be parsed into at least two distinct subpopulations, one of which is enriched in neurokinin-3 receptor (NK3R), the receptor for neurokinin B (NKB), encoded by the Tac2 gene. This tachykininergic ligand-receptor system has been implicated in reproduction, fear memory, and stress in other brain regions, but NKB interactions with LHAMCH neurons are poorly understood. We first identified how LHAMCH subpopulations may be distinguished anatomically and electrophysiologically. To dissect functional connectivity between NKB-expressing neurons and LHAMCH neurons, we used Cre-dependent retrograde and anterograde viral tracing in male Tac2-Cre mice and identified Tac2/EYFP+ neurons in the bed nucleus of the stria terminalis and central nucleus of the amygdala, the central extended amygdala, as major sources of NKB input onto LHAMCH neurons. In addition to innervating the LHA, these limbic forebrain NKB neurons also project to midbrain and brainstem targets. Finally, using a dual-virus approach, we found that optogenetic activation of these inputs in slices evokes GABA release onto a subset of LHAMCH neurons but lacked specificity for the NK3R+ subpopulation. Overall, these data define parallel tachykininergic/GABAergic limbic forebrain projections that are positioned to modulate multiple nodes of homeostatic and behavioral control.SIGNIFICANCE STATEMENT The LHA orchestrates fundamental behavioral states in the mammalian hypothalamus, including arousal, energy balance, memory, stress, and motivated behavior. The neuropeptide MCH defines one prominent population of LHA neurons, with multiple roles in the regulation of homeostatic behavior. Outstanding questions remain concerning the upstream inputs that control MCH neurons. We sought to define neurochemically distinct pathways in the mouse brain that may communicate with specific MCH neuron subpopulations using viral-based retrograde and anterograde neural pathway tracing and optogenetics in brain slices. Here, we identify a specific neuropeptide-defined forebrain circuit that makes functional synaptic connections with MCH neuron subpopulations. This work lays the foundation for further manipulating molecularly distinct neural circuits that modulate innate behavioral states.
Assuntos
Núcleo Central da Amígdala/citologia , Região Hipotalâmica Lateral/citologia , Vias Neurais/citologia , Neurônios/citologia , Animais , Hormônios Hipotalâmicos/metabolismo , Masculino , Melaninas/metabolismo , Camundongos , Camundongos Transgênicos , Vias Neurais/metabolismo , Neurocinina B/metabolismo , Neurônios/metabolismo , Hormônios Hipofisários/metabolismoRESUMO
Glutamatergic neurons in the retrotrapezoid nucleus (RTN) function as respiratory chemoreceptors by regulating breathing in response to tissue CO2/H+. The RTN and greater parafacial region may also function as a chemosensing network composed of CO2/H+-sensitive excitatory and inhibitory synaptic interactions. In the context of disease, we showed that loss of inhibitory neural activity in a mouse model of Dravet syndrome disinhibited RTN chemoreceptors and destabilized breathing (Kuo et al., 2019). Despite this, contributions of parafacial inhibitory neurons to control of breathing are unknown, and synaptic properties of RTN neurons have not been characterized. Here, we show the parafacial region contains a limited diversity of inhibitory neurons including somatostatin (Sst)-, parvalbumin (Pvalb)-, and cholecystokinin (Cck)-expressing neurons. Of these, Sst-expressing interneurons appear uniquely inhibited by CO2/H+. We also show RTN chemoreceptors receive inhibitory input that is withdrawn in a CO2/H+-dependent manner, and chemogenetic suppression of Sst+ parafacial neurons, but not Pvalb+ or Cck+ neurons, increases baseline breathing. These results suggest Sst-expressing parafacial neurons contribute to RTN chemoreception and respiratory activity.
Assuntos
Dióxido de Carbono/metabolismo , Células Quimiorreceptoras/metabolismo , Epilepsias Mioclônicas/metabolismo , Hidrogênio/metabolismo , Núcleos Intralaminares do Tálamo/metabolismo , Pulmão/inervação , Respiração , Somatostatina/metabolismo , Animais , Modelos Animais de Doenças , Epilepsias Mioclônicas/genética , Epilepsias Mioclônicas/fisiopatologia , Feminino , Ácido Glutâmico/metabolismo , Núcleos Intralaminares do Tálamo/fisiopatologia , Masculino , Camundongos da Linhagem 129 , Camundongos Endogâmicos C57BL , Camundongos Transgênicos , Canal de Sódio Disparado por Voltagem NAV1.1/genética , Canal de Sódio Disparado por Voltagem NAV1.1/metabolismo , Inibição Neural , Somatostatina/genética , Transmissão SinápticaRESUMO
The ventral posterior hypothalamus (VPH) is an anatomically complex brain region implicated in arousal, reproduction, energy balance, and memory processing. However, neuronal cell type diversity within the VPH is poorly understood, an impediment to deconstructing the roles of distinct VPH circuits in physiology and behavior. To address this question, we employed a droplet-based single-cell RNA sequencing (scRNA-seq) approach to systematically classify molecularly distinct cell populations in the mouse VPH. Analysis of >16,000 single cells revealed 20 neuronal and 18 non-neuronal cell populations, defined by suites of discriminatory markers. We validated differentially expressed genes in selected neuronal populations through fluorescence in situ hybridization (FISH). Focusing on the mammillary bodies (MB), we discovered transcriptionally-distinct clusters that exhibit neuroanatomical parcellation within MB subdivisions and topographic projections to the thalamus. This single-cell transcriptomic atlas of VPH cell types provides a resource for interrogating the circuit-level mechanisms underlying the diverse functions of VPH circuits.
Assuntos
Hipotálamo Posterior/citologia , Animais , Feminino , Perfilação da Expressão Gênica , Hipotálamo Posterior/anatomia & histologia , Hipotálamo Posterior/fisiologia , Masculino , Camundongos , Camundongos Endogâmicos C57BL , RNA/genética , Análise de Sequência de RNA , Análise de Célula ÚnicaRESUMO
BACKGROUND: Cancer of unknown primary (CUP), representing approximately 3-5% of all malignancies, is defined as metastatic cancer where a primary site of origin cannot be found despite a standard diagnostic workup. Because knowledge of a patient's primary cancer remains fundamental to their treatment, CUP patients are significantly disadvantaged and most have a poor survival outcome. Developing robust and accessible diagnostic methods for resolving cancer tissue of origin, therefore, has significant value for CUP patients. METHODS: We developed an RNA-based classifier called CUP-AI-Dx that utilizes a 1D Inception convolutional neural network (1D-Inception) model to infer a tumor's primary tissue of origin. CUP-AI-Dx was trained using the transcriptional profiles of 18,217 primary tumours representing 32 cancer types from The Cancer Genome Atlas project (TCGA) and International Cancer Genome Consortium (ICGC). Gene expression data was ordered by gene chromosomal coordinates as input to the 1D-CNN model, and the model utilizes multiple convolutional kernels with different configurations simultaneously to improve generality. The model was optimized through extensive hyperparameter tuning, including different max-pooling layers and dropout settings. For 11 tumour types, we also developed a random forest model that can classify the tumour's molecular subtype according to prior TCGA studies. The optimised CUP-AI-Dx tissue of origin classifier was tested on 394 metastatic samples from 11 tumour types from TCGA and 92 formalin-fixed paraffin-embedded (FFPE) samples representing 18 cancer types from two clinical laboratories. The CUP-AI-Dx molecular subtype was also independently tested on independent ovarian and breast cancer microarray datasets FINDINGS: CUP-AI-Dx identifies the primary site with an overall top-1-accuracy of 98.54% in cross-validation and 96.70% on a test dataset. When applied to two independent clinical-grade RNA-seq datasets generated from two different institutes from the US and Australia, our model predicted the primary site with a top-1-accuracy of 86.96% and 72.46% respectively. INTERPRETATION: The CUP-AI-Dx predicts tumour primary site and molecular subtype with high accuracy and therefore can be used to assist the diagnostic work-up of cancers of unknown primary or uncertain origin using a common and accessible genomics platform. FUNDING: NIH R35 GM133562, NCI P30 CA034196, Victorian Cancer Agency Australia.
Assuntos
Inteligência Artificial , Biomarcadores Tumorais/genética , Biologia Computacional/métodos , Neoplasias Primárias Desconhecidas/diagnóstico , Neoplasias Primárias Desconhecidas/genética , RNA , Software , Algoritmos , Biologia Computacional/normas , Bases de Dados Genéticas , Genômica/métodos , Humanos , Aprendizado de Máquina , Metástase Neoplásica/diagnóstico , Metástase Neoplásica/genética , Redes Neurais de Computação , Reprodutibilidade dos Testes , Fluxo de TrabalhoRESUMO
Patients with systemic lupus erythematosus (SLE) display a complex blood transcriptome whose cellular origin is poorly resolved. Using single-cell RNA sequencing, we profiled ~276,000 peripheral blood mononuclear cells from 33 children with SLE with different degrees of disease activity and 11 matched controls. Increased expression of interferon-stimulated genes (ISGs) distinguished cells from children with SLE from healthy control cells. The high ISG expression signature (ISGhi) derived from a small number of transcriptionally defined subpopulations within major cell types, including monocytes, CD4+ and CD8+ T cells, natural killer cells, conventional and plasmacytoid dendritic cells, B cells and especially plasma cells. Expansion of unique subpopulations enriched in ISGs and/or in monogenic lupus-associated genes classified patients with the highest disease activity. Profiling of ~82,000 single peripheral blood mononuclear cells from adults with SLE confirmed the expansion of similar subpopulations in patients with the highest disease activity. This study lays the groundwork for resolving the origin of the SLE transcriptional signatures and the disease heterogeneity towards precision medicine applications.
Assuntos
Leucócitos Mononucleares/fisiologia , Lúpus Eritematoso Sistêmico/genética , Análise de Célula Única/métodos , Adolescente , Adulto , Células Cultivadas , Criança , Estudos de Coortes , Progressão da Doença , Feminino , Perfilação da Expressão Gênica , Humanos , Interferons/genética , Masculino , Análise de Sequência de RNA , Índice de Gravidade de Doença , TranscriptomaRESUMO
Cross-linking of high-affinity immunoglobulin E (IgE) results in the life-threatening allergic reaction anaphylaxis. Yet the cellular mechanisms that induce B cells to produce IgE in response to allergens remain poorly understood. T follicular helper (TFH) cells direct the affinity and isotype of antibodies produced by B cells. Although TFH cell-derived interleukin-4 (IL-4) is necessary for IgE production, it is not sufficient. We report a rare population of IL-13-producing TFH cells present in mice and humans with IgE to allergens, but not when allergen-specific IgE was absent or only low-affinity. These "TFH13" cells have an unusual cytokine profile (IL-13hiIL-4hiIL-5hiIL-21lo) and coexpress the transcription factors BCL6 and GATA3. TFH13 cells are required for production of high- but not low-affinity IgE and subsequent allergen-induced anaphylaxis. Blocking TFH13 cells may represent an alternative therapeutic target to ameliorate anaphylaxis.
Assuntos
Anafilaxia/imunologia , Imunoglobulina E/imunologia , Interleucina-13/metabolismo , Linfócitos T Auxiliares-Indutores/imunologia , Adolescente , Animais , Criança , Fator de Transcrição GATA3/metabolismo , Fatores de Troca do Nucleotídeo Guanina/genética , Fatores de Troca do Nucleotídeo Guanina/metabolismo , Humanos , Interleucina-13/genética , Camundongos , Camundongos Endogâmicos C57BL , Camundongos Mutantes , Proteínas Proto-Oncogênicas c-bcl-6/metabolismoRESUMO
Long-term maintenance of spermatogenesis in mammals is supported by GDNF, an essential growth factor required for spermatogonial stem cell (SSC) self-renewal. Exploiting a transgenic GDNF overexpression model, which expands and normalizes the pool of undifferentiated spermatogonia between Plzf +/+ and Plzf lu/lu mice, we used RNAseq to identify a rare subpopulation of cells that express EOMES, a T-box transcription factor. Lineage tracing and busulfan challenge show that these are SSCs that contribute to steady state spermatogenesis as well as regeneration following chemical injury. EOMES+ SSCs have a lower proliferation index in wild-type than in Plzf lu/lu mice, suggesting that PLZF regulates their proliferative activity and that EOMES+ SSCs are lost through proliferative exhaustion in Plzf lu/lu mice. Single cell RNA sequencing of EOMES+ cells from Plzf +/+ and Plzf lu/lu mice support the conclusion that SSCs are hierarchical yet heterogeneous.
Assuntos
Proteína com Dedos de Zinco da Leucemia Promielocítica/genética , Espermatogênese/genética , Espermatogônias/citologia , Proteínas com Domínio T/genética , Animais , Diferenciação Celular/genética , Proliferação de Células/genética , Autorrenovação Celular/genética , Regulação da Expressão Gênica no Desenvolvimento/genética , Masculino , Camundongos , RNA-Seq , Espermatogônias/metabolismo , Células-Tronco/citologia , Células-Tronco/metabolismo , Testículo/crescimento & desenvolvimentoRESUMO
Cancer-associated fibroblasts (CAF) are major players in the progression and drug resistance of pancreatic ductal adenocarcinoma (PDAC). CAFs constitute a diverse cell population consisting of several recently described subtypes, although the extent of CAF heterogeneity has remained undefined. Here we use single-cell RNA sequencing to thoroughly characterize the neoplastic and tumor microenvironment content of human and mouse PDAC tumors. We corroborate the presence of myofibroblastic CAFs and inflammatory CAFs and define their unique gene signatures in vivo. Moreover, we describe a new population of CAFs that express MHC class II and CD74, but do not express classic costimulatory molecules. We term this cell population "antigen-presenting CAFs" and find that they activate CD4+ T cells in an antigen-specific fashion in a model system, confirming their putative immune-modulatory capacity. Our cross-species analysis paves the way for investigating distinct functions of CAF subtypes in PDAC immunity and progression. SIGNIFICANCE: Appreciating the full spectrum of fibroblast heterogeneity in pancreatic ductal adenocarcinoma is crucial to developing therapies that specifically target tumor-promoting CAFs. This work identifies MHC class II-expressing CAFs with a capacity to present antigens to CD4+ T cells, and potentially to modulate the immune response in pancreatic tumors.See related commentary by Belle and DeNardo, p. 1001.This article is highlighted in the In This Issue feature, p. 983.
Assuntos
Apresentação de Antígeno/imunologia , Fibroblastos Associados a Câncer/imunologia , Fibroblastos Associados a Câncer/metabolismo , Carcinoma Ductal Pancreático/etiologia , Carcinoma Ductal Pancreático/metabolismo , Neoplasias Pancreáticas/etiologia , Neoplasias Pancreáticas/metabolismo , Animais , Fibroblastos Associados a Câncer/patologia , Carcinoma Ductal Pancreático/patologia , Imunofluorescência , Antígenos de Histocompatibilidade Classe II/genética , Antígenos de Histocompatibilidade Classe II/imunologia , Humanos , Camundongos , Modelos Biológicos , Neoplasias Pancreáticas/patologia , Análise de Célula Única , Microambiente Tumoral/imunologiaRESUMO
The protein kinase catalytic domain is one of the most abundant domains across all branches of life. Although kinases share a common core function of phosphoryl-transfer, they also have wide functional diversity and play varied roles in cell signaling networks, and for this reason are implicated in a number of human diseases. This functional diversity is primarily achieved through sequence variation, and uncovering the sequence-function relationships for the kinase family is a major challenge. In this study we use a statistical inference technique inspired by statistical physics, which builds a coevolutionary "Potts" Hamiltonian model of sequence variation in a protein family. We show how this model has sufficient power to predict the probability of specific subsequences in the highly diverged kinase family, which we verify by comparing the model's predictions with experimental observations in the Uniprot database. We show that the pairwise (residue-residue) interaction terms of the statistical model are necessary and sufficient to capture higher-than-pairwise mutation patterns of natural kinase sequences. We observe that previously identified functional sets of residues have much stronger correlated interaction scores than are typical.
Assuntos
Evolução Molecular , Proteínas Quinases/química , Proteínas Quinases/metabolismo , Motivos de Aminoácidos , Método de Monte Carlo , Mutação , Probabilidade , Proteínas Quinases/genéticaRESUMO
Understanding the complex mutation patterns that give rise to drug resistant viral strains provides a foundation for developing more effective treatment strategies for HIV/AIDS. Multiple sequence alignments of drug-experienced HIV-1 protease sequences contain networks of many pair correlations which can be used to build a (Potts) Hamiltonian model of these mutation patterns. Using this Hamiltonian model, we translate HIV-1 protease sequence covariation data into quantitative predictions for the probability of observing specific mutation patterns which are in agreement with the observed sequence statistics. We find that the statistical energies of the Potts model are correlated with the fitness of individual proteins containing therapy-associated mutations as estimated by in vitro measurements of protein stability and viral infectivity. We show that the penalty for acquiring primary resistance mutations depends on the epistatic interactions with the sequence background. Primary mutations which lead to drug resistance can become highly advantageous (or entrenched) by the complex mutation patterns which arise in response to drug therapy despite being destabilizing in the wildtype background. Anticipating epistatic effects is important for the design of future protease inhibitor therapies.
Assuntos
Farmacorresistência Viral/genética , Protease de HIV/genética , Sequência de Aminoácidos , Simulação por Computador , Epistasia Genética/genética , Infecções por HIV , Protease de HIV/metabolismo , HIV-1/genética , Humanos , Modelos Moleculares , Mutação , Alinhamento de SequênciaRESUMO
Potts Hamiltonian models of protein sequence co-variation are statistical models constructed from the pair correlations observed in a multiple sequence alignment (MSA) of a protein family. These models are powerful because they capture higher order correlations induced by mutations evolving under constraints and help quantify the connections between protein sequence, structure, and function maintained through evolution. We review recent work with Potts models to predict protein structure and sequence-dependent conformational free energy landscapes, to survey protein fitness landscapes and to explore the effects of epistasis on fitness. We also comment on the numerical methods used to infer these models for each application.
Assuntos
Evolução Molecular , Modelos Moleculares , Proteínas/genética , Proteínas/metabolismo , Epistasia Genética , Proteínas/química , TermodinâmicaRESUMO
We describe binding free energy calculations in the D3R Grand Challenge 2015 for blind prediction of the binding affinities of 180 ligands to Hsp90. The present D3R challenge was built around experimental datasets involving Heat shock protein (Hsp) 90, an ATP-dependent molecular chaperone which is an important anticancer drug target. The Hsp90 ATP binding site is known to be a challenging target for accurate calculations of ligand binding affinities because of the ligand-dependent conformational changes in the binding site, the presence of ordered waters and the broad chemical diversity of ligands that can bind at this site. Our primary focus here is to distinguish binders from nonbinders. Large scale absolute binding free energy calculations that cover over 3000 protein-ligand complexes were performed using the BEDAM method starting from docked structures generated by Glide docking. Although the ligand dataset in this study resembles an intermediate to late stage lead optimization project while the BEDAM method is mainly developed for early stage virtual screening of hit molecules, the BEDAM binding free energy scoring has resulted in a moderate enrichment of ligand screening against this challenging drug target. Results show that, using a statistical mechanics based free energy method like BEDAM starting from docked poses offers better enrichment than classical docking scoring functions and rescoring methods like Prime MM-GBSA for the Hsp90 data set in this blind challenge. Importantly, among the three methods tested here, only the mean value of the BEDAM binding free energy scores is able to separate the large group of binders from the small group of nonbinders with a gap of 2.4 kcal/mol. None of the three methods that we have tested provided accurate ranking of the affinities of the 147 active compounds. We discuss the possible sources of errors in the binding free energy calculations. The study suggests that BEDAM can be used strategically to discriminate binders from nonbinders in virtual screening and to more accurately predict the ligand binding modes prior to the more computationally expensive FEP calculations of binding affinity.